实验四MATLAB在方程求解和级数中的应用.pptVIP

实验四MATLAB在方程求解 和级数中的应用 ;一、利用MATLAB进行级数运算的方法和技能;1.常数项级数的求和与审敛;例1：求级数 ， 的和。 解：利用MATLAB函数symsum设计如下程序：;本例是收敛的情况，如果发散，则求得的和为inf，因此，本方法就可以同时用来解决求和问题和收敛性问题。 例2：求级数 ， 的和。 解：MATLAB程序如下：;运行结果为： I3 =1/6*sin(x)*pi^2 I4 =log(1+x) 从这个例子可以看出，symsum()这个函数不但可以处理常数项级数，也可以处理函数项级数。 ;2.函数的泰勒展开;在MATLAB中，用于幂级数展开的函数为taylor()，其具体格式为： taylor(function,n,x,a) function是待展开的函数表达式，n为展开项数，缺省时展开至5次幂，即6项，x是function中的变量，a为函数的展开点，缺省为0，即麦克劳林展开。;例3：将函数sin（x）展开为 x 的幂级数，分别展开至5次和20次。 解：MATLAB程序为： clear syms x f=sin(x); taylor(f) taylor(f,20);ans=x-1/6*x^3+1/120*x^5- 1/5040*x^7+1/362880*x^9-1x^11+1/6227020800*x^13-1/00*x^15+1/355687428096000*x^17-1/1216451*x^19 ;例4：将函数(1+x)m展开为x 的幂级数，为任意常数。展开至4次幂。 解：MATLAB程序为： clear syms x m f=(1+x)^m; taylor(f,5);例5：将函数 展开为 的幂级数。 解：MATLAB程序为： clear syms x f=1/(x^2+5*x-3); taylor(f,5,x,2) pretty(ans);结果为： ans =29/121-9/121*x+70/1331*(x-2)^2-531/14641*(x-2)^3+4009/161051*(x-2)^4 29 70 2 531 3 4009 4 ----- - 9/121 x + ------ (x - 2) - ------- (x - 2) + ---------- (x - 2) 121 1331 14641 161051;Nanjing University of Posts and Telecommunications;二、线性方程组、非线性方程、非线性方程组的求解。;函数null(A)则用来解线性方程组AX=O的基础解系，实际是求系数矩阵A的零空间，在null函数中可加入参数’r’，表示有理基。通过求系数矩阵的秩和增广矩阵的秩，可以判定方程组是否有解，以及是否需要求基础解系。另外，还可以用函数fzero来求解非线性方程。用法与fsolve类似。 ;例1：求解方程 的MATLAB程序为：X=solve(‘x^2-x-6=0’,’x’) 结果为：X=3, -2 例2：求解方程组 的程序为：; 例3：求解方程组 的程序为： clear format rat A=[5, 0, 4, 2;1, -1, 2, 1;4, 1, 2,0;1,1,1,1]； B=[3;1;1;0]； X=A\B; 例4：求方程组 的通解的程序为： clear format rat A=[1,2,2,1;2,1,-2,-2;1,-1,-4,-3] C=null(A,r) ％求出矩阵A的解空间的有理基。 结果如下： ;;求出的通解为： X = [ 2*k1+5/3*k2] [ -2*k1-4/3*k2] [ k1] [ k2]; 例5：求方程组 的通解的程序为： clear format rat A