太原理工大学数字信号分析与处理.docVIP

数字信号分析与处理 1、通过MATLAB编程找到下述信号的频率谱。（自行给出，每位同学不同；附MATLAB程序以及结果图）。(试题等级：容易) matlab 程序如下 %%%%%原始信号 %%%%%%%% x=@(t) 0.25+2.*sin(2.*pi.*1000.*t)+sin(2.*pi.*3000.*t)+1.5.*sin(2.*pi.*4000.*t) +0.5.*sin(2.*pi.*2000.*t) F=16000； T=1/F %%% 采样频率为 16000hz 采样周期为 8.3333e-05,这里取采样点32个%%%% a=[] syms p； p=2； while p33 t=(p-1).*T，a(p)=x(t)， p=p+1 end a(1)=x(0)； %%%%% 一维矩阵a即为采样后的采样点x(0)x(1)x(2)x(3)......x(15)% A=[] %%%%定义DFT变换后元素X(n)所称的矩阵% k=1 while k33 A(k)=0,p=1， while p33 A(k)=A(k)+a(p).*exp(-1.*i.*2.*pi.*k.*(p-1)/32) p=p+1 end k=k+1 end %%%% 此时经过DFT变换的采样点已全部赋值给矩阵A[]%%%%%% %%%%%作图%%%%% k=[0:1:31]X=A(k+1)； XR=real(X)； XI=imag(X)； XAMP=abs(X) XPHA=angle(X)；stem(k,XAMP)；xlabel(n),ylabel(X的幅值) figure stem(k,XPHA)， xlabel(n),ylabel(X的相位) 幅频图 相频图  2.已知输入信号为, 若采样点个数有限, 要求绘出① x(n)的DFT图，此时; ②若①中不能很好地区分频谱，采用什么办法来提高谱分辨率？编写MATLAB程序并实现。（试题等级：中等）  （1）Matlab程序 当数据点数为8时： n=[1:1:8]; x=cos(0.49.*pi.*n)+2.*cos(0.51.*pi.*n); X=fft(x,8); stem(n,X); XAMP=abs(X); XPHA=angle(X); stem(n,XAMP); xlabel(n),ylabel(X的幅值); figure; stem(n,XPHA); xlabel(n),ylabel(X的相位); 幅频图 相频图 （2）上面的DFT图不能完全反应频谱特性，可增加数据点数，以增加采样分辨率，降数据点数改为128. 可得到较为清晰的频谱图，matlab程序如下： n=[1:1:128]; x=cos(0.49.*pi.*n)+2.*cos(0.51.*pi.*n); X=fft(x,128); stem(n,X); XAMP=abs(X); XPHA=angle(X); stem(n,XAMP); xlabel(n),ylabel(X的幅值); figure; stem(n,XPHA); xlabel(n),ylabel(X的相位) 幅频图 相频图 3、已知，。① 编写MATLAB程序实现两个输入信号的线性卷积；② 编写MATLAB程序实现两个输入信号的循环卷积，要求循环卷积的结果等于①。（要求附上MATLAB程序和结果图）(试题等级：容易)  (1)、计算线性卷积;X= Matlab 程序如下; xn=0:3; x=[1 2 2 1]; yn=0:3; y=[-1 2 2 -1]; zn=0:(3+3); z=conv(x,y); stem(zn,z);grid;axis([-1 7 -2 7]); xlabel(n);ylabel(z(n)); title(线性卷积g(n)); （2）、计算循环卷，matlab程序如下： %%先计算线性卷积%%% xn=0:3; x=[1 2 2 1]; yn=0:3; y=[-1 2 2 -1]; zn=0:(3+3); z=conv(x,y); %循环卷积%%%%% N=length(zn); %为了得到与线性卷积相同的结果，取周期为线性卷积序列的长度 k=1;%取周期延拓后的主序列 nt=min(zn):k*N+min(zn)-1; t=zeros(1,k*N); for i=0:k*N-1 t(i+1)=z(mod(i,N)+1); end ,stem(nt,t,.);grid;axis([-1 7 -2 7]); xlabel(n);ylabel(t(