Matlab和化工数值计算-第3讲矩阵操作和线性方程组求解.pptVIP

第三讲 矩阵操作与线性方程组求解 隋志军 化工学院软件应用教科组 2006-10 本章知识要点 本章的所要解决的典型问题 病态线性方程组 线性方程组的求解方法 迭代解法： 雅可比（Jacobian）迭代 高斯-赛德尔（Gauss- Siedel）迭代 松弛（SOR）迭代 共轭梯度（CG）迭代 在Matlab线性方程组的解法归结为矩阵的除法命令 矩阵的生成 一维数组的寻访和赋值 逻辑1标识 矩阵的特殊操作 矩阵的基本性质 左除和右除的比较 用以下程序生成条件数很大的方程 randn(state,0); A=gallery(randsvd,100,2e13,2);%产生条件数为2e13的100阶矩阵 x=ones(100,1); %生成真解 b=A*x; %用A和x生成b 比较左除法和求逆法的求解所用的时间和相对残差 本章开始问题的求解 线性方程组的迭代解法 * 数值计算 线性方程组的直接解法 线性方程组的迭代解法 MATLAB 矩阵的生成与操作 矩阵运算函数 数学模型： 线性方程组在化工模型中的作用 多组分体系的物料衡算 各种化合物的物理化学性质 稳态动力学计算 微分方程的差分方程求解 x=1,y=1 x=10,y=-2 x=2,y=-2 x=0.991,y=-0.487 直接求解方法： 高