原帖由 咕噜噜 于 2007-5-3 15:19 发表
我这样理解你的程序,x是一个矩阵,t, sig1, sig2是三个未知矩阵,他们三个分别为矩阵x的第3*i+1,3*i+2,3*i+3个值组成?
假如是这样的话,那肯定是你的x矩阵没有赋值。我似乎觉得你给的这一段不完整, ...
fs=1/0.008
N=length(s);
for i=0:N/3-1
t(i+1)=x(3*i+1);
sig1(i+1)=x(3*i+2);
sig2(i+1)=x(3*i+3);
end
figure(1);
subplot(2,2,1);
plot(t,sig1,'LineWidth',2);
xlabel('时间 t/s');
ylabel('电压 V/mV');
sig=hilbert(sig1);
sig=sig-mean(sig);
[ifr,tt]=instfreq(sig',2:length(sig)-1);
subplot(2,2,2);
plot(tt/fs,(1-ifr)*fs/2, 'LineWidth',2);
xlabel('时间 t/s');
ylabel('瞬时频率 f/Hz');
上面是所有的程序。我也觉得是x没赋值,但我不知道应该怎么写才算正确的赋值 |