第四节数值计算.pptVIP

目前三十页\总数五十一页\编于八点 目前三十一页\总数五十一页\编于八点 目前三十二页\总数五十一页\编于八点 §4.5 矩阵运算 det(A) A的行列式 eig(A) A的特征值 norm(A,1) A的1范数 norm(A) A的2范数 norm(A,inf) A的无穷范数 norm(A,‘fro’) A的F范数 rank(A) A的秩 inv(A) A的逆矩阵 trace(A) A的迹：矩阵对角线元素之和 目前三十三页\总数五十一页\编于八点 线性方程组的解法 它可以表示为 ,其中 线性方程组 目前三十四页\总数五十一页\编于八点 利用左除运算符left division operation x=A\B 若方程是xC=d形式，就用“右除”，即x=d/C。 采用逆阵法x=inv(A)*B 但该法只适用于当A非奇异、维数较低且条件 较好时，使用范围有限。 上述方程组的解为： 目前三十五页\总数五十一页\编于八点 例：分别采用逆阵法和左除法求解。 已知 ，求 ，使 解1：逆阵法 A=[1 0 1; 2 1 0; -3 2 -5]; b=[1 2 -1]’; X=inv(A)*b 解2：直接法 tic X=A\b toc 所以，左除法速度快，而且精度高。 Elapsed time is 36.672815 seconds. Elapsed time is 3.346355 seconds. 目前三十六页\总数五十一页\编于八点 第四节数值计算演示文稿 目前一页\总数五十一页\编于八点 （优选）第四节数值计算 目前二页\总数五十一页\编于八点 本节内容 掌握统计分析中常用函数的应用 掌握复数运算中常用函数的格式和使用 掌握多项式运算和矩阵运算 掌握曲线拟合和数据插值的应用 了解matlab在微积分方面的应用 目前三页\总数五十一页\编于八点 §4.1 复数运算 1.复数的表示形式 其中：c是复数，a为实部，b为虚部，i是 虚数单位 其中：z是矢量长度， 是弧度 另一种表示形式： 目前四页\总数五十一页\编于八点 两种表示形式的换算关系 实轴 虚轴 目前五页\总数五十一页\编于八点 常用复数计算函数 conj(x)：求复数x的共轭数 real： 求复数的实部 imag：求复数的虚部 abs：求复数的模值 angle：求复数的相位角 isreal：判断变量是否为实数，是则返回1 目前六页\总数五十一页\编于八点 §4.2 数据统计分析 求最大值max和最小值min 可查阅datafun函数库中相应的函数。 doc datafun 调用格式如下： y=max(x)：返回向量x的最大值并存入y； 若x中包含复数元素，则按模取最大值 [y,I ]=max(x)：返回向量x的最大值并存入y； 最大值的序号存入I 目前七页\总数五十一页\编于八点 Y=max(A)：返回矩阵A的每一列上的最大值 [Y,U ]=max(A)：返回值Y和U均为行向量，Y表 示A的每列的最大值，U表示每列最大值的行号 Y=max(A,[ ],dim)：dim取1时，表示求每一列 上的最大值；dim取2时，表示求矩阵A每一行 上的最大值 求最小值的函数是min，其用法和max完全相同 目前八页\总数五十一页\编于八点 例：求3阶魔方矩阵中各列和各行元素的最大值， 并求整个矩阵的最大值和最小值。 程序如下： x=magic(3) max(x) %求矩阵x中各列元素中的最大值 max(x,[ ],2) %求x中各行元素中的最大值 max(max(x)) %求x的最大值9 min(min(x)) %求矩阵x的最小值1 x = 8 1 6 3 5 7 4 9 2 ans= 8 9 7 ans= 8 7 9 目前九页\总数五十一页\编于八点 数据统计常用函数 max 求矩阵的最大元素 min 求最小元素 mean 求平均值 median 求中值 sort 各列元素按升序排列 prod 各元素的积 sortrows 各行元素按升排列 std 求矩阵的标准方差 sum 求矩阵各列元素的和 目前十页\总数五十一页\编于八点 例:分别求该数组x中各列和各行元素的和。 x=[4 5 6; 1 4 8] x = 4 5 6