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

您所在位置:网站首页 > 海量文档

 > 高等教育 > 实验设计

7e3e25dd1884ce9b7725a69a698d5758.gif

Matlab数字信号处理实验报告.doc26页

本文档一共被下载:2377011b61454b208b34ad77b4c31088.gif次,您可全文免费在线阅读后下载本文档。

072ad097814a33d1ad18e33ca14a5a9e.png

9ec2792a743ebdf4575ee09a03c275ea.png

7af312d1be2cfc1e59f15880bb7a7e7a.png

5f207613e329ca76dc87c76d3ff77e64.png

adf7aa3071c520b8a2063089f549261c.png

下载提示

1.本站不保证该用户上传的文档完整性,不预览、不比对内容而直接下载产生的反悔问题本站不予受理。

2.该文档所得收入(下载+内容+预览三)归上传者、原创者。

3.登录后可充值,立即自动返金币,充值渠道很便利

数字信号处理实验报告

基础实验篇

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

实验原理

利用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定义为

反变换为

??有限长序列的DFT是其Z变换在单位圆上的等距采样,或者说是序列Fourier变换的等距采样,因此可以用于序列的谱分析。

??FFT是为了减少DFT运算次数的一种快速算法。它是对变换式进行一次次分解,使其成为若干较短序列的组合,从而减少运算量。常用的FFT是以2为基数的,其长度 。它的效率高,程序简单,使用非常方便,当要变换的序列长度不等于2的整数次方时,为了使用以2为基数的FFT,可以用末位补零的方法,使其长度延长至2的整数次方。

??用FFT可以实现两个序列的圆周卷积。在一定的条件下,可以使圆周卷积等于线性卷积。一般情况,设两个序列的长度分别为N1和N2,要使圆周卷积等于线性卷积的充要条件是FFT的长度

N≥N1+N2

对于长度不足N的两个序列,分别将他们补零延长到N。

二、实验目的

加深理解离散傅立叶变换及快速傅立叶变换概念;

学会应用FFT对典型信号进行频谱分析的方法;

研究如何利用FFT程序分析确定性时间连续信号;

熟悉应用FFT实现两个序列的

发表评论

请自觉遵守互联网相关的政策法规,严禁发布色情、暴力、反动的言论。

用户名:

验证码:

c9f2bced460b0329ba0aadbbc3f0fc71.png

匿名?

发表评论

你可能感兴趣的:(matlab,及数字信号实验报告)