系统响应及系统稳定性.doc

实验一及课堂作业
实验一：系统响应及系统稳定性
一、 实验原理与方法
1、在时域求系统响应的方法有两种：第一种是通过解差分方程求得系统输出；第二种是已知系统的单位脉冲响应，通过求输入信号和系统单位脉冲响应的线性卷积求得系统输出。
2、检验系统的稳定性，其方法是在输入端加入单位阶跃序列，观察输出波形，如果波形稳定在一个常数值（包括零）上，系统稳定，否则不稳定。
3、系统的频域特性包括传输函数/特性（系统单位脉冲响应的傅里叶变换——幅频、相频）、系统函数/特性（系统单位脉冲响应的Z变换）、零极点分布等。分析系统的频域特性是为了知晓系统对不同频率的输入信号所产生的响应结果，因为零、极点分布对系统的频域特性有影响，通过控制系统函数的零、极点分布就可以设计出不同特性需求的系统。
二、实验内容
1、编制程序，包括产生输入信号、单位脉冲响应序列的子程序，用filter函数或conv函数求解系统输出响应的主程序。程序中要有绘制信号波形的功能。
2、给定一个低通滤波器的差分方程为

输入信号

（1）分别求出和)的系统响应，并画出其波形。
（2）求出系统的单位脉冲响应，画出其波形。程序见附录1.1、实验结果见图1.1。
3、给定系统的单位脉冲响应为 

用线性卷积法求)分别对系统和)的输出响应并画出波形。程序见附录1.2、实验结果见图1.2。
4、给定一谐振器的差分方程为 

令，谐振器的谐振频率为。 
（1）用实验方法检查系统是否稳定，输入信号为时，画出系统输出波形。
（2）给定输入信号为 

求出系统的输出响应，并画出其波形。程序见附录1.3、实验结果见图1.3。
三、 实验结果和分析、讨论及结论
1、实验结果：

图1.1 依据差分方程求取系统脉冲响应和输出响应

实验分析、讨论及结论：
（a） 中25个点数和程序所写一致。Filter函数实现线性常系数差分方程的递推求解，调用格式如下：Y=[filter(B,A,x)] ***计算系统对输入信号x的零状态响应输出信号向量Y，B、A是差分方程的系数向量。即B=[a1,a2……am] A=[b1,b2……bn]
2、实验结果：

图1.2 线性卷积求取二个不同系统输出响应

实验分析、讨论及结论：
（d）（f）单位脉冲响应点数与程序要求一致；（e）（g）卷积点数满足M+N-1的要求，图形也满足要求。
Conv函数用于计算两个有限长序列的卷积；C=conv（A,B）计算两个有限长序列向量A和B的卷积

3、 实验结果：

图1.3 依据差分方程求取系统脉冲响应和输出响应

