版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

samhx

金虫 (小有名气)

应助: 1 (幼儿园)
金币: 660.4
散金: 5
红花: 1
帖子: 88
在线: 115.4小时
虫号: 704586
注册: 2009-02-20
专业: 金属材料表面科学与工程

[求助] 求助：用matlab解方程组，如下所示，请牛人指教！

求F1、F2、F3，各位大牛，请指教！

1.png

回复此楼

1楼2013-07-21 21:24:32

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

专家经验: +1059
应助: 1712 (讲师)
贵宾: 8.888
金币: 68120.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
jjdg: 金币+1, 感谢参与 2013-07-22 00:49:15
samhx: 金币+5, ★★★很有帮助, 太羡慕牛人了，大赞 2013-07-25 21:12:57

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

function f=nonl(x)
F1=x(1);
F2=x(2);
F3=x(3);
f(1)=F1^(2/3)+3.75*F2^(2/3)-2.75*F3^(2/3);
f(2)=F2-3.75*F1+196.8;
f(3)=F3+2.75*F1-448.9;

结果：
x =

80.6164 105.5114 227.2049

不知道这样的解是否符合实际要求，不行的话，可以试着调一下初值。

赞一下(1人)

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

2楼2013-07-21 22:49:32

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zyxzyh921

禁虫 (初入文坛)

★ ★ ★ ★ ★
感谢参与，应助指数 +1
samhx: 金币+5, ★★★很有帮助, 非常羡慕啊，搞学术的好人 2013-07-25 21:13:29

本帖内容被屏蔽

3楼2013-07-22 10:19:52

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

专家经验: +1059
应助: 1712 (讲师)
贵宾: 8.888
金币: 68120.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

引用回帖:

3楼: Originally posted by zyxzyh921 at 2013-07-22 10:19:52
直接写：
=solve('x1^2/3+3.75*x2^2/3-2.75*x3^2/3=0','x2=3.75*x1-196.8','x3=-2.75*x1+448.9')
得到结果：
x1 =

93.993030742840714792347444854146
-132.08373283012724230657894390538

x2 =
...

方程组中的第一个方程，貌似写错了吧：
“x1^2/3+3.75*x2^2/3-2.75*x3^2/3=0”
按照原方程，应该这样输入吧：
“x1^(2/3)+3.75*x2^(2/3)-2.75*x3^(2/3)=0”

赞一下(2人)

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

4楼2013-07-22 10:37:26

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