基于二分思想的圆弧计算工具的设计与实现.docVIP

基于二分思想的圆弧计算工具的设计与实现 　　摘要：二分法在超越方程的代数运算求解中起着重要的作用。该文以在开发系统过程中如何实现“圆弧计算”功能之一的“已知弧长和弦长，求拱高和半径（或夹角）” 为例，分析了几个参数之间的数学关系，最终用二分法求出近似值。 　　关键词：超越方程；二分法 　　中图分类号：TP312文献标识码：A文章编号：1009-3044(2011)30-7443-02 　　Design and Implementation of Arc Computational Tools on Dichotomy Application of Dichotomy in Arc Computing 　　ZHAO Zhen-qi 　　(Wuxi Machinery and Electron Higher Professional and Technical School, Wuxi 214028, China) 　　Abstract: Dichotomy in the transcendental equations in solving the algebraic operations, plays an important role in. This paper in the system development process how to implement arc computing function of known arc and chord length, and arch height and radius (or angle) as an example, analyzes the mathematical relationship between several parameters, finally use dichotomy to seek out the approximate value. 　　Key words: transcendental equation; dichotomy 　　在人类用智慧架设的无数座从未知通向已知的金桥中，方程的求解理所当然是其中璀璨的一座。在中学阶段通过代数运算能求出一元一次方程，一元二次方程及三角方程等基本初等方程的解，而指数方程、对数方程等超越方程和五次以上的高次代数方程不能用代数运算的求解。但其数值解法随着现代计算技术的发展得到了广泛的应用，如二分法等。[1] 　　笔者在开发玻璃深加工管理系统过程中，按照客户需求，开发弯钢玻璃计算辅助工具，其中一个“圆弧计算”，要求实现四个基本功能：1）已知弦长和拱高，求弧长；2）已知半径和弦长，求弧长和拱高；3）已知半径和拱高，求弧长和弦长；4）已知弧长和弦长，求拱高和半径(或夹角)。前三个是求解三角函数，属于初等数学范畴，最后一个要求解超越方程，属于高等数学范畴。下面将分析如何得到此超越方程，确定二分法解决方案，并最终用C#语言具体实现其功能。 　　1 数学原理 　　已知：L为玻璃的初始长度，w为弦长，求解：弯曲后玻璃弧形拱高和半径（或夹角）。 　　数学分析过程： 　　令夹角的一半为θ，θ=L/2r，sin(θ)=w/2r，由此可推出公式w=2rsin(L/2r)， 　　公式两边再同时除以L，则w/L=2r/Lsin(L/2r)， 　　令x=L/2r，即x为弧长一半所对应的圆心角（x=θ），则题目转化为公式w/L=sin(x)/x， 　　设f(x)= sin(x)/x，此函数显然当x∈(0,π)时是单调递减的，因此可以采用二分求解，得到x的近似值。 　　最后，确定拱高h与x的关系表达式，从而求出h值，如图1弧形几何作图所示，在弦长一半BC（即w/2）与所求高度AC（即h）所围成的直角三角形△ABC中，由正弦定理得：h/sin(z)=(w/2)/sin(y)，其中y=Π/2-x/2，z =x/2（OA=OB，所以△OAB为等腰三角形，∠OAB=∠OBA，∠OAB=(Π-∠AOB)/2=(Π-x)/2=Π/2-x/2=y，∠ABC=Π/2-∠BAC=Π/2-(Π-x)/2=x/2=z，公式化简为：h=(1-cos(x))*w/(2*sin(x))。 　　显然，根据勾股定理，r^2=(r-h)^2+(w/2)^2，得到r=(h*h+w*w/4)/(2*h)。 　　2 解决方案 　　超越方程一般没有解析解，而只有数值解或近似解，只有特殊的超越方程才可以求出解析解来。常用的求解超越方程的近似解法有牛顿切线法、幂级数解法等。那么程序设计又应该采用什么方法呢？ 　　解题思路：推出公式后，就可以用二分枚举区间内对应的角度求解，分析： 　　这里设s=0，ed=Math.PI；开始二分，当ed?C