EDA365电子论坛网

标题: 我也有个程序出不来仿真结果,大神能不能帮助帮助? [打印本页]

作者: lahhse    时间: 2022-8-9 13:15
标题: 我也有个程序出不来仿真结果,大神能不能帮助帮助?
N = 33;! `8 l; n' K- M6 ^
H =[ones(1,9),0.4,zeros(1,15),0.4,ones(1,9)]; %确定抽样点的幅度
* V1 z6 B) Z& ~7 t5 Jk = 0:N-1;
0 V+ m# r, t9 N; e, O  D. z( WA =exp(-j*pi*k*(N-1)/N); %抽样点相位大小
: ?9 `4 H4 R! i- O2 G& VHK =H.*A; %求抽样点的H(k)
& O) @  U: s& Z  Mhn =ifft(HK,N); %求出FIR的单位冲激响应h(n)  p3 g8 n1 E3 G
freqz(hn,1,256); %画出幅频相频曲线
5 ]2 k9 X4 Y- V2 k  A1 _figure(2);
3 t3 E; ^, o6 s! u! Bstem(real(hn),'.'); %绘制单位冲激响应的实部, T  \* Q; D' i; K7 E" h
line([0,35],[0,0]);xlabel('n');ylabel('Real(h(n))');  q9 F6 }& P  l

作者: Blah    时间: 2022-8-9 14:20
你这个程序可这样改:
2 S: l& B, M4 c% d) v; f  zN = 35;2 |# ~: n5 @' r, N! H
H =[ones(1,9),0.4,zeros(1,15),0.4,ones(1,9)]; %确定抽样点的幅度) A5 |; z6 O/ m  }% p
k = 0:N-1;" \1 B5 ?5 m% a; n! A6 j' L8 d; ^
A =exp(-1i*pi*k*(N-1)/N); %抽样点相位大小1 T2 \% B' m7 y1 y+ G, C" f
HK =H.*A; %求抽样点的H(k)
) c- `/ L0 W; i# yhn =ifft(HK,N); %求出FIR的单位冲激响应h(n)6 G1 q; P) }. _& D$ ]8 Y/ x5 S; k: X
freqz(hn,1,256); %画出幅频相频曲线
9 |' S8 e  R3 S& R1 Cfigure(2);" X. K' o2 k7 q: u! T# o
stem(real(hn),'.'); %绘制单位冲激响应的实部1 i4 l$ a, Y; S8 D) J& d
line([0,35],[0,0]);xlabel('n');ylabel('Real(h(n))');5 g- @5 p+ @" H5 h

作者: land    时间: 2022-8-9 16:39
可以试一下楼上的方法7 w. r5 D8 A  x/ ~, Z5 V) b( Y  e

作者: Memory00    时间: 2022-8-9 17:09
再看看别人是怎么说的" n( H2 J, N  E% c; W$ c1 s





欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/) Powered by Discuz! X3.2