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

《机电控制系统仿真实验报告》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、word完美格式科技学院ScienceandTechnologyCollege上机调试报告课程名称：机电控制系统仿真指导教师：专业班级：姓名：学号：成绩：2012年12月精心整理学习帮手word完美格式第一次上机：MATLAB语言基础1.1问题描述本次上机共进行了3个问题的调试，分别是：①命令窗口中矩阵的定义与常用矩阵的生成；②矩阵的运算；③二维图形的绘制④多项式根的求解与反求⑤循环语句的应用1.2程序编制对应于问题①的程序如下：a=ones(3,2)b=ones(3)v=[1,2,3,4]a=diag(v)a

2、=rand(3,2)a=1;b=2;c=3;x=[5bc;a*ba+cc/b]a=linspace(1,10,10)a=eye(3,2)b=eye(3)a=zeros(3,2)b=zeros(3)对应于问题②的程序如下：a=[123;456;235];b=inv(a)c=det(a)d=rank(a)a=[123;345];[m,n]=size(a)length(A)=max(size(A))a=[123;456]'a1=[123;456].'b=[1+2i2-7i]'b2=[1+2i2-7i].'对应于问题③

3、的程序如下：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])t=0:0.01:2*pi;y

4、=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)精心整理学习帮手word完美格式对应于问题④的程序如下：p=[1-12025116];a=roots(p)r=[11.7473,2.7028,-1.2251+1.4672i,-1.

5、2251-1.4672i];p=poly(r)a=[123];b=[12];c=conv(a,b)对应于问题⑤的程序如下：A=zeros(9);fori=1:1:9forj=1:1:9a(i,j)=i*j;ifi

6、251-1.4672ip=1.0000-11.9999-0.001525.0012116.0023c=1476对应于问题⑤，运行程序后得到结果如下：精心整理学习帮手word完美格式第二次上机：MATLAB的符号运算2.1问题描述本次上机共进行了5个问题的调试，分别是：①命令窗口中符号矩阵的定义；②命令窗口中符号矩阵的修改；③计算不定积分④拉氏变换和反拉氏变换⑤方程求解⑥用matlab表示传递函数2.2程序编制对应于问题①的程序如下：A=sym('[a,2*b;3*a,0]')B=['[a,2*b]';'[3*a

7、,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')%指定求解变量为b对应于问题⑤的程序如下：num=[1,1];den=[1,2,1];%传递函数分

8、子分母系数(降幂排列)G=tf(num,den)2.3运行结果对应于问题①，运行程序后得到结果如下：B=[a,2*b][3*a,0]A=[a,2*b][3*a,0]对应于问题②，运行程序后得到结果如下：A1=[a,2*(c)][3*a,4*b]A1=[a,2*(c)][3*a,0]精心整理学习帮手word完美格式对应于问题③，运行程序后得到结果如下：m=1/y*exp(-x*y)对应于