[理学]MATLAB教学.ppt

Z = null(A) 矩阵A的零空间标准正交基Z = null(A,r)零空间的有理数基 8、左除运算计算非齐次方程组的特解 例7：求解方程组 的通解 A=[2,4,-1,4,16; -3,-6,2,-6,-23; 3,6,-4,6,19; 1,2,5,2,19]; [R, s] =rref ([A b]) x=null(A,r) b=[-2;7;-23;43]; x0=A\b % 用矩阵左除运算求得方程组特解x0 9、orth把向量组正交规范化 ， ， 例8：设向量组 求由这三个向量生成的子空间V的一个标准正交基。 a1=[1;2;3];a2=[-1;1;2];a3=[5;1;0]; A=[a1,a2,a3]; P=orth(A) 10、eig计算矩阵的特征值和特征向量 例9：计算 的特征值和特征向量 A=[1 2 3; 2 1 3; 1 1 2]; [ V, D ] = eig(A) 11、eig把二次型标准化 例10： 用正交变换法将二次型，化为标准形 A=[1,0,0;0,2,2;0,2,2]; [V,D]=eig(A); % 其中矩阵V即为所求正交矩阵 12、eig判断二次型的正定性 例11： 判断矩阵的正定性 A=[1,1,-1;1,2,-1;-1,-1,5]; eig(A) * 化学方程的配平 实 际 应 用 * 要使方程配平，x1 , x2 , x3 , x4 必须满足： 将所有项移到左端，并写成矩阵相乘的形式，就有： * 对矩阵A进行行阶梯变换 A=[3,0,1,0;8,0,0,2;0,2,2,1] U0?rref(A) 得到 注意四个列对应于四个变量的系数，有 * 此处可取 x4=4 ，则 x1 , x2 , x3 均有整数解。因而配平后的化学方程为： 应该在MATLAB中先规定取有理分式格式。即先键入format rat，结果为： 对于比较复杂的反应过程，为了便于得到最小整数的解，在解化学配平的线性方程组时， * 混凝土配料中的应用 超强型A 通用型B 长寿型C 水泥c 20 18 12 水w 10 10 10 沙s 20 25 15 石g 10 5 15 灰f 0 2 8 1. 假如某客户要求的混凝土的五种成分为16,10,21,9,4，试问A,B,C三种类型应各占多少比例？如果客户总共需要5000公斤混凝土，则三种类型各占多重？ * 2. 如果客户要求的成分为16,12,19,9,4，则这种材料能用A,B,C三种类型配成吗？为什么？ 解：从数学上看，这三种基本类型的混凝土相当于三个基向量 vA , vB , vC ，待配混凝土相当于两个合成向量 w1 , w2 ，其数值如下： * w1 , w2是否是向量 vA , vB , vC 的线性组合？ 或w1 , w2是否在vA , vB , vC 张成的向量空间内？ 解的思路是先分析三个基向量所张成的空间是几维?把 w1,w2加进去后维数是否没有增加？这就可以知道w1,w2是不是在vA，vB，vC张成的向量空间内。因为系统的阶数已达到5，手工解太费时了，必须藉助于计算机，程序ea581如下： * * 为了计算 w1所需的混合比例，问题归结为解下列线性代数方程： 为此，可对增广矩阵[v,w1]作行阶梯变换，求U0=rref([v,w1])，得到 * 所以三种基本混凝土应按8%，56%，36%的比例调配，对于5000公斤混凝土，三种基本类型混凝土的用量分别为800，2800，1400公斤。 例6.1 在钢板热传导的研究中，常常用节点温度来描述钢板温度的分布。假设图6.1中钢板已经达到稳态温度分布，上下、左右四个边界的温度值如图所示，而表示钢板内部四个节点的温度。若忽略垂直于该截面方向的热交换，那么内部某节点的温度值可以近似地等于与它相邻四个节点温度的算术平均值，如 。请计算该钢板的温度分布。 解：根据已知条件可以得到以下线性方程组： 化简为标准的矩阵形式如下： ℃， ℃， ℃， ℃。 得到方程组的解为： 西安电子科技大学 电子工程学院 工程线性代数（MATLAB） 教材：《 线性代数机算与应用 》 西安电子科技大学 MATLAB让您站在巨人的肩膀上, 高瞻远瞩 ! MATLAB软件简介 MathWorks公司 MATLAB?是MathWorks公司于1982年推出的一套高性能的数值计算和可视化软件, 它集数值分析、矩阵运算、信号处理和图形显示于一体, 构成了一个方便的, 界面友好的用户环境。 在这个环境下, 对所要求解的问题, 用户只需简单地列