乐视网估值仅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)