机电控制系统仿真-实验报告

对应于问题④的程序如下：
p=[1 -12 0 25 116];a=roots(p) r=[11.7473, 2.7028, -1.2251 + 1.4672i, -1.2251 - 1.4672i]; p=poly(r) a=[1 2 3] ; b=[1 2];c=conv(a,b)
对应于问题⑤的程序如下：
length(A)=max(size(A))
对应于问题③的程序如下：
t=0:0.01:2*pi; y=sin(t); y1=sin(t+0.5*pi); y2=cos(t); y3=cos(t+0.5*pi); figure(1) subplot(1,2,1) plot(t,y) xlabel('time/s') ylabel('amplitude/1') axis([0,6.28,-1.5,1.5]) subplot(1,2,2) plot(t,y1)
A=zeros(9); for i=1:1:9
for j=1:1:9 a(i,j)=i*j; if i<j a(i,j)=0 end
end end
1.3 运行结果 对应于问题①，运行程序后得到结果如下：
对应于问题②，运行程序后得到结果如下： 对应于问题③，运行程序后得到结果如下：
对应于问题④，运行程序后得到结果如下： a= 11.7473 2.7028 -1.2251 + 1.4672i -1.2251 - 1.4672i p = 1.0000 -11.9999 -0.0015 25.0012 c= 1 4 7 6
1.2 程序编制
对应于问题①的程序如下：
a=1; b=2; c=3; x=[5 b c; a*b a+c c/b] a=linspace(1,10,10) a= eye(3,2) b= eye(3) a=zeros(3,2) b= zeros(3)
a=ones(3,2) b=ones(3) v=[1,2,3,4] a=diag(v) a=rand(3,2)
xlabel('time/s') ylabel('amplitude/1') axis([0,6.28,-1.5,1.5]) figure(2) subplot(2,1,1) plot(t,y2) xlabel('time/s') ylabel('amplitude/1') axis([0,6.28,-1.5,1.5]) subplot(2,1,2) plot(t,y3) xlabel('time/s') ylabel('amplitude/1') axis([0,6.28,-1.5,1.5])
对应于问题②的程序如下：
a=[1 2 3;4 5 6]' a=[1 2 3; 4 5 6; 2 3 5]; b=inv(a)
a1=[1 2 3;4 5 6].' c=det(a)
b=[1+2i 2-7iቤተ መጻሕፍቲ ባይዱ'
d=rank(a)
b2=[1+2i 2-7i].' a=[1 2 3;3 4 5]; [m,n]=size(a)
f='a*x^2+b*x+c'; solve(f) %方程f = ax2+bx+c 求解
f='a*x^2+b*x+c'; solve(f, 'b') %指定求解变量为b 对应于问题⑤的程序如下：
num=[1,1];den=[1,2,1]; %传递函数分子分母系数(降幂排列) G=tf(num,den) 2.3 运行结果 对应于问题①，运行程序后得到结果如下：
S 1 S2 2S 1
2.2 程序编制
对应于问题①的程序如下：
A = sym('[a , 2*b ; 3*a , 0]') B =['[ a,2*b]'; '[3*a, 0]'] 对应于问题②的程序如下： A1=subs(A, 'b', 'c') A1(2,2)='4*b' 对应于问题③的程序如下： m=int(int('x*exp(-x*y)','x'),'y') 对应于问题④的程序如下： n=laplace(m) k=ilaplace(n) 对应于问题⑤的程序如下：
对应于问题⑤，运行程序后得到结果如下：
116.0023
第二次上机：MATLAB 的符号运算
2.1 问题描述
本次上机共进行了 5 个问题的调试，分别是： ① 命令窗口中符号矩阵的定义；
② 命令窗口中符号矩阵的修改；
③ 计算不定积分
④ 拉氏变换和反拉氏变换
⑤ 方程求解
⑥
用 matlab 表示传递函数 G(s)
A=
B=
[ a, 2*b]
[ a,2*b]
[ 3*a, 0]
[3*a, 0]
对应于问题②，运行程序后得到结果如下：
A1 =
A1 =
[ a, 2*(c)] [ a, 2*(c)]
[ 3*a, 0] [ 3*a, 4*b]
对应于问题③，运行程序后得到结果如下：
m=
1/y*exp(-x*y) 对应于问题④，运行程序后得到结果如下：
n=
k=
1/y/(s+y) 1/y*exp(-y*t) 对应于问题⑤，运行程序后得到结果如下：
ans =
1/2/a*(-b+(b^2-4*a*c)^(1/2))
1/2/a*(-b-(b^2-4*a*c)^(1/2)) 对应于问题⑥，运行程序后得到结果如下：
Transfer function:
s+1
-------------
s^2 + 2 s + 1
第三次上机：控制系统的 MATLAB 分析
3.1 问题描述
本次上机共进行了 8 个问题的调试，分别是： ① 系统传递函数的求解； ② 带反馈的传递函数系统阶跃响应和单位脉冲响应求解； ③ 阻尼比变化时的阶跃响应求解 ④ 系数变化时的系统阶跃响应求解 ⑤ 输入为任意函数的系统响应求解 ⑥ 判定系统稳定性 ⑦ 典型环节的波德（Bode）图绘制 ⑧ 奈奎斯特图的绘制
科技学院
Science and Technology College
上机调试报告
课程名称： 机电控制系统仿真 指导教师： 专业班级： 姓 名： 学 号： 成 绩：
2012 年 12 月
第一次上机：MATLAB 语言基础
1.1 问题描述
本次上机共进行了 3 个问题的调试，分别是： ① 命令窗口中矩阵的定义与常用矩阵的生成； ② 矩阵的运算； ③ 二维图形的绘制 ④ 多项式根的求解与反求 ⑤ 循环语句的应用