| 查看: 1106 | 回复: 3 | ||
[求助]
matlab解方程组求助 已有1人参与
|
|
已知圆心不在坐标原点的圆上两点相对圆心的坐标，两点到原点的连线分别与圆心到原点连线的夹角也已知，求圆心坐标。 用matlab写的程序，能够得到正确结果，但是改变两个点的位置，比如将theta改为theta = [0, pi/3]，就得不到正确的解，用folve命令也尝试过，得到的结论一样，就大神指教 x0 = 15; %圆心（未知） y0 = 18; %圆心（未知） r = 3; %圆半径（已知） theta = [0, pi/6]; %圆上两点相对圆心的极角（已知） %theta = [0, pi/3]; x = x0 +r*cos(theta); %圆上两点的坐标（未知） y = y0 +r*sin(theta); %圆上两点的坐标（未知） phi0 = atan(y0/x0); %圆心相对原点的极角（已知） phi = atan(y./x); %圆上两点相对原点的极角（已知） syms x0_c y0_c x_c = x0_c +r*cos(theta); y_c = y0_c +r*sin(theta); eq1 = phi0 - phi(1) - acos((x_c(1)*x0_c+y_c(1)*y0_c)/(sqrt(x_c(1)^2+y_c(1)^2)*sqrt(x0_c^2+y0_c^2))); eq2 = phi0 - phi(2) - acos((x_c(2)*x0_c+y_c(2)*y0_c)/(sqrt(x_c(2)^2+y_c(2)^2)*sqrt(x0_c^2+y0_c^2))); [x0_c, y0_c] = solve(eval(eq1), eval(eq2)); x0_c = double(x0_c); y0_c = double(y0_c); |
回复此楼» 猜你喜欢
求助中医院临床中医中标的青年标书学习，谢谢
已经有6人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
物理学I论文润色/翻译怎么收费?
已经有279人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
独孤神宇
版主 (知名作家)
- 应助: 490 (硕士)
- 贵宾: 0.008
- 金币: 31030.3
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5603
- 在线: 856.8小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
2楼2018-05-21 16:46:40
3楼2018-05-21 16:57:54
独孤神宇
版主 (知名作家)
- 应助: 490 (硕士)
- 贵宾: 0.008
- 金币: 31030.3
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5603
- 在线: 856.8小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
4楼2018-05-21 20:12:39
50