实验分析、讨论及结论：
在系统的输入端加入单位阶跃序列，如果系统的输出趋近一个常数（包括零），就可以断定系统是稳定的，（h）中的输出显然趋近于零，所以系统是稳定的。
在（i）中，谐振器具有对某个频率进行谐振的性质，本实验中的谐振器的谐振频率是0.4 rad,因此稳定波形为sin(0.4n)。
四、 思考题
1、如果输入信号为无限长序列，系统的单位脉冲响应是有限长序列，可否用线性卷积法求系统的响应? 如何求？
答：如果输入信号为无限长序列，系统的单位脉冲响应是有限长序列，可用分段线性卷积法求系统的响应。具体步骤是对输入信号序列分段；‚求单位脉冲响应h（n）与各段的卷积；ƒ将各段卷积结果相加。具体实现方法有第三章P91介绍的重叠相加法和重叠保留法。
2、如果信号经过低通滤波器，把信号的高频分量滤掉，时域信号会有何变化，用前面第一个实验结果进行分析说明。
答：如果信号经过低通滤波器，则信号的高频分量将被过滤掉，时域信号的剧烈变化变的平滑，在有阶跃处附近产生过渡带。由实验内容1的结果图可见，经过系统低通滤波使得输入信号和的阶跃变化变的缓慢上升与下降。
五、 总结与心得体会
实验总结即在实验原理中说明的两点：
1、在时域求系统响应的方法有两种：第一种是通过解差分方程求得系统输出；第二种是已知系统的单位脉冲响应，通过求输入信号和系统单位脉冲响应的线性卷积求得系统输出。
2、检验系统的稳定性，其方法是在输入端加入单位阶跃序列，观察输出波形，如果波形稳定在一个常数值（包括零）上，系统稳定，否则不稳定。
实验的心得体会见下：
在此次试验中，通过课堂所留的三个例子，温习了关于MATLAB软件的操作及应用，基本使用方法和它的运行环境。又进一步地通过实验加深了对MATLAB软件的了解，体会到了MATLAB具有完备的图形处理功能，实现计算结果和编程的可视化等功能。通过做实验的过程以及实验分析的结果，了解并学会了filter函数和conv函数的基本用法，前者可计算知道输入信号的前提下求解输出响应的序列，后者则可以通过输入信号和系统的单位脉冲响应的线性卷积，求出系统的响应。
当然，在实验过程中，也遇到了一些问题，比如课堂作业的第三个问题里面，由于程序里面缺少“；”，导致少了一个结果图，通过检查并修改程序，解决了问题。总得来说，实验还是比较圆满的。
通过这次的实验。极大地提升了自己对于程序编辑的熟练度，增加了对于书本里面知识点的应用，更深一层的加深了对MATLAB软件的使用。这对自己以后的实验积累了丰富的经验。
六、附件：MATLAB原程序清单
1.1调用conv函数计算卷积
close all;clear all
A=[1,-0.9];B=[0.05,0.05]; %系统差分方程系数向量B和A
x1n=[1 1 1 1 1 1 1 1 zeros(1,40)]; %产生信号x1(n)=R8(n)，用zeros用来加点的个数
x2n=ones(1,60); %产生信号x2(n)=u(n)
hn=impz(B,A,40); %求系统单位脉冲响应h(n)
subplot(3,1,1);stem(hn); %调用函数stem绘图
title('(a) 系统单位脉冲响应h(n)');
y1n=filter(B,A,x1n); %求系统对x1(n)的响应y1(n)
subplot(3,1,2);stem(y1n);
title('(b) 系统对R8(n)的响应y1(n)');
y2n=filter(B,A,x2n); %求系统对x2(n)的响应y2(n)
subplot(3,1,3);stem(y2n);
title('(c) 系统对u(n)的响应y2(n)');

1.2利用线性卷积求取二个不同系统输出响应
close all;clear all
x1n=[ones(1,8)]; %产生信号x1(n)=R8(n)
h1n=[ones(1,10) zeros(1,10)];
h2n=[1 2.5 2.5 1 zeros(1,10)];
y21n=conv(h1n,x1n);
y22n=conv(h2n,x1n);
subplot(2,2,1);stem(h1n);
title('(d) 系统单位脉冲响应h1(n)');
subplot(2,2,3);stem(y21n);
title('(e) h1(n)与R8(n)的卷积y21(n)');
subplot(2,2,2);stem(h2n);
title('(f) 系统单位脉冲响应h2(n)');
subplot(2,2,4);stem(y22n);
title('(g) h2(n)与R8(n)的卷积y22(n)');

1.3 调用filter解差分方程求取系统单位脉冲响应和输出响应
close all;clear all
un=ones(1,300); %产生信号u(n)
n=0:299;
xsin=sin(0.014*n)+sin(0.4*n); %产生正弦信号
A=[1,-1.8237,0.9801];B=[1/100.49,0,-1/100.49]; %系统差分方程系数向量B和A
y1n=filter(B,A,un); %谐振器对u(n)的响应y31(n)
y2n=filter(B,A,xsin); %谐振器对u(n)的响应y31(n)
subplot(2,1,1);stem(y1n);
title('(h) 谐振器对u(n)的响应y31(n)');
subplot(2,1,2);stem(y2n);
title('(i) 谐振器对正弦信号的响应y32(n)');

课堂作业
1、绘出的频谱。
解：实验程序：
1.1
n=200;
stept=2*pi/n;
w=stept:stept:2*pi;
y=sin(2.5*w)./sin(0.5*w);
plot(w, y,w,zeros(size(w)));
axis([stept 2*pi -2 6]);
ylabel('y=sin(2.5*pi)/sin(0.5*pi)');
xlabel('w=0~2*pi');
grid on;

1.2
n=200;
stept=2*pi/n;
w=stept:stept:2*pi;
y=sin(2.5*w)./sin(0.5*w);
plot(w,abs(y),w,zeros(size(w)));
axis([stept 2*pi 0 6]);
ylabel('y=sin(2.5*pi)/sin(0.5*pi)');
xlabel('w=0~2*pi');
grid on;

实验结果：

图1.1 图1.2

