TA的每日心情 | 开心 2023-6-2 15:15 |
---|
签到天数: 1 天 [LV.1]初来乍到
|
你这个程序可这样改:1 M4 U' a: x' r: r7 [: R- F
N = 35;" r# F' t2 J4 l* K* ^! h
H =[ones(1,9),0.4,zeros(1,15),0.4,ones(1,9)]; %确定抽样点的幅度0 ?( B7 X! H. d
k = 0:N-1;
' e! ]. {1 x5 x7 u3 R. |9 W zA =exp(-1i*pi*k*(N-1)/N); %抽样点相位大小+ f2 t6 E6 G$ Y n1 R* i
HK =H.*A; %求抽样点的H(k)! _6 d" g6 }* |
hn =ifft(HK,N); %求出FIR的单位冲激响应h(n)' c; m4 A' ~: n8 d7 D' x
freqz(hn,1,256); %画出幅频相频曲线# r, T8 s3 r8 D3 I9 O& l. J j
figure(2);* D+ a' D& f! c; O( Z: v1 f% g/ i
stem(real(hn),'.'); %绘制单位冲激响应的实部# |- p, ]* J! o# f9 m2 z4 w
line([0,35],[0,0]);xlabel('n');ylabel('Real(h(n))');) P+ O- N, G/ ~3 N- e
|
|