【求助】关于STFT 和WVD的仿真
跪求啊 谢谢 着急啊这是我在网上下的两个程序，运行后貌似少了信号源还是 什么 ？求高手看下，求指导。
STFT部分：三个怎么设
function = STFT(D,M,fm);
% y: input signal
% D: the data matrix: D(:,1)=t; D(:,2)=y;
% M: the lenght of the FT window
% 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=*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 =Wigner(D,fm)
谢谢
1.不喜欢重覆帖
2.信号源应该使用者输入
3.不是有输入的注解吗? 回复 2 # ChaChing 的帖子
我是看到没人回我才发的
我自己试了一个不行，所以请教下，你能给点什么建议吗？ 回复 3 # magic0elf 的帖子
求助完整格式：出错代码和出错提示 回复 4 # ChaChing 的帖子
报错
??? Input argument "D" is undefined.
Error in ==> STFT at 11
t=D(:,1);y=D(:,2);
感觉这是一个函数 是要自己编程调用的直接是不能运行的 回复 5 # magic0elf 的帖子
Ref: 4F, 常见的程序出错问题整理 http://forum.vibunion.com/thread-46001-1-1.html
from http://forum.vibunion.com/home-space-uid-63979-do-blog-id-18250.html
页:
[1]