机电控制系统仿真_实验报告.doc

科 技 学 院 Science and Technology College 上机调试报告 课程名称： 机电控制系统仿真 指导教师： 专业班级： 姓 名： 学 号： 成 绩： 2012年12月 第一次上机：MATLAB语言基础 对应于问题②的程序如下： 对应于问题③的程序如下： 对应于问题④的程序如下： 对应于问题⑤的程序如下： 1.3 运行结果 对应于问题①，运行程序后得到结果如下： 对应于问题②，运行程序后得到结果如下： 对应于问题③，运行程序后得到结果如下： 对应于问题④，运行程序后得到结果如下： a = 11.7473 2.7028 -1.2251 + 1.4672i -1.2251 - 1.4672i p = 1.0000 -11.9999 -0.0015 25.0012 116.0023 c = 1 4 7 6 对应于问题⑤，运行程序后得到结果如下： 第次上机：MATLAB传递函数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) 对应于问题⑤的程序如下： f=a*x^2+b*x+c; solve(f) %方程f = ax2+bx+c 求解 f=a*x^2+b*x+c; solve(f,b) %指定求解变量为bG=tf(num,den) 2.3 运行结果 对应于问题①，运行程序后得到结果如下： 对应于问题②，运行程序后得到结果如下： 对应于问题③，运行程序后得到结果如下： 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）图绘制 奈奎斯特图的绘制 3.2 程序编制 对应于问题①的程序如下： 对应于问题②的程序如下： 对应于问题③的程序如下： 对应于问题④的程序如下： 对应于问题⑤的程序如下： 对应于问题⑥的程序如下： 对应于问题⑦的程序如下： 对应于问题⑧的程序如下： 3.3 运行结果 对应于问题①，运行程序后得到结果如下： 对应于问题②，运行程序后得到结果如下： 对应于问题③，运行程序后得到结果如下： 对应于问题④，运行程序后得到结果如下： 对应于问题⑤，运行程序后得到结果如下： 对应于问题⑥，运行程序后得到结果如下： 对应于问题⑦，运行程序后得到结果如下： 对应于问题⑧，运行程序后得到结果如下： 第次上机：模型建立如下： “Start Simulation”之后，双击“Scope”的下图： “Command”验证程序及结果如下： a=ones(3,2) b=ones(3) v=[1,2,3,4] a=diag(v) a=rand(3,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=[1 2 3; 4 5 6; 2 3 5]; b=inv(a) c=det(a) d=rank(a) a=[1 2 3;3 4 5]; [m,n]=size(a) length(A)=max(size(A)) a=[1 2 3;4 5 6] a1=[1 2 3;4 5 6]. b=[1+2i 2-7i] b2=[1+2i 2-7i]. xlabel(