应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 三元二次的三角方程如何用MATLAB求解？？
111/2返回列表12下一页
查看: 4235  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

cyclonelee

至尊木虫 (著名写手)

[求助] 三元二次的三角方程如何用MATLAB求解？？

三个方程如下：

r1=(-k*sin(a)*sin(b)+z*tan(b)+m)^2+(z-k*sin(a)*cos(b))^2
r2=(-k*sin(a)*sin(b)+z*tan(b)-m)^2+(z-k*sin(a)*cos(b))^2
r3=(k*sin(a)*sin(b)+z*tan(b))^2+(2*k*cos(a)-n)^2+(k*sin(a)*cos(b)+z)^2

a,b,z 为待求未知量，其余为已知量

求高手帮忙！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
既然认准一条道路，何必去打听要走多久。
1楼 2013-08-22 17:57:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

信彼南山

木虫 (著名写手)
你要解析解么
如果要数值解的话，为什么不把已知量给出呢
赞一下 回复此楼
2楼2013-08-22 18:44:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cyclonelee

至尊木虫 (著名写手)
引用回帖:
2楼: Originally posted by 信彼南山 at 2013-08-22 18:44:10
你要解析解么
如果要数值解的话，为什么不把已知量给出呢

已知量其实可以任取的,如：
r1=120 r2=110 r3=100
m=40 n=30 k=10
赞一下 回复此楼
既然认准一条道路，何必去打听要走多久。
3楼2013-08-22 22:06:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

星月夜dh

金虫 (小有名气)

【答案】应助回帖

★ ★ ★
感谢参与，应助指数 +1
cyclonelee: 金币+2, ★有帮助, 多谢！ 2013-08-22 23:02:34
xiegangmai: 金币+1, 谢谢参与！ 2013-08-25 00:39:32
这个有些类似于公式推倒的样子，相当于一种符号计算，最近在一本书上正好看到过，找一本软件学习的书吧，上面应该都会有的。。。
赞一下 回复此楼
4楼2013-08-22 22:45:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
cyclonelee: 金币+8, ★★★★★最佳答案, 非常感谢！！！！我先试试去！多谢 2013-08-22 23:03:28
xiegangmai: 金币+1, 谢谢参与！ 2013-08-25 00:39:37
数值解可参考：
http://muchong.com/bbs/viewthread.php?tid=6150808&authorid=1122189
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
5楼2013-08-22 22:48:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cyclonelee

至尊木虫 (著名写手)
引用回帖:
4楼: Originally posted by 星月夜dh at 2013-08-22 22:45:13
这个有些类似于公式推倒的样子，相当于一种符号计算，最近在一本书上正好看到过，找一本软件学习的书吧，上面应该都会有的。。。

能否告知书名
回复此楼
既然认准一条道路，何必去打听要走多久。
6楼2013-08-22 22:55:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cyclonelee

至尊木虫 (著名写手)
引用回帖:
5楼: Originally posted by 月只蓝 at 2013-08-22 22:48:47
数值解可参考：
http://muchong.com/bbs/viewthread.php?tid=6150808&authorid=1122189

试做了一下：
程序如下：

function solvenonl
clear all;clc
x0=100*[1 1 1];
x=fsolve(@nonl,x0)


function f=nonl(x)
a=x(1);
b=x(2);
z=x(3);
r1=120;
r2=110;
r3=90;
m=40;
n=30;
k=10;
f(1)=(-k*sin(a)*sin(b)+z*tan(b)+m)^2+(z-k*sin(a)*cos(b))^2-r1;
f(2)=(-k*sin(a)*sin(b)+z*tan(b)-m)^2+(z-k*sin(a)*cos(b))^2-r2;
f(3)=(k*sin(a)*sin(b)+z*tan(b))^2+(2*k*cos(a)-n)^2+(k*sin(a)*cos(b)+z)^2-r3;

结果：

Solver stopped prematurely.

fsolve stopped because it exceeded the function evaluation limit,
options.MaxFunEvals = 300 (the default value).


