线性代数matlab实验指导

《线性代数matlab实验指导》由会员分享，可在线阅读，更多相关《线性代数matlab实验指导（82页珍藏版）》请在装配图网上搜索。

1、第1章 矩阵与行列式【矩阵与行列式简介】在计算机日益发展的今天，线性代数起着越来越重要的作用。线性代数起源于解线性方程组的问题，而利用矩阵来求解线性方程组的Gauss消元法至今仍是十分有效的计算机求解线性方程组的方法。矩阵是数学研究和应用的一个重要工具，利用矩阵的运算及初等变换可以解决求解线性方程组等问题。特殊的矩阵方阵的数字特征之一是方阵的行列式，使用行列式可以描述方阵的一些重要的性质。通过计算行列式可求逆矩阵，n个 第1章 矩阵与行列式未知量n个方程的线性方程组的惟一解等问题。向量也是研究矩阵的有力工具，可通过向量组的秩来定义矩阵的秩。向量与矩阵、行列式都是线性代数的重要基本概念，它们是建

2、立线性方程组的解的构造理论与系统求解方法的三个基本工具。 第1章 矩阵与行列式验证性实验实验一 矩阵的运算【实验目的】1理解矩阵、逆矩阵的概念2掌握矩阵的线性运算、乘法、转置、逆、方阵的幂的运算【实验要求】理解矩阵赋值命令、符号变量说明syms、加法+、乘法*、转置、逆矩阵inv、方阵的幂等命令 第1章 矩阵与行列式【实验内容】1已知下列矩阵：（1） ， ； （2） ， 计算 ， ， ， ， ， ， 321 212 113A 101 012 111B dc baA baB 11 BA AB A6 cA A 1A 5A 第1章 矩阵与行列式【实验过程】1（1） A=3 1 1;2 1 2;1 2

3、 3; B=1 1 -1;2 -1 0;1 0 1; C=A+B运行结果：C = 4 2 0 4 0 2 2 2 4 第1章 矩阵与行列式 AB=A*B运行结果：AB = 6 2 -2 6 1 0 8 -1 2 D=6*A运行结果：D = 18 6 6 12 6 12 6 12 18 第1章 矩阵与行列式 sym c; cA=c*A运行结果：cA = 3*c, c, c 2*c, c, 2*c c, 2*c, 3*c F=A运行结果：F = 3 2 1 1 1 2 1 2 3 第1章 矩阵与行列式 G=inv(A)运行结果：G = 1/4 1/4 -1/4 1 -2 1 -3/4 5/4 -1

4、/4 H=A5运行结果：H = 1492 1006 1460 1558 1069 1558 1914 1331 1946 第1章 矩阵与行列式（2） A=sym(a b;c d); B=sym(1 a;1 b); C=A+B运行结果：C = a+1, b+a c+1, d+b AB=A*B运行结果：AB = b+a, a2+b2 c+d, c*a+d*b 第1章 矩阵与行列式 D=6*A运行结果：D = 6*a, 6*b 6*c, 6*d syms c; cA=c*A运行结果：cA = c*a, c*b c2, c*d 第1章 矩阵与行列式 F=A运行结果：F = conj(a), conj(

5、c) conj(b), conj(d) % conj为复数共轭即 G=inv(A)运行结果：G = d/(a*d-c*b), -b/(a*d-c*b) -c/(a*d-c*b), a/(a*d-c*b) 即 db caA cbad acbad c cbad bcbad dA 1 第1章 矩阵与行列式实验二 矩阵的初等变换【实验目的】1理解矩阵初等变换的概念 2掌握矩阵的初等变换及用初等变换求矩阵的逆矩阵【实验要求】掌握矩阵的表示、符号变量说明syms、逆矩阵inv等命令 【实验内容】1已知矩阵 ，求对矩阵实施如下的初等变换后所得矩阵。矩阵的第2行乘以m；矩阵的第3列的n倍加到第1列上去；矩阵的

