24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 3865  |  回复: 10

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 为待求未知量,其余为已知量

求高手帮忙!
回复此楼

» 猜你喜欢

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

既然认准一条道路,何必去打听要走多久。
已阅   回复此楼   关注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
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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿211,335分,0856,求调剂院校和导师 +15 倾____萧 2026-03-27 16/800 2026-04-02 22:50 by JourneyLucky
[考研] 296求调剂 +4 sdhu 2026-04-02 4/200 2026-04-02 21:29 by baoball
[考研] 一志愿武汉理工0856,初试334 +3 26考研材料 2026-04-02 3/150 2026-04-02 21:22 by dongzh2009
[基金申请] esi高被引论文是不是能对中标有所加分和帮助呢 +4 redcom 2026-04-01 5/250 2026-04-02 20:48 by redcom
[考研] 085602化工求调剂(331分) +9 111@127 2026-03-30 9/450 2026-04-02 20:00 by dick_runner
[考研] 275学硕081000服从调剂到其他专业,保不住本专业了 +7 一只小小水牛 2026-04-02 8/400 2026-04-02 14:23 by alice-2022
[考研] 一志愿厦门大学材料工程专硕354找调剂!!! +7 贝呗钡钡 2026-03-30 7/350 2026-04-02 12:19 by 1753564080
[考研] 085601一志愿中山大学深圳材料工程330求调剂 +8 pipiver 2026-03-30 8/400 2026-04-02 12:01 by ms629
[考研] 求调剂,一志愿 南京航空航天大学 ,080500材料科学与工程学硕,总分289分 +11 @taotao 2026-03-29 11/550 2026-04-02 10:04 by realme321
[考研] 一志愿北京科技,085601总分305求调剂 +9 半生瓜! 2026-04-01 11/550 2026-04-02 08:28 by Wang200018
[考研] 一志愿南昌大学324求调剂 +12 hanamiko 2026-03-27 12/600 2026-04-01 22:34 by 运气yunqi
[考研] 食品学硕362求调剂 +3 xuanxianxian 2026-04-01 3/150 2026-04-01 21:05 by 啊李999
[考研] 332求调剂 +8 Lyy930824@ 2026-03-29 8/400 2026-04-01 18:40 by 千钧澄玉宇啊
[考研] 材料专硕306英一数二 +7 z1z2z3879 2026-03-31 7/350 2026-04-01 14:50 by ZXlzxl0425
[考研] 254材料与化工求调剂 +3 翰冬林楠 2026-03-30 4/200 2026-03-31 17:53 by yishunmin
[考研] 求调剂 生物学 377分 +6 zzll03 2026-03-31 6/300 2026-03-31 17:33 by 唐沐儿
[考研] 334求调剂 +7 Trying] 2026-03-31 7/350 2026-03-31 12:33 by 无际的草原
[考研] 328求调剂 +8 嗯滴的基本都 2026-03-27 8/400 2026-03-30 17:20 by Wang200018
[考研] 071010 323 分求调剂 +3 Baekzhy 2026-03-27 3/150 2026-03-30 14:24 by andresqi
[考研] 298调剂 +3 jiyingjie123 2026-03-27 3/150 2026-03-27 11:57 by wxiongid
信息提示
请填处理意见