武汉工程大学matlab实验三离散时间系统的时域分析Word文档下载推荐.docx
- 文档编号:20798862
- 上传时间:2023-01-25
- 格式:DOCX
- 页数:10
- 大小:18.60KB
武汉工程大学matlab实验三离散时间系统的时域分析Word文档下载推荐.docx
《武汉工程大学matlab实验三离散时间系统的时域分析Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《武汉工程大学matlab实验三离散时间系统的时域分析Word文档下载推荐.docx(10页珍藏版)》请在冰豆网上搜索。
m)
3.系统的因果性
系统的因果性即系统的可实现性。
如果系统时刻的输出取决于时刻及时刻以前的输入,而和时刻以后的输入无关,则该系统是可实现的,是因果系统。
系统具有因果性的充分必要条件为h(n)?
0,n?
4.系统的稳定性
稳定系统是指有界输入产生有界输出(BIBO)的系统。
如果对于输入序列,存在一个不变的正有限值,对于所有值满足|x(n)|?
M?
?
则称该输入序列是有界的。
稳定性要求对于每个有界输入存在一个不变的正有限值,对于所有值,输出序列满足|y(n)|?
K?
系统稳定的充分必要条件是系统的单位取样响应绝对可和,用公式表示n?
|h(n)|?
5.系统的冲激响应
设系统输入x(n)?
(n),系统输出y(n)的初始状态为零,这时系统输出用h(n)表示,即h(n)?
T[?
(n)],则称h(n)为系统的单位脉冲响应。
对于任意输入信号x(n),系统输出y(n)?
T[
?
x(m)?
(n?
m)]为利用系统满足m?
叠加原理得y(n)?
x(m)T[?
m)]
m?
m?
利用系统时不变性质得到y(n)?
x(m)h(n?
m)?
x(n)?
h(n)?
上式的运算关系称为卷积运算。
6.卷积的性质
1)交换律y(n)?
h(n)?
x(n)
2)结合律
y(n)?
h1(n)?
h2(n)?
[x(n)?
h1(n)]?
h2(n)
h2(n)]?
h1(n)
[h1(n)?
h2(n)]
3)分配律x(n)?
x(n)?
四、实验内容与步骤
1.离散时间系统的仿真
1y(n)?
1)M点因果滑动平滑系统的仿真,时域表达为MM?
1k?
0?
x(n?
k)
通过上述时域平滑系统可实现由若干个正弦信号之和所组成的信号中滤出高频分量。
据此,可以理解M点因果滑动平滑系统。
n=0:
100;
s1=sin(2*pi*0.34*n);
s2=sin(2*pi*0.02*n);
x=s1+s2;
M=input(&
#39;
desiredlengthofthefilter=&
);
num=ones(1,M);
y=filter(num,1,x)/M;
clf;
subplot(2,2,1);
stem(n,s1);
axis([0,100,-2,2]);
xlabel(&
n&
ylabel(&
s1(n)&
title(&
正弦序列1&
subplot(2,2,2);
stem(n,s2);
s2(n)&
正弦序列2&
subplot(2,2,3);
stem(n,x);
x(n)&
输入信号&
subplot(2,2,4);
stem(n,y);
y(n)&
输出信号&
desiredlengthofthefilter=6
2)线性与非线性离散时间系统的仿真。
简单的非线性系统实例:
y(n)?
1)x(n?
1)n=0:
x=cos(2*pi*0.1*n);
x1=[x00];
x2=[0x0];
x3=[00x];
y=x2.*x2-x1.*x3;
y=y(1:
101);
subplot(2,1,1);
)subplot(2,1,2);
)
2
简单的线性系统实例:
0.4y(n?
1)?
0.75y(n?
2)?
2.2403x(n)?
2.4908x(n?
2.2403x(n?
2)n=0:
50;
p=3;
q=-4;
x1=cos(2*pi*0.2*n);
x2=cos(2*pi*0.5*n);
x=p*x1+q*x2;
b=[2.24032.49082.2403];
a=[1-0.40.75];
y1=filter(b,a,x1);
y2=filter(b,a,x2);
y=filter(b,a,x);
yt=p*y1+q*y2;
d=y-yt;
subplot(3,1,1)
ylabel(&
y1(n)&
T[a1x1(n)+a2x2(n)]&
);
subplot(3,1,2);
stem(n,yt);
y2(n)&
;
a1T[x1(n)]+a2T[x2(n)]&
subplot(3,1,3);
stem(n,d);
y1(n)-y2(n)&
d&
);
3)时变与时不变系统的仿真。
时不变系统实例:
60;
m=8;
p=2;
q=-5;
x=p*cos(2*pi*0.1*n)+q*cos(2*pi*0.4*n);
xd=[zeros(1,m)x];
a=[1-0.40.75];
yd=filter(b,a,xd);
d=y-yd(1+m:
61+m);
subplot(3,1,1);
信号x(n)的响应y(n)&
grid;
stem(n,yd(1:
61));
y(n-m)&
信号x(n-m)的响应y(n-m)&
波幅&
y(n-m)-y(n-m)&
时变系统实例:
n=[-10:
10];
nx(n)?
1)
x1=sin(0.1*pi*n);
stem(n,x1);
x1&
x2=sin(0.1*pi*(n-1));
stem(n,x2);
x2&
y=n.*x1+x2;
stem(n,y)
y&
2.线性时不变系统仿真
1)冲激响应的计算
用MATLAB语言编程实现线性时不变系统的冲激响应计算。
线性时不变系统实例:
y
(n)?
2x(n)?
3x(n?
2x(n?
2)
m=50;
x=[1zeros(1,m-1)];
b=[232];
K=0:
1:
m-1;
stem(K,y);
title(&
冲激响应&
h(n)&
2)在实际应用中高阶因果线性时不变系统可以用低阶因果线性时不变系统级联得到,这可简化系统的设计与实现。
例如,对于三阶线性时不变系统
8y(n)?
10y(n?
6y(n?
3)?
16x(n?
40x(n?
3)可以用一个一阶和一个二阶系统级联实现。
第一级y1(n)?
0.25y1(n?
第二级y
2(n)?
y2(n?
0.5y2(n?
y1(n?
2y1(n?
用MATLAB语言编程验证系统的级联。
x=[1zeros(1,20)];
20;
a=[8-106-1];
b=[016-4016];
a1=[1-0.25];
b1=[2-1];
a2=[1-10.5];
b2=[01-2];
y1=filter(b1,a1,x);
y2=filter(b2,a2,y1);
d=y-y2;
信号通过高阶LTI系统的响应&
subplot(3,1,2);
stem(n,y2);
信号通过低阶LTI系统级联的响应&
subplot(3,1,3);
y(n)-y2(n)&
3.线性时不变系统的稳定性
若一个线性时不变系统的冲激响应是绝对可和,则此系统就是BIBO的稳定系统。
由此,无限冲激响应线性时不变系统稳定的必要条件是,随着输入序列点的增加,冲激响应衰减到零。
用MATLAB语言编程计算一个IIR线性时不变系统冲激响应的绝对值的和,验证稳定特性。
b=[0121];
a=[1-0.5-0.0050.3];
zplane(b,a);
%由分子分母多项式的系数画出零-极点分布图
Re&
jIm&
axis([-22-11]);
h=impz(b,a);
stem(h);
单位脉冲响应&
k&
[H,w]=freqz(b,a);
subplot(3,1,3)
plot(w/pi,abs(H));
频率\omega&
频率响应&
4.滤波概念实验
通过具体的时间系统理解信号滤波概念。
如:
系统1y(n)?
0.5x(n)?
0.27x(n?
0.77x(n?
系统2
0.53y(n?
0.46y(n?
0.45x(n)?
0.5x(n?
0.45x(n?
20?
n200?
)?
)0?
n?
99对于输入信号256256
实现各系统的滤波输出结果。
b1=[0.50.270.77];
a1=[1];
b2=[0.450.50.45];
a2=[1-0.530.46];
99;
xn=cos((20*pi*n)/256)+cos((200*pi*n)/256);
hn1=filter(b1,a1,xn);
hn2=filter(b2,a2,xn);
n1=0:
length(hn1)-1;
n2=0:
length(hn2)-1;
figure;
stem(n1,hn1);
subplot(2,1,2);
h1(n)&
系统1的冲激响应&
subplot(2,1,1);
subplot(2,1,2);
stem(n2,hn2);
h2(n)&
系统2的冲激响应&
五、实验扩展与思考
1.线性与非线性系统在信号输入/输出上有何不同?
时变与时不变系统又有何不同呢?
答:
由线性系统与非线性系统的仿真图可知,线性系统的输入信号与输出信号满足齐次性和可加性,而非线性系统的的输入与输出不满足齐次性和可加性。
2.冲激响应的计算实验中,就此系统计算它的阶跃响应,并与冲激响应比较,理解他们之间的关系。
Matlab产生阶跃响应的程序
x=[1ones(1,m-1)];
阶跃响应&
xlabel(&
冲激响应求和就是阶跃响应
3.系统级联实验中,三阶线性时不变系统若改用并联实现,又该如何进行?
3)11.2y1(n)?
0.25x(n?
4.8y2(n)?
5.6y2(n?
16x(n)?
y1(n)?
y2(n)
Matlab实现系统并联的程序
b1=[11.2];
b2=[4.8-5.6];
y2=filter(b2,a2,x);
yt=-16*x+y1+y2;
yt(n)&
信号通过低阶LTI系统并联的响应&
y(n)-yt(n)&
4.滤波概念实验中,两个系统的输出有何不同,为什么?
若改用最低频率为0、最高频率为0.5,长度为301的扫频正弦序列,则结果又将如何?
答:
系统1的冲激响应与输入信号仿真图一致,没有起到滤波的作用,序列值显得杂乱无章;
系统2的冲激响应与输入信号相比,系统2滤掉了部分序列值,使序列值变得平滑有序。
结果不变,系统输出特性与系统本身有关,与输入无关。
六、总结实验心得体会。
心得体会:
这次实验给我体会最深的就是系统级联与并联的转化,一开始对系统级联的概念不是特别理解,虽然通过老师给的分级系统完成了matlab程序的仿真,但是在进行系统并联的转化时,自己却无法直接得出并联的分级系统。
后来通过对课本的复习与深刻理解,完成了并联系统的转化,使自己受益匪浅。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 武汉 工程 大学 matlab 实验 离散 时间 系统 时域 分析