云南大学信息学院信号系统实验二参考Word文件下载.docx
《云南大学信息学院信号系统实验二参考Word文件下载.docx》由会员分享,可在线阅读,更多相关《云南大学信息学院信号系统实验二参考Word文件下载.docx(16页珍藏版)》请在冰豆网上搜索。
(1)LTI连续时间系统的响应
(2)系统的冲激响应和阶跃响应
(3)卷积积分与卷积和
2、在编写Matlab程序时用到的函数及主要功能
(1)lism
连续系统在指定的任意时间范围内系统响应的时域波形及输出信号的时域波形。
(2)impulse
求连续系统冲激响应并绘制时域波形。
(3)step
求连续系统阶跃响应并绘制其时域波形图。
(4)filiter
对指定时间范围内的输入序列产生的响应的数值解。
(5)conv
该函数用于求两个有限时间区间非零的离散时间序列卷积
实验记录(程序、数据、波形等)
8.1
(1)y'
'
(t)+3y'
(t)+3y(t)=2f'
(t)+6f'
(t),f(t)=u(t)
MATLAB源程序如下:
a=[1,3,2];
b=[2,6];
sys=tf(b,a);
t=0:
0.2:
10;
f=exp(0*t);
lsim(sys,f,t)
(3)y'
(t)+2y'
(t)+2y(t)=f'
MATLAB源程序如下:
a=[1,2,2];
b=[1,0];
8.3
(2)y'
(t)+6y'
(t)+8y(t)=f'
(t)
a=[1,6,8];
y1=step(b,a,0:
1:
10)
y2=impulse(b,a,0:
subplot(2,1,1)
step(b,a,10)
subplot(2,1,2)
impulse(b,a,10)
(4)y'
(t)+4y(t)=2f(t)
a=[1,0,4];
b=[2];
y1=
0
0.7081
0.8268
0.0199
0.5728
0.9195
0.0781
0.4316
0.9788
0.1698
0.2960
y2=
0.9093
-0.7568
-0.2794
0.9894
-0.5440
-0.5366
0.9906
-0.2879
-0.7510
0.9129
8.4
(2)y(n)+1/2y(n-1)=x(n)+2x(n-1),x(n)=2cos(nπ/3)u(n)
a=[1,0.5];
b=[1,2];
n=0:
20;
x=2.*cos(n.*pi./3);
y=filter(b,a,x)
stem(n,x,'
filled'
)
stem(n,y,'
y=
Columns1through10
2.00004.0000-1.0000-3.5000-3.25000.62503.68753.1563-0.5781-3.7109
Columns11through20
-3.14450.57233.71393.1431-0.5715-3.7142-3.14290.57143.71433.1429
Column21
-0.5714
8.6
(2)y(n)-y(n-2)=x(n)
a=[1,0,-1];
b=[1];
n=1:
y=impz(b,a,20)
z=dstep(b,a,20)
stem(n,z,'
1
z=
2
3
4
5
6
7
8
9
有单位序列响应的时域波形可见,该系统不是稳定系统。
(4)y(n)-4y(n-1)+8y(n-2)=x(n)
a=[1,-4,8];
axis([0,20,-50000,50000])
-64
-256
-512
4096
16384
32768
-262144
-1048576
-2097152
16777216
67108864
134217728
13
-51
-307
-819
3277
19661
52429
-209715
-1258291
-3355443
13421773
80530637
8.7
x1=[0,1,2,1,0,0];
x2=[0,1,1,1,1,0,0];
n1=-2:
3;
n2=-1:
5;
x=conv(x1,x2);
ns=n1
(1)+n2
(1);
leg=length(x1)+length(x2)-2;
n=ns:
(ns+leg);
)
8.9
(1)
建立M文件
functionf=gggfconv(f1,f2,t1,t2)
f=conv(f1,f2);
f=f*0.01;
ts=t1
(1)+t2
(1);
leg=length(t1)+length(t2)-2;
t=ts:
0.01:
(ts+leg*0.01);
subplot(2,2,1)
plot(t1,f1)
axis([min(t1),max(t1),min(f1)-min(f1)*0.2,max(f1)+max(f1)*0.2])
title('
f1(t)'
xlabel('
t'
subplot(2,2,2)
plot(t2,f2)
axis([min(t2),max(t2),min(f2)-min(f2)*0.2,max(f2)+max(f2)*0.2])
f2(t)'
subplot(2,2,3)
plot(t,f);
axis([min(t),max(t),min(f)-min(f)*0.2,max(f)+max(f)*0.2])
p=get(gca,'
position'
);
p(3)=2.4*p(3);
set(gca,'
p);
f(t)=f1(t)*f2(t)'
t1=0:
2;
t2=1:
t4=-2:
f1=Heaviside(t1)-Heaviside(t1-2);
f2=Heaviside(t2-1)-Heaviside(t2-3);
f4=(0.5.*t4+1).*(t4<
0)+(-0.5.*t4+1).*(t4>
=0);
gggfconv(f1,f2,t1,t2)
gggfconv(f1,f4,t1,t4)
六、实验分析及实验小结
通过本次试验,我对MATLAB有了更深入的了解,试验中也遇到了些问题,有的程序编了又编才运行出来,实验过后了解了连续时间系统的时域分析。