乐视网估值仅3.91元:3.3 方程求解

来源:百度文库 编辑:中财网 时间:2024/07/04 18:31:59
 

3.3.1代数方程

solve(f)             解符号方程式f

solve(f1,…,fn)      解由f1,…,fn组成的代数方程组

例:

>>solve(f)     %求一元二次方程的根

ans =

[ 1/2/a*(-b+(b^2-4*c*a)^(1/2))]

[ 1/2/a*(-b-(b^2-4*c*a)^(1/2))]

>>solve(f,a)     %指定要求解的变量是a

ans =

-(b*x+c)/x^2

 

>>eq1 = 'x+y+z=10';

>>eq2 = 'x-y+z=0';

>>eq3 = '2*x-y-z=-4';

>>[x,y,z]=solve(eq1,eq2,eq3)     % 解三个联立方程式

x =

2

y =

5

z =

3 

3.3.2 常微分方程

    dsolve('equation','condition')

    其中equation代表常微分方程式,condition为初始条件,如果初始条件没有给出,则给出通解形式。在equation中,用字母D来表示求微分,其后的数字表示几重微分,默认所有这些变量都是对自变量t求导。

如:

>>dsolve( 'Dy=5' )             %求微分方程y’=5的通解

ans =

5*t+C1

>>dsolve( ' D2y=1+Dy ' )      %求微分方程y''=1+y'的通解

ans =

-t+C1+C2*exp(t)

>>[x,y]=dsolve(‘Dx=y+x,Dy=2*x’,’x(0)=0,y(0)=1’)     %求微分方程组的解,加初始条件

x =

1/3*exp(2*t)-1/3*exp(-t)

y =

2/3*exp(-t)+1/3*exp(2*t)