LMD.zip_lmd_lmd 分解_sleepmnm_信号分解_信号的分解

clc clear all %close all %%% adams 仿真压缩机数据 % temp='F:\lianganjiasudu.txt'; % %temp='F:\jiechuli.txt'; % %temp='F:\quzhoujiasudu.txt'; % yy= importdata(temp); % xx=yy.data; % x=xx(50:1857,1); % 轴承正常状态数据 %load('f:\98.mat') %x=X098_DE_time(1:1024); %ts=1/fs; %t=0:r-1; %%%压缩机测试数据 % temp='F:\yasuojiceshishuju\zhoucheng1\t1n1\BLACKBOX\BD0001.DAT'; % % [filename , pathname ]= uigetfile ( '*.dat','dat'); %temp=strcat(pathname ,filename); %fid=fopen(temp); %a=fread(fid,'short'); %fclose(fid); [filename , pathname ]= uigetfile ( '*.*','当前文件名'); fid=fopen(strcat(pathname ,filename),'rb'); a=freaad(fid,'short'); a1=(a(:,1)); n=1; x=a1(73+1024*n:1097+1024*n-1,1); [r,s]=size(x); fs=r/0.3; t=0:r-1; y=x'; figure;set(gcf,'outerposition',get(0,'screensize')); plot(x); %-------------------------- for i=1:s [pf,a,si]=lmd2(y(i,:)); [M,N] = size(pf); figure;set(gcf,'outerposition',get(0,'screensize')); for k1 = 1:M subplot(M,1,k1), plot(t,pf(k1,:)); title('LMD分解结果'); xlabel('Time'); end figure;set(gcf,'outerposition',get(0,'screensize')); for k1 = 1:M-1 subplot(M-1,1,k1), plot(t,si(k1,:)); end figure;set(gcf,'outerposition',get(0,'screensize')); for k1 = 1:M-1 subplot(M-1,1,k1), plot(t,a(k1,:)); end figure;set(gcf,'outerposition',get(0,'screensize')); for k1 = 1:M-1 subplot(M-1,1,k1), plot(t,abs(fft(a(k1,:)))); end ff=shunpin(si).*fs;%实际频率要乘以采样率 figure;set(gcf,'outerposition',get(0,'screensize')); for k1 = 1:M-1 subplot(M-1,1,k1), plot(t,ff(k1,:)); end end