短时傅立叶窗函数问题

tongxiaoxia

%采样频率
fs=25000;
%实验信号
n=length(A0); t=1:n; s=A0;
%时域波形
figure(1); plot(0:1/fs:(n-1)/fs,s); xlabel('时间 t/n'); ylabel('幅值 A/v');
%设置窗函数1
h1=window(165,'hamming');
%计算短时傅立叶变换
s=hilbert(s); [tfr,t,f]=tfrstft(s,1:500,500,h1);
%时域表示
figure(2); contour(t,2*f(1:length(A0)/2),abs(tfr)); xlabel('时间 t'); ylabel('频率 f');
%设置窗函数2
h2=window(17,'hamming');
%计算短时傅立叶变换
s=hilbert(s); [tfr,t,f]=tfrstft(s,1:500,500,h2);
%时域表示
figure(3); contour(t,2*f(1:length(A0)/2),abs(tfr)); xlabel('时间 t'); ylabel('频率 f');

运行结果为：
??? Error using ==> feval
Argument must contain a string or function_handle.
Error in ==> window at 52
    w = feval(wname,N,varargin{:});
Error in ==> duanshi at 18
h1=window(165,'hamming');

这问题如何解决

[ 本帖最后由 ChaChing 于 2009-3-27 13:26 编辑 ]

songzy41

是否应把window函数写成:
h1=window(@hamming,165);