非线性方程求解..ppt

数值计算的误差 算法设计的原则 非线性方程求解 引言 引言 在945.36kPa（9.33atm）、300.2K时，容器中充以2mol氮气，试求容器体积。已知此状态下氮气的P-V-T关系符合范德华方程，其范德华常数为a＝4.17atm?L/mol2， b＝0.0371L/mol 数学模型：范德华方程变形可得关于V的非线性方程 非线性方程 非线性方程包括：高次代数方程、超越方程(具有未知量的对数函数、指数函数、三角函数、反三角函数等的方程)及其它们的组合 与线性方程相比，非线性方程求解问题无论从理论上还是从计算公式上都要复杂得多 对于高次代数方程，当次数4时，则没有通解公式可用，对于超越方程既不知有几个根，也没有同样的求解方式。实际上，对于n≥3代数方程以及超越方程都采用数值方法求近似根。 MATLAB求解非线性方程函数 多项式求根函数roots 多项式的表达式约定如下： 对于多项式，用以下行向量表示： 这样就把多项式问题转化为向量问题 函数roots r = roots(P)，用于求解多项式的根 其中，行向量P的元素是多项式的系数，按多项式次数降序排列 如果P中含有n+1个元素，则多项式为n次 roots可以获得多项式的所有根 例题： 求方程 的根 例题： 在945.36kPa（9.33atm）、300.2K时，容器中充以2mol氮气，试求容器体积。已知此状态下氮气的P-V-T关系符合范德华方程，其范德华常数为a＝4.17atm?L/mol2，b＝0.0371L/mol。 数学模型：范德华方程变形可得 这是关于V的三次方程，可以由roots 非线性方程求解函数fzero fzero 对于一般的含单个未知数的超越方程，可以采用fzero函数求解 fzero函数结合使用二分法、割线法和可逆二次内插法 函数fzero [x,fval,exitflag,output] = fzero(fun,x0,options, p1, p2, ...) 此函数的作用求函数fun在x0附件的零值点，x0是标量 x 所求解 fval 函数在解x处的值 exitflag 程序结束情况: 0，程序收敛于解；0，程序没有收敛；＝0，计算达到了最大次数 output 是一个结构体，提供程序运行的信息；output.iterations，迭代次数；output.functions，函数fun的计算次数；output.algorithm，使用的算法 options 选项，可用optimset函数设定选项的新值 fun可以是函数句柄或匿名函数。 P1,P2是传递的参数 说明 1) 采用句柄函数定义函数 function Cha2demo1 x=fzero(@fun,1) function y=fun(x) y=x^3-2*sin(x); 2) 初值的选择对于解有影响，不同的初值可能获得不同的解 可以根据感兴趣的解的区间确定初值范围 可以作出函数在一定范围内的曲线，直观的确定解的大致范围 3) fzero不能获得多项式的多重根，尤其是复数根。而roots函数求解，则可获得所有根 例题7： 计算以下方程的根 1) 求sinx在3附近的零点; 2) 求cosx在[1,2]范围内的零点； 3) 4) 例题： 在945.36kPa（9.33atm）、300.2K时，容器中充以2mol氮气，试求容器体积。已知此状态下氮气的P-V-T关系符合范德华方程，其范德华常数为a＝4.17atm?L/mol2，b＝0.0371L/mol。 数学模型：范德华方程为： function Cha2demo3 P = 9.33; % atm T = 300.2; % K n = 2; % mol R = 0.08206; V0 = n*R*T/P % 5.2807 [V,fval] = fzero(@PVTeq,V0) % ------------------------------------------- function f = PVTeq(V) P = 9.33; % atm T = 300.2; % K n = 2; % mol a = 4.17; b = 0.0371; R = 0.08206; f = (P + a*n^2/V^2) * (V-n*b) - n*R*T; 函数fsolve 与fzero函数只能求解单个方程的根不同，fsolve函数可求解非线性方程组的解。其算法采用的是最小二乘法。 调用格式： [x,fval,exitflag,output,jacobian]