请教高手啊!主程序改怎么写？？？

sxd1168

谐波信号的 HHT 分析
   设间谐波信号为
                     2
            s(t) =  ∑ Ai sin(2πmi ft)           (13)
                    i=1
式中 A1=1V，A2=0.3V，m1=1，m2=5.5，f=50Hz。
采样频率为3200Hz。
要对这个信号进行HHT分析，我用的是下面三个自定义程序，请教高手给怎么写主程序调用呢？？

sxd1168

我用的上面几个附件里有人说还是有小错误，希望高手指点一下，我在写毕业论文，好恼火啊。。。

pf317

这是我在论坛里找的，我现在也用的是这个，你用试试吧，里面的几个程序要下工具箱才能用
t=1:500;
t=t*1/400;
x1=sin(2*pi*20*t);
x2=0.4*sin(2*pi*40*t+140);
f=400;
y=linspace(1,2,500);
z=x1+x2+y
%在后面添加一个趋势向可以把三个分量完全分解出来
imf=emd(z);
emd_visu(z,1:length(z),imf,1)
figure (4)
dt=1/f;
h = nspab(imf',500,0,500,dt);
h(1:2,:)=0;


ms=mspc(h,t)

subplot(224)
plot((1:length(ms)),ms)
xlabel('频率(Hz)','FontSize',12);
ylabel('幅值','FontSize',12);
subplot(211)
surf(h(1:200,1:end-20))
shading interp
xlabel('时间(点数)','FontSize',12);
Ylabel('频率(Hz)','FontSize',12);
zlabel('幅值','FontSize',12);
title('a','FontSize',12)


view([-75,25])
yt=subplot(223)
imagesc(h(1:200,:))
Ylabel('频率(Hz)','FontSize',12);
xlabel('时间(点数)','FontSize',12);
set(yt,'ydir','nor')
title('b','FontSize',12);title('c','FontSize',12);