Matlab上机作业部分参考答案.pptVIP

Matlab 上机课作业 吴梅红 2012.11.26 冠庄脾蝉瓦只蹲腆耗瑚驼炯肇剔铺横欢刽秧蘑簇撵赦懦苗来讹滇贯裹磺胶Matlab上机作业部分参考答案Matlab上机作业部分参考答案 请将下面给出的矩阵A和B输入到MATLAB环境中，并将它们转换成符号矩阵。 对矩阵A、B两个矩阵进行分析，判定它们是否为奇异矩阵，得出矩阵的秩、行列式、迹和逆矩阵，检验得出的逆矩阵是否正确。同时，给出A和B矩阵的特征多项式、特征值与特征向量，并对它们进行LU分解。 拦左金彼腋习舱识币俄饱奶凑莎钦豌悍棵疼逸厅绝丙赦剥澡寡该脱须阁医Matlab上机作业部分参考答案Matlab上机作业部分参考答案 参考答案 A=[5,7,6,5,1,6,5; 2,3,1,0,0,1,4; 6,4,2,0,6,4,4; 3,9,6,3,6,6,2;10,7,6,0,0,7,7; 7,2,4,4,0,7,7; 4,8,6,7,2,1,7]; A=sym(A); rank(A) det(A) trace(A) B=inv(A); A*B eig(A) p=poly(A) [L U]=lu(A) p=sym2poly(p) 末灶督娜视敞素镰争地甥苗埃瘸笑琉处历察配孪炽喇睡翅夜秃收券协丝融Matlab上机作业部分参考答案Matlab上机作业部分参考答案 参考答案 2. 求方程组的根 解： a=[1 4 -3;2 5 -1;1 6 1]; b=[2;11;12]; x=a\b x = 4 1 2 x1+4x2-3x3=2 2x1+5x2-x3=11 x1+6x2+x3=12 撤舟肺铂廊喘议津教超繁达狂仅啄悲游缺鼎鬃孺揖队损样系瓜厉诞箍蛋怠Matlab上机作业部分参考答案Matlab上机作业部分参考答案 3. 解方程组: [答案] a=[2 9 0;3 4 11;2 2 6]; b=[13 6 6]; x=a\b x = 7.4000 -0.2000 -1.4000 奎耗动袒梦奶茧笋断肠软耘翅保演滤瑞詹唤枢瓜玩柞抢旅掇敞恐彬嘉脆嘛Matlab上机作业部分参考答案Matlab上机作业部分参考答案 4. 试判定下面的线性代数方程是否有解。 解：由秩判定矩阵可以得出如下结果。 A=[16,2,3,13; 5,11,10,8; 9,7,6,12; 4,14,15,1];B=[1; 3; 4; 7]; [rank(A), rank([A B])] ans = 3 4 由得出的结果看，A, [A;B] 两个矩阵的秩不同，故方程是矛盾方程，没有解。 哎椅蚁邱星鞋巢凯炽铺庙恳旅等安钦滨贴倡名暇燕汗抿卤颗逮境示跋警衍Matlab上机作业部分参考答案Matlab上机作业部分参考答案 5. 试求下面齐次方程的基础解系 压潮滤秉官拾走撼雍疹喜隶乓宏除擒麓选毒溃瞪缎汲更鳖醛涪演肖渗犯痈Matlab上机作业部分参考答案Matlab上机作业部分参考答案 【求解】可以将方程写成矩阵形式，得出的两列向量为方程的基础解系。 A=[6,1,4,-7,-3; -2,-7,-8,6,0; -4,5,1,-6,8; -34,36,9,-21,49; -26,-12,-27,27,17]; A=sym(A); rank(A) ans = 3 null(A) ans = [ 191/34, 95/17] [ 0, 1] [ 1, 0] [ 109/34, 103/34] [ 173/34, 151/34] 挨岁液敝吟嗣冰土冀秃啃朗境翻释阳砌暮蒙搓悉讣教闰技润馁蚁筛肋蝉悠Matlab上机作业部分参考答案Matlab上机作业部分参考答案 6. 求解方程组的通解 搽鞋撂殷氟美滑歹强涯腰菌薄煞笆条殿乖上琉伶酥从感敌虫惭篇莉成汹溢Matlab上机作业部分参考答案Matlab上机作业部分参考答案 A=[1 2 4 6 -3 2;2 4 -4 5 1 -5;3 6 2 0 5 -9;2 3 0 4 0 1;0 -4 -5 2 1 4 ;5 5 -3 6 6 -4 ]; b=[4 3 -1 8 -5 2]; B=[A b]; C=rref(B) ? A=[1 2 4 6 -3 2;2 4 -4 5 1 -5;3 6 2 0 5 -9;2 3 0 4 0 1;0 -4 -5 2 1 4 ;5 5 -3 6 6 -4 ]; b=[4 3 -1 8 -5 2]; B=[A b]; n=6; R_A=rank(A) R_B=rank(B) format rat if R_A==R_BR_A==n X=A\b elseif