2.1)MATLAB及数值运算.pptVIP

* 5.多项式微积分 matlab提供了polyder函数多项式的微分。 命令格式： polyder(p): 求p的微分 例：a=[1 2 3 4 5]; poly2str(a,x) ans = x^4 + 2 x^3 + 3 x^2 + 4 x + 5 b=polyder(a) b = 4 6 6 4 poly2str(b,x) ans =4 x^3 + 6 x^2 + 6 x + 4 * polyint求多项式函数的不定积分： 命令格式： p=polyint(a): 求a的不定积分，常数项为0 例：a=[1 2 3 4 5]; poly2str(a,x) ans = x^4 + 2 x^3 + 3 x^2 + 4 x + 5 b=polyint(a) b = 0.2000 0.5000 1.0000 2.0000 5.0000 0 poly2str(b,x) ans = 0.2 x^5 + 0.5 x^4 + x^3 + 2 x^2 + 5 x * 七、代数方程组求解 matlab中有两种除运算左除和右除。 对于方程ax=b，a 为an×m矩阵，有三种情 况： ? 当n=m时，此方程成为“恰定”方程 ? 当nm时，此方程成为“超定”方程 ? 当nm时，此方程成为“欠定”方程 matlab定义的除运算可以很方便地解上 述三种方程 * 1.恰定方程组的解 方程ax=b(a为非奇异) x=a-1 b 矩阵逆 两种解: x=inv(a)?b — 采用求逆运算解方程 x=a\b — 采用左除运算解方程 * 方程ax=b a=[1 2;2 3];b=[8;13]; ?x=inv(a)*b ? x=a\b x = x = 2.00 2.00 3.00 3.00 = a x = b 例: x1+2x2=8 2x1+3x2=13 * 2.超定方程组的解 方程 ax=b ,mn时此时不存在唯一解。 方程解 (a a)x=a b x=(a a)-1 a b —— 求逆法 x=a\b —— matlab用最小二乘法找一 个准确地基本解。 定理 :当RT·R可逆时，以上超定方程组存在最小二乘解. * 例: x1+2x2=1 2x1+3x2=2 3x1+4x2=3 a=[1 2;2 3;3 4];b=[1;2;3]; 解1 x=a\b 解2 x=inv(a?a) ? a ? b x = x = 1.00 1.00 0 0.00 = a x = b * 3.欠定方程组的解 当方程数少于未知量个数时,即不定 情况,有无穷多个解存在。 matlab可求出两个解： 用除法求的解x是具有最多零元素的解 是具有最小长度或范数的解，这个解是基于伪逆pinv求得的。 * 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