方程求根与解常微分方程.PDF

第七章 方程求根与解常微分方程 7.1实验目的 了解微分方程的通解、特解和近似解的概念。熟悉方程求根和常微分方程解的概念，熟悉Mathematica软 件的方程求根和求常微分方程解的命令，掌握用数学软件处理方程求根和常微分方程解的有关问题． 7.2实验准备 7.2.1数学概念 1． 微分方程 2． 微分方程的通解、特解 7.2.2数学软件命令 1. Solve[eqn, x] 功能：求多项式方程eqn的所有根，当多项式方程的次数n≤4时，给出eqn所有根的准确形式, 当n4时， 不一定能求出所有的根, 此时，命令输出形式为 {ToRules[Roots[eqn, x ]]} n次多项式方程的一般形式为: 2 n a +ax+ ax +…+a x = 0 0 1 2 n 式中a ,a, a,…,a 为常数。 0 1 2 n 2. Solve[{eqn1, eqn2, …, eqnk}, {x1, x2,…, xk}] 功能：求多项式方程组{eqn1, eqn2, …, eqnk}的所有根, 当其中每个多项式方程的次数n≤4时, 给出所 有根的准确形式, 否则，不一定能求出所有的根, 此时，命令输出形式为{ToRules[Roots[{eqn1, eqn2, …, eqnk}, {x1, x2,…, xk} ]]} 。 3. NSolve[eqn, x] 功能：求多项式方程eqn的所有根的近似形式。 4. NSolve[{eqn1, eqn2, …, eqnk}, {x1, x2,…, xk}] 功能：求多项式方程组{eqn1, eqn2, …, eqnk}所有根的近似形式。 5. FindRoot[eqn, {x, x0}] 功能：求方程eqn的在初值x0附近的一个近似根。 6. FindRoot[{eqn1,eqn2, ... }, {x, x0}, {y, y0}, ... ] 功能：求方程组{eqn1, eqn2, …}在初值(x0,y0,…)附近的一个近似根。 7.DSolve[eqn，y[x]，x] 功能：求常微分方程eqn的通解y=f(x,C). 8.DSolve[{eqn，初始条件}，y[x]，x] 功能：求常微分方程eqn满足初始条件的特解. 9.DSolve[{eqn1,eqn2, ... }, {y1[x],y2[x], ... }, x] 功能：求常微分方程组{eqn1,eqn2, ... }的解. 10.NDSolve[ eqns, y, {x, xmin, xmax}] 功能：求出自变量范围为[xmin, xmax]且满足给定常微分方程及初值条件eqns的未知函数y的数值解. 得 到的解是以 {{未知函数名-InterpolatingFunction[range, ]}} 的形式给出的,其中的InterpolatingFunction[range, ]是所求的插值函数表示的数值解, range就是 所求数值解的自变量范围。 11.NDSolve[ eqns, {y1,y2, ... }, {x, xmin, xmax}] 功能：求出自变量范围为[xmin, xmax]且满足给定常微分方程及初值条件eqns的未知函数y1,y2, ... 的 数值解。 7.3实验任务 7.3.1基础实验： 熟悉数学软件命令操作。 1. 求下列代数方程的根 3 2 1）设 x -4x +9x - 10 = 0 ，求方程的所有根。 2 2）设x -ax - 4b=0，求方程的所有根，a，b 为常数。 ììììx+3y=0 x+3y=0 ïïïïïïxx++33yy==00 ïï 3）求方程组 íííí 2 2 的所有根。 2 2 ïïïï 22 22