资源描述：
《数字信号处理西电.docx》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、数字信号处理上机第一次实验实验一：设给定模拟信号，的单位是ms。(1)利用MATLAB绘制出其时域波形和频谱图（傅里叶变换），估计其等效带宽（忽略谱分量降低到峰值的3%以下的频谱）。(2)用两个不同的采样频率对给定的进行采样。。。比较两种采样率下的信号频谱，并解释。实验一MATLAB程序：（1）clc;fs=5000;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','
2、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)yy=abs(abs(sp)-3);[aa,freqind]=min(yy);(freqind-SPL/2)*fs/SPLclc;fs=1000;ts=1/fs;N=1000;t=(
3、-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=-fs/2:fs/SPL:fs/2-fs/SPL;plot(freqb,abs(sp))xlabel('频率')ylabel('频谱幅度')set(gca,'fontweight','b','fonts
4、ize',12)yy=abs(abs(sp)-3);[aa,freqind]=min(yy);(freqind-SPL/2)*fs/SPL实验三：设，，编写MATLAB程序，计算：（１）５点圆周卷积；（２）６点圆周卷积；（３）线性卷积；（４）画出的，和时间轴对齐。a=[1,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(
5、n)');axis([06015])subplot(312)stem(y2);gridontitle('六点圆周卷积y2(n)');xlabel('n'),ylabel('y2(n)');axis([06015])subplot(313)stem(y3);gridontitle('线性卷积y3(n)');xlabel('n'),ylabel('y3(n)');axis([06015])实验四：给定因果系统：（１）求系统函数并画出零极点示意图。（２）画出系统的幅频特性和相频特性。（３）求脉冲响应并画序列图。提示：在ＭＡＴＬＡＢ中，zplane(b,a
6、)函数可画零极点图；Freqz(b,a,N)可给出范围内均匀间隔的点频率响应的复振幅；Impz(b,a,N)可求的逆变换（即脉冲响应）。clca=[1,0]b=[1,-0.9]figure(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
7、,phi);gridon;title('相频特性');xlabel('f/Hz'),ylabel('W(f)');数字信号处理第二次实验1.给定模拟信号，对其进行采样，用DFT（FFT）进行信号频谱分析。（１）确定最小采样频率和最小采样点数。（２）若以秒进行采样，至少需要取多少采样点？（３）用DFT的点数画出信号的点DFT的幅度谱，讨论幅度谱结果。（４）分别为和，能否分辨出信号的所有频率分量。（５）在（３）和（４）的条件下做补0FFT，分析结果。（６）在不满足最小采样点数的情况下做补0DFT，观察是否可以分辨出两个频率分量。程序如下：clearcl
8、oseallclc%(1)确定最小采样频率和最小采样点数w1=4*pi;w2=8*pi;f1=w1/(2*pi);f2=w
