用Matlab及建模实践改造工科线代数课程.doc

用Matlab及建模实践改造工科线性代数课程实验作业姓 名： 学 号： 所在学院： 自动化 班 级： 指导教师： 2009年5月一、基础篇：写出用Matlab求解下列问题的过程（在命令窗口的执行过程或M文件均可）。（每题10分）1、利用函数rand和函数round构造一个的随机正整数矩阵A和B，验证一下等式是否成立：（1）AB=BA；（2）（A+B）（A-B）=（3）2、已知，其中 用求逆矩阵和矩阵右除两种方法求矩阵X。3、已知向量组，求它的最大无关组，并用该最大无关组来线性表示其他向量。4、求下面非齐次线性方程组的通解。5、设，计算的特征值与特征向量.6、计算正交的变量替换，化二次型为标准型.7、利用主子式法判定二次型 的正定性.二、应用篇（任选一题或根据自己的专业自行编制与线性代数相关的题目，分析解决问题，写出所用模型及Matlab求解过程）（30分）1、已知：某调料公司用7中成分来制造多种调味品，以下表格列出了6种调味品A、B、C、D、E、F每包所需各成分的量（以盎司为单位）ABCDEF红辣椒31.54.57.594.5姜 黄240816胡 椒120423欧莳萝120413大蒜粉0.510221.5盐0.510221.5丁香油0.250.50210.75问题：(1)、一个顾客为了避免购买全部6种调味制品，它可以只购买其中一部分并用它们配制出其余几种调味制品，为了能配制出其余几种调味品，这位顾客必须购买的最少调味品的种类是多少？写出所需最少的调味品的集合。（2）、由（1）中得到的最小调味品集合是否唯一？能否找到一个最小调味品集合？（3）利用在（1）中找到的最小调味品的集合，按下列成分配制一种新的调味品，并写出需要的调味品的包数。红辣椒：18，姜黄：18，胡椒：9，欧莳萝：9，大蒜粉：4.8，盐4.5，丁香油3.25。（4）另一顾客希望按下列成分表配制一种调味品能，它要购买的自小调味品集合是什么？红辣椒：12，姜黄：14，胡椒：7，欧莳萝：7，大蒜粉：35，盐35，丁香油175。 (5) 在上述问题的解答中，用到了哪些知识点？请列出所用知识点。2、某农场饲养的动物所能达到的最大年龄为15岁，将其分为三个年龄组：第一组，0-5岁；第二组，6-10岁；第三组，11-15岁，动物从第二年龄组起开始繁殖后代，经过长期统计，第二年龄组的动物在其年龄段平均繁殖4个后代，第三组在其年龄段平均繁殖3个后代，第一年龄组和第二年龄组的动物能顺利进入下一个年龄组的存活率分别是1/2和1/4。假设农场现有三个年龄段的动物各1000头，问15 年后农场饲养的动物总数及农场三个段的动物各将达到多少头？指出15年间，动物总增长多少头及总增长率。 答案：1.（1），M文件：A=round(rand(5)*10);B=round(rand(5)*10);E=A*B;D=B*A;if E=D fprintf(等式成立)else fprintf(等式不成立)end;运行结果：等式不成立（2）M文件：A=round(rand(5)*10);B=round(rand(5)*10);C=A+B;D=A-B;E=A2-B2;F=C*D;if F=E fprintf(等式成立)else fprintf(等式不成立)end;运行结果：等式不成立（3）M文件：A=round(rand(5)*10);B=round(rand(5)*10);C=(A*B);D=B*A;if C=D fprintf(等式成立)else fprintf(等式不成立)end;运行结果：等式成立2.用逆矩阵求解 M文件：A=5 3 1;1 -3 -2;-5 2 1;B=-8 3 0;-5 9 0;-2 15 0;C=inv(A);format rat;X=B*C运行结果：X = 1 2 3 4 5 6 7 8 9 矩阵右除求解 M文件：A=5 3 1;1 -3 -2;-5 2 1;B=-8 3 0;-5 9 0;-2 15 0;X=B/A 运行结果：X = 1 2 3 4 5 6 7 8 9 3.M文件：a1=1 1 0 2 2;a2=3 4 0 8 3;a3=2 3 0 6 1;a4=9 3 2 1 2;a5=6 -2 2 -9 2;A=a1 a2 a3 a4 a5;R,jb=rref(A);A(:,jb) 运行结果：ans = 1 3 9 1 4 3 0 0 2 2 8 1 2 3 2 由此可看出，a1,a2,a4为极大线性无关组。a3=-a1+a2a5=3a1-2a2+a44.M文件：A=2 4 -1 4 16;-3 -64 2 -6 -23;3 6 -4 6 19;1 2 5 2 19;b=-2 7 -23 43;B=A b;RA=rank(A);RB=rank(B);if RA=RB&RA=5 X=Abelse RA=RB&RA In e01 at 9X0 = 0 -0.0000 7.3333 0 0.3333D = -2 -9 0 0 0 -2 1 0 0 1 X = -2*k1-9*k2 -3354844070832297/40564819207303340847894502572032 22/3-2*k2 k1 1/3+k2 5.M文件：A=1 3 1 2;3 -1 1 3;1 1 2 5;2 3 5 2;eig(A);V,D=eig(A)V*D*inv(V);a1=V(:,1)a2=V(:,2)a3=V(:,3)a4=V(:,4) 运行结果：e02V = -0.2430 -0.5449 -0.7154 0.3636 0.6643 0.5003 -0.4248 0.3578 0.4054 -0.4947 0.5246 0.5619 -0.5791 0.4562 0.1802 0.6512D = -4.1026 0 0 0 0 -2.5208 0 0 0 0 1.5443 0 0 0 0 9.0791a1 = -0.2430 0.6643 0.4054 -0.5791a2 = -0.5449 0.5003 -0.4947 0.4562a3 = -0.7154 -0.4248 0.5246 0.1802a4 = 0.3636 0.3578 0.5619 0.65126.M文件：A=2 1 0 -1;1 1 -1 0;0 -1 1 1;-1 0 1 2;syms x1 x2 x3 x4;X=x1 x2 x3 x4;f=X*A*X;P,D=eig(A);syms y1 y2 y3 y4;Y=y1;y2;y3;y4;X=P*Yf=Y*D*Y 运行结果：e03 X = -6696817237882403/18014398509481984*y1+209275538683825/562949953421312*y2-5417838953669963/9007199254740992*y3-2708919476834979/4503599627370496*y4 5417838953669959/9007199254740992*y1-5417838953669963/9007199254740992*y2-3348408618941201/9007199254740992*y3-1674204309470601/4503599627370496*y4 1354459738417489/2251799813685248*y1+2708919476834979/4503599627370496*y2-3348408618941205/9007199254740992*y3+3348408618941197/9007199254740992*y4 -6696817237882399/18014398509481984*y1-1674204309470601/4503599627370496*y2-5417838953669955/9007199254740992*y3+677229869208745/1125899906842624*y4 f = -347922205179541/562949953421312*conj(y1)*y1+777977701663083/562949953421312*conj(y2)*y2+910872158600853/562949953421312*conj(y3)*y3+4073544130886953/1125899906842624*conj(y4)*y4 7.M文件：A=2 1 1;1 2 1;1 1 1;c=1;for i=1:3 fprintf(第%d阶主子式为,i) B=A (1:i,1:i) fprintf(第%d阶主子式的值为,i) det(B) if (det(B)由运行结果可看出，该顾客最少应买4种调味品。所需购买的调味品种类的集合为：A，B，D，E（2）不唯一（3）新的调味品的成分组合与（1）中的4种调味品构成非齐次线性方程组求解，若有解，则可以求出所需4种调味品的包数，否则不能调配。M文件：a=3 1.5 7.5 9;b=2 4 8 1;c=1 2 4 2;d=1 2 4 1;e=0.5 1 2 2;f=0.5 1 2 2;g=0.25 0.5 2 1;A=a;b;c;d;e;f;g;m=18 18 9 9 4.8 4.5 3.25;B=A m;RA=rank(A);RB=rank(B);if (RA=RB&RA=n) X=Amelse if(RA=RB&RA（4）把新的调味品与这6种调味品构成新的矩阵，求其极大线性无关组即可。M文件：A=3 2 1 1 0.5 0.5 0.25;B=1.5 4 2 2 1 1 0.5;C=4.5 0 0 0 0 0 0;D=7.5 8 4 4 2 2 2;E=9 1 2 1 2 2 1;F=4.5 6 3 3 1.5 1.5 0.75;G=12 14 7 7 35 35 175;H=A B C D E F G;R,jb=rref(H)H(:,jb)运行结果：e07R = 1 0 2 0 0 1 0 0 1 -1 0 0 1 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0jb = 1 2 4 5 7ans = 3.0000 1.5000 7.5000 9.0000 12.0000 2.0000 4.0000 8.0000 1.0000 14.0000 1.0000 2.0000 4.0000 2.