我在用EEMD处理凯斯西储大学的轴承故障数据，具体代码如下

close all;

clear all;

clc;

load('C:\Users\Administrator\Desktop\123\130.mat');%导入数据

fs=12000;%采样率

N=10240;%采样点数

n=0:N-1;

t=n/fs;%信号时间序列

X=X130_DE_time(1:N);%装载驱动计数端的内圈故障数据

y=X;%信号幅值序列

k_in=kurtosis(y);%峭度系数，正常轴承为3左右

figure;%画原始信号时域和频域图

subplot(211);plot(t,y);title('原始信号时域波形');

subplot(212);hua_fft(y,fs,1);title('原始信号频谱');

figure;%原始信号的包络谱

subplot(211);hua_baoluo(y,fs,1);title('原始信号包络谱');

subplot(212);hua_baoluo(y,fs,1);title('原始信号部分频段包络谱');

subplot(212);hua_baoluo(y,fs,1);axis([0 400 0 0.7]);title('原始信号部分频段包络谱');

NE=50;

Nstd=0.01;

imf=eemd(y,0.01,50);%经验模态分解

figure;%画出所有imf及残值

%画出每个IMF分量的辐频曲线

[a,b]=size(imf);

for i=1:N-1;

subplot(a+1,1,i);

plot(t,imf(i,:))

ylabel(['imf',int2str(i)])

end

%画出最后一个剩余分量residual的辐频曲线

subplot(a,1,a);

plot(t,imf(a,:))

ylabel(['r'])

figure;%前五个IMF分量

subplot(511);plot(t,imf(1,:));title('IMF1时域波形图');

subplot(512);plot(t,imf(2,:));title('IMF2时域波形图');

subplot(513);plot(t,imf(3,:));title('IMF3时域波形图');

subplot(514);plot(t,imf(4,:));title('IMF4时域波形图');

subplot(515);plot(t,imf(5,:));title('IMF5时域波形图');

figure;%前五个IMF分量频谱

subplot(511);hua_fft(imf(1,:),fs,1);title('IMF1频谱');

subplot(512);hua_fft(imf(2,:),fs,1);title('IMF2频谱');

subplot(513);hua_fft(imf(3,:),fs,1);title('IMF3频谱');

subplot(514);hua_fft(imf(4,:),fs,1);title('IMF4频谱');

subplot(515);hua_fft(imf(5,:),fs,1);title('IMF5频谱');

figure;%前五个IMF分量选择频段内的包络谱

xf1=0;%需要查看的包络谱频率段起点频率

xf2=1000;%需要查看的包络谱频率段终止频率

subplot(511);hua_baoluo(imf(1,:),fs,1,xf1,xf2);title('IMF1包络谱');

subplot(512);hua_baoluo(imf(2,:),fs,1,xf1,xf2);title('IMF2包络谱');

subplot(513);hua_baoluo(imf(3,:),fs,1,xf1,xf2);title('IMF3包络谱');

subplot(514);hua_baoluo(imf(4,:),fs,1,xf1,xf2);title('IMF4包络谱');

subplot(515);hua_baoluo(imf(5,:),fs,1,xf1,xf2);title('IMF5包络谱');

出现了错误使用 plot

矢量长度必须相同。

出错 innerguzhangeemd (line 30)

plot(t,imf(i,:))

请问该如何处理