本科生实验报告
实验课程 《MATLAB仿真技术》
学院名称 核技术与自动化工程学院 专业名称 核工程与核技术 学生姓名 杨洲 学生学号201506020219
指导教师 杨剑波
实验地点 6C802 实验成绩
二〇 一六 年 十 月 二〇一六 年 十一 月
1.对于典型二阶系统
?G?s??22s?2??ns??n 编写程序并绘制出无阻尼自然振荡?n=8，阻尼比?分别为0.2,0.4,…,1.0,2.0时系统的单位阶跃响应曲线 程序:
%典型二阶系统的单位阶跃响应曲线 clc,clear;
Wn=8; %无阻尼自然振荡 t=linspace(0,5);
for a=0.2:0.2:0.8 %（0<阻尼比a<1,欠阻尼） Wd=Wn.*sqrt(1-a.^2);
c=1-exp(-a.*Wn.*t).*(cos(Wd.*t)+a/sqrt(1-a.^2).*sin(Wd.*t)); plot(t,c); hold on end
hold on
a=1.0; %（阻尼比a=1,临界阻尼） c=1-(1+Wn.*t).*exp(-Wn.*t); plot(t,c); hold on
a=2.0; %(阻尼比a=2.0>1,过阻尼） s1=-a.*Wn+sqrt(a.^2-1).*Wn; s2=-a.*Wn-sqrt(a.^2-1).*Wn;
c=1-1/(2.*sqrt(a.^2-1)).*(exp(s1.*t)/(a-sqrt(a.^2-1))-exp(s2.*t)/(a+sqrt(a.^2-1))); plot(t,c);
xlabel('t'),ylabel('c(t)'); title('单位阶跃响应曲线')
2n绘制图形如下：
2.已知某负反馈系统的开环传递函数为
G(s)?0.05s?0.045(s2?1.8s?0.9)(s2?5s?6)
编写程序并绘制系统的根轨迹。 程序：
%负反馈系统的根轨迹 clc,clear;
%G=(0.05*s+0.045)/((s.^2-1.8*s+0.9)*(s.^2+5*s+6) num=[0.05 0.045];
den=[1 3.2 -2.1 -6.3 5.4]; rlocus(num,den); title('Root locus'）
图形如下：
3.如下40个gamma数据，绘制成2D图所示。
（1）编写程序绘制出如图所示的2D图；
（2）编写程序求出峰位和半高宽，并显示。
程序：
%gamma数据处理 clc,clear;
x=[1:40];
y=[0,0,14,45,81,124,172,225,283,342,401,457,506,... 546,574,587,585,566,533,486, ...
430,368,305,246,196,158,135,127,135,155,... 184,217,248,272,285,283,266,236,193,142]; plot(x,y); %gamma数据转化为2D图 [a,b]=max(y);
fprintf('峰位是: x=%d\\n',b);%求峰位 m1=b;
for x=b:30
if abs(a/2-y(x))0 m1=x; end end
m2=m1-1; n1=1;
for x=1:b
if abs(a/2-y(x))0 n1=x; end end
n2=n1+1;
fprintf('半高宽是： %d 图形如下：