[信号与系统利用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= +num2str(N)]) axis([-2 2 -0.2 1.2]); （N=31） （N=7） 4-33 function y=sf1(t,w); y=(t=-1t=1).*(1-abs(t)).*exp(-j*w*t); w=linspace(-6*pi,6*pi,512); N=length(w);X=zeros(1,N); for k=1:N X(k)=quadl(sf1,-1,1,[],[],w(k)); end figure(1); plot(w,real(X)); xlabel(\ omega); ylabel(X(j\ omega)); figure(2); plot(w,real(X)-sinc(w/2/pi).^2); xlabel(\ omega); title(计算误差); （计算误差） （三角波信号的近似频谱） 例4-34 function y=sf2(t) y=2*sinc(t).*sinc(t); fm=linspace(0,5,256); N=length(fm);E=zeros(1,N); for k=1:N E(k)=quadl(sf2,0,fm(k)); end plot(fm,E); xlabel(Hz); ylabel(E);（矩形信号在频域的能量E随fm变化的曲线） 例5-15 w=linspace(0,2*pi,200); b=[1]; a=[1 2 2 1]; H=freqs(b,a,w); subplot(2,1,1); plot(w,abs(H)); set(gca,xtick,[0 1 2 3 4 5]); set(gca,ytick,[0 0.4 0.707 1]);grid; xlabel(\omega); ylabel(|H(j\omega)|) subplot(2,1,2); plot(w,angle(H)); set(gca,xtick,[0 1 2 3 4 5]);grid; xlabel(\omega); ylabel(\phi(\omega));  例5-16 RC=0.04; t=linspace(-2,2,1024); w1=5;w2=100; H1=j*w1/(j*w1+1/RC); H2=j*w2/(j*w2+1/RC); x=cos(5*t)+cos(100*t); y=abs(H1)*cos(w1*t+angle(H1))+abs(H2)*cos(w2*t+angle(H2)); subplot(2,1,1); plot(t,x); ylabel(x(t)); xlabel(时间（秒）); subplot(2,1,2); plot(t,y); ylabel(y(t)); xlabel(时间（秒）);  例5-17 b=[1,1];a=[1,-1,0.5]; [H,w]=freqz(b,a); plot(w,abs(H)) xlabel(Frequency(rad)); ylabel(Magnitude); title(Magnitude respnse); 例5-18 Fm=10;Fc=100;Fs=1000; N=1000;k=0:N-1; t=k/Fs; x=sin(2