MATLAb实验报告第六次.docVIP

MATLAB数学实验报告 姓名：李犁 学号： 1304120724 指导老师：易昆南 专业班级：统计1203 目录 1 用三种方法求解恰定方程组Ax=b····3 2 判断矩阵是否可对角化·············5 3 互付工资问题·····················6 4 交通流量问题·····················8 学号 1304120724 班级 统计1203 姓名 李犁 指导教师 易昆南 实验题目 判断矩阵是否可对角化 评 分 1、设计（实习）目的： 通过matlab判断矩阵是否可对角化 2、实验内容： 给出判断矩阵是否可对角化的程序,若可对角化返回1,否则返回0; 3.详细设计： function y=trigle(A); y=1;c=size(A); if c(1)~=c(2) y=0; return; end e=eig(A);n=length(A); while 1 if isempty(e) return; end d=e(1); f=sum(abs(e-d)10*eps); g=n-rank(A-d*eye(n)); if f~=g y=0; return; end e(find(abs(e-d)10*eps))=[ ]; end 4：实验结果： A=[4 -3 1 2;5 -8 5 4;6 -12 8 5;1 -3 2 2] trigle(A) ans = 0 学号 1304120724 班级 统计1203 姓名 李犁 指导教师 易昆南 实验题目 互付工资问题 评 分 设计（实习）目的： 了解MATLAB在实际问题中的应用,特别是最优化算法的应用 2、实验内容： 现有一个木工，一个电工和一个油漆工，三人相互同意彼此装修他们自己的房．在 装修之前，他们达成了如下协议：(1)每人总共工作10天(包括给自己家干活在内)；(2) 每人的日工资根据一般的市价在60一80元之间；(3)每人的日工资数应使得每人的总 收入与总支出相等．表12. 1是他们协商后制定出的的工作天数的分配方案，如何计 算出他 们每人应得的工资？ 3.详细设计： 设木工, 电工, 油漆工的日工资分别为x, y, z元, 则由下表 各家应付工资和各人应得收入 在谁家 工人 木工 电工 油漆工 各家应付工资 木工家 2x 1y 6z 2x + y + 6z 电工家 4x 5y 1z 4x + 5y + z 油漆工家 4x 4y 3z 4x + 4y + 3z 各人应得收入 10x 10y 10z 可得 , 即 在Matlab命令窗口输入以下命令 A = [-8,1,6;4,-5,1;4,4,-7]; x = null(A,’r’); format rat, x’ 4：实验结果： Matlab执行后得 ans = 31/36 8/9 1 可见上述齐次线性方程组的通解为x = k(31/36, 8/9, 1)T. 因而根据“每人的日工资一般的市价在60~80元之间”可知 60 (k k k ( 80, 即 ( k ( 80. 也就是说, 木工, 电工, 油漆工的日工资分别为k元, k元, k元, 其中( k ( 80. 为了简便起见, 可取k = 72, 于是木工, 电工, 油漆工的日工资分别为62元, 64元, 72元. 学号 1304120724 班级 统计1203 姓名 李犁 指导教师 易昆南 实验题目 交通流量问题 评 分 1、设计（实习）目的： 用matlab解决实际问题 学会构造m函数文件 2、实验内容： 图给出了某城市部分单行街道的交通流量(每小时过车数).假设： (1)全部流入网络的流量等于全部流出网络的流量； (2)全部流入一个节点的流量等于全部流出此节点的流量。 试建立数学模型确定该交通网络未知部分的具体流量。 3.详细设计： 由网络流量假设，所给问题满足如下线性方程组 ? ? ? ? ? ? ? ? 由 Matlab 软件的编辑器构造 M 函数文件： jt. m 如下： A5 =［ 0 1 -1 1 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 -1 1 0 0 0 1 1 0 0 0 0