6、第1行与第2行交换。1） syms m;A=sym(a b c d;e f g h;i j k l);A(2,:)=m*A(2,:)第1章 矩阵与行列式 lkji hgfe dcbaA 第1章 矩阵与行列式运行结果： A = a, b, c, d m*e, m*f, m*g, m*h i, j, k, l2） syms n;A=sym(a b c d;e f g h;i j k l);A(:,1)=A(:,1)+n*A(:,3)运行结果： A = a+n*c, b, c, d e+n*g, f, g, h i+n*k, j, k, l 第1章 矩阵与行列式3） A=sym(a b c d;e

7、f g h;i j k l);A(2,1,:)=A(1,2,:)运行结果：A = e, f, g, h a, b, c, d i, j, k, l 第1章 矩阵与行列式2已知矩阵 ，提取矩阵的第2、3、4行与第3、4列的元素构成矩阵B A=1 2 3 4;5 6 7 8;9 10 11 12;13 14 15 16;B=A(2:4,3:4)运行结果：B = 7 8 11 12 15 16 16151413 1211109 8765 3321A 3已知 ， ， 且 ，求 A=1 0 1;-1 1 1;2 -1 1;B=1 1; 0 1;-1 0;X=inv(A)*B运行结果：X = 3 1 5

8、2 -2 0第1章 矩阵与行列式 112 111 101A 01 10 11B 33 22 11 yx yx yxX BAX BAX 1 实验三 Gauss消元法【实验目的】掌握解线性方程组的Gauss消元法【实验要求】掌握矩阵赋值命令、初等变换相关命令、简化矩阵为阶梯形式rref等命令【实验内容】1用Gauss消元法解线性方程组：（1） ； 第1章 矩阵与行列式 922 1332 1032 82 321 321 321 321 xxx xxx xxx xxx 【实验过程】1（1）解法一：Gauss消元法A=1 2 1 8;1 2 3 10;2 3 1 13;1 2 2 9 ;A(2,:)=A

9、(2,:)-A(1,:);A(3,:)=A(3,:)-2*A(1,:);A(4,:)=A(4,:)-A(1,:)运行结果：A = 1 2 1 8 0 0 2 2 0 -1 -1 -3 0 0 1 1 A(2,3,:)=A(3,2,:)运行结果：A = 1 2 1 8 0 -1 -1 -3 0 0 2 2 0 0 1 1 第1章 矩阵与行列式 A(2,:)=(-1)*A(2,:);A(3,:)=1/2*A(3,:)运行结果：A = 1 2 1 8 0 1 1 3 0 0 1 1 0 0 1 1 A(4,:)=A(4,:)-A(3,:);A(1,:)=A(1,:)-A(3,:);A(2,:)=A(

10、2,:)-A(3,:)运行结果：A = 1 2 0 7 0 1 0 2 0 0 1 1 0 0 0 0 第1章 矩阵与行列式 A(1,:)=A(1,:)-2*A(2,:)运行结果：A = 1 0 0 3 0 1 0 2 0 0 1 1 0 0 0 0由上可知，方程组有惟一解解法二： A=1 2 1 8;1 2 3 10;2 3 1 13;1 2 2 9;A=rref(A)运行结果：A = 1 0 0 3 0 1 0 2 0 0 1 1 0 0 0 由上可知，结果同解法一。 第1章 矩阵与行列式 实验四 行列式及应用【实验目的】1. 了解行列式的概念，掌握行列式的性质2掌握行列式的计算方法3掌握

11、Gramer法则求解线性方程组【实验要求】掌握计算行列式det、解线性方程组solve、生成Vandermonde行列式vander等命令【实验内容】1计算下列行列式的值：（1） ；（2） ；10782 5513 71391 3152 abbbb aabbb bbabb bbbab bbbba 第1章 矩阵与行列式 第1章 矩阵与行列式（1） A=-2 5 -1 3;1 -9 13 7;3 -1 5 -5;2 8 -7 -10;det(A)运行结果：ans = 312（2） A=sym(a b b b b;b a b b b;b b a b b;b b b a b;b b b b a);det

