第7章 非线性方程求解.pptVIP

第7章 非线性方程求解 在科学技术中，许多问题常归结为解一元函数方程f(x)=0根的数值问题。方程按f(x)是多项式或超越函数分别称为代数方程或超越方程。 7.1 求解非线性方程的MATLAB符号法 Matlab中设有求出方程f(x)=0解析解或精确解的符号命令solve，由他得出的符号量结果，可以转换成任意位有效数字的数值解。该命令的使用格式为 solve (s1, s2,…,sn, v1, v2,…, vn) solve (s1, s2,…,sn, v1,v2,…,vn) [z1,z2,…,zn]= solve (s1, s2,…,sn, v1, …, vn) 1）输入参量s1,s2,…,sn为待解方程组f(x)=0或函数f(x)的字符、符号表达式，或者代表他们的变量名。待解方程可以是任意线性、非线性或超越方程。 2）输入参量’v1’,’v2’,…,’vn’是与方程对应的未知量，它的数目必须与方程的数目相等；若有输出变量名z1,z2,…,zn，且与方程数相等，则输入变量‘v1’,’v2’,…,’vn’可以缺省。 3）输出参量z1,z2,…,zn是指定的输入变量名，方程解的结果分别赋值给它们。但是，赋值顺序并不是输入参量’v1’,’v2’,…,’vn的排序，而是按未知量名在字母表的排序输出，求解方程组时，这些输出参量不可省略，而且必须根方程数相等。否则只输出方程解的维数。 4）当方程组不存在解析解或精确解时，该指令输出方程的数字形式符号量解。 5）解析解表达式太冗长或含有不熟悉的特殊函数时，可用vpa命令转换成数值解。 例7.1求解方程组 解：在matlab命令窗口中输入 a=x^2+sqrt(5)*x=-1; b=x+3*z^2=4;c=y*z+1=0; [u v w]=solve(a,b,c) 回车得到 u = -1/2-1/2*5^(1/2) -1/2-1/2*5^(1/2) 1/2-1/2*5^(1/2) 1/2-1/2*5^(1/2) v = 1/76*(54+6*5^(1/2))^(1/2)*(-9+5^(1/2)) -1/76*(54+6*5^(1/2))^(1/2)*(-9+5^(1/2)) 1/44*(42+6*5^(1/2))^(1/2)*(-7+5^(1/2)) -1/44*(42+6*5^(1/2))^(1/2)*(-7+5^(1/2)) w = 1/6*(54+6*5^(1/2))^(1/2) -1/6*(54+6*5^(1/2))^(1/2) 1/6*(42+6*5^(1/2))^(1/2) -1/6*(42+6*5^(1/2))^(1/2) 再输入 vpa(u,6),vpa(v,6),vpa(w,6) ans = -1.61804 -1.61804 -.618040 -.618040 ans = -.730749 .730749 -.805994 .805994 ans = 1.36846 -1.36846 1.24071 -1.24071 例7.2分别求出一元二次方程ax2+bx+5=0和psin(x)=r三角方程的根 解 在matlab命令窗口中输入 x=solve(a*x^2+b*x+5) 回车得到 x = 1/2/a*(-b+(b^2-20*a)^(1/2)) 1/2/a*(-b-(b^2-20*a)^(1/2)) 即 再输入 x=solve(p*sin(x)=r) 回车得到 x = asin(r/p) 即x=arcsin(r/p) 7.2非线性方程的求解方法 本部分的几种数值方法，在不加说明的情况下，都是在某个区间内求的方程的一个根，其理论依据为：如果f(a)f(b)0,则在区间[a,b]上方程f(x)=0至少存在一个实根。 7.2.1二分法 2.二分法的matlab实现 function root=HalfInterval(f,a,b,eps) if(nargin==3) eps=1.0e-4; end f1=subs(sym(f),findsym(sym(f)),a); f2=subs(sym(f),findsym(sym(f)),b); if(f1==0) root=a; end if(f2==0) root=b; end if(f1*f20) disp(两端点函数值乘积大于0!); return; else root=FindRoots(f,a,b,eps); end function r=FindRoots(f,a,b,eps) f_1=subs(sym(f),findsym(sym(f)),a); f_2=subs(sy