|
把程序改了一下,其中有些小错误已改,可以执行了。同时窗函数的调用,我用的是7.1版本的。
sig1=fmlin(128,0,0.4);
sig2=fmlin(128,0.1,0.45);
sig=sig1+sig2;
figure(1)
plot(real(sig),'LineWidth',2);
hold on;
plot(imag(sig),'ro--');
legend('实部','虚部');
xlabel('时间 t');
ylabel('幅值 A');
h1=window(@gausswin,23);
%h1=window(23,'gauss');
[tfr,t,f]=tfrsp(sig,1:128,128,h1);
figure(2)
contour(t,f,tfr);
xlabel('时间 t');
ylabel('频率 f');
h2=window(@gausswin,63);
%h2=window(63,'gauss');
[tfr,t,f]=tfrsp(sig,1:128,128,h2);
figure(3)
contour(t,f,tfr);
xlabel('时间 t');
ylabel('频率 f'); |
评分
-
1
查看全部评分
-
|