数字信号处理上机实验报告.pdf

数字信号处理上机实验报告 实验一 熟悉 MATLAB环境 一、实验目的 1、 熟悉 MATLAB的主要操作命令。 2、 学会简单的矩阵输入和数据读写。 3、 掌握简单的绘图命令。 4、 用 MATLAB编程并学会创建函数。 5、 观察离散系统的频率响应。 二、实验容 认真阅读本章附录， 在 MATLAB环境下重新做一遍附录中的例子， 体会各条命令的含义。 在熟悉 MATLAB基本命令的基础上，完成以下实验。 上机实验容： 1、 数组的加减乘除和乘方运算， 输入 A 1 2 3 4 ，B 3 4 5 6 ，求 C A B ， D A B ， E A. B ， F A./ B ， G A.^ B ，并用 stem 语句画出 A、B、C、D、E、 F、G。 程序： >> A=[1 2 3 4];B=[3 4 5 6]; C=A+B; 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,'.') 2、 用 MATLAB实现下列序列。 a) x(n) 0.8n 0 n 15 b) ( ) (0.2 3j ) n 0 15 x n e n c) x(n) 3cos(0.125 n 0.2 ) 2 sin( 0.25 n 0.1 ) 0 n 15 程序： A) clear;clc; n=[0:15]; x1=0.8.^n; subplot(3,1,1),stem(x1) title('x1=0.8^n') xlabel('n'); ylabel('x1'); 2 B) 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); subplot(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'); 3 3、 绘出下列时间常数的图形，对 x 轴， y 轴以及图形上方均须加上适当的标注： a) x(t ) sin( 2 t ) 0 t 10s b) x(t ) cos(100 t )sin( t ) 0 t 4s >> m=0:0.01:10; n=0:0.01:4; x1t=sin(2*pi*m); x2t=cos(100*pi*n).*sin(pi*n); subplot