数字信号处理课后习题Matlab作业..docxVIP

第 PAGE \* MERGEFORMAT 数字信号处理MATLAB习题 数字信号处理MATLAB习题 M1-1 已知，，，以抽样频率对上述三个信号进行抽样。在同一张图上画出，和及抽样点，对所得结果进行讨论。 解： 从以上两幅图中均可看出，三个余弦函数的周期虽然不同，但它们抽样后相应抽样点所对应的值都相同。那么这样还原回原先的函数就变成相同的，实际上是不一样的。这是抽样频率太小的原因，我们应该增大抽样频率才能真实还原。如下图：f=50Hz 程序代码 f=10; t=-0.2:0.001:0.2; g1=cos(6.*pi.*t); g2=cos(14.*pi.*t); g3=cos(26.*pi.*t); k=-0.2:1/f:0.2; h1=cos(6.*pi.*k); h2=cos(14.*pi.*k); h3=cos(26.*pi.*k); % subplot(3,1,1); % plot(k,h1,r.,t,g1,r); % xlabel(t); % ylabel(g1(t)); % subplot(3,1,2); % plot(k,h2,g.,t,g2,g); % xlabel(t); % ylabel(g2(t)); % subplot(3,1,3); % plot(k,h3,b.,t,g3,b); % xlabel(t); % ylabel(g3(t)); plot(t,g1,r,t,g2,g,t,g3,b,k,h1,r.,k,h2,g.,k,h3,b.) xlabel(t); ylabel(g(t)); legend(g1(t),g2(t),g3(t)); M2-1 利用DFT的性质，编写一MATLAB程序，计算下列序列的循环卷积。 (1) g[k]={1,-3,4,2,0,-2,},h[k]={3,0,1,-1,2,1}； (2) x[k]=cos(k/2),y[k]=3k,k=0,1,2,3,4,5。 解： (1) 循环卷积结果 6.0000 -3.0000 17.0000 -2.0000 7.0000 -13.0000 程序代码 g=[1 -3 4 2 0 -2]; h=[3 0 1 -1 2 1]; l=length(g); L=2*l-1; GE=fft(g,L); HE=fft(h,L); y1=ifft(GE.*HE); for n=1:l if n+l=L y2(n)=y1(n)+y1(n+l); else y2(n)=y1(n); end end y2 stem(0:l-1,y2) xlabel(k) ylabel(y(k)) title(循环卷积) (2) 循环卷积结果 -71.0000 -213.0000 89.0000 267.0000 73.0000 219.0000 程序代码 k=0:5; x=cos(pi.*k./2); y=3.^k; l=length(x); L=2*l-1; GE=fft(x,L); HE=fft(y,L); y1=ifft(GE.*HE); for n=1:l if n+l=L y2(n)=y1(n)+y1(n+l); else y2(n)=y1(n); end end y2 stem(0:l-1,y2) xlabel(k) ylabel(y’(k)) title(循环卷积) M2-2 已知序列 (1)计算序列DTFT的表达式,并画出N=10时，的曲线。 (2)编写一MATLAB程序，利用fft函数，计算N=10时，序列x[k]的DTFT在的抽样值。利用hold函数，将抽样点画在的曲线上。 解： (1) 程序代码 N=10; k=-N:N; x=cos(k.*pi./(2*N)); W=linspace(-pi,pi,512); X=zeros(1,length(W)); for k=-N:N X1=x(k+N+1).*exp(-j.*W.*k); X=X+X1; end plot(W,abs(X)) xlabel(W); ylabel(abs(X)); (2) 程序代码 N=10; k=-N:N; x=cos(k.*pi./(2*N)); X_21=fft(x,21); L=-10:10; W=linspace(-pi,pi,1024); X=zeros(1,length(W)); for k=-N:N X1=x(k+N+1).*exp(-j.*W.*k); X=X+X1; end plot(W,abs(X));hold on; plot(2*pi*L/21,fftshift(a