资源描述：
《数字信号处理上机实验报告.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、数字信号处理上机实验报告实验一熟悉MATLAB环境一、实验目的1、熟悉MATLAB的主要操作命令。2、学会简单的矩阵输入和数据读写。3、掌握简单的绘图命令。4、用MATLAB编程并学会创建函数。5、观察离散系统的频率响应。二、实验内容认真阅读本章附录，在MATLAB环境下重新做一遍附录中的例子，体会各条命令的含义。在熟悉MATLAB基本命令的基础上，完成以下实验。上机实验内容：1、数组的加减乘除和乘方运算，输入，，求，，，，，并用stem语句画出A、B、C、D、E、F、G。程序：>>A=[1234];B=[3456];C=A+B
2、;D=A-B;E=A.*B;F=A./B;G=A.^B;subplot(2,4,1);stem(A,'.');subplot(2,4,2);stem(B,'.');subplot(2,4,3);stem(C,'.');subplot(2,4,4);stem(D,'.');subplot(2,4,5);stem(E,'.');subplot(2,4,6);stem(F,'.');subplot(2,4,7);stem(G,'.')161、用MATLAB实现下列序列。a)b)c)程序：A)clear;clc;n=[0:15];x1=
3、0.8.^n;subplot(3,1,1),stem(x1)title('x1=0.8^n')xlabel('n');ylabel('x1');16B)clear;clc;n=[0:15];x2=exp((0.2+3j)*n);subplot(3,1,1),stem(x2)title('x2=exp((0.2+3j)*n)')xlabel('n');ylabel('x2');C)clear;clc;n=[0:15];x3=3*cos(0.125*pi*n+0.2*pi)+2*sin(0.25*pi*n+0.1*pi);subpl
4、ot(3,1,1),stem(x3)title('x3=3*cos(0.125*pi*n+0.2*pi)+2*sin(0.25*pi*n+0.1*pi)')xlabel('n');ylabel('x3');161、绘出下列时间常数的图形，对轴，轴以及图形上方均须加上适当的标注：a)b)>>m=0:0.01:10;n=0:0.01:4;x1t=sin(2*pi*m);x2t=cos(100*pi*n).*sin(pi*n);subplot(2,1,1);plot(m,x1t);subplot(2,1,2);plot(n,x2t);
5、2、给定一因果系统H(z)=(1+)/()，求出并绘制H(z)的幅频响应与相频响应。>>b=[1,sqrt(2),1];a=[1,-0.67,0.9];[h,w]=freqz(b,a);am=20*log10(abs(h));%am=20*log10(abs(h))为幅频响应取dBsubplot(2,1,1);16plot(w,abs(h));xlabel('w');ylabel('幅频响应');title('系统响应')ph=angle(h);subplot(2,1,2);plot(w,ph);xlabel('w');ylab
6、el('相频响应');1、计算序列和序列的离散卷积，并作图表示卷积结果。>>a=[8-2-123];b=[23-1-3];c=conv(a,b);M=length(c)-1;n=0:1:M;stem(n,c);xlabel('n');ylabel('幅度');161、求以下差分方程所描述系统的单位脉冲响应：程序：>>N=50;a=[1-2];b=[10.1-0.06];x=[1zeros(1,N-1)];k=0:1:N-1;y=filter(a,b,x);stem(k,y);xlabel('n');ylabel('幅度');16
7、实验二快速傅里叶变换（FFT）及其应用一、实验目的1、在理论学习的基础上，通过本实验，加深对FFT的理解，熟悉MATLAB中的有关函数。2、应用FFT对典型信号进行频谱分析。3、了解应用FFT进行信号频谱分析过程中可能出现的问题，以便在实际中正确应用FFT。4、应用FFT实现序列的线性卷积和相关。二、实验内容实验中用到的信号序列a)高斯序列b)衰减正弦序列16a)三角波序列b)反三角波序列上机实验内容：（1）观察三角波和反三角波序列的时域和幅频特性，用点FFT分析信号序列和的幅频特性，观察两者的序列形状和频谱曲线有什么异同？绘出
8、两序列及其幅频特性曲线。程序：n=[0:3];k=[1:8];Xc(n+1)=n;Xc(n+5)=4-n;Xd(n+1)=4-n;Xd(n+5)=n;三角波特性subplot(2,2,1);plot(k-1,Xc);xlabel('n');ylabel('时域特
