24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1210  |  回复: 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

铁杆木虫 (著名写手)

ส็็็

【答案】应助回帖

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

纠正关于 半径 的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的回帖
查看全部 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的回帖

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 通信工程求调剂!!! +6 zlb770521 2026-04-14 6/300 2026-04-15 20:00 by 学员JpLReM
[考研] 211本科材料化工求调剂 +19 YHLAH 2026-04-11 23/1150 2026-04-14 22:25 by fenglj492
[考研] 调剂 +12 月@163.com 2026-04-11 12/600 2026-04-14 15:37 by zs92450
[考研] 一志愿沪9,326求生物学调剂 +10 刘墨墨 2026-04-13 10/500 2026-04-14 15:16 by zs92450
[考研] 考研调剂 +13 长弓傲 2026-04-13 14/700 2026-04-14 14:44 by zs92450
[考研] 300分求调剂 (085501机械专硕,本科扬大) +9 xu@841019 2026-04-11 10/500 2026-04-14 08:48 by 木木mumu~
[考研] 一志愿西交机械专硕求调剂 +9 求上岸的小王 2026-04-10 9/450 2026-04-13 16:08 by jiangguiquan11
[考研] 一志愿浙大生物325分求调剂 +9 zysheng 2026-04-12 9/450 2026-04-12 22:31 by yuyin1233
[硕博家园] 新一代电子信息294求调剂 不挑学校 +7 Ytyt11 2026-04-09 8/400 2026-04-12 16:57 by ajpv风雷
[考研] 307求调剂 +10 tzq94092 2026-04-10 10/500 2026-04-12 08:18 by wise999
[考研] 一志愿厦大0856,306求调剂 +15 Bblinging 2026-04-11 15/750 2026-04-11 22:53 by 314126402
[考研] 求调剂 +6 archer.. 2026-04-09 8/400 2026-04-11 10:55 by zhq0425
[考研] 288求调剂 +15 代fish 2026-04-09 16/800 2026-04-11 10:26 by wwj2530616
[考研] 22408 327分求调剂 +4 韵风kon 2026-04-10 4/200 2026-04-11 09:51 by 猪会飞
[考研] 282,求调剂 +12 jggshjkkm 2026-04-09 14/700 2026-04-11 09:39 by 猪会飞
[考研] 085402通信工程调剂,有4项学科竞赛国奖(电赛国二),硕士研究生调剂自荐信。 +5 m永o不v言o弃m 2026-04-09 5/250 2026-04-11 09:33 by zhq0425
[考研] 346,工科0854求调剂,专硕 +7 moser233 2026-04-10 8/400 2026-04-11 08:52 by 猪会飞
[考研] 一志愿京区985,085401,与本科专业一致,电子信息工程, +4 阳光开朗的男孩 2026-04-10 4/200 2026-04-10 18:27 by shenrf
[考研] 一志愿华东师范生物学326分,求调剂 +8 刘墨墨 2026-04-09 8/400 2026-04-10 12:00 by pengliang8036
[考研] 083200 初试305分 求调剂 暂不考虑跨专业 +15 Claireyyyy 2026-04-09 15/750 2026-04-09 16:11 by zhuimr
信息提示
请填处理意见