MATLAB上机实验报告MATLAB求解线性方程组和矩阵的初等计算.doc

计算机语言类课程实验报告 课程名称MATLAB语言及应用实验机房计算机公修实验室三院系电气信息工程学院专业电气工程及其自动化班级学号姓名机器号实验日期任课教师实验学时2实验成绩一．实验名称：MATLAB求解线性方程组和矩阵的初等计算二．实验目的和要求 1、掌握利用MATLAB程序编辑器编写应用程序的方法； 2、掌握MATLAB求解线性方程组的方法； 3、掌握MATLAB进行矩阵的初等计算的方法三．实验内容 教材（《MATLAB及其在理工课程中的应用指南，陈怀琛，西安电子科技大学出版社》）P93-1,2四．实验设计方案（实验步骤或开发过程） 1、设a矩阵为各个方程的变量的系数，s为各个变量的列向量，b为等式右边的自然数的列向量，利用s=a\b,即可求出该方程组的解。 2、利用转置和乘法，以及元素乘法分别算出C1,C2,C3，再通过求逆矩阵inv（）这个函数求出C1,C2,C3的逆矩阵。 五．实验中存在问题及解决办法六．实验结果五．实验中存在问题及解决办法六．实验结果五．实验中存在问题及解决办法六．实验结果五．实验中存在问题及解决办法六．实验结果 五．实验中存在问题及解决办法 在第一题中当使用右除时，运行出现错误，只有使用左除，这与矩阵除法的定义有关。 六．实验结果 1、该方程组的解为：s=[-1.4841 -0.6816 0.5337 -1.2429] 2、C1= [ 19 -82 30 d1= [ 0.0062 0.0400 -0.0106 12 27 3 -0.0046 0.0169 0.0030 -38 54 29 ] 0.0168 0.0209 0.0150] C2 = [ -15 16 -24 36 d2= 1.0e+015 * 63 -17 93 -105 [-0.9553 -0.2391 -0.1997 0.2700 22 6 117 -60 0.9667 0.2420 0.2021 -0.2732 19 46 84 -10] -0.4473 -0.1120 -0.0935 0.1264 -1.1259 -0.2818 -0.2353 0.3182] C3= [ 5 16 24 -26 d3= 不能求C3的逆矩阵 -18 -12 -15 72 -2 -21 108 -56] 七．附录（源程序清单） 第一题：求该方程组的解 2012.09.17 clc clear all a=[3 4 -7 -12;5 -7 4 2;1 0 8 -5;-6 5 -2 10] b=[4;-3;9;-8] s=a\b 第二题：求出C1,C2,C3,以及它们的逆矩阵d1,d2,d3 2012.09.17 clc clear all A=[1 4 8 13;-3 6 -5 -9;2 -7 -12 -8] B=[5 4 3 -2;6 -2 3 -8;-1 3 -9 7] C1=A*B C2=A*B C3=A.*B d1=inv(C1) d2=inv(C2) [m,n]=size(C3); if m==n d3=inv(C3) else disp 不能求C3的逆矩阵 end