N1=2;x1=10;n=01N1-1; k1=01N1-1;wn=exp-j2piN1;nk=nk1;wnk=wn^nk;xk=x1wnk;figure subplot211;stem1234x1x1;grid;xlabel时间;ylabel幅度;title 时间域波形; subplot212;stemk12piN1absxk;grid;xlabel频率弧度;ylabel幅度;title 系统的

% N1表示采样点数 % x1为输入信号 % n为时间序列 % k1为频率序列 % wn为旋转因子 % nk为矩阵，表示n与k1的乘积 % wnk为旋转因子的nk次幂 % xk为傅里叶变换结果，即系统的幅频响应 % subplot为画图函数，第一个参数表示图的行数，第二个参数表示图的列数，第三个参数表示当前图的位置 % stem为绘制离散信号的函数，第一个参数为横坐标，第二个参数为纵坐标，'*'表示绘制的点的形状 % xlabel为横坐标标签，ylabel为纵坐标标签，title为图的标题 % abs为求复数的模，即幅度