N=3;x=111; nx=0N-1;K=1000;dw=2piK;k=floor0K-05;X=xexp- jdwnxk; figuresubplot211;stemnxx;grid;xlabel时间;ylabel幅度;title 时间域波形;subplot212;plotkdwabsX;grid;xlabel频率弧度;ylabel幅度;title 系统的幅频响应;将以上代码进行注释

%定义参数 N=3; %信号长度 x=[1,1,1]; %信号序列 nx=0:N-1; %信号的时间序列 K=1000; %FFT点数 dw=2*pi/K; %频率分辨率

%计算FFT k=floor(0:(K-0.5)); %频率序列 X=xexp(- jdw*nx'*k); %计算频域波形

%绘图 figure subplot(2,1,1); stem(nx,x);grid; xlabel('时间'); ylabel('幅度'); title ('时间域波形');

subplot(2,1,2); plot(k*dw,abs(X));grid; xlabel('频率（弧度）'); ylabel('幅度'); title ('系统的幅频响应');