资源描述：
《中南大学自动化胡杨系统仿真实验报告完整版.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、中南大学系统仿真实验报告指导老师：实验者：学号：专业班级：完成时间：4949实验一MATLAB中矩阵与多项式的基本运算基本命令训练：1．eye(m)取n=3，程序如下：>>eye(3)ans=100010001结论：eye(n)用于产生n×n维的单位矩阵，在这里n取3，故产生3×3维单位矩阵。2．one(n)、ones(m,n)对ones(n)取n=5，对ones(m,n)取m=3，n=5，程序如下：>>ones(5)ans=1111111111111111111111111>>ones(3,5)ans=111111111111111结论：ones(n)用于产生n×n维的全
2、1矩阵，在这里n取5，故产生5行5列全1矩阵。ones(m,n)用于产生m×49n维的全1矩阵，在这里m取3，n取5，故产生3行5列的全1矩阵。1．zeros(m,n)取m=3，n=2，程序如下：>>zeros(3,2)ans=000000结论：zeros(m,n)用于产生m×n维全0矩阵，在这里m取3，n取2，故产生3行2列全0矩阵。4．rand(m,n)取m=3，n=4，程序如下：>>rand(3,4)ans=0.95010.48600.45650.44470.23110.89130.01850.61540.60680.76210.82140.7919结论：rand(m
3、,n)用于产生m×n维平均分布的随机矩阵，在这里m取3，n取4，故产生了3行4列的随机矩阵5．diag(v)先创建3×3的魔方矩阵v，在进行diag(v)运算，程序如下：>>v=magic(3)diag(v)v=816357492ans=49852结论：diag(v)用于得到矩阵v的对角元素6．AB、A/B、inv(A)*B、B*inv(A)先创建A、B两个矩阵，在进行运算，程序如下：>>A=[1,2;3,4];>>B=[5,6;7,8];>>a=ABb=A/Bc=inv(A)*Bd=B*inv(A)a=-3-445b=3.0000-2.00002.0000-1.000
4、0c=-3.0000-4.00004.00005.000049d=-1.00002.0000-2.00003.0000结论：’/’表示矩阵右除，’’表示矩阵左除，inv(A)表示求A的逆矩阵，由实验结果可知，矩阵左除与右除结果不一样，矩阵左乘与右乘结果也不一样，AB是求AX=B的解，A/B是求XB=A的解。所以编程求解的时候要注意区分他们的区别。7、roots(p)>>symsx;>>a=3*x.^3+2*x+5;>>p=[3,0,2,5]>>roots(p)p=3025ans=0.5000+1.1902i0.5000-1.1902i-1.0000结论：roots(p)
5、函数用于求多项式的根，以向量形式输入多项式的系数，对应降幂排列，然后调用函数，即可求得对应多项式的根。8、poly>>A=[1,2;3,4];>>poly(A)ans=491.0000-5.0000-2.0000结论：poly(A)用于求矩阵A的特征多项式的系数9．conv、deconv>>A=[1,2];>>B=[3,4];>>a=conv(A,B)b=deconv(A,B)a=3108b=0.3333结论：使用conv函数对多项式进行乘法运算，其使用格式为c=conv(a,b),其中a和b为两个多项式的系数向量，c为相乘所生成的多项式的系数向量。使用deconv(a,b
6、)完成除法运算。10．A*B与A.*B的区别>>A=[1,2];>>B=[5,6]';>>a=A*BA=[1,2];B=[5,6];b=A.*Ba=17b=512结论：A.*B称为“点乘”、“位乘“，即为两个行列数相同的矩阵，对应位置一一相乘，得到的结果依位置对应到结果矩阵中，而A*B为矩阵乘法，要求前者A的列数与后者B行数对应。4911．who与whos的使用>>A=[1,2;3,4];>>whowhosYourvariablesare:ANameSizeBytesClassAttributesA2x232double结论：who给出变量的名称清单；而whos给出所有变量
7、的详细信息。12．disp、size(a)、length(a)的使用>>a='helloworld';>>disp(a)a=[1,2,3,4];B=size(a)C=length(a)helloworldB=14C=4结论：disp函数的作用是直接将内容输出在Matlab命令窗口中，size(a)表示矩阵每个维度的长度，length(a)表示矩阵a的最大的长度。49实验二MATLAB绘图命令基本命令训练1．plot2．loglog3．semilogx4．semilogy5．polar6．title7．xlabel8
