资源描述：
《数字信号处理西电》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、-数字信号处理上机第一次实验实验一：设给定模拟信号xate1000t，t的单位是ms。(1)利用MATLAB绘制出其时域波形和频谱图（傅里叶变换），估计其等效带宽（忽略谱分量降低到峰值的3%以下的频谱）。(2)用两个不同的采样频率对给定的xat进行采样。15000样本秒采样xat得到x1n。○以fs画出x1n及其频谱X1ej。21000样本秒采样xat得到x1n。○以fs画出x1n及其频谱X1ej。比较两种采样率下的信号频谱，并解释。实验一MATLAB程序：（1）○1clc;fs=5000;ts=1/fs;N=1000;t=(-N:N)*ts;s
2、=exp(-abs(t));plot(t,s,'linewidth',1.5)xlabel('时间')ylabel('幅度')set(gca,'fontweight','b','fontsize',12)SPL=N*100;figuresp=fftshift(fft(s,SPL));sp=sp/max(sp)*100;freqb=-fs/2:fs/SPL:fs/2-fs/SPL;plot(freqb,abs(sp))xlabel('频率')ylabel('频谱幅度')--set(gca,'fontweight','b','fontsize',12
3、)yy=abs(abs(sp)-3);[aa,freqind]=min(yy);(freqind-SPL/2)*fs/SPL--○2--clc;--fs=1000;ts=1/fs;N=1000;t=(-N:N)*ts;s=exp(-abs(t));plot(t,s,'linewidth',1.5)xlabel('时间')ylabel('幅度')set(gca,'fontweight','b','fontsize',12)SPL=N*100;figuresp=fftshift(fft(s,SPL));sp=sp/max(sp)*100;freqb=
4、-fs/2:fs/SPL:fs/2-fs/SPL;plot(freqb,abs(sp))xlabel('频率')ylabel('频谱幅度')set(gca,'fontweight','b','fontsize',12)yy=abs(abs(sp)-3);[aa,freqind]=min(yy);(freqind-SPL/2)*fs/SPL--实验三：设xn1,2,2，x2n1,2,3,4，编写MATLAB程序，计算：1（１）５点圆周卷积y1n；（２）６点圆周卷积y2n；（３）线性卷积y3n；（４）画出的y1n，y2n和y3n时间轴对齐。a=[1,
5、2,2];b=[1,2,3,4];y1=cconv(a,b,5)y2=cconv(a,b,6)y3=conv(a,b)figure(1);subplot(311)stem(y1);gridontitle('五点圆周卷积y1(n)');--xlabel('n'),ylabel('y1(n)');axis([06015])subplot(312)--stem(y2);gridontitle('六点圆周卷积y2(n)');xlabel('n'),ylabel('y2(n)');axis([06015])subplot(313)stem(y3);grid
6、ontitle('线性卷积y3(n)');xlabel('n'),ylabel('y3(n)');axis([06015])--实验四：给定因果系统：yn0.9yn1xn（１）求系统函数Hz并画出零极点示意图。（２）画出系统的幅频特性Hej和相频特性。（３）求脉冲响应hn并画序列图。--提示：在ＭＡＴＬＡＢ中，zplane(b,a)函数可画零极点图；Freqz(b,a,N)可给出0,范围内--均匀间隔的N点频率响应的复振幅；Impz(b,a,N)可求Hz的逆变换（即脉冲响应）。--clca=[1,0]b=[1,-0.9]figu
7、re(1)zplane(b,a);title('零极点分布图')w=[-3*pi:0.01:3*pi];[h,phi]=freqz(b,a,w);figure(2);subplot(2,1,1);plot(w,abs(h));gridon;title('幅频特性');xlabel('f/Hz'),ylabel('H(f)');subplot(2,1,2);plot(w,phi);gridon;title('相频特性');xlabel('f/Hz'),ylabel('W(f)');----数字信号处理第二次实验1.给定模拟信号xt2sin4t5co
8、s8t，对其进行采样，用DFT（FFT）进行信号频谱分析。（１）确定最小采样频率和最小采样点数。（２）若以t0.01nn0:N1秒进行采
