南邮matlab上机实验二,南邮MATLAB上机实验二.doc

南邮MATLAB上机实验二

实验二 MATLAB的符号运算MATLAB中，符号计算虽以数值计算的补充身份出现，但涉及符号计算的指令使用、运算符操作、计算结果可视化、程序编制以及在线帮助系统都是十分完整、便捷的。

一、实验目的

1、熟悉符号对象和符号表达式；

2、掌握符号表达式和符号函数的操作；

3、学习利用符号代数方程的求解 。

二、实验基本知识

1、基本运算指令 ；

2. 掌握符号运算的应用。

三、上机练习

1. (1)通过符号计算求的导数。(2)然后根据此结果，求和。

>> clear all

syms t

y=abs(sin(t))

df=diff(y)

df1=limit(df,t,0,'left')

df2=subs(df,'t',sym(pi/2))

y =

abs(sin(t))

df =

sign(sin(t))*cos(t)

df1 =

-1

df2 =

0

2.计算二重积分。

>>clear all

syms x y

f=x^2+y^2

Fint=(int(int(f,y,1,x^2),x,1,2))

f =

x^2 + y^2

Fint =

1006/105

3. 已知

(1)求;(2)求.

>> clear all

syms x y

f=x^2+y^2;

dfdy2=diff(f,y,2)

dfdy2 =

2

>> intf=int(f,x,-1,1)

intf =

2*y^2 + 2/3

4.解方程:

>> clear all

syms x

eq=solve(3*x^2+4*x-6)

eq =

-2/3+1/3*22^(1/2)

-2/3-1/3*22^(1/2)

5.已知方程:

求通解;(2)求满足的特解.

>>clear all

>> y=dsolve('D3y+3*Dy+2*y-6=0','x')

y =

3+C1*exp((1+2^(1/2))^(1/3)*(-1-(1+2^(1/2))^(1/3)+(1+2^(1/2))^(1/3)*2^(1/2))*x)+C2*exp(-1/2*(1+2^(1/2))^(1/3)*(-1-(1+2^(1/2))^(1/3)+(1+2^(1/2))^(1/3)*2^(1/2))*x)*sin(1/2*3^(1/2)*(1+2^(1/2))^(1/3)*(1-(1+2^(1/2))^(1/3)+(1+2^(1/2))^(1/3)*2^(1/2))*x)+C3*exp(-1/2*(1+2^(1/2))^(1/3)*(-1-(1+2^(1/2))^(1/3)+(1+2^(1/2))^(1/3)*2^(1/2))*x)*cos(1/2*3^(1/2)*(1+2^(1/2))^(1/3)*(1-(1+2^(1/2))^(1/3)+(1+2^(1/2))^(1/3)*2^(1/2))*x)

>> y=dsolve('D3y+3*Dy+2*y-6=0','y(0)=3,D2y(0)=1,Dy(0)=2','x')

y =

3+(-1/12*(1+2^(1/2))^(1/3)*2^(1/2)+1/6*(1+2^(1/2))^(2/3)+1/3*(1+2^(1/2))^(1/3)-1/4*(1+2^(1/2))^(2/3)*2^(1/2))*exp((1+2^(1/2))^(1/3)*(-1-(1+2^(1/2))^(1/3)+(1+2^(1/2))^(1/3)*2^(1/2))*x)+(1/3*3^(1/2)*(1+2^(1/2))^(1/3)-1/12*3^(1/2)*(1+2^(1/2))^(1/3)*2^(1/2)-1/6*3^(1/2)*(1+2^(1/2))^(2/3)+1/4*3^(1/2)*(1+2^(1/2))^(2/3)*2^(1/2))*exp(-1/2*(1+2^(1/2))^(1/3)*(-1-(1+2^(1/2))^(1/3)+(1+2^(1/2))^(1/3)*2^(1/2))*x)*sin(1/2*3^(1/2)*(1+2^(1/2))^(1/3)*(1-(1+2^(1/2))^(1/3)+(1+2^(1/2))^(1/3)*2^(1/2))*x)+(1/12*(1+2^(1/2))^(1/3)*2^(1/2)-1/6*(1+2^(1/2))^(2/3)-1/3*(1+2^(1/2))^(1/3)+1/4*(1+2^(1/2))^(2/3)*2^(1/2))*exp(-1/2*(1+2^(1/2))^(1/3)*(-1-(1+2^(1/2))^(1/3)+(1+2^(1/2))^(1/3)*2^(1/2))*x)*cos(1/2*3^(1/2)*(1+2^(1/2))^(1/3)*(1-(1+2^(1/2