TA的每日心情 | 开心 2023-6-2 15:15 |
---|
签到天数: 1 天 [LV.1]初来乍到
|
你这个程序可这样改:
% @" J$ L7 F' {9 X7 PN = 35;4 u1 _. B: Z' e8 ], H0 D7 i
H =[ones(1,9),0.4,zeros(1,15),0.4,ones(1,9)]; %确定抽样点的幅度+ j: {8 W' Z( w! O# a
k = 0:N-1;
: v( J+ C- ]: m" a& x5 s# zA =exp(-1i*pi*k*(N-1)/N); %抽样点相位大小
6 g+ u0 g3 ]' _* H9 W$ c1 [! pHK =H.*A; %求抽样点的H(k)* r( ]7 I- T+ t
hn =ifft(HK,N); %求出FIR的单位冲激响应h(n)3 w; w# K& t2 f% v9 K/ G- m7 o
freqz(hn,1,256); %画出幅频相频曲线
( v( l; p' n! M7 l m' xfigure(2);
8 j$ j& o* \9 k$ q, K E8 r4 ?stem(real(hn),'.'); %绘制单位冲激响应的实部7 E8 e0 T; G$ i, O* ~
line([0,35],[0,0]);xlabel('n');ylabel('Real(h(n))');
/ t( [* w& m6 G) m( r |
|