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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » matlab求解一个 三元 非线性方程组啊
51/1返回列表
查看: 2078  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

pl03818948

新虫 (初入文坛)

[求助] matlab求解一个 三元 非线性方程组啊 已有1人参与

急求急求
f(1)=0.1724239171*sin(x1/180*pi)+cow(x1/180*pi)*0.9832199726*cos((15*12.7-x2-60)/180*pi)-x3/(x3^2+1.149625826^2)^0.5;
f(2)=0.1724239171*sin(x1/180*pi)+cow(x1/180*pi)*0.9832199726*cos((15*12.75-x2-60)/180*pi)-x3/(x3^2+1.149625826^2)^0.5;
f(3)=0.1724239171*sin(x1/180*pi)+cow(x1/180*pi)*0.9832199726*cos((15*12.8-x2-60)/180*pi)-x3/(x3^2+1.149625826^2)^0.5;
知道很简单，但是这里我不了解，不胜感激
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-09-11 23:01:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pl03818948

新虫 (初入文坛)
问题是求解过程中总是需要提供大概估值，但是我根本无从下手
赞一下 回复此楼
2楼2015-09-12 02:31:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
pl03818948: 金币+10, ★★★★★最佳答案 2015-09-15 01:12:17
代码如下：
CODE:
function as
format long
clear;clc
u0=[0 0 0];
[u fval]=fsolve(@eqs,u0);
[u ;fval]'


function f=eqs(u)
x1=u(1);
x2=u(2);
x3=u(3);

f(1)=0.1724239171*sin(x1/180*pi)+cos(x1/180*pi)*0.9832199726*cos((15*12.7-x2-60)/180*pi)-x3/(x3^2+1.149625826^2)^0.5;
f(2)=0.1724239171*sin(x1/180*pi)+cos(x1/180*pi)*0.9832199726*cos((15*12.75-x2-60)/180*pi)-x3/(x3^2+1.149625826^2)^0.5;
f(3)=0.1724239171*sin(x1/180*pi)+cos(x1/180*pi)*0.9832199726*cos((15*12.8-x2-60)/180*pi)-x3/(x3^2+1.149625826^2)^0.5;

x                              方程残差
-89.999999999960266  -0.000000000068213
   0.006138053265927  -0.000000000068220
  -0.201236946236002  -0.000000000068227

sin cos是周期性函数，存在多解，上述结果只是其中一个。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
3楼2015-09-12 07:52:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sinyfiny

新虫 (初入文坛)
引用回帖:
3楼: Originally posted by 月只蓝 at 2015-09-12 07:52:01
代码如下：
function as
format long
clear;clc
u0=;
=fsolve(@eqs,u0);
'


function f=eqs(u)
x1=u(1);
x2=u(2);
x3=u(3);

f(1)=0.1724239171*sin(x1/180*pi)+cos(x1/180*pi)*0.9832199726*cos((1 ...

能不能帮我也看看，我没得权限，折腾一晚上了
赞一下(1人) 回复此楼
4楼2015-09-12 08:10:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

引用回帖:
4楼: Originally posted by sinyfiny at 2015-09-12 08:10:54
能不能帮我也看看，我没得权限，折腾一晚上了...

发个帖子到计算模拟大区--计算模拟板块。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
5楼2015-09-12 08:27:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 pl03818948 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