用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;
ifE==D
fprintf('等式成立')
else
fprintf('等式不成立')
end;
运行结果：
等式不成立
（2）M文件：
A=round(rand(5)*10);
B=round(rand(5)*10);
C=A+B;
D=A-B;
E=A^2-B^2;
F=C*D;
ifF==E
fprintf('等式成立')
else
fprintf('等式不成立')
end;
运行结果：
等式不成立
（3）M文件：
A=round(rand(5)*10);
B=round(rand(5)*10);
C=(A*B)';
D=B'*A';
ifC==D
fprintf('等式成立')
else
fprintf('等式不成立')
end;
运行结果：
等式成立
2.用逆矩阵求解
M文件：
A=[531;1-3-2;-521];
B=[-830;-590;-2150];
C=inv(A);
formatrat;
X=B*C
运行结果：
X=

123
456
789

>>
矩阵右除求解
M文件：
A=[531;1-3-2;-521];
B=[-830;-590;-2150];
X=B/A
运行结果：
X=

123
456
789

>>
3.M文件：
a1=[11022];
a2=[34083];
a3=[23061];
a4=[93212];
a5=[6-22-92];
A=[a1'a2'a3'a4'a5'];
[R,jb]=rref(A);
A(:,jb)
运行结果：
ans=

139
143
002
281
232
由此可看出，a1,a2,a4为极大线性无关组。
a3=-a1+a2
a5=3a1-2a2+a4
4.M文件：
A=[24-