第五章_用Matlab进行现代科学运算.pptVIP

第五章 用Matlab进行现代科学计算 5.1 解析解与数值解 在实际的工程技术中，一般通过数值解法来获得问题的解。至少有两个原因： 1、解析解不存在 如圆周率的值就没有解析解，超越积分也没有解析解。在这种情况下，就必须采用数值解技术。 2、解析解存在但不实用 如求解n元一次代数方程组的问题，从理论上讲，总可以把多元一次方程组简化成解析可解的形式。 然而当n较大时，需要的基本运算次数是非常惊人的，对于求解实际问题所需要的计算成本是根本无法接受的，只能采用数值解法。 5.2 数值线性代数问题及求解 5.2.1 特殊矩阵的MATLAB输入 零矩阵： zeros(m,n); (m,n) 定义矩阵维数 幺矩阵： ones(m,n); 单位矩阵： eye(m,n); 随机元素矩阵： rand(m,n); [0,1]上均匀分布 randn(m,n); 正态分布 对角阵： diag(v); v为对角元素组成的向量 5.2.2 矩阵的特征参数运算 矩阵的行列式： det(A) 矩阵的迹： trace(A)，即对角线元素之和 矩阵的秩： rank(A)，即线性无关的列数或行数 矩阵的范数 norm(A,p) 1 1-范数，A的最大列和，m