第七讲 常微分方程及级数.docVIP

第七讲 实验目的： １．学习用matlab求解微分方程命令dsolve. ２．学习matlab泰勒级数展开命令. ３．巩固幂级数的收敛半径、和等概念. 实验内容： １．学习matlab命令. matlab求解微分方程命令dsolve，调用格式为： dsolve(‘微分方程’)给出微分方程的解析解，表示为t的函数. dsolve(‘微分方程’，‘初始条件’)给出微分方程初值问题的解，表示为t的函数. dsolve(‘微分方程’，‘变量x’)给出微分方程的解析解，表示为x的函数. dsolve(‘微分方程’，‘初始条件’，‘变量x’)给出微分方程初值问题的解，表示为x的函数. 求已知函数的taylor展开式taylor命令，调用格式为： taylor(函数f(x)) f(x)的5次taylor多项式. taylor(函数f(x)，n) f(x)的n-1次taylor多项式. taylor(函数f(x)，a) f(x)在a点的taylor多项式. 求级数的和命令symsum调用格式为： symsum(S,n)，求 symsum(S,k,m,n),求 matlab求极限命令limit调用格式为： limit(函数f(x)，变量x，自变量的趋向值) ２．求解一阶微分方程. 微分方程在输入时，应输入Dy，应输入D2y等，D应大写． 　　例１：求微分方程的通解． 　　解：输入命令： dsolve(Dy+2*x*y=x*exp(-x^2)) 结果为 ans = 1/2*(1+2*exp(-2*x*t)*C1*exp(x^2))/exp(x^2) 系统默认的自变量是t，显然系统把x当作常数，把y当作t的函数求解．输入命令： dsolve(Dy+2*x*y=x*exp(-x^2),x) 得正确结果： ans = 1/2*(x^2+2*C1)/exp(x^2) 　　例２：求微分方程在初始条件下的特解． 　　解：输入命令： dsolve(x*Dy+y-exp(x)=0,y(1)=2*exp(1),x) 得结果为： ans = 1/x*(exp(x)+exp(1)) 　　例３：求微分方程在初始条件下的特解． 　　解：输入命令： dsolve((x^2-1)*Dy+2*x*y-cos(x)=0,y(0)=1,x) 得结果为 ans = 1/(x^2-1)*(sin(x)-1) ３．求解二阶微分方程． 　　例４：求的通解． 　　解：输入命令： dsolve(D2y+3*Dy+exp(x)=0,x) 得结果： ans = -1/4*exp(x)+C1+C2*exp(-3*x) 　　例５：求解微分方程． 　　解：输入命令： dsolve(D2y-exp(2*y)*Dy=0,x) 得结果： ans = 1/2*log(-2*C1/(-1+exp(2*x*C1+2*C2*C1)))+x*C1+C2*C1 ４．taylor展开式． 　　例６：求函数y=cosx在x=0点处的5阶taylor展开式及在处的6阶taylor展开式． 　　解：输入命令： syms x; taylor(cos(x)) 得结果： ans = 1-1/2*x^2+1/24*x^4 输入命令： taylor(cos(x),pi/3,7) 得结果： ans = 1/2-1/2*3^(1/2)*(x-1/3*pi)-1/4*(x-1/3*pi)^2+1/12*3^(1/2)*(x- 1/3*pi)^3+1/48*(x-1/3*pi)^4-1/240*3^(1/2)*(x-1/3*pi)^5-1/1440*(x-1/3*pi)^6 ５．级数求和． 　　例７：求． 　　解：输入命令： syms n;symsum(1/2^n,1,inf) 得结果： ans = 1 　　例８：求幂级数的和函数． 　　解：输入命令： symsum(x^n/(n*2^n),n,1,inf) 得结果 ans = -log(1-1/2*x) 　　例９：求幂级数的和