x =

  101.4392  100.5481   98.0448

问下，为什么只有一组解，要是想保留多组解该怎么处理？
上面的英文是不是说明结果太多，溢出了？
赞一下 回复此楼
既然认准一条道路，何必去打听要走多久。
7楼2013-08-22 23:24:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

星月夜dh

金虫 (小有名气)
引用回帖:
6楼: Originally posted by cyclonelee at 2013-08-22 22:55:28
能否告知书名...

额。。。前段时间看的，叫什么Matlab有入门到精通还是啥的，具体忘了，好多书上都有这方面的内容的，你找到书后可以先看看目录的啊。。。
赞一下 回复此楼
8楼2013-08-22 23:36:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
cyclonelee: 金币+5, ★★★很有帮助, 多谢！ 2013-08-23 09:19:47
xiegangmai: 金币+2, 谢谢参与！ 2013-08-25 00:39:47
引用回帖:
7楼: Originally posted by cyclonelee at 2013-08-22 23:24:37
试做了一下：
程序如下：

function solvenonl
clear all;clc
x0=100*;
x=fsolve(@nonl,x0)


function f=nonl(x)
a=x(1);
b=x(2);
z=x(3);
r1=120;
r2=110;
r3=90;
m=40;
n=30;
k=10;
f(1)=(- ...

第三行
x0=100*[1 1 1]; 这里，你先给未知数一个大概的初值，原来这里的x0只适合那个帖子中的情形。
初值肯定得调的，而且初值给在不同的数值区间，多组解也能算出来。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
9楼2013-08-22 23:47:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cyclonelee

至尊木虫 (著名写手)
引用回帖:
9楼: Originally posted by 月只蓝 at 2013-08-22 23:47:58
第三行
x0=100*; 这里，你先给未知数一个大概的初值，原来这里的x0只适合那个帖子中的情形。
初值肯定得调的，而且初值给在不同的数值区间，多组解也能算出来。...

第三行初值的问题我明白了！
但是多组解还是没有搞懂！我这个求解的是三角函数中的角度值，是不是可以这样理解，在0~2PI设置不同的初值来获得不同的解？但是，初值的选取怎么选呢？如何可以选择到所有合理的初值我还是没想出来，望解答！多谢，辛苦了！
赞一下 回复此楼
既然认准一条道路，何必去打听要走多久。
10楼2013-08-23 09:23:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cyclonelee 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 6/300 2026-09-14 22:37 by vZfe6xYu34yj
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-14 22:25 by vZfe6xYu34yj
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-14 22:24 by vZfe6xYu34yj
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-14 22:01 by vZfe6xYu34yj
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 3n8v2C8RimXI 2026-09-13 5/250 2026-09-14 22:00 by vZfe6xYu34yj
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +4 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:33 by vZfe6xYu34yj
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:25 by vZfe6xYu34yj
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:24 by vZfe6xYu34yj
[公派出国] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 QUjhNVAcOSff 2026-09-13 7/350 2026-09-14 21:21 by vZfe6xYu34yj
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-14 18:01 by Aj1rhIDL5ixY
[考博] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-14 17:49 by Aj1rhIDL5ixY
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-14 17:48 by Aj1rhIDL5ixY
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +5 QUjhNVAcOSff 2026-09-13 6/300 2026-09-14 15:21 by Aj1rhIDL5ixY
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 15:09 by Aj1rhIDL5ixY
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 14:57 by Aj1rhIDL5ixY
[基金申请] 求教各位大神：2026教育部人文社科青年基金项目何时公示呀？ +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[找工作] 浙江师范大学是怎么坑我的 +6 recruit123 2026-09-09 6/300 2026-09-14 12:00 by 萨仁其其格
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
[基金申请] 面上有专家说实验设备不是我们单位的 +6 zj_muchong 2026-09-08 8/400 2026-09-11 06:56 by suifengdao
信息提示
关闭
请填处理意见
关闭