资源描述：
《matlab语言矩阵运算》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、一、实验程序1、下列运算是否合法，为什么？如合法，结果是多少？(1)result1=a'>>a=[1,2,3;4,5,6];>>b=[2,4,-1;1,3,5];>>c=[1;0;-2];>>d=[1,4,7;8,5,2;3,6,0];>>a'ans=142536(2)result2=a*b>>a*b???Errorusing==>mtimesInnermatrixdimensionsmustagree.(3)result3=a+b>>a+bans=3625811(4)result4=b*d>>b*d
2、ans=312222404913(5)result5=[b;c']*d>>[b;c']*dans=312222404913-5-87(6)result6=a.*b>>a.*bans=28-341530(7)result7=a./b>>a./bans=0.50000.5000-3.00004.00001.66671.2000(8)result8=a.*c>>a.*c???Errorusing==>timesMatrixdimensionsmustagree.(9)result9=a.b>>a.ban
3、s=2.00002.0000-0.33330.25000.60000.8333(10)result10=a.^2>>a.^2ans=149162536(11)result11=a^2>>a^2???Errorusing==>mpowerMatrixmustbesquare.(12)result11=2.^a>>2.^aans=2481632642、用MATLAB求下面的的方程组(1)>>a=[7,2,1,-2;9,15,3,-2;-2,-2,11,5;1,3,2,13];>>b=[4;7;-1;0];
4、>>a^(-1)*bans=0.49790.14450.0629-0.0813>>abans=0.49790.14450.0629-0.0813(2)>>c=[1,1,1,0;1,2,1,-1;2,-1,0,-3;3,3,5,-6];>>d=[1;8;3;5];>>c^(-1)*dans=1.00005.0000-5.0000-2.0000>>cdans=1.00005.0000-5.0000-2.0000>>symsxyzw;>>eq1=x+y+z-1;>>eq2=x+2*y+z-w-8;>>e
5、q3=2*x-y-3*w-3;>>eq4=3*x+3*y+5*z-6*w-5;>>[x,y,z,w]=solve(eq1,eq2,eq3,eq4)x=-2y=1z=5w=-53、已知(1)求矩阵A的秩(rank)>>A=[7,2,1,-2;9,15,3,-2;-2,-2,11,5;1,3,2,13];>>rank(A)ans=4(2)求矩阵A的行列式(determinant)>>det(A)ans=12568(3)求矩阵A的逆(inverse)>>A^(-1)ans=0.1744-0.0303-0.0
6、1250.0270-0.10500.0789-0.01210.00060.00830.01730.0911-0.03110.0095-0.0185-0.01030.0795(4)求矩阵A的特征值及特征向量(eigenvalueandeigenvector)>>eig(A)ans=4.855412.6460+1.8333i12.6460-1.8333i15.85264、关系运算与逻辑运算已知a=20,b=-2,c=0,d=1(1)r1=a>b>>symsabcd;>>a=20;>>b=-2;>>c=0;
7、>>d=1;>>a>bans=1(2)r2=a>b&c>d>>(a>b)&(c>d)ans=0(3)r3=a==b*(-10)>>a==b*(-10)ans=1(4)r4=~b
8、c>>~b
9、cans=0一、实验结果及结果分析（结果见实验程序部分）1、下列运算是否合法，为什么？如合法，结果是多少？(1)result1=a'(2)result2=a*b(3)result3=a+b(4)result4=b*d(5)result5=[b;c']*d(6)result6=a.*b(7)result7=a./b(
10、8)result8=a.*c(9)result9=a.b(10)result10=a.^2(11)result11=a^2(12)result11=2.^a结果分析：矩阵的基本运算操作如下+矩阵加法/矩阵的左除-矩阵减法’矩阵的转置*矩阵乘法.*矩阵对应元素相乘^矩阵的幂./矩阵对应元素左除矩阵的右除.矩阵对应元素右除.^矩阵对应元素的幂a．矩阵的乘法一定要注意维数的对应，前一矩阵的列数要等于后一矩阵的行数；b．X=AB等价于A*X=B,X=B/
