TA的每日心情 | 开心 2023-6-2 15:15 |
|---|
签到天数: 1 天 [LV.1]初来乍到
|
你这个程序可这样改: N. Z# |- G B3 F
N = 35;
4 n, W" {3 N3 P2 r1 t4 r$ P gH =[ones(1,9),0.4,zeros(1,15),0.4,ones(1,9)]; %确定抽样点的幅度
* p8 W+ U% [; w6 T$ n4 A, \k = 0:N-1;
2 g/ x5 ?2 @ c7 t. {! sA =exp(-1i*pi*k*(N-1)/N); %抽样点相位大小
0 l! V) Y- Z0 ?& Q5 B- T+ HHK =H.*A; %求抽样点的H(k)0 V) O) x+ p- F! i0 a8 C1 a6 R9 R
hn =ifft(HK,N); %求出FIR的单位冲激响应h(n)* @6 o( T% f7 j8 {1 s" i5 G0 O% P1 q
freqz(hn,1,256); %画出幅频相频曲线- Q6 J' n3 V3 o" T$ o% {
figure(2);
+ f" X( w* i1 H/ N" qstem(real(hn),'.'); %绘制单位冲激响应的实部
2 d# l" [7 t4 Eline([0,35],[0,0]);xlabel('n');ylabel('Real(h(n))');
! a$ [* O4 a0 u% \* y |
|