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

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

实验一 熟悉MATLAB环境

一、实验目的

1、 熟悉MATLAB的主要操作命令。 2、 学会简单的矩阵输入和数据读写。 3、 掌握简单的绘图命令。

4、 用MATLAB编程并学会创建函数。 5、 观察离散系统的频率响应。

二、实验内容

认真阅读本章附录，在MATLAB环境下重新做一遍附录中的例子，体会各条命令的含义。在熟悉MATLAB基本命令的基础上，完成以下实验。

上机实验内容：

1、 数组的加减乘除和乘方运算，输入A??求C?A?B，1 2 3 4?，B??3 4 5 6?，

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) x(n)?e(0.2?3j)n 0?n?15

?n?0.2?)?2sin(0.25?n?0.1?) 0?n?15 c) x(n)?3cos(0.125程序： 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轴以及图形上方均须加上适当的标注：

0?t?10s a) x(t)?sin(2?t) 100?t)sin(?t) 0?t?4s b) x(t)?cos(>> 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);

4、 给定一因果系统H(z)=(1+2z?z)/(1-0.67z的幅频响应与相频响应。

>> 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))为幅频响应取dB subplot(2,1,1);

-1-2?1?z?2)，求出并绘制H(z)

4

plot(w,abs(h));xlabel('w');ylabel('幅频响应');title('系统响应') ph=angle(h); subplot(2,1,2);

plot(w,ph);xlabel('w');ylabel('相频响应');

5、 计算序列?8 -2 -1 2 3?和序列?2 3 -1 -3?的离散卷积，并作图表示卷积结果。

>> a=[8 -2 -1 2 3];

b=[2 3 -1 -3]; c=conv(a,b); M=length(c)-1; n=0:1:M; stem(n,c);

xlabel('n');ylabel('幅度');

5

本文来源：https://www.bwwdw.com/article/yl2v.html