matlab上机教学

1、 线性代数线性代数上机上机 作业格式： 1.姓名，班级，学号N= A= , B= .要求：要求： 写成201165083张三.txt, 201165083张三.doc文件 N= 201165083共9位a= 后两位 83b=第4-5位 86c=第6-7位 d=第4,8位e=第1,8位f=第5,9位g=第4,9位h=第5,7位 60上机界面上机界面MATLAB中基本代数运算符中基本代数运算符 运算运算 符号符号 举例举例加法，加法，a+b + 5+3减法，减法，a-b - 5-3乘法，乘法，ab * 5*3 除法，除法，ab / or 48/4=448=12乘幂，乘幂，ab 52=25变量及数组

2、输入变量及数组输入 2：矩阵方式输入（：矩阵方式输入（请记下红色部分！请记下红色部分！）a=1,2,3;4,5,6;7,8,0 %矩阵输入矩阵输入 （a为为3阶方阵）阶方阵）b=366;804;351 %列矩阵输入列矩阵输入c=366;804;351 %行矩阵（转置）输入行矩阵（转置）输入det(a) %方阵行列式方阵行列式inv(a) %方阵的逆方阵的逆x=ab %ax=b方程组的解方程组的解y=inv(a)*b %与与x相同相同r=rank(A), %矩阵的秩矩阵的秩disp(a,b,x) %显示矩阵显示矩阵变量及数组输入变量及数组输入 3：矩阵的简单运算：矩阵的简单运算c=inv(a)

3、%方阵的逆阵方阵的逆阵y=c*b %矩阵乘积矩阵乘积d=a b;disp(d) %矩阵拼接矩阵拼接d=a; disp(d) %矩阵转置矩阵转置g=2*a+3 %常数乘矩阵，各元素加常数乘矩阵，各元素加3p=eye(3) %3阶单位矩阵阶单位矩阵y=a.*p %两矩阵对应元素乘积两矩阵对应元素乘积zeros(3) %3阶零矩阵阶零矩阵常见任务常见任务 矩阵的赋值和其加、减、乘、除（求逆）命令； 矩阵化为最简行阶梯型的计算命令；U0,ip=rref(A) 多元线性方程组MATLAB求解的几种方法；x=inv(A)*b, U=rref(A) 行列式的几种计算机求解方法；D=det(A),L,U=lu

