24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 1201  |  回复: 7

jiyude

木虫 (正式写手)

[求助] 圆的半径和圆心角 已有2人参与

已知弧长36.7米,拱高10.4米的圆弧,求圆的半径和圆心角(最好附带Matlab程序)。怎么限制x是锐角?

clear all;
syms x y
[x,y]=solve('2*x*y-36.7=0','y*(1-cos(x))-10.4=0','x','y');
solution=[x,y]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

努力努力
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็

【答案】应助回帖

感谢参与,应助指数 +1
CODE:
s=36.7;h=10.4;  r=fsolve(@(r)2*h^2+1-4*h*r-cos(s/r),4.1)

ส็็็็็็็็็็็็็็็็็็็็
2楼2015-07-15 06:30:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็

1. 这种问题不太可能得到符号解,只能数值求解。
2. 锐角也无法限制,只能通过添加 r 与 h 关系的约束确定是优弧劣弧,比如 r > h >0 , 0<h<2r
3. 综合看,你的问题更适合用带约束条件的非线性最优化来解决,或者求单个根的数值方法。

具体条件根据自己的要求琢磨
ส็็็็็็็็็็็็็็็็็็็็
3楼2015-07-15 06:37:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็

【答案】应助回帖

前面的公式写错了,结果会错误。

纠正关于 半径 的0值函数为:



绘出它的图像可以看出有三个根适用:

圆的半径和圆心角

如果用Matlab求数值解,则带信赖域的牛顿法初值很重要,决定解
代码仍旧类似

初值为1,则结果是
CODE:
s=36.7;h=10.4; x0=1; r=fsolve(@(r)2*h^2+r^2-4*h*r-r^2*cos(s/r),x0); disp(num2str(r, '% 16.9f'))

得到 5.323271244

初值10
CODE:
s=36.7;h=10.4; x0=10;r=fsolve(@(r)2*h^2+r^2-4*h*r-r^2*cos(s/r),x0);disp(num2str(r,'%16.9f'))

得到
11.183055764

初值18
CODE:
s=36.7;h=10.4; x0=18;r=fsolve(@(r)2*h^2+r^2-4*h*r-r^2*cos(s/r),x0);disp(num2str(r,'%16.9f'))

得到
13.999546228
ส็็็็็็็็็็็็็็็็็็็็
4楼2015-07-15 07:13:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer


【答案】应助回帖

感谢参与,应助指数 +1
设半径为r,对应的圆心角为θ,由题意有:
r*θ=36.7                                          (1)
r*[1-Cos(θ/2)]=10.4                           (2)
先假设0<θ<π/2,
由泰勒公式,先只取三项:
1-Cos(θ/2)≈1/8*θ^2-1/384*θ^4         (3)
(1)、(3)代入方程(2):
θ*[θ^3-64*θ+108.8174]=0
由θ≠0,故:
θ^3-64*θ+108.8174=0
由函数知识可分析得到:它在-8/sqrt(3)有极大值305.8863,在8/sqrt(3)有极小值-88.2515。而在这个区间里函数单调下降,故有一个且只有一个根。8=0时函数值为108.8174;为π/2时为12.1622,故不存在此区域的根。
5楼2015-07-15 13:29:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็

引用回帖:
4楼: Originally posted by cooooldog at 2015-07-15 07:13:10
前面的公式写错了,结果会错误。

纠正关于 半径 r的0值函数为:

f(r)=r^2+r^2 \left(-\cos \left(\dfrac{36.7}{r}\right)\right)-41.6 r+216.32

绘出它的图像可以看出有三个根适用:



如果用Matlab求 ...

太小的根应该去掉; 因为导致s/r>2Pi显然不合理

看上去两个近似解

而 s/h <大约5.36 时 不存在锐角的可能, 反求也只能找到临界比值的近似解.
ส็็็็็็็็็็็็็็็็็็็็
6楼2015-07-16 06:25:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็

引用回帖:
6楼: Originally posted by cooooldog at 2015-07-16 06:25:10
太小的根应该去掉; 因为导致s/r>2Pi显然不合理

看上去两个近似解

而 s/h <大约5.36 时 不存在锐角的可能, 反求也只能找到临界比值的近似解....

应该是弧长与拱高之比小于则不会有锐角
ส็็็็็็็็็็็็็็็็็็็็
7楼2015-07-16 06:34:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
jiyude: 金币+10 2015-07-18 09:23:59
引用回帖:
7楼: Originally posted by cooooldog at 2015-07-16 06:34:42
应该是弧长与拱高之比小于\dfrac{\pi }{2-\sqrt{2}}\approx 5.363则不会有锐角...

应该是只有唯一解;

三个解中, 第一个不满足角度条件

第二个得到角度为优弧,高度为劣弧,也不是正确解

第三个是唯一正确的, 作图出来是:

圆的半径和圆心角-1
ส็็็็็็็็็็็็็็็็็็็็
8楼2015-07-16 10:04:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jiyude 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂学校 +3 不会吃肉 2026-04-13 3/150 2026-04-13 23:13 by 木鱼.
[考研] 290求调剂 +16 luoziheng 2026-04-10 16/800 2026-04-13 20:34 by biomen
[考研] 复试调剂 +17 积极向上; 2026-04-10 19/950 2026-04-13 19:53 by 闲人终南山
[考研] 368化学求调剂 +14 wwwwabcde 2026-04-07 15/750 2026-04-13 08:36 by lhj2009
[考研] 0831一轮调剂失败求助 +10 小熊睿睿_s 2026-04-11 10/500 2026-04-12 22:43 by 长弓傲
[考研] 本科南方医科大学 一志愿985 药学学硕284分 求调剂 +5 弱水听文 2026-04-09 5/250 2026-04-12 13:16 by shengxi123
[考研] 307求调剂 +10 tzq94092 2026-04-10 10/500 2026-04-12 08:18 by wise999
[考研] 085404 293求调剂 +9 勇远库爱314 2026-04-08 9/450 2026-04-12 02:24 by 秋豆菜芽
[考研] 调剂 +10 只叙离别辞 2026-04-09 12/600 2026-04-11 20:57 by 逆水乘风
[考研] 343求调剂 +9 王国帅 2026-04-10 9/450 2026-04-11 20:31 by dongdian1
[考研] 296求调剂 +14 汪!?! 2026-04-08 15/750 2026-04-11 20:28 by dongdian1
[考研] 296求调剂 +6 汪!?! 2026-04-09 6/300 2026-04-11 11:25 by zhq0425
[考研] 085800 能源动力求调剂 +6 阿biu啊啊啊啊啊 2026-04-10 6/300 2026-04-10 15:03 by hemengdong
[考研] 机械还有还有名额吗?太难了 +6 笑笑袁 2026-04-10 6/300 2026-04-10 11:54 by 高维春
[考研] 一志愿中南大学物理学,英一66,求调剂 +4 长烟旖旎 2026-04-08 5/250 2026-04-10 10:31 by 颖果儿
[考研] 265求调剂 +4 风说她早忘了 2026-04-07 4/200 2026-04-09 13:59 by only周
[考研] 086000生物与医药调剂 +7 awwwwwooooo 2026-04-09 7/350 2026-04-09 13:31 by 北极159263
[考研] 求调剂 +13 柒luck 2026-04-07 13/650 2026-04-08 22:46 by 猪会飞
[考研] 331求调剂 +5 luoxin0706. 2026-04-08 5/250 2026-04-08 22:15 by zhouyuwinner
[考研] 电子信息346 +4 zuoshaodian 2026-04-08 4/200 2026-04-08 11:54 by zzucheup
信息提示
请填处理意见