[])matlab的数值运算.pptVIP

* x1+2x2+3x3=1 2x1+3x2+4x3=2 a=[1 2 3;2 3 4];b=[1;2]; x=a\b x=pinv(a)?b x = x = 1.00 0.83 0 0.33 0 -0.17 = a x = b * 七、数据分析与插值函数 max —— 各列最大值 mean —— 各列平均值 sum —— 各列求和 std —— 各列标准差 var —— 各列方差 sort —— 各列递增排序 * 八、拟合与插值 1. 多项式拟合 x0=0:0.1:1； y0=[-.447 1.978 3.11 5.25 5.02 4.66 4.01 4.58 3.45 5.35 9.22]; p=polyfit(x0,y0,3) p = 56.6915 -87.1174 40.0070 -0.9043 xx=0:0.01:1;yy=polyval(p,xx); plot(xx,yy,-b,x0,y0,or) * 2.插值 插值的定义——是对某些集合给定的数据点之间函数的估值方法。 当不能很快地求出所需中间点的函数时，插值是一个非常有价值的工具。 Matlab提供了一维、二维、 三次样条等许多插值选择 * intep1 —— interp2 —— spline —— 利用已知点确定未知点 粗糙—— 精确 集合大的—— 简化的 插值函数 * 本节介绍了matlab语言的数值运算 功能，通过学习应该掌握： 如何创建矩阵、修改矩阵 符号的用法 矩阵及数组运算 多项式运算 线性方程组与微分运算 小 结 * Task 例1 数组与矩阵乘法的对比示例 P34 矩阵乘法 例2 中矩阵的除法与幂运算 P34第3节 例3 P39 矩阵分解，掌握多种分解方法，包括特征值分解、奇异值分解、LU分解、cholesky分解、QR分解。 例4 P61 基于数值运算的多项式创建、求值、求根乘除法、微积分 例5 代数方程组求解 分三种情况，用左除和求逆运算来求解x，判断解是否相同？并进行分析。 * 思考问题 为什么矩阵在数值运算中那么重要? 说明曲线拟合的最小二乘法 * 由线性代数我们知道A非奇异时,A的行列式不为0,此时方程的解是唯一的。 在实际应用中,除法解方程的速度要比求逆法快2.5倍精确度更高,明显优于求逆法,所以推荐尽量使用除运算,少用逆运算. * * a^0.5 ans = 0.4498 + 0.7623i 0.5526 + 0.2068i 0.6555 -0.3487i 1.0185 + 0.0842i 1.2515 + 0.0228i 1.4844 - 0.0385i 1.5873 - 0.5940i 1.9503 - 0.1611i 2.3134 + 0.2717i * inv —— 矩阵求逆 det —— 行列式的值 eig —— 矩阵的特征值 diag —— 对角矩阵 ’ —— 矩阵转置 sqrt —— 矩阵开方 5. 矩阵的其它运算 * 6. 矩阵的一些特殊操作 矩阵的变维 a=[1:12];b=reshape(a,3,4) c=zeros(3,4);c(:)=a(:) 矩阵的变向 rot90:旋转; fliplr:左右翻; flipud:上下翻 矩阵的抽取 diag:抽取主对角线;tril: 抽取主下三角; triu:抽取主上三角 矩阵的扩展 * 关系运算 关系符号 意义 = = == ~= 小于 小于或等于 大于 大于或等于 等于 不等于 * 逻辑运算 逻辑符号 意义 ＆ ｜ ～ 与（AND） 或（OR） 非（NOT） * 关系函数和逻辑函数 函数名 功能 函数名 功能 all 是否为全1矩阵 isinf 是否无穷大 any 找非零元素 isnan 是否非值 exist 存在性与类别 issparse 是否稀疏 find 找非零元素 isstr 是否字串 isempty 是否为空 isglobal 是否全局 isfinite 是否有限 xor(x,y) 异或运算 * 数组运算指元素对元素的算术运算， 与通