24小时热门版块排行榜    

查看: 1601  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

周宁_fish

新虫 (初入文坛)

[求助] matlab解方程组求助攻 已有2人参与

方程组如下,求解x  y  z
f1=0.00124*y+0.114094*z-0.1*y+0.1*y*exp(-900*x)=0
f2=0.00549*y+0.648922*z-0.1*y+0.1*y*exp(-1800*x)=0
f3=0.00607*y+0.750408*z-0.1*y+0.1*y*exp(-2700*x)=0
f4=0.00821*y+1.251459471*z-0.1*y+0.1*y*exp(-3600*x)=0
还有f5...
用matlab的fsolve函数求解的时候3个未知数只能编辑三个方程,否则会提示Warning: Trust-region-dogleg algorithm of FSOLVE cannot handle non-square systems; using Levenberg-Marquardt algorithm instead. 用前三个方程来解方程,x的解和所选的初值完全相同。请问这是什么原因,程序如下请问解出来的解对不对。
function q=myfun(p)
x=p(1);
y=p(2);
z=p(3);
q(1)=0.00124*y+0.114094*z-0.1*y+0.1*y*exp(-900*x);
q(2)=0.00549*y+0.648922*z-0.1*y+0.1*y*exp(-1800*x);
q(3)=0.00607*y+0.750408*z-0.1*y+0.1*y*exp(-2700*x);
x = fsolve('myfun',[10,100000,3000]',optimset('Largescale','off'))

在网上搜的 Levenberg-Marquardt 方法来解这个方程组,程序如下
f()={"0.00124*y+0.114094*z-0.1*y+0.1*y*exp(-900*x)","0.00549*y+0.648922*z-0.1*y+0.1*y*exp(-1800*x)","0.00607*y+0.750408*z-0.1*y+0.1*y*exp(-2700*x)"}
x()={"x","y","z"}
x0()={10,100000,3000}
d()={0,0,0}
erro=1e-8
Loopn=3000
nonlin_LM(f,x,x0,d,erro,Loopn)
运行时显示Error: File: Untitled4.m Line: 1 Column: 2
An indexing expression on the left side of an assignment must have at least one subscript.

我对matlab不懂,请大神帮忙解一下这个方程组,望附上程序,谢谢!
回复此楼

» 猜你喜欢

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

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

戴钢盔的猪头

木虫 (知名作家)

【答案】应助回帖

哦不好意思,之前看的不够仔细,实际上你的方程组应该有无数组解。爱渴死可以取任意值,歪和贼只要都是零就能满足了。还是那句话,你要定义三个变量的可行取值范围才有可能计算出唯一解。

[ 发自手机版 http://muchong.com/3g ]
5楼2014-02-18 09:43:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
方程没问题吗?如果个方程右边均为0,那就只有解:
y: 0
z: 0
x: 0
2楼2014-02-17 14:49:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

周宁_fish

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by dingd at 2014-02-17 14:49:14
方程没问题吗?如果个方程右边均为0,那就只有解:
y: 0
z: 0
x: 0

方程没有问题    怎么会只有0解   方程组的系数矩阵肯定不是0,那肯定有非零解的啊
3楼2014-02-17 15:54:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

戴钢盔的猪头

木虫 (知名作家)

【答案】应助回帖

感谢参与,应助指数 +1
你的问题没表述清楚吧?首先二楼说的全为零的解肯定是可行解,其次你求解超定问题,需要定义成求最小二乘解,然后你需要对变量做不等式约束,否则只能找到全零解。信赖域折线法可以求解超定最小二乘问题的,但是你的问题提的不对。

[ 发自手机版 http://muchong.com/3g ]
4楼2014-02-18 01:29:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 281求调剂(0805) +4 烟汐忆海 2026-03-16 11/550 2026-03-18 11:57 by djl2006
[考研] 307求调剂 +3 冷笙123 2026-03-17 3/150 2026-03-18 09:55 by macy2011
[考研] 材料专硕306英一数二 +9 z1z2z3879 2026-03-16 12/600 2026-03-18 08:53 by zhukairuo
[考研] 341求调剂 +5 捣蛋猪猪 2026-03-11 7/350 2026-03-17 19:09 by 捣蛋猪猪
[考研] 085600材料与化工求调剂 +5 绪幸与子 2026-03-17 5/250 2026-03-17 16:40 by laoshidan
[考研] 289求调剂 +6 步川酷紫123 2026-03-11 6/300 2026-03-17 10:23 by Sammy2
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[考研] 0703化学调剂 +6 妮妮ninicgb 2026-03-15 9/450 2026-03-16 16:40 by houyaoxu
[考研] 085600材料与化工 求调剂 +13 enenenhui 2026-03-13 14/700 2026-03-16 15:19 by 了了了了。。
[考研] 294求调剂 +3 Zys010410@ 2026-03-13 4/200 2026-03-15 10:59 by zhq0425
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 281求调剂 +9 Koxui 2026-03-12 11/550 2026-03-13 20:50 by Koxui
[考研] 【考研调剂求收留】 +3 Ceciilia 2026-03-11 3/150 2026-03-13 20:18 by JourneyLucky
[硕博家园] 085600 260分求调剂 +3 天空还下雨么 2026-03-13 5/250 2026-03-13 18:46 by 天空还下雨么
[考研] 314求调剂 +7 无懈可击的巨人 2026-03-12 7/350 2026-03-13 15:40 by JourneyLucky
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 333求调剂 +3 152697 2026-03-12 4/200 2026-03-13 07:08 by Iveryant
信息提示
请填处理意见