TA的每日心情 | 开心 2023-6-2 15:15 |
|---|
签到天数: 1 天 [LV.1]初来乍到
|
你这个程序可这样改:" f- o6 m6 X t7 {+ ^. k% Y- e
N = 35;) W3 |7 e3 Q! V" m/ H0 u( i
H =[ones(1,9),0.4,zeros(1,15),0.4,ones(1,9)]; %确定抽样点的幅度
( Q0 {' D8 H; F/ E' e1 k8 W4 q9 Mk = 0:N-1;
- w/ F7 V2 ?! g* V8 wA =exp(-1i*pi*k*(N-1)/N); %抽样点相位大小5 L: j |, Z5 E: V7 D6 r
HK =H.*A; %求抽样点的H(k)
2 |$ D; H4 G$ lhn =ifft(HK,N); %求出FIR的单位冲激响应h(n); {& O% ^+ O: T
freqz(hn,1,256); %画出幅频相频曲线. c! A; B2 `6 s) V, S% s
figure(2);( L9 T; w3 |( L: q& M! [' y
stem(real(hn),'.'); %绘制单位冲激响应的实部# \ F- |/ g' R+ |
line([0,35],[0,0]);xlabel('n');ylabel('Real(h(n))');
0 B2 F8 J# b7 C. s/ L# ^ |
|