- 1、本文档共5页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
数字信号处理——实验二
武汉工程大学 电气信息学院 通信工程
红烧大白兔
实验目的
在时域中仿真离散时间系统,进而理解离散时间系统对输入信号或延时信号进行简单运算处理,生成具有所需特性的输出信号的方法。
仿真并理解线性与非线性、时变与时不变等离散时间系统。
掌握线性时不变系统的冲激响应的计算并用计算机仿真实现。
仿真并理解线性时不变系统的级联、验证线性时不变系统的稳定特性。
实验设备
计算机,MATLAB语言环境
实验基础理论
系统的线性
线性性质表现为系统满足线性叠加原理:若某一输入是由N个信号的加权和组成的,则输出就是由系统对这N个信号中每一个的响应的相应加权和组成的。设x1(n)和 x2(n)分别作为系统的输入序列,其输出分别用y1(n)和y2(n)表示,即
Y1(n)=T[x1(n)], y2(n)=T[x2(n)] 若满足T[a1x1(n)+a2x2(n)]=a1y1(n)+a2y2(n)
则该系统服从线性叠加原理,或者称为该系统为线性系统。
系统的时不变特性
若系统的变换关系不随时间变化而变化,或者说系统的输出随输入的移位而相应移位但形状不变,则称该系统为时不变系统。对于时不变系统,
若y(n)=T[x(n)] 则T[x(n-m)]=y(n-m)
系统的因果性
系统的因果性既系统的可实现性。如果系统n时刻的输出取决于n时刻及n时刻以前的输入,而和以后的输入无关,则该系统是可实现的,是因果系统。系统具有因果性的充分必要条件是h(n)=0,n0
系统的稳定性
稳定系统是指有界输入产生有界输出(BIBO)的系统。如果对于输入序列x(n),存在一个不变的正有限值M,对于所有n值满足|x(n)|≤M∞则称该输入序列是有界的。稳定性要求对于每个有界输入存在一个不变的正有限值K,对于所有n值,输出序列y(n)满足 |y(n)|≤K∞ 系统稳定的充分必要条件是系统的单位取样响应绝对可和,用公式表示为
系统的冲激响应
设系统输入x(n)=δ(n),系统输出y(n)的初始状态为零,这时系统输出用h(n)表式,即h(n)=T[δ(n)]则称h(n)为系统的单位脉冲响应。也就是说,单位脉冲响应h(n)是系统对δ(n)的零状态响应,它表征了系统的时域特性。两个序列的卷积是一个序列与另一个序列翻转后逐次移位乘积之和,故称为离散卷积,也称为两序列的线性卷积。
6、卷积的性质
交换律、结合律、分配率
实验内容与步骤
离散时间系统的仿真
M点因果滑动平滑系统的仿真
n=0:100;
s1=cos(2*pi*0.05*n);
s2=cos(2*pi*0.47*n);
x=s1+s2;
M=input(desired length of the filter =);
desired length of the filter =6
num=ones(1,M);
y=filter(num,1,x)/M;
clf;
subplot(2,2,1);
plot(n,s1);
axis([0,100,-2,2]);
xlabel(Time index n);ylabel(Amolitude);
title(Signal #1);
subplot(2,2,2);
plot(n,s2);
axis([0,100,-2,2]);
xlabel(Time index n);ylabel(Amolitude);
title(Signal #2);
subplot(2,2,3);
plot(n,x);
axis([0,100,-2,2]);
xlabel(Time index n);ylabel(Amolitude);
title(Input Signal);
subplot(2,2,4);
plot(n,y);
axis([0,100,-2,2]);
xlabel(Time index n);ylabel(Amolitude);
title(Ouput Signal);
线性与非线性离散时间系统的仿真
n=0:200;
x=cos(2*pi*0.05*n);
x1=[x 0 0];
x2=[0 x 0];
x3=[0 0 x];
y=x2.*x2-x1.*x3;
y=y(2:202);
subplot(2,1,1);
plot(n,x)
xlabel(Time index n);ylabel(Amolitude);
title(Input Signal)
subplot(2,1,2)
plot(n,y)
xlabel(Time index n);ylabel(Amplitude);
title(Output Signal)
线性离散时间系统的仿真
n=0:40;
a=2;b=-3;
x1=cos(2*pi*0.1*n);
文档评论(0)