我也有个程序出不来仿真结果，大神能不能帮助帮助？

lahhse

N = 33;
1 Y' M4 H7 x' r4 T( X8 nH =[ones(1,9),0.4,zeros(1,15),0.4,ones(1,9)]; %确定抽样点的幅度
k = 0:N-1;
: _8 n: H) B. d, x' fA =exp(-j*pi*k*(N-1)/N); %抽样点相位大小
HK =H.*A; %求抽样点的H(k)
hn =ifft(HK,N); %求出FIR的单位冲激响应h(n)
freqz(hn,1,256); %画出幅频相频曲线
2 l! W2 }( l! |% m# g5 G: wfigure(2);
stem(real(hn),'.'); %绘制单位冲激响应的实部
line([0,35],[0,0]);xlabel('n');ylabel('Real(h(n))');

Blah

你这个程序可这样改：
N = 35;
/ W/ X4 H# c) i2 v8 S5 aH =[ones(1,9),0.4,zeros(1,15),0.4,ones(1,9)]; %确定抽样点的幅度
k = 0:N-1;
A =exp(-1i*pi*k*(N-1)/N); %抽样点相位大小
+ h# [  Q3 B$ _; A& L; X# K) @HK =H.*A; %求抽样点的H(k)
hn =ifft(HK,N); %求出FIR的单位冲激响应h(n)
7 v+ O: q0 g9 n- ]4 m4 x, N5 y1 xfreqz(hn,1,256); %画出幅频相频曲线
figure(2);
stem(real(hn),'.'); %绘制单位冲激响应的实部
/ K8 n  E: ?% p6 K5 tline([0,35],[0,0]);xlabel('n');ylabel('Real(h(n))');

land

可以试一下楼上的方法

Memory00

再看看别人是怎么说的