matlab加窗函数减少频谱泄露的问题

e3gyp2002

clc;
clear all;
* P2 p) R. W/ @6 dfs=2000;
dt=1./fs;
+ f% M% }0 e0 @7 `# {N=512;
n=0:1:N-1;
0 ^9 b; z! V# Dt=n*dt;
' b/ l- i1 f/ R. l; }Wd=hamming(N);
; r8 X- d; U) {( i, c8 x! J1 n+ a: of=n*fs./N;
s=2*cos(2*30*t);
sc=s.*Wd;
, ~) s# M7 _: R$ _figure(1)
7 r3 J. L6 C/ asubplot(211)
8 Z; U6 N5 d* @$ t+ z) l6 T! [4 Eplot(t,s)
) H- d9 Q3 w. r  usubplot(212)
plot(t,sc)
s1=fft(s);
s2=fft(sc);
figure(2)
subplot(211)
9 ^9 H+ Y( [" ^5 Lplot(f-fs/2,abs([s1(N/2+1:N) s1(1:N/2)]));
axis([-400,400,0,500]);
subplot(212)
+ k" o) }4 R0 B6 |7 Q1 X9 nplot(f-fs/2,abs([s2(N/2+1:N) s2(1:N/2)]));
, @/ M1 J9 \' E% e2 {axis([-400,400,0,500]);
9 ]& J3 a1 l, q! Z# x# z
( e1 A4 l4 g% y" X/ [) }' U! d7 |出现错误：
??? Error using ==> unknown
Matrix dimensions must agree.

Error in ==> Untitled at 22
sc=s.*Wd;
5 U! e7 {6 s$ c) y; V
; O: \' d$ |2 K- o3 x# C) C4 k我是初学者。请问信号与窗函数是怎么乘的？是不是产生窗函数的方法错了？
% F/ u1 j1 ]# X  E谢谢！-%-40740-%-

norika

恩，可以查找一下附有程序代码的文献，我用matlab写过这方面的程序，不难

qizi1209

0000xiexie