24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 1199  |  回复: 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

铁杆木虫 (著名写手)

ส็็็

引用回帖:
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的回帖
查看全部 8 个回答

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求助调剂,跨调 +14 X十甫寸Y 2026-04-11 15/750 2026-04-14 00:36 by 王珺璞
[考研] 求调剂 +12 璃茉一定上岸 2026-04-10 13/650 2026-04-14 00:08 by Equinoxhua
[考研] 279学硕食品专业求调剂院校 20+6 孤独的狼爱吃羊 2026-04-12 23/1150 2026-04-13 20:42 by biomen
[考研] 266调剂 +10 daya sun 2026-04-07 11/550 2026-04-13 10:12 by fenglj492
[考研] 344 材料专业 求调剂211 无地域要求 +8 hualkop 2026-04-11 8/400 2026-04-12 22:24 by fqwang
[找工作] 山东高校教师考核超级无底线,员工过不下去啦 +4 qut2026 2026-04-09 9/450 2026-04-12 00:54 by qut2026
[考研] 277 数一104,学硕,求调剂 +21 瓶子PZ 2026-04-09 23/1150 2026-04-11 23:12 by labixiaoqiao
[考研] 086003调剂求助 +21 苏弋万 2026-04-09 22/1100 2026-04-11 20:25 by dongdian1
[考研] 085410-273求调剂 +6 X1999 2026-04-10 6/300 2026-04-11 10:32 by Delta2012
[考研] 工科273调剂 +6 X1999 2026-04-09 7/350 2026-04-11 10:23 by zhq0425
[考研] 广东省 085601 329分求调剂 +14 Eddieddd 2026-04-10 14/700 2026-04-11 09:58 by bljnqdcc
[考研] 0854调剂 +4 长弓傲 2026-04-09 4/200 2026-04-11 09:18 by 猪会飞
[考研] 302分求调剂 +9 凡语祈愿 2026-04-08 10/500 2026-04-10 23:26 by 314126402
[考研] 一志愿沪9,326生物学求相关专业调剂 +4 刘墨墨 2026-04-09 4/200 2026-04-10 12:07 by pengliang8036
[考研] 292求调剂 +9 笑笑袁 2026-04-09 9/450 2026-04-10 10:05 by LHGeng
[考研] 本科211 工科085400 280分求调剂 可跨专业 +3 LZH(等待调剂中 2026-04-09 3/150 2026-04-09 21:29 by wutongshun
[考研] 材料专硕初试分332一志愿西北工业大学, +12 故人?? 2026-04-09 12/600 2026-04-09 18:34 by Ccclqqq
[考研] 348求调剂 +3 candyyyi 2026-04-09 3/150 2026-04-09 17:20 by 段伟艳
[论文投稿] 求助文献原文 10+3 18500821399 2026-04-08 3/150 2026-04-09 16:56 by 北京莱茵润色
[考研] 22408 266求调剂 +11 masss11222 2026-04-07 14/700 2026-04-08 11:06 by yulian1987
信息提示
请填处理意见