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

e3gyp2002

clc;
clear all;
3 |$ }4 o  Q+ [& p* R# j1 e2 Bfs=2000;
  o% Q3 c- d% ?, @) ~2 tdt=1./fs;
N=512;
7 d' ^; I5 ?( g& l2 _8 ^; e! Dn=0:1:N-1;
t=n*dt;
Wd=hamming(N);
f=n*fs./N;
8 z' S9 w( v: m% x3 d# u0 ~s=2*cos(2*30*t);
sc=s.*Wd;
( X+ I7 \# y6 Q9 c8 s, F& o& ]figure(1)
subplot(211)
9 \0 K/ h# R4 Z/ _plot(t,s)
. u0 |2 N9 r2 G" K' zsubplot(212)
4 A( ~0 w) l: u6 [5 v9 D: Gplot(t,sc)
% S1 D6 M! [8 o: [& @+ As1=fft(s);
9 i" e* k, [$ Us2=fft(sc);
figure(2)
subplot(211)
. v2 [2 q$ t0 M6 t6 `0 L: [  [4 Q8 dplot(f-fs/2,abs([s1(N/2+1:N) s1(1:N/2)]));
' ^# l9 \6 A: J* G  ^axis([-400,400,0,500]);
subplot(212)
plot(f-fs/2,abs([s2(N/2+1:N) s2(1:N/2)]));
( e$ ?' w: g5 ^% V4 f* e/ [8 eaxis([-400,400,0,500]);
! M5 o0 P3 L1 A2 d
出现错误：
0 }1 o6 H1 I; R& I, j0 Q??? Error using ==> unknown
* B: |! c0 s0 W3 v7 o- H2 q4 {Matrix dimensions must agree.
+ U7 S8 }$ U; R0 r8 s
  _) S, k$ \5 fError in ==> Untitled at 22
sc=s.*Wd;
" K" s9 C4 r3 ~" f, T5 G, r  C+ i
) X3 y* j& [) I2 Y. [我是初学者。请问信号与窗函数是怎么乘的？是不是产生窗函数的方法错了？
谢谢！-%-40740-%-

norika

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

qizi1209

0000xiexie