MATLAB中如何对一个正弦函数矩阵（行为时间采样点，列为频率逐渐增大）的逐行求EMD

Question

fs=400;%采样率
ts=1/fs;
N=500;%采样点数（100倍）
M=0;%采样数据段的起始位置
n=M:N-1;x=zeros(100,500);
t=n/fs;%信号时间序列
for f0=1:100
x=sin(f0*2*pi*t);
y(f0,:)=x;
end
plot(t,y);
for k=1:100;
imf=emd(y(k,:));
r=size(imf,1);%提取imf行数
for i=2:(r+1);
eng_error=sum(abs(imf(i,:)))/N;
sum=sum+eng_error;
end
end
我要对每个频率的正弦函数进行EMD分解，在对得出的IMF求误差，得到一个误差关于频率的图像。但是每个频率EMD分解的IMF个数也不确定，不知道怎么 用矩阵了，求大侠帮忙！

Answer 1

追问

这是个发病毒的片子，fuck uuuu