【求助】关于STFT 和WVD的仿真

magic0elf

跪求啊   谢谢   着急啊
这是我在网上下的两个程序，运行后貌似少了信号源还是 什么 ？求高手看下，求指导。
STFT部分：  三个怎么设
function [T,F,SP] = STFT(D,M,fm);

% y: input signal
% D: the data matrix: D(:,1)=t; D(:,2)=y;
% M: the lenght of the FT window


% [0 fm] The Frequency range for the signal
% SP: the spectogram (short time fourier transform) of D

%fm % the maximum freq

t=D(:,1);y=D(:,2);
N=length(y);

W=hanning(M);
ye=zeros(N+M,1);
ye(M/2+1:M/2+N)=y;
%yea=hilbert(ye);
dt=t(2)-t(1);% time resolution
fs=1/dt; % the sampling frequency

df=fs/M; %frequency resolution
L=round(fm/df); % the index of the max freq.
f=[0:L-1]*df;

i=1;
while i<=N
   
   s=ye(i:i+M-1).*W;
   sa=hilbert(s);
   
   S=fft(sa);
   B(:,i)=S(1:L);
   i=i+1;% increase by one (next sample)
end

SP = (1/N)*abs(B.*B);
SP = SP/max(max(SP)); % normalize it
T  = t;
F  = f;



% mesh(T,F,SP)
% contour(T,F,SP)



还有一个是function [T,F,Z]=Wigner(D,fm)  



谢谢

ChaChing

1.不喜欢重覆帖
2.信号源应该使用者输入
3.不是有输入的注解吗?

magic0elf

回复 2 # ChaChing 的帖子

我是看到没人回我才发的        
我自己试了一个不行，所以请教下，你能给点什么建议吗？

ChaChing

回复 3 # magic0elf 的帖子

求助完整格式：出错代码和出错提示

magic0elf

回复 4 # ChaChing 的帖子

报错
??? Input argument "D" is undefined.

Error in ==> STFT at 11
t=D(:,1);y=D(:,2);

taosky

感觉这是一个函数 是要自己编程调用的  直接是不能运行的

ChaChing

回复 5 # magic0elf 的帖子

Ref: 4F, 常见的程序出错问题整理 http://forum.vibunion.com/thread-46001-1-1.html
from http://forum.vibunion.com/home-s ... -blog-id-18250.html