2、 输入，单位脉冲响应，求输出序列，
解：实验程序：
N=5;
M=6;
L=N+M-1;
x=[1,2,3,4,5];
h=[6,2,3,6,4,2];
y=conv(x,h);
nx=0:N-1;
nh=0:M-1;
ny=0:L-1;
subplot(2 3 1);stem(nx,x,'.k');
xlabel('n');ylabel('x(n)');grid on;
subplot(2 3 2);
stem(nh,h,'.k');
xlabel('n');ylabel('h(n)');grid on;
subplot(2 3 3);stem(ny,y,'.k');
xlabel('n');ylabel('y(n)');grid on;

实验结果：
图2.1 图2.2 图2.3
3、分析频谱（a，b，c保证幅频特性的最大值为1）。
，，

解：实验程序：

p=0.8;r=0.85;alpha=pi/4;N=25;
b1=[1,1];a1=[1 -p];a=(1-p)/2;b1=b1*a;
b2=[1,-1];a2=[1 -p];b=(1+p)/2;b2=b2*b;
b3=[1 0 -1];
a3=[1 -2*r*cos(alpha) r*r];
c1=exp(j*2*alpha);
c=abs(1-c1)/((1-r)*abs(1-r*c1));
b3=b3/c;h1=impz(b1,a1,N);
subplot(331);stem(h1,'.k');subplot(332);zplane(b1,a1)
[H1,P]=freqz(b1,a1,256,'whole',1);
subplot(333);plot(P,abs(H1));grid on;
h2=impz(b2,a2,N);subplot(334);stem(h2,'.k')
hold on;plot(zeros(size(h2)));subplot(335);zplane(b2,a2)
[H2,P]=freqz(b2,a2,256,'whole',1);
subplot(336);plot(P,abs(H2));grid on;
h3=impz(b3,a3,N);subplot(337);stem(h3,'.k')
hold on;plot(zeros(size(h3)));subplot(338);zplane(b3,a3)
[H3,P]=freqz(b3,a3,256,'whole',1);
subplot(339);plot(P,abs(H3));grid on;

实验结果：
图3.1 图3.2 图3.3
图3.4 图3.5 图3.6
图3.7 图3.8 图3.9

其中专业理论知识内容包括：保安理论知识、消防业务知识、职业道德、法律常识、保安礼仪、救护知识。作技能训练内容包括：岗位操作指引、勤务技能、消防技能、军事技能。
二．培训的及要求培训目的
安全生产目标责任书
为了进一步落实安全生产责任制，做到“责、权、利”相结合，根据我公司2015年度安全生产目标的内容，现与财务部签订如下安全生产目标：
一、目标值：
1、全年人身死亡事故为零，重伤事故为零，轻伤人数为零。
2、现金安全保管，不发生盗窃事故。
3、每月足额提取安全生产费用，保障安全生产投入资金的到位。
4、安全培训合格率为100%。
二、本单位安全工作上必须做到以下内容：
1、对本单位的安全生产负直接领导责任，必须模范遵守公司的各项安全管理制度，不发布与公司安全管理制度相抵触的指令，严格履行本人的安全职责，确保安全责任制在本单位全面落实，并全力支持安全工作。
2、保证公司各项安全管理制度和管理办法在本单位内全面实施，并自觉接受公司安全部门的监督和管理。
3、在确保安全的前提下组织生产，始终把安全工作放在首位，当“安全与交货期、质量”发生矛盾时，坚持安全第一的原则。
4、参加生产碰头会时，首先汇报本单位的安全生产情况和安全问题落实情况；在安排本单位生产任务时，必须安排安全工作内容，并写入记录。
5、在公司及政府的安全检查中杜绝各类违章现象。
6、组织本部门积极参加安全检查，做到有检查、有整改，记录全。
7、以身作则，不违章指挥、不违章操作。对发现的各类违章现象负有查禁的责任，同时要予以查处。
8、虚心接受员工提出的问题，杜绝不接受或盲目指挥；
9、发生事故，应立即报告主管领导，按照“四不放过”的原则召开事故分析会，提出整改措施和对责任者的处理意见，并填写事故登记表，严禁隐瞒不报或降低对责任者的处罚标准。
10、必须按规定对单位员工进行培训和新员工上岗教育；
11、严格执行公司安全生产十六项禁令，保证本单位所有人员不违章作业。
三、 安全奖惩：
1、对于全年实现安全目标的按照公司生产现场管理规定和工作说明书进行考核奖励；对于未实现安全目标的按照公司规定进行处罚。
2、每月接受主管领导指派人员对安全生产责任状的落