弦截法非线性方程.docxVIP

一、科学计算的算法及其举例应用和利用 MATLAB自带函 数实现 科学计算法：弦截法非线性方程 实例：采用弦截法求方程lgx+x = 2在区间［1, 4］上的一个根。 1?算法说明 弦截法的算法过程如下： 过两点(a,f(a)),(b,f(b))作一直线，它与x轴有一个交点，记为x1.. 如果f(a)f(x1)0,过两点(a,f(a)),(x1,f(x1))作一直线，它与x轴的交点记为 x2, 否则过两点(b,f(b)),(x1,f(x1))作一直线，它与x轴的交点记为x2 ; 如此下去，直到|xn-xn-1| §就可以认为xn为f(x)=0在区间［a,b］上的一个 根。 xk-1-a?? x=a-f(a),f(a)f(x)0kk-1 ?? f(xk-1)-f(a) ?? (4) Xk 的递推公式为：??且 xk-1-b?xk=b-f(b),f(a)f(xk-1)0 ??? f(xk-1)-f(b) ?? x仁a-b-af(a) f(b)-f(a) 在MATLAB中编程实现的弦截法的函数为：Secant. 功能：用弦截法求函数在某个区间的一个零点。 调用格式：root=Seca nt(f,a,b,eps). 其中，f为函数名； a为区间左端点； b为区间左端点； eps为根的精度； root为求出的函数零点。 2?流程图 弦截法流程图 实例流程图： 开始 输入条件： sqrt(x)+log(x)-2 积分区间：(1,4) 调用Secan进行运算 输出结果 r =1.8773 结束 3.源程序代码 function root=Secant(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 tol=1; fa=subs(sym(f),findsym(sym(f)),a); fb=subs(sym(f),findsym(sym(f)),b); root=a-(b- a)*fa/(fb-fa); while(toleps) r仁root; fx=subs(sy m(f),fin dsym(sy m(f) ),r1); s=fx*fa; if(s==0) root=r1; else if(s0) root=b-(r1-b)*fb/(fx-fb); else root=a-(r1-a)*fa/(fx-fa); end end tol=abs(root-r1); end end 实例源代码： r=Seca nt(sqrt(x)+log(x)-2,1,4) 运行结果： r = 1.8773 二、分析单自由度阻尼系统的阻尼系数对其固有振动模态的 影响 1.流程图 引」 2?源程序代码: 2 3 4 5 6 E 9 10 11 12 14 15 16 IF 18 】9 20 21 22 WTl=10h x0=l;张初贻垃置 v0=0; 歸初始速度 tf=2H 1 - ：ll for j=l:10 leta(j)s0. J*j;鴛阻居至魏 wd(j)=vti*^rt11-teta(j)⑵； A=sqrt ( (wn^xO^zet at j +v0 2+(x0*wd(j)) *2)/wd(j) ; * F^ - a=it ari2 (urd( j) *x0, vO+zeta1； j) rwn+s J)； 訝二工先 t=O:tf/1000：tf x(j, :?-A*eip(-let 1(j)*wnat).?占in(wd(jO *t+a) end subplot (1, 2,1＞站结制耳八二维圈徵 for j=l：2：10 plot(t?3t (jF ：)} text (0, 3a j, 151 )j [ zetap=, nu?2str (zeta( j))]) hold o：. end ^label C t* ),ylabel(*x* )| grid Mubpl忒⑴ 2,2)昵hguwi m；； 3?运行结果: BHK Ficure 1 三.Appollo卫星的运行轨迹 1.流程图 2?源程序代码： 函数文件： function dp=appollo(t? p) - u~l/82.45; - in^l-u; - r l=sqrt ( (p (1) -Hi) 2+p (3) * 2) j - r2=sqrt ((p(I)+ni) 2+p(3) 2): - dp= [p (2) 2*p(4)-^j(i)-j