资源描述:
《matlab控制系统计算机辅助设计》由会员上传分享,免费在线阅读,更多相关内容在学术论文-天天文库。
1、实验目录实验一:Matlab环境熟悉与基本运算(设计型)实验二:Matlab语言程序设计(设计型)实验三:控制系统模型的建立(设计型)实验四:Simulink仿真入门(验证型)实验五:控制系统时域仿真分析(设计型)实验六:Simulink环境下时域仿真实验七:控制系统根轨迹仿真分析实验八:控制系统频域仿真分析(设计型)实验名称:Matlab环境熟悉与基本运算(设计型)1、矩阵运算(1)矩阵的乘法A=[12;34];B=[55;78];y=A^2*By=105115229251(2)矩阵除法A=[123;456;789];B=[100;020;003];y1=AB警告
2、:矩阵接近奇异值,或者缩放错误。结果可能不准确。RCOND=1.541976e-18。y1=1.0e+16*-0.45041.8014-1.35110.9007-3.60292.7022-0.45041.8014-1.3511y2=A/By2=1.00001.00001.00004.00002.50002.00007.00004.00003.0000(3)矩阵的转置及共轭转置A=[5+i,2-i,1;6*i,4,9-i];y1=A.'y1=5.0000+1.0000i0.0000+6.0000i2.0000-1.0000i4.0000+0.0000i1.0000+0.
3、0000i9.0000-1.0000iy2=A'y2=5.0000-1.0000i0.0000-6.0000i2.0000+1.0000i4.0000+0.0000i1.0000+0.0000i9.0000+1.0000i(4)使用冒号选出指定元素A=[123;456;789];y1=A(1:2,3)y1=36y2=A(2:3,:)y2=456789(5)复数矩阵的生成:symsa1a2a3a4b1b2b3b4A=[a1+b1*ja2+b2*j;a3*exp(b3*j)a4*exp(b4*j)]A=[a1+b1*1i,a2+b2*1i][a3*exp(b3*1i),a
4、4*exp(b4*1i)]a1=3;a2=-2;a3=9;a4=23;b1=5;b2=3;b3=6;b4=33;A=[a1+b1*ja2+b2*j;a3*exp(b3*j)a4*exp(b4*j)]A=3.0000+5.0000i-2.0000+3.0000i8.6415-2.5147i-0.3054+22.9980i2、多项式p=1024roots(p)%用roots函数求多项式的根ans=0.5898+1.7445i0.5898-1.7445i-1.1795+0.0000i①A=[1.2350.9;51.756;3901;1234];P=poly(A);poly2
5、sym(P)ans=x^4-(69*x^3)/10-(3863*x^2)/50-(8613*x)/100+12091/20②polyval(P,20)ans=7.2778e+04③polyvalm(P,A)ans=1.0e-11*-0.4093-0.4849-0.3876-0.4519-0.5002-0.8072-0.6004-0.5684-0.4704-0.6196-0.6480-0.5230-0.3297-0.4455-0.3595-0.36383、基本绘图命令(1)t=0:0.01:2*pi;y=cos(t);plot(t,y)运行结果如右图所示:(2)t=0:
6、0.01:2*pi;y1=cos(t-0.25);y2=sin(t-0.5);plot(t,y1,t,y2)运行结果如右图所示:4、基本绘图控制t=0:0.01:4*pi;x1=10*sin(t);plot(t,x1,'-.b*')axis([-2020-2020])gridontitle('正弦曲线x1=10*sin(t)图')xlabel('t=0:0.01:4*pi')ylabel('x1')text(pi/2,10,'最大值')运行结果如右图所示:实验名称:Matlab语言程序设计(设计型)1、编写命令文件:计算1+2+…+n<5000时的最大n值;解:建立M
7、文件,文件名称记为two1,编辑程序内容如下所示:s=0;n=1;while(s<5000),s=s+n;n=n+1;end,n保存,并在命令窗口中输入two1,运行结果如下所示:two1n=1012、编写函数文件:解:①用for循环结构编写程序:建立函数文件,文件名称记为findsum1,编辑程序内容如下所示:function[s1]=findsum1()%findsum1用for循环编写的程序%求2的0到30次幂的和,s表示和,i表示执行次数s1=0;fori=0:30,s1=s1+2^i;Endend保存,并在命令窗口中输入[s1]=finds