牛顿法解流程图
- 格式:doc
- 大小:75.00 KB
- 文档页数:1
应用牛顿运动定律解题的方法和步骤Document serial number【LGGKGB-LGG98YT-LGGT8CB-LGUT-§3.4应用牛顿运动定律解题的方法和步骤应用牛顿运动定律的基本方法是隔离法,再配合正交坐标运用分量形式求解。
解题的基本步骤如下:(1)选取隔离体,即确定研究对象一般在求某力时,就以此力的受力体为研究对象,在求某物体的运动情况时,就以此物体为研究对象。
有几个物体相互作用,要求它们之间的相互作用力,则必须将相互作用的物体隔离开来,取其中一物体作研究对象。
有时,某些力不能直接用受力体作研究对象求出,这时可以考虑选取施力物体作为研究对象,如求人在变速运动的升降机内地板的压力,因为地板受力较为复杂,故采用人作为研究对象为好。
在选取隔离体时,采用整体法还是隔离法要灵活运用。
如图3-4-1要求质量分别为M 和m 的两物体组成的系统的加速度a ,有两种方法,一种是将两物体隔离,得方程为 另—种方法是将整个系统作为研究对象,得方程为 显然,如果只求系统的加速度,则第二种方法好;如果还要求绳的张力,则需采用前一种方法。
(2)分析物体受力情况:分析物体受力是解动力学问题的一个关键,必须牢牢掌握。
①一般顺序:在一般情况下,分析物体受力的顺序是先场力,如重力、电场力等,再弹力,如压力、张力等,然后是摩擦力。
并配合作物体的受力示意图。
大小和方向不受其它力和物体运动状态影响的力叫主动力,如重力、库仑力;大小和主向与主动力和物体运动状态有密切联系的力叫被动力或约束力,如支持力、摩擦力。
这m图3-4-1就决定了分析受力的顺序。
如物体在地球附近不论是静止还是加速运动,它受的重力总是不变的;放在水平桌面上的物体对桌面的压力就与它们在竖直方向上有无加速度有关,而滑动摩擦力总是与压力成正比。
②关于合力与分力:分析物体受力时,只在合力或两个分力中取其一,不能同时取而说它受到三个力的作用。
一般情况下选取合力,如物体在斜面上受到重力,一般不说它受到下滑力和垂直面的两个力。
1:二分法流程图:二分法基本思路:一般地,对于函数f(x),如果存在实数c,当x=c 时,若f(c)=0,那么把x=c 叫做函数f(x)的零点。
解方程即要求f(x)的所有零点。
假定f(x)在区间(x ,y )上连续先找到a 、b 属于区间(x ,y ),使f(a),f(b)异号,说明在区间(a,b)内一定有零点,然后求f[(a+b)/2],现在假设f(a)<0,f(b)>0,a<b① 如果f[(a+b)/2]=0,该点就是零点,如果f[(a+b)/2]<0,则在区间((a+b)/2,b)内有零点,(a+b)/2>=a ,从①开始继续使用② 中点函数值判断。
如果f[(a+b)/2]>0,则在区间(a,(a+b)/2)内有零点,(a+b)/2<=b ,从①开始继续使用 中点函数值判断。
这样就可以不断接近零点。
通过每次把f(x)的零点所在小区间收缩一半的方法,使区间的两个端点逐步迫近函数的零点,以求得零点的近似值,这种方法叫做二分法。
从以上可以看出,每次运算后,区间长度减少一半,是线形收敛。
另外,二分法不能计算复根和重根。
二分法步骤:用二分法求方程()0f x =的根*x 的近似值k x 的步骤① 若对于a b <有()()0f a f b <,则在(,)a b 内()0f x =至少有一个根。
② 取,a b 的中点12a b x +=计算1()f x ③ 若1()0f x =则1x 是()0f x =的根,停止计算,运行后输出结果*1x x =若1()()0f a f x <则在1(,)a x 内()0f x =至少有一个根。
取111,a a b x ==;若1()()0f a f x >,则取111,a x b b ==;④ 若12k k b a ε-≤(ε为预先给定的要求精度)退出计算,运行后输出结果*2k ka b x +≈,反之,返回步骤1,重复步骤1,2,3二分法Mtalab 程序syms x;fun=input('(输入函数形式)fx=');a=input('(输入二分法下限)a=');b=input('(输入二分法上限)b=');d=input('输入误差限d=')%二分法求根%f=inline(x^2-4*x+4);%修改需要求解的inline函数的函数体f=inline(fun);%修改需要求解的inline函数的函数体e=b-a; k=0 ;while e>dc=(a+b)/2;if f(a)*f(c)<0b=c;elseif f(a)*f(c)>0a=c;elsea=c;b=cende=e/2; k=k+1;endx=(a+b)/2;x%x为答案k%k为次数2,牛顿法及流程图:方程f(x)=0的根就是曲线y=f(x)与x轴交点的横坐标x*,当初始近似值x0选取后,过( x0,f(x0))作切线,其切线方程为:y- f(x0)=f′(x0)(x-x0)它与x轴交点的横坐标为x一般地,设是x*的第n次近似值,过( x,f(x))作y=f(x)的切线,其切线与x轴交点的横坐标为:x = - 即用切线与x轴交点的横坐标近似代曲线与x轴交点的横坐标,如图牛顿法正因为有此明显的几何意义,所以也叫切线法。
《MATLAB 程序设计实践》课程考核1、编程实现以下科学计算算法,并举一例应用之。
(参考书籍《精通MALAB科学计算》,王正林等著,电子工业出版社,2009年)牛顿-科茨(Newton-Cotes)公式算法:Step 1:判断type类型,1转Step 2,2转Step 3,3转Step 4;否则输出值为0;Step 2:计算科茨公式:Step 3:计算牛顿-科茨六点公式:Step 4:计算牛顿-科茨七点公式:“牛顿-科茨(Newton-Cotes)公式”function I = NewtonCotes(f,a,b,type)%type = 1 科茨公式%type = 2 牛顿-科茨六点公式%type = 3 牛顿-科茨七点公式I=0;switch typecase 1,I=((b-a)/90)*(7*subs(sym(f),findsym(sym(f)),a)+...32*subs(sym(f),findsym(sym(f)),(3*a+b)/4)+...12*subs(sym(f),findsym(sym(f)),(a+b)/2)+...32*subs(sym(f),findsym(sym(f)),(a+3*b)/4)+7*subs(sym(f),findsym(sym(f)),b));case 2,I=((b-a)/288)*(19*subs(sym(f),findsym(sym(f)),a)+...75*subs(sym(f),findsym(sym(f)),(4*a+b)/5)+...50*subs(sym(f),findsym(sym(f)),(3*a+2*b)/5)+...50*subs(sym(f),findsym(sym(f)),(2*a+3*b)/5)+...75*subs(sym(f),findsym(sym(f)),(a+4*b)/5)+19*subs(sym(f),findsym(sym(f)),b));case 3,I=((b-a)/840)*(41*subs(sym(f),findsym(sym(f)),a)+...216*subs(sym(f),findsym(sym(f)),(5*a+b)/6)+...27*subs(sym(f),findsym(sym(f)),(2*a+b)/3)+...272*subs(sym(f),findsym(sym(f)),(a+b)/2)+...27*subs(sym(f),findsym(sym(f)),(a+2*b)/3)+...216*subs(sym(f),findsym(sym(f)),(a+5*b)/6)+41*subs(sym(f),findsym(sym(f)),b));end实例:计算积分dxx100sin。