TA的每日心情 | 开心 2023-6-2 15:15 |
|---|
签到天数: 1 天 [LV.1]初来乍到
|
你这个程序可这样改:' q: g( P" t) u. A
N = 35;
: M/ Y4 U$ a# {9 B+ vH =[ones(1,9),0.4,zeros(1,15),0.4,ones(1,9)]; %确定抽样点的幅度
, `( U2 d. P8 \* A* ]k = 0:N-1;2 {: f, m2 V% z# b2 m
A =exp(-1i*pi*k*(N-1)/N); %抽样点相位大小
1 T; n4 j) B7 s3 Z+ {, _HK =H.*A; %求抽样点的H(k)4 u3 X1 n/ R" r* d% s* p3 z# X
hn =ifft(HK,N); %求出FIR的单位冲激响应h(n)
- E$ z: M' X: Q* J4 |6 pfreqz(hn,1,256); %画出幅频相频曲线
8 {1 [6 j( A; C# J& Z( Wfigure(2);
, ?% ^6 Y' B" L6 Qstem(real(hn),'.'); %绘制单位冲激响应的实部
7 C# H& A* u4 X) n$ O2 S8 |/ Bline([0,35],[0,0]);xlabel('n');ylabel('Real(h(n))'); ? k3 r' x' k# s) ?
|
|