12、(A)运行结果：ans =a5-10*a3*b2+20*a2*b3-15*a*b4+4*b5即行列式的值为 5432235 4152010 babbabaa 2用Gramer法则解线性方程组 A=2 1 -5 1;1 4 -7 6;1 -3 0 -6;0 2 -1 2;A1=8 1 -5 1;0 4 -7 6;9 -3 0 -6;-5 2 -1 2;A2=2 8 -5 1;1 0 -7 6;1 9 0 -6;0 -5 -1 2;A3=2 1 8 1;1 4 0 6;1 -3 9 -6;0 2 -5 2;A4=2 1 -5 8;1 4 -7 0;1 -3 0 9;0 2 -1 -5;a=det(

13、A);a1=det(A1);a2=det(A2);a3=det(A3);a4=det(A4);X=a1/a,a2/a,a3/a,a4/a运行结果：X = 3 -4 -1 1 即得方程组的解为 ， ， ， 第1章 矩阵与行列式 522 963 0674 852 432 421 4321 4321 xxx xxx xxxx xxxx 31 x 42 x 13 x 14 x 实验五 向量【实验目的】理解向量、向量的线性组合与线性表示、向量组的线性相关与线性无关的概念掌握向量组线性相关、线性无关的有关性质及判别法理解向量组的极大线性无关组和向量组的秩的概念会求向量组的极大线性无关组和秩5掌握矩阵秩的求

