线性代数_matlab实验指导要点详解.pptVIP

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 第4章 二次型 AE(3,:)=AE(3,:)-2*AE(2,:); AE(:,3)=AE(:,3)-2*AE(:,2); AE(4,:)=AE(4,:)+2*AE(2,:); 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章 二次型 得 即正交变换 将原二次型化为标准形 设计性实验1 房屋装修的工资问题 【实验目的】 1．理解矩阵特征值概念 2．能根据实际问题，建立模型然后使用Matlab相关命令求解 【实验要求】 掌握求解特征值的eig命令、生成对角矩阵的diag命令等 【实验内容】 有三个技术个人分别是木工、电工和管道工，他们准备合作装修自己的新房子。在装修之前约定：每人总共工作20天（包括在自己家）；每人每日的工资平均为100元；每人的日工资应使得每人的总收入和总支出等。需要计算每人的日工资分别是多少，以确定他们的工作日交换是否平衡，如果不平衡，将由谁买单。一个初步的工作日分配方案如下 表3-1 工作日分配方案 工作日 工种 木工 电工 管道工 木工家 4 2 12 电工家 8 10 2 管道工家 8 8 6 【实验方案】 设木工、电工和管道工的日工资分别为：，，。由总收入和总支出相等的约定，建立线性议程组 整理，得 显然问题与矩阵特征值问题有联系，由于矩阵 是正矩阵且每列元素之和均为20，所以20是该矩阵的牲值，于是 就是属于特征值 的特征向量。按约定总工作量决定总工资应该为6000元，则应该有 【实验过程】 MATLAB程序如下 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 Index about eigvalu=20:=”后输入索引值1。 程序继续运行，得出最后计算结果为 daily = 93.94 96.97 109.09 pay = 375.76 193.94 1309.09