12、MATLAB程序设计与应用刘卫国（第三版）课后实验十二：级数与方程符号求解

目录

 一、

 二、

 三、

 四、

 五、

---

一、

级数符号求和。

  ---------------------------------------示例代码---------------------------------------------

syms n x;
s1=symsum(1/(2*n-1),n,1,10)
s2=symsum(n*2+x^(n-1),n,1,inf)
s3=symsum(n^2/5^n,n,1,inf)

---------------------------------------运行结果---------------------------------------------

 二、

将lnx在x=1处按5次多项式展开为泰勒级数。

 ---------------------------------------示例代码---------------------------------------------

x=sym('x');
f=log(x);
taylor(f,x,1,'Order',6)

---------------------------------------运行结果---------------------------------------------

 三、

求下列方程的符号解。

  ---------------------------------------示例代码---------------------------------------------

syms x y;
x1=solve(log(x+1)-5/(1+sin(x))==2)
x2=solve(x^2+9*sqrt(x+1)-1)
x3=vpasolve(3*x*exp(x)+5*sin(x)-78.5)
[x4,y4]=solve([sqrt(x^2+y^2)-100',3*x+5*y-8])

---------------------------------------运行结果---------------------------------------------

 四、

求微分方程初值问题的符号解,并与数值解进行比较

---------------------------------------示例代码---------------------------------------------

dsolve('D2y+4*Dy+29*y','y(0)=0','Dy(0)=15','x')

function yy=test12_4_2(x,y)
yy=[-4*y(2)-29*y(1);y(2)];

[x,y]=ode23(@test12_4_2,[1.0e-9,20],[0;15])

 五、

求微分方程组的通解。

 ---------------------------------------示例代码---------------------------------------------

[x,y,z]=dsolve('Dx-2*x+3*y-3*z','Dy-4*x+5*y-3*z','Dz-4*x+4*y-2*z','t')

---------------------------------------运行结果---------------------------------------------