第三章MATLAB的数值计算资料.ppt

第三章 MATLAB的数值计算 第三章 MATLAB的数值计算 §3-1-3利用M文件建立矩阵 §3-1-4矩阵元素的求和 §3-1-5 矩阵的下标 §3-1-6 矩阵的转置 §3-1-7 矩阵的对角元素 §3-1-8 矩阵的最大值 §3-2 矩阵的基本运算§3-2-1 矩阵的加减运算 §3-2-2 矩阵的乘积 §3-2-3 矩阵的冒号运算符 §3-3 矩阵的乘方和指数运算 §3-3-2 矩阵的指数运算 §3-4 矩阵的除法 §3-5 行列式及矩阵的逆 §3-5-1矩阵的秩 §3-6 特征值分解 §3-7求解代数方程组 §3-8 多项式§3-8-1多项式的表示和处理多项式的函数 §3-8-2 多项式的根和系数 §3-8-3 多项式的值和特征多项式 §3-8-4 多项式的乘法和除法 §3-8-5 多项式的导数 §3-9 多项式的计算 §3-10 多项式的部分分式 §3-10-2 函数residue还可以用来对多项式的部分分式作逆运算 作业：求下列线性方程组的解A*X=b。用三种方法实现 A=[2,7,3,1;3,5,2,2;9,4,1,7];b=[6,4,2]; 2求解非线性方程组用solve 函数 solve 函数用于求解非线性代数方程的根。 调用格式： g=solve(eq) 缺省情形下，采用系统默认的未知数。 g=solve(eq1,eq2,……eqn) %解方程组 例： syms a b c x y u v %用syms表示引出的为符号 [x,y]=solve(x^2+x*y+y=3,x^2-4*x+3=0)%表示解的是 二元一次方程组 x = [ 1] [ 3] y = [ 1] [ -3/2] solve(‘a*x^2+b*x+c=0’) % 即solve(‘a*x^2+b*x+c’)，求的是 %x的 值， ans = [ 1/2/a*(-b+(b^2-4*a*c)^(1/2))] [ 1/2/a*(-b-(b^2-4*a*c)^(1/2))] 在上面的solve函数中，没有指定未知数，由系统默认确定为x。 我们还可以自己来指定函数中的未知参数，如下例中指定b为 未知数。 solve(a*x^2+b*x+c,b) %后面的b表示指定b为未知数 ans = -(a*x^2+c)/x 在上面的方程中都确认为多项式的值为0，当然也可以给多项式附其他值，如： solve(p*sin(x)=r) ans = asin(r/p) %asin表示反正弦 例： 解方程。 命令如下： x=solve(1/(x+2)+4*x/(x^2-4)=1+2/(x-2),x) %解方程(1) f=sym(x-(x^3-4*x-7)^(1/3)=1); x=solve(f) %解方程(2)，等于solve(x-(x^3-4*x-7)^(1/3)=1); x=solve(2*sin(3*x-pi/4)=1) %解方程(3) x=solve(x+x*exp(x)-10,x) %解方程(4)。仅标出方程的左端 x = 1 x = 3 x = 5/36*pi x = 1.6335061701558463841931651789789 例：求方程组的解。 命令如下： [x y]=solve(1/x^3+1/y^3=28,1/x+1/y=4,x,y) %解方程组(1) x = [ 1] [ 1/3] y = [ 1/3] [ 1] [x v]=solve(x^2+y^2-5,2*x^2-3*x*y-2*y^2) %解方程组(4) x = [ -1] [ 1] [ 2] [ -2] v = [ 2] [ -2] [ 1] [ -1] 对于多项式：p(x)=a0xn+a1xn-1+a2xn-2+……an-1x+an 可以用行向量p=[a0 a1 a2………an-1 an]来表示。 例如：多项式p(x)=x3+5x2-6 可以用p=[1,5,0,-6]来表示。 注意：对于系数为0的项，必须用0来填充。 行向量中各元素按照多项式项的次数从高到低排列。 在MATLAB中，是使用行向量来表示多项式的系数。 MATLAB提供了一些处理多项式的基本函数，如表3-1所示： 表3-1 MATLAB中的多项式函数