EDA365电子论坛网
标题:
我也有个程序出不来仿真结果，大神能不能帮助帮助？
[打印本页]
作者:
lahhse
时间:
2022-8-9 13:15
标题:
我也有个程序出不来仿真结果，大神能不能帮助帮助？
N = 33;
# \# i, I0 f9 z' f8 C8 z$ C
H =[ones(1,9),0.4,zeros(1,15),0.4,ones(1,9)]; %确定抽样点的幅度
3 L$ @/ n L( z0 q6 v9 p$ H
k = 0:N-1;
' v# r' l8 I4 _* L$ D5 A
A =exp(-j*pi*k*(N-1)/N); %抽样点相位大小
8 H" ^" z8 D# g, j
HK =H.*A; %求抽样点的H(k)
4 c* `1 R* f; D
hn =ifft(HK,N); %求出FIR的单位冲激响应h(n)
* f" z4 u- c8 g
freqz(hn,1,256); %画出幅频相频曲线
# M9 O# O+ h" ^7 K: q
figure(2);
# p1 A8 r2 D, V. O8 b
stem(real(hn),'.'); %绘制单位冲激响应的实部
; J5 A+ q9 p1 B
line([0,35],[0,0]);xlabel('n');ylabel('Real(h(n))');
) i& R, R2 }* v! ]. h+ M. C" r
作者:
Blah
时间:
2022-8-9 14:20
你这个程序可这样改：
" ]9 h9 u3 U# ?
N = 35;
& `) W9 N% l( f+ X( |( a1 e6 {
H =[ones(1,9),0.4,zeros(1,15),0.4,ones(1,9)]; %确定抽样点的幅度
: ]* x! ?7 K* z# d" g( P4 p
k = 0:N-1;
# b0 q+ F4 A. V$ `
A =exp(-1i*pi*k*(N-1)/N); %抽样点相位大小
4 J1 Y# ~6 h3 f: M
HK =H.*A; %求抽样点的H(k)
! `5 K' J3 V, C
hn =ifft(HK,N); %求出FIR的单位冲激响应h(n)
) v8 v8 S+ H, `1 u& v5 B& z
freqz(hn,1,256); %画出幅频相频曲线
2 T: J% q* J2 z( W1 \5 m' k
figure(2);
) I1 J' y% w+ k
stem(real(hn),'.'); %绘制单位冲激响应的实部
- G0 m( }" |% M7 G8 a
line([0,35],[0,0]);xlabel('n');ylabel('Real(h(n))');
6 y9 q) m) @: K3 ^& b
作者:
land
时间:
2022-8-9 16:39
可以试一下楼上的方法
4 Y9 U+ I( I0 A2 `
作者:
Memory00
时间:
2022-8-9 17:09
再看看别人是怎么说的
, i( d/ Z* ^9 D9 L8 l
欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/)
Powered by Discuz! X3.2