4、(A);D=prod(diag(L) n个m维向量组的相关性及其秩的计算方法和命令; r=rank(A),U=rref(A) 求欠定线性方程组的基础解系及超定方程解的MATLAB命令；xb=null(A) 矩阵的特征方程、特征根和特征向量的计算命令；f=poly(A);P,D=eig(A) 化二次型为标准型的MATLAB命令；yTDy=xTAx; 其中y=P-1x, 例例1 求下列矩阵列向量组的一个最大无关组. 解：解：在Matlab中输入： a=1,-2,-1,0,2;-2,4,2,6,-6;2,-1,0,2,3;3,3,3,3,4; b=rref(a) 求得: 所以 是一个极大无关组,且向

5、量组的最大无关组向量组的最大无关组3b = 1.0000 0 0.3333 0 1.7778 0 1.0000 0.6667 0 -0.1111 0 0 0 1.0000 -0.3333 0 0 0 0 0Warning: Could not get change notification handle for local C:Program FilesMATLABR2007bwork.Performance degradation may occur due to on-disk directory change checking. a=1,-2,-1,0,2;-2,4,2,6,-6;2,-

6、1,0,2,3;3,3,3,3,4; b=rref(a) %不打；则计算机将显示b的值b = 1.0000 0 0.3333 0 1.7778 0 1.0000 0.6667 0 -0.1111 0 0 0 1.0000 -0.3333 0 0 0 0 0 解：解： 在Matlab中输入：112311232242,.3333baaa baaa4故、解线性方程组、解线性方程组 直接解法直接解法 利用左除运算符的直接解法利用左除运算符的直接解法 对于线性方程组对于线性方程组Ax=b，可以利用左除运，可以利用左除运算符算符“”求解：求解： x=Ab例例1 用直接解法求解下列线性方程组用直接解法求解下

7、列线性方程组.命令如下：命令如下：A=2,1,-5,1;1,-5,0,7;0,2,1,-1;1,6,-1,-4;b=13,-9,6,0;x=Ab5 A=2,1,-5,1;1,-5,0,7;0,2,1,-1;1,6,-1,-4;b=13,-9,6,0; x=Ab %不打；则计算机将显示x的值x = -66.5556 25.6667 -18.7778 26.5556 例子三、求秩例子三、求秩 A=2,1,-5,1;1,-5,0,7;0,2,1,-1;1,6,-1,-4; r=rank(A); r= % = 计算机不显示r的值? r= |Error: Expression or statement

8、is incomplete or incorrect. rank(A) %不打；则计算机将显示rank(A)的值ans = 4 rr = 4 三.Help四、上机作业四、上机作业123412152217abcdAefgh1234121522175716deabBN= 200865083共9位a= 后两位 83b=第4-5位 86c=第6-7位 d=第4,8位e=第1,8位f=第5,9位g=第4,9位h=第5,7位 6012bfc 1. A,B列向量组的一个最大无关组列向量组的一个最大无关组.2. A*B+B*A 3. Ax=b,Bx=b的解的解4. A,B的行列式，逆，秩的行列式，逆，秩5.

9、A-B,A+B,A*B的行列式，逆，秩的行列式，逆，秩6. 见后面一页见后面一页习题习题11页页 第第3题的题的(1)-(3)24页第页第1题题 38页第页第1题题(5)-(7)49页第页第1题题(1)(2)49页第页第2题题(1)-(3) 126页第页第1题题(1)-(2)N= 200865083共9位1= 最后1位 32=最后第2位 8例如：15： 变成36 121223344,0.1, 0.2,0.5, 0.7, 0.9, 0.9512152217abcdbAfcefgh 选做：选做：Ax=b的解在下列不同的取值时变化如何？的解在下列不同的取值时变化如何？作业格式： 1.姓名，班级，学号

10、N= A= , B= .要求：要求： 写成200865083张三.txt, 200865083张三.doc文件 N= 200865083共9位a= 后两位 83b=第4-5位 86c=第6-7位 d=第4,8位e=第1,8位f=第5,9位g=第4,9位h=第5,7位 60上机出现的问题上机出现的问题1 . a=1,2,3;4,5,6;7,8,0；a=(1 2 3,4 5 6, 7 8 0)2. inv(a)Inf Inf Inf Inf3. x=Ab；x=b/A , b为行向量，为行向量， Inv(A)*b通过M文件创建矩阵 当矩阵的规模较大时，直接输入法就力不从心，出现差错也不容易修改。因此

11、可以使用M文件生成矩阵。 方法是：建立一个M文件，其内容是生成矩阵的命令，在MATLAB的命令窗口中输入此文件名，即将矩阵调入工作空间（写入内存）。 用建立M文件的方式生成矩阵 (1) 建立M文件mydata.m内容如下 %生成矩阵 A=1，2，3；4，5，6；7，8，9 (2) 运行M文件mydata.m mydata 则生成矩阵A。 求基础解系：Z = null(A) 标准正交基 for the null space of A Z = null(A,r) rational basis for the null space obtained from the reduced row eche

12、lon form A = 1,2, 3; 1, 2, 3; 1, 2, 3; ZR = null(A,r) ZR = -2 -3 1 0 0 1 B=null(A)B = 0.96362411165943 0 -0.14824986333222 -0.83205029433784 -0.22237479499833 0.55470019622523Warning: Could not get change notification handle for local C:Program FilesMATLABR2007bwork.Performance degradation may occur

13、 due to on-disk directory change checking. a=1,-2,-1,0,2;-2,4,2,6,-6;2,-1,0,2,3;3,3,3,3,4; b=rref(a) %不打；则计算机将显示b的值b = 1.0000 0 0.3333 0 1.7778 0 1.0000 0.6667 0 -0.1111 0 0 0 1.0000 -0.3333 0 0 0 0 0 上机作业上机作业123412152217abcdAefgh1234121522175716deabBN= 201165083共9位a= 后两位 83b=第4-5位 86c=第6-7位 d=第4,8位e=第1,8位f=第5,9位g=第4,9位h=第5,7位 6012bfc A,B特征值，特征向量特征值，特征向量.A的的20次方次方121223344,0.1, 0.2