MATLAB的数值计算(二)线性方程求解
一、代数方程及方程组的求解
MATLAB提供了求解方程的一些函数,包括代数方程和微分方程。
利用solve函数命令求解代数方程,其调用格式为:
[x1,x2,...,xn]=solve('eqn1','eqn2','eqn')
表示对n个未知变量x1,x2,...,xn求解n个方程。
例:一元二次方程,求x的解。
例:求代数方程组 的解。
二、微分方程及方程组的求解
利用dsolve函数命令求解微分方程,其调用格式为:
[x1,x2,...,xn]=dsolve('eqn1','eqn2',...,'eqn')
其中,输入量eqn可以是微分方程,也可以是求解微分方程的初始条件,还可以是独立变量。一般前面的是微分方程,中间的是初始条件,最后的是独立变量生命。
在表达微分方程时,n阶导数表示为Dny,例如,
,
,
,分别表示为Dy,Dx,D2y,D3y。
例:求解。
例:求解一阶微分方程组的解,初始条件为x(0)=0,y(0)=1。