14、法【实验要求】掌握简化矩阵为阶梯形式rref、计算行列式det、计算矩阵的秩rank等命令【实验内容】设向量： ， ， ， ，问b能否 由 线性表示？ 21011a 52432a 90413a 1445b321 , aaa 第1章 矩阵与行列式 第1章 矩阵与行列式 A=-1 3 1;0 4 4;1 -2 0;2 5 9;b=5;4;-4;1;B=A,b;r=rank(A),rank(B)运行结果：r = 1 2由上可知 ，故方程组有解。2)()( BrAr 2求向量 在基 ， ， 下的坐标.即求满足方程 的解。 A1=1;1;0; A2=1;0;1; A3=0;1;1; A=A1,A2,A3

15、; b=3;-5;9; X=inv(A)*b 输出X = -5.5000 8.5000 0.5000 953 953 1011a 1103a第1章 矩阵与行列式332211 axaxax 第2章 线性方程组【线性方程组简介】线性方程组的求解问题促进了线性代数的产生和发展，利用矩阵、行列式和向量这三个基本工具可较好的解决线性方程组的求解问题。利用解向量所构成的基础解系可方便的描述解空间的基本特征及写出通解，从而较好地描述了线性方程组解的结构问题。 第2章 线性方程组验证性实验实验一 线性方程组【实验目的】理解齐次线性方程组的基础解系、通解及解空间的概念掌握齐次线性方程组的基础解系和通解的求法3理

16、解非齐次线性方程组解的结构及通解的概念【实验要求】掌握分数数据格式format rat、求基础解系null、简化矩阵为阶梯形式rref、解方程组solve等命令 第2章 线性方程组【实验内容】1.求齐次线性方程组 的基础解系及通解。 02 0423 0431 4321 4321 xxx xxxx xxxx 第2章 线性方程组【实验过程】1解法一： format ratA=1 1 1 1;1 3 2 4;2 0 1 -1 ;B=rref(A)运行结果：B = 1 0 1/2 -1/2 0 1 1/2 3/2 0 0 0 0 第2章 线性方程组由上可知，方程组有解 ，其中 ， 是自由未知量。故得方

17、程组的基础解系为 ， 通解为 ，其中 为任意常数。 432 431 2321 2121 xxx xxx3x 4x 012121 1 10232122211 kk 21,kk 第2章 线性方程组解法二： format rat A=1 1 1 1;1 3 2 4;2 0 1 -1 ; B=null(A,r)运行结果：B = -1/2 1/2 -1/2 -3/2 1 0 0 1 syms k1 k2 X=k1*B(:,1)+k2*B(:,2) 第2章 线性方程组运行结果：X = -1/2*k1+1/2*k2 -1/2*k1-3/2*k2 k1 k2即原方程组的通解为 ，其中 为任意常数。 10232

18、1012121 21 kkX 21,kk 第2章 线性方程组2.求方程组 的基础解系及通解。3求方程组 的基础解系及通解。 99129 771210 66321 321 321 xxx xxx xxx 393423 26222 132 54321 54321 54321 xxxxx xxxxx xxxxx 第2章 线性方程组2解法一：A=1 1 1;-10 12 1;1 -9 12;b=66;77;99;r=rank(A),rank(A,b)运行结果：r = 3 3 即系数矩阵的秩等于增广矩阵的秩3，且等于未知量的个数，故原方程组有惟一解。 X=inv(A)*b % X=Ab 运行结果：X =

19、 21 22 23 第2章 线性方程组解法二： syms x1 x2 x3;f1=x1+x2+x3-66;f2=-10*x1+12*x2+x3-77;f3=x1-9*x2+12*x3-99; x1 x2 x3=solve(f1,f2,f3,x1,x2,x3)运行结果： x1 =21x2 =22x3 = 23 第2章 线性方程组3解法一：A=1 1 -2 1 3;2 -1 2 2 6;3 2 -4 -3 -9;b=1;2;3;rA=rank(A)运行结果：rA = 3 rAb=rank(A,b)运行结果： rAb = 3 即系数矩阵的秩等于增广矩阵的秩3，故原方程组有解。 第2章 线性方程组 x

20、0=Ab运行结果：x0 = 1 0 0 0 0 即原线性方程组的一个特解 000010 第2章 线性方程组B=rref(A)运行结果：B = 1 0 0 0 0 0 1 -2 0 0 0 0 0 1 3 由上可知，原方程组的导出组的解为 ，即可得其导出组的基础解系为 ， 故原方程组的通解为 ，其中 为任意常数。 54 32 1 320 xx xxx 001201 130002 22110 kk 21,kk 第2章 线性方程组解法二：A=1 1 -2 1 3;2 -1 2 2 6;3 2 -4 -3 -9;b=1;2;3;X=Ab运行结果：X = 1 0 0 0 0 第2章 线性方程组 B=nu

21、ll(A,r)运行结果：B = 0 0 2 0 1 0 0 -3 0 1故原方程组的通解为 ，其中 为任意常数。 130000012000001 21 kk 21,kk 【矩阵的特征值与特征向量简介】 矩阵的特征值与特征向量是矩阵的数字特征，利用矩阵的特征值与特征向量可判断矩阵的相似、解决矩阵对角化及实对称矩阵正交化等问题，促进了矩阵理论的进一步发展及应用。 第3章 矩阵的特征值与特征向量 验证性实验 矩阵的特征值与特征向量【实验目的】理解矩阵的特征值和特征向量的概念会求矩阵的特征值和特征向量掌握将矩阵化为相似对角矩阵的方法【实验要求】掌握求矩阵的特征多项式poly、求矩阵的特征值和特征向量e

22、ig、矩阵的范数norm、值空间正交化orth、单位阵eye等命令第3章 矩阵的特征值与特征向量 【实验内容】1、设 ，求矩阵A的特征多项式和特征值。第3章 矩阵的特征值与特征向量 310 810 001A 1 A=1 0 0;0 1 8;0 1 3; poly(A)运行结果：ans = 1 -5 -1 5 即矩阵A的特征多项式为 lamda=eig(A)运行结果：lamda = 5 -1 1即矩阵A的特征值为 ， ， 55 23 xxx 51 12 13 第3章 矩阵的特征值与特征向量 3设矩阵 ，求正交矩阵T，使得 为对角矩阵。第3章 矩阵的特征值与特征向量 7311 3711 1173

23、1137A ATT 3解法一： A=7 -3 -1 1;-3 7 1 -1;-1 1 7 -3;1 -1 -3 7 ; kesai,lamda=eig(A)运行结果：kesai = -0.0000 0.7071 0.5000 -0.5000 -0.0000 0.7071 -0.5000 0.5000 0.7071 -0.0000 0.5000 0.5000 0.7071 0 -0.5000 -0.5000lamda = 4.0000 0 0 0 0 4.0000 0 0 0 0 8.0000 0 0 0 0 12.0000即所求正交矩阵为 5000.05000.007071.0 5000.05

24、000.00000.07071.0 5000.05000.07071.00000.0 5000.05000.07071.00000.0T 第3章 矩阵的特征值与特征向量 kesai*A*kesai运行结果：ans = 4 * * * * 4 * * 0 * 8 * 0 * * 12 即经验证有 norm(kesai*kesai-eye(4)运行结果：ans = 9.7171e-016由上可知，所求正交矩阵精度很高。 12*0 *8*0 *4* *4 ATT第3章 矩阵的特征值与特征向量 解法二： A=7 -3 -1 1;-3 7 1 -1;-1 1 7 -3;1 -1 -3 7 ;T=orth

25、(A)运行结果：T = -0.5000 0.5000 -0.7071 0 0.5000 -0.5000 -0.7071 -0.0000 0.5000 0.5000 0.0000 0.7071 -0.5000 -0.5000 -0.0000 0.7071 norm(T*T-eye(4)运行结果：ans = 7.6679e-016由上可知，所求正交矩阵精度很高。第3章 矩阵的特征值与特征向量 实验二 矩阵的三角分解【实验目的】1.理解矩阵的三角分解（又称为LU分解）2.掌握 函数的两种调用方法【实验要求】掌握Matlab软件中有关矩阵LU分解的命令【实验内容】分别用两种方法调用MATLAB中的 函

26、数，实现矩阵LU分解问题。第3章 矩阵的特征值与特征向量 lu lu 【实验方案】 矩阵的三角分解又称为LU分解，它的目的是将一个矩阵分解成一个下三角矩阵L和一个上三角矩阵U的乘积，亦即A=LU,其中L和U矩阵可以分别写成 111 2121 nn lllL nnnnuuu uuuU 222 11211 第3章 矩阵的特征值与特征向量 【实验过程】 （1）求出三角分解矩阵。 1 116 2 3 13; 5 11 10 8; 9 7 6 12; 4 14 15 1 ;,AL U lu A 1 1 0 0 05 16 83 108 1 09 16 47 108 1 11 4 1 0 0L 1 16

27、2 3 130 27 2 57 4 9 40 0 17 9 17 30 0 0 *U 第3章 矩阵的特征值与特征向量 可见，这样得出的 矩阵并非下三角矩阵，这是因为再分解过程中采用了主元素交换的方法。现在考虑 函数的另一中调用方法。1L lu第3章 矩阵的特征值与特征向量 , , L U P lu A1 0 0 01 4 1 0 05 16 83 108 1 09 16 47 108 1 1L 16 2 3 130 27 2 57 4 9 40 0 17 9 17 30 0 0 *U 1 0 0 00 0 0 1 0 1 0 00 0 1 0P 注意，这里得出的P矩阵不是一个单位矩阵，而是单位

28、矩阵的置换矩阵。结合得出的 矩阵可以看出，P矩阵的 ，表明需要将 矩阵的第4行换到第2行， 表明需要将 的第2行换至第3行，将原来第3行换至第4行，这样就可以得出一个真正的下三角矩阵L了。将L,P,U代入并检验，可以精确地还原A矩阵。1L 2.4 1p 1L 3,2 4,3 1p p 1L第3章 矩阵的特征值与特征向量 16 2 3 135 11 10 89 7 6 124 14 15 1inv p L Uans 第3章 矩阵的特征值与特征向量 第4章 二次型【二次型简介】 非线性问题广泛存在于各个科学技术领域，而某些非线性问题在一定的条件下可以转化为线性问题来进行研究。方法之一是通过矩阵的方

29、法将二次型化为标准形，具体包括合同变化法和正交变换法。 第4章 二次型验证性实验二次型及标准形【实验目的】掌握二次型及其矩阵表示了解二次型秩、二次型的标准形的概念会用正交变换等方法化二次型为标准形理解正定二次型、正定矩阵的概念，并掌握其判别法【实验要求】掌握分数数据格式format rat、计算矩阵的秩rank、求矩阵的特征值和特征向量eig、单位阵eye等命令 第4章 二次型【实验内容】1求二次型 的矩阵和二次型的秩。2用合同变换将二次型 化为标准形。yzxzxyzyxf 22332 222 434232413121 222222 xxxxxxxxxxxxf 第4章 二次型【实验过程】1 f

30、ormat ratA=1 -3/2 -1;-3/2 2 1;-1 1 3运行结果：A = 1 -3/2 -1 -3/2 2 1 -1 1 3 rA=rank(A)运行结果：rA = 3 第4章 二次型2 format rat A=0 1 1 -1;1 0 -1 1;1 -1 0 1;-1 1 1 0; E=eye(4); AE=A,E运行结果：AE =0 1 1 -1 1 0 0 0 1 0 -1 1 0 1 0 0 1 -1 0 1 0 0 1 0 -1 1 1 0 0 0 0 1 AE(1,:)=AE(1,:)+AE(2,:); AE(:,1)=AE(:,1)+AE(:,2) 第4章 二次

31、型运行结果：AE =2 1 0 0 1 1 0 0 1 0 -1 1 0 1 0 0 0 -1 0 1 0 0 1 0 0 1 1 0 0 0 0 1 AE(2,:)=AE(2,:)-1/2*AE(1,:); AE(:,2)=AE(:,2)-1/2*AE(:,1)运行结果:AE =2 0 0 0 1 1 0 0 0 -1/2 -1 1 -1/2 1/2 0 0 0 -1 0 1 0 0 1 0 0 1 1 0 0 0 0 1 第4章 二次型 AE(3,:)=AE(3,:)-2*AE(2,:); AE(:,3)=AE(:,3)-2*AE(:,2); AE(4,:)=AE(4,:)+2*AE(2,

32、:); AE(:,4)=AE(:,4)+2*AE(:,2)运行结果：AE =2 0 0 0 1 1 0 0 0 -1/2 0 0 -1/2 1/2 0 0 0 0 2 -1 1 -1 1 0 0 0 -1 2 -1 1 0 1 第4章 二次型得 即正交变换 将原二次型化为标准形 008660.05000.0 5774.05774.02887.05000.0 7887.02113.02887.05000.0 2113.07887.02887.05000.0T TYX 242322213 yyyyf 设计性实验1 房屋装修的工资问题 【实验目的】 1理解矩阵特征值概念 2能根据实际问题，建立模型然

33、后使用Matlab相关命令求解 【实验要求】 掌握求解特征值的eig命令、生成对角矩阵的diag命令等 【实验内容】 有三个技术个人分别是木工、电工和管道工，他们准备合作装修自己的新房子。在装修之前约定：每人总共工作20天（包括在自己家）；每人每日的工资平均为100元；每人的日工资应使得每人的总收入和总支出等。需要计算每人的日工资分别是多少，以确定他们的工作日交换是否平衡，如果不平衡，将由谁买单。一个初步的工作日分配方案如下 表3-1 工作日分配方案 工作日 工种木工电工管道工木工家4 2 1 2电工家8 1 0 2管道工家8 8 6 【实验方案】 设木工、电工和管道工的日工资分别为：，。由总

34、收入和总支出相等的约定，建立线性议程组 整理，得1x 2x 3x1 2 3 11 2 3 21 2 3 34 2 12 208 10 2 208 8 6 20 x x x xx x x xx x x x 1 12 23 34 2 128 10 2 208 8 6 x xx xx x 显然问题与矩阵特征值问题有联系，由于矩阵 是正矩阵且每列元素之和均为20，所以20是该矩阵的牲值，于是 就是属于特征值 的特征向量。按约定总工作量决定总工资应该为6000元，则应该有 4 2 128 10 28 8 6A 1 2 3 , , Tx x x 20 1 2 3 6000 x x x 【实验过程】 MAT

35、LAB程序如下 A=4,2,12;8,10,2;8,8,6; P,D=eig(A); disp(diag(D) II=input(input Index about eigvalu=20:=); if II=0,error(problem have no solution),end alpha=P(:,II); R=alpha./sum(alpha); format bank daily=300*R pay=A*diag(daily) 运行结果：在MATLAB命令窗口中运行程序，屏幕将显示出A的三个特征值 20.00 -2.00 2.00 由于第一个特征值恰好为20，在提示符“input In

36、dex about eigvalu=20:=”后输入索引值1。 程序继续运行，得出最后计算结果为 daily = 93.94 96.97 109.09 pay = 375.76 193.94 1309.09 751.52 969.70 218.18 751.52 775.76 654.55 每人的日工资由变量daily的数据给出。 结果表明：表3-2 日工资列表 最后的二维数组给出了二维数组，表明付款明细账，行表示支付，列表示收取。显然第一行相加等于第一列相加，第二行相加等于第二列相加，第三行相加等于第三列相加。 表3-3 工资支付收取方案 工种木工电工管道工日工资9 3 .9 4 9 6 .

37、9 7 1 0 9 .0 9 支付收取木工电工管道工木工3 7 5 .7 6 1 9 3 .9 4 1 3 0 9 .0 9电工7 5 1 .5 2 9 6 9 .7 0 2 1 8 .1 8管道工7 5 1 .5 2 7 7 5 .7 6 6 5 4 .5 5 设计性实验2 卷烟叶组配方设计【实验目的】掌握线性方程组的各种解法。能根据实际问题，使用Matlab建立相应的线性方程组并求解。【实验要求】1掌握几种线性方程组（定解方程组、不定方程组、超定方程组、奇异方程组、符号方程组）的解法。2能用Matlab求解不同类型线性方程组的方法。 【实验内容】 如何提高卷烟抽吸时的感官质量，以及如何降低

38、烟气中的有害成分始终是卷烟制造工业的重中之中。卷烟的叶组配方，即卷烟中的混合烟丝，是由多种单料烟叶按照某种特定的百分比例组合而成的，其化学成分含量（包括总糖、总碱、氯、氮、磷、氧化钾的含量）与其感官质量指标（包括光泽、香气、谐调、杂气、刺激性、余味）和烟气化学成分含量（包括焦油量、CO量、烟气烟碱量）之间存在着一定的映射关系，也就是说特定化学成分的叶组配方对应着其特定的感官质量和烟气化学成分，叶组配方化学成分的含量从另一个角度反映了其感官质量和烟气化学成分含量。因此，在对叶组配方进行设计时，通常要求在确定叶组配方化学成分含量的前提下来确定进入叶组的各 种单料烟叶的百分比例，这样既保证了卷烟叶组

39、的感官质量，又确保了其烟气化学成分含量不会太高。本实验要求设计出根据叶组配方化学成分含量要求确定各种单料烟叶百分比例的数学模型，并对模型求解，给出问题的结果。 配方设计师根据将要生产的卷烟的抽吸风格的需要选择了12种单料烟叶进入叶组，其化学成分含量与叶组所要求的化学成分含量如表3-10所示，要求根据叶组所要求的化学成分含量确定出各种单料烟叶在叶组中所含的百分比例。 【实验方案】 该问题的目的是要在确定叶组化学成分的前提下求出每种单料烟叶在叶组中所占的百分比例，而叶组的某种化学成分是由各种单料烟叶相对应的化学成分按照其百分比例组合而成的，并且各种单料烟叶的百分比例之和应该为100%，因此，我们可

40、据此列出线性方程组，求解该线性方程组即可求得每种单料烟叶的百分比例。 设编号为i的单料烟叶在该叶组中所占的百分比例为 ，即编号依次为1，2，12的单料烟叶在叶组中所占的百分比例分别为 ， ， ，根据前面表中给出的数据可列出下面的线性方程组 ix1x 2x 12x 求解该线性方程组即可求得12种单料烟叶的百分比例。 观察前面所列出的方程组，未知数的个数大于方程组的个数，该线性方程组是不定方程组，有多个解，可利用线性代数中求解不定线性方程组的方法，求出该方程组的特解与通解。 1 18.243.222.226.289.125.237.236.210.211.229.295.101.2 23.019.

41、020.021.026.020.021.020.027.025.026.022.020.0 05.239.221.206.279.122.297.136.299.114.286.184.105.2 23.045.026.027.036.007.013.013.030.022.022.021.024.0 43.248.241.347.207.200.359.291.276.135.277.107.296.2 26.2544.2236.1976.2484.1972.1867.2213.2431.2850.2564.2649.3275.31 121110987654321 1211109876543

42、21 121110987654321 121110987654321 121110987654321 121110987654321 121110987654321 xxxxxxxxxxxx xxxxxxxxxxxx xxxxxxxxxxxx xxxxxxxxxxxx xxxxxxxxxxxx xxxxxxxxxxxx xxxxxxxxxxxx 【实验过程】 clear all clc %输入方程组的系数矩阵A=31.75,32.49,26.64,25.5,28.31,24.13,22.67,18.72,19.48,24.76,19.36,22.44;2.96,2.07,1.77, 2.35,

43、1.76,2.91,2.59,3,2.07,2.47,3.41,2.48;0.24,0.21,0.22,0.22,0.3 ,0.13,0.13,0.07,0.36,0.27,0.26,0.45;2.05,1.84,1.86,2.14,1.99,2.36,1.97,2.22,1.79,2.06,2.21,2.39;0.2,0.22,0.26,0.25,0.27,0.2,0.21,0.2,0.26,0.21,0.2,0.19;2.01,1.95,2.29,2.11,2.1,2.36,2.37,2.25,1.89,2.26,2.22,2.43;1,1,1,1,1,1,1,1,1,1,1,1; B=2

44、5.26;2.43;0.23;2.05;0.23;2.18;1; X0=AB %求方程组的一个特解运行结果：X0 = 0 0.1339 0 0 0.3995 0 0.1740 0.0834 -0.0372 0 0.2310 0.0153 X=null(A) %求原方程组对应的齐次线性方程组的基础解系运行结果：X = 0.0639 -0.0745 -0.1521 -0.5375 -0.1014 0.0263 0.0253 -0.0824 0.5424 0.0405 0.4675 0.1258 -0.2093 0.0146 -0.3859 -0.4757 -0.2368 0.0695 0.1650

45、 -0.5181 0.1032 -0.2302 0.1131 -0.0731 0.6750 -0.0780 0.7395 -0.0729 -0.0438 0.1515 -0.5769 -0.1814 -0.1848 -0.0828 0.2353 0.3931 -0.3908 0.0010 -0.1888 -0.0240 -0.1066 0.3212 -0.1129 -0.2211 -0.0168 0.0058 0.0857 0.8900 -0.0913 -0.1117 0.1865 -0.0251 -0.0056 0.5313 0.1327 -0.0090 -0.1586 -0.2537 -0

46、.0149 -0.0771 所以，原方程组的通解为 0153.0 2310.00 0372.00834.0 1740.00 3995.000 1339.000771.01327.0 1117.0 0168.0 0240.02353.0 1515.0 6750.0 5181.0 3859.00405.0 1014.00149.05313.0 0913.0 2211.0 1888.0 0828.0 0438.0 0731.01650.0 0146.0 5424.0 5375.02537.0 0056.08900.0 1129.00010.0 1848.0 0729.01131.0 0695.0 2

47、093.0 0824.0 1521.01586.0 0251.00857.0 3212.0 3908.0 1814.07395.0 2302.0 2368.01258.0 0253.0 0745.00090.01865.0 0058.0 1066.0 3931.0 5769.0 0780.01037.0 4757.04675.0 0263.0 0639.0 54321 12111098765 4321 kkkkkxxxxxxxxxxxx 如果要求出原方程组的一个特定解，即求出叶组配方中每种单料烟叶的具体百分比例，就将 , , , , 取某特定的值，代入上式进行计算就可以了，但一定要确保计算出的

48、特定解的所有元素都要大于0且小于1，因为每种单料烟叶在叶组中所占的百分比例必定是大于0且小于1的。如取 =-0.001， =0.15， =0.001， =-0.01， =-0.2，计算出的结果为x1=0.0143，x2=0.1241，x3=0.0952，x4=0.0670，x5=0.2308，x6=0.0811，x7=0.1010，x8=0.0311，x9=0.0165，x10=0.0370，x11=0.1952，x12=0.0068，即表示编号为1，2，12的单料烟叶在叶组中所占的百分比例分别为1.43%，12.14%，9.52%，6.7%，23.08%，8.11%，10.1%，3.11%，1.65%，3.7%，19.52%，0.68%。然后，配方设计师就可按照每种单料烟叶的百分比例制作叶组配方来进行抽吸评定，如果抽吸评定的结果不符要求，则将 , , , , 重新取值再进行计算，直到找到满意且符合生产要求的叶组配方为止。1k 2k 3k 4k 5k1k 2k3k 4k 5k 1k 2k 3k 4k 5k