想做个分析谐波的时频程序.大家给看下

caolizhong

我想做个想做个分析谐波的时频程序,
大家给看下有什么错误,这程序产生不了时频图;
或都哪个大虾给个程序啊
s=0:0.001:0.6;
x=sin(2*pi*50*s)+sin(2*pi*120*s);
y=x(1:256);
sig=y.';
%y需要转置
figure(1)
plot(1000*s(1:256),sig)
title('含噪信号')
xlabel('时间/us')
tfr=tfrstft(sig);
figure(2)
imagesc(abs(tfr));


另外这个程序可以用来做时频处理,可是我把原始信号换成自己的后,也就是我给的第一段程序.运行不了,我对                fmconst(128,0.2)函数不明白.请大家给指点下,0.2是归一化频率,那它的实际频率是多少

sig=[fmconst(128,0.2);fmconst(128,0.4)];
plot(real(sig))
tfr=tfrstft(sig);
figure
subplot(211);imagesc(abs(tfr));
subplot(212);imagesc(angle(tfr));

gzslbs2003

想求归一化频率,你用AR模型试一下