信号与系统利用MATLAB分析信号频谱作业(第四章).docVIP

4-31 计算周期矩形序列的频谱。 程序如下： %Program4_1计算周期矩形序列的频谱 N=32;M=4;%定义周期矩形序列的参数 x=[ones(1,M+1),zeros(1,N-2*M-1),ones(1,M)];%产生周期矩形序列 X=fft(x);%计算DFS系数 m=0:N-1; stem(m,real(X));%画出频谱X的实部 title(X[m]的实部); xlabel(m); figure; stem(m,imag(X));%画出频谱X的虚部 title(X[m]的虚部); xlabel(m); xr=ifft(X);%由X计算x以重建序列x[k] figure; stem(m,real(xr));%画出重建序列xr[k] xlabel(k); title(重建的x[k]); (矩形序列的频谱实部)  (矩形序列的频谱虚部) (由DFS系数重建的x[k]:) 4-32 t=-2:0.001:2; N=input(N=); c0=0.5; xN=c0*ones(1,length(t)); for n=1:2:N xN=xN+cos(pi*n*t)*sinc(n/2); end plot(t,xN);title([N= +nu