TA的每日心情 | 开心 2023-6-2 15:15 |
---|
签到天数: 1 天 [LV.1]初来乍到
|
你这个程序可这样改:1 B% O2 w+ m# U, l$ N
N = 35;" x/ E7 D+ @, n/ k' e# I
H =[ones(1,9),0.4,zeros(1,15),0.4,ones(1,9)]; %确定抽样点的幅度* C7 `5 t2 {1 s2 R* ^
k = 0:N-1;
) i* D# p$ |. }5 O8 nA =exp(-1i*pi*k*(N-1)/N); %抽样点相位大小 d3 M$ q. m: C, X
HK =H.*A; %求抽样点的H(k)# C8 W, C! O% ]$ B% J, V( r6 u
hn =ifft(HK,N); %求出FIR的单位冲激响应h(n)8 A% g0 t& O' A- S$ D/ M
freqz(hn,1,256); %画出幅频相频曲线) P. V# E) Q( V' ]4 v
figure(2);
6 V; w- m/ ~1 fstem(real(hn),'.'); %绘制单位冲激响应的实部& ^- w8 ?# ~1 y6 h) B
line([0,35],[0,0]);xlabel('n');ylabel('Real(h(n))');6 A3 v8 T% k2 B) Y1 n2 u, v
|
|