24小时热门版块排行榜    

查看: 1652  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

dearzyh

新虫 (初入文坛)

[求助] 关于弓形面积求圆半径

已知弓形面积以及其对应弦长,怎么求这个圆的半径?我列方程出来时个超越方程,解不出来,求教!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer


【答案】应助回帖

设弓形的弦长为a,面积为S,所对的圆心角为θ。
    则S=θ/(2*π)*π*r2=θ*r2/2
    而  Sin(θ/2)=(a/2)/r=a/(2*r)  ,θ=2*arcsin[a/(2*r)]
    故有:S= r2  * arcsin[a/(2*r)] 。
此方程中,S和a均已知,因此其为关于未知数r的一元非线性方程,先试着手工计算给出一个近似值,然后可用欧拉割线法或牛顿切线法求出其近似数值解。具体可参考数学分析中一元微分学的应用相关章节。例如:若S=10000mm2, a=100mm.
令 f(r)=S-r^2*arcsin[a/(2*r)] , f(r)=10000-r2*arcsin(50/r)
则 df/dr=-2*r*arcsin[a/(2*r)] -r2/sqrt{1-[a/(2*r)]2}*[-a/(2*r2)]
           =a/sqrt{4-[a/r]2}-2*r*arcsin[a/(2*r)]
d2f/dr2=-2*arcsin[a/(2*r)] +a/sqrt[r2-(a/2)2]-{a/{r*sqrt[4-(a/r)2]}}3
r=150时,f(150)=10000-22500*arcsin(1/3)=2354>0
r=200时,f(200)=10000-40000*arcsin(1/4)=-107.2<0
故知:       150<r<200
  因为d2f/dr2(r=200)=-6.2*10-3<0,与f(200)同号,故以r=200为初值进行迭代。迭代公式:
Rn= Rn-1 +f(Rn-1)/df/dr(R=Rn-1)
  r1=200-107.2/49.43=197.83128,  r2=197.83128+0.020037/49.41868=197.83169
由于abs(r2-r1)/r1=2.05*10-6,可以认为r的值已经收敛。实际上,f(197.83169)=-0.0403。
9楼2013-11-06 20:08:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer


【答案】应助回帖

非常抱歉,刚才的应答是从word中复制粘贴的,由于此处上下标不分,导致里面的r^2全变成了r2,请楼主阅读时注意。若觉得可取,可告知信箱,我可发word文档。
10楼2013-11-06 20:14:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dearzyh 的主题更新
信息提示
请填处理意见