matlab数字信号处理实验报告,Matlab数字信号处理实验报告.pdf

数字信号处理实验报告

基础实验篇

实验一 离散时间系统及离散卷积

一、实验原理

利用Matlab软件计算出系统函数的零极点分布、单位脉冲响应和系统

频率响应等的图像并于笔算结果进行比较，找出异同。编译合适程序能计

算取值范围不同的离散卷积。

二、实验目的

(1)熟悉MATLAB软件的使用方法。

(2)熟悉系统函数的零极点分布、单位脉冲响应和系统频率响应等概念。

(3)利用MATLAB绘制系统函数的零极点分布图、系统频率响应和单位脉冲

响应。

三、实验步骤

(1)自编并调试实验程序，并且，给实验程序加注释；

(2)按照实验内容完成笔算结果；

(3)验证计算程序的正确性，记录实验结果。

(4)至少要求一个除参考实例以外的实验结果，在实验报告中，要描述清楚

实验结果对应的系统，并对实验结果进行解释说明。

四、实验源程序及实验结果

实验1-1 运行结果

xlabel ('n');

ylabel ('h (n)');

figure (2)

[z,p,g] tf2zp (b,a);

zplane (z,p)

title ('零极点');

function

[x,n] chongji (n1,n2)

n [n1:n2];

x [n 0];

function shiyan1()

a [1,-1,0.9];

b 1;

x chongji (-20,120);

n -20:120;

h filter (b,a,x);

figure (1)

stem (n,h);

title ('冲击响应');

实验1-2 运行结果

b [0.0181,0.0543,

0.0543,0.0181];

a [1.000,-1.76,

1.1829,-0.2781];

w pi*freqspace (500);

H freqz(b,a,w);

MH abs(H);

AH angle (H);

subplot (2,1,1);

plot (w/pi,MH);

grid;

axis([0,1,0,1]);

xlabel ('w (pi)');

ylabel ('|H |');

title ('幅度、相位响应');

subplot (2,1,2);

plot (w/pi,AH);

grid;

xlabel ('w (pi)');

ylabel ('angle (H)');

实验1-3 运行结果

n 0:30;

%输入x (n)和冲激响应h (n)

x zeros(1,length (n));

h zeros(1,length (n));

x ([find ((n> 0)&(n< 4))]) 1;

h ([find ((n> 0)&(n< 8))]) 0.

5;

figure (1)

subplot (3,1,1);

stem (n,x);

axis([0,30,0,2]);

title ('输入序列');

xlabel ('n');

ylabel ('x (n)');

subplot (3,1,2);

stem (n,h);

axis([0,30,0,2]);

title ('冲激响应序列');

xlabel ('n');

ylabel ('h (n)');

%输出响应

y conv (x,h);

subplot (3,1,3);

n 0:length (y)-1;

stem (n,y);

title ('输出响应');

xlabel ('n');

ylabel ('y (n)');

实验二 离散傅立叶变换与快速傅立叶变换

一、 实验原理

对有限长序列使用离散Fouier变换(DFT)可以很好的反映序列的频谱特性，

而且易于用快速算法在计算机上实现，当序列x(n)的长度为N 时，它的DFT定

义为

N1