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

登录注册

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

24小时热门版块排行榜

返回列表

qq642867921

银虫 (小有名气)

应助: 2 (幼儿园)
金币: 273.1
红花: 4
帖子: 146
在线: 30.5小时
虫号: 3860444
注册: 2015-05-10
性别: GG
专业: 钢铁冶金

[求助] matlab解多元方程组程序，以前我算出来了结果，现在怎么有问题，计算不出来已有2人参与

求大神，非常感谢！！

程序如下：
function f=CaOSiO2MgOAl2O3FeOCr2O329(N)
global number;
number=number+1;
a1=0.03;a2=0.12;a3=0.3;b1=0.35;b2=0.05;b3=0.15;%a1=Cr2O3,a2=Al2O3,a3=SiO2,b1=CaO,b2=MgO,b3=FeO;
f(1)=N(1)+N(2)+N(3)+N(4)+N(5)+N(6)+N(7)+N(8)+N(9)+N(10)+N(11)+N(12)+N(13)+N(14)+N(15)+N(16)+N(17)+N(18)+N(19)+N(20)+N(21)+N(22)+N(23)+N(24)+N(25)+N(26)+N(27)+N(28)+N(29)-1;
f(2)=a1*(0.5*N(1)+N(7)+2*N(8)+3*N(9)+N(10)+N(11)+3*N(12)+N(13)+12*N(14)+N(22)+2*N(24)+N(25)+3*N(26)+N(27)+2*N(28)+N(29))-b1*(N(4)+N(21)+N(22)+N(23));
f(3)=a1*(0.5*N(2)+N(15)+2*N(16)+N(17)+N(23)+N(26)+N(27)+N(28)+N(29))-b2*(N(4)+N(21)+N(22)+N(23));
f(4)=a1*(0.5*N(3)+N(19)+2*N(20)+N(21))-b3*(N(4)+N(21)+N(22)+N(23));
f(5)=a1*(N(5)+N(10)+2*N(11)+N(12)+6*N(13)+7*N(14)+N(17)+3*N(18)+N(19)+N(24)+N(25))-a2*(N(4)+N(21)+N(22)+N(23));
f(6)=a1*(N(6)+N(7)+N(8)+N(9)+N(15)+N(16)+2*N(18)+N(20)+N(24)+2*N(25)+2*N(26)+N(27)+2*N(28)+2*N(29))-a3*(N(4)+N(21)+N(22)+N(23));
f(7)=N(7)-659.1868765*N(1)*N(6);
f(8)=N(8)-48840.50342*N(1)^2*N(6);
f(9)=N(9)-6410.811308*N(1)^3*N(6);
f(10)=N(10)-30.59284903*N(1)*N(5);
f(11)=N(11)-62.92544394*N(1)*N(5)^2;
f(12)=N(12)-43.76471764*N(1)^3*N(5);
f(13)=N(13)-262.9506479*N(1)*N(5)^6;
f(14)=N(14)-55613987772*N(1)^12*N(5)^7;
f(15)=N(15)-6.723640959*N(2)*N(6);
f(16)=N(16)-44.56734705*N(2)^2*N(6);
f(17)=N(17)-12.59153449*N(2)*N(5);
f(18)=N(18)-14.10214903*N(5)^3*N(6)^2;
f(19)=N(19)-0.0000927114*N(3)*N(5);
f(20)=N(20)-4.209805425*N(3)^2*N(6);
f(21)=N(21)-6.753781389*N(3)*N(4);
f(22)=N(22)-5324.185271*N(1)*N(4);
f(23)=N(23)-6.684150735*N(2)*N(4);
f(24)=N(24)-170498.3476*N(1)^2*N(5)*N(6);
f(25)=N(25)-59100.22934*N(1)*N(5)*N(6)^2;
f(26)=N(26)-31859855.35*N(1)^3*N(2)*N(6)^2;
f(27)=N(27)-1907.157919*N(1)*N(2)*N(6);
f(28)=N(28)-0.053804824*N(1)^2*N(2)*N(6)^2;
f(29)=N(29)-3593.935175*N(1)*N(2)*N(6)^2;

N0=[0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1];
N=fsolve(@myfun5,N0)

1.png

matlab解多元方程组程序，以前我算出来了结果，现在怎么有问题，计算不出来-1

360反馈意见截图18241017123144138.png

回复此楼

» 猜你喜欢

材料与化工专硕调剂已经有6人回复
086500 325 求调剂已经有3人回复
材料080500调剂求收留已经有7人回复
一志愿南昌大学，327分，材料与化工085600 已经有6人回复
环境工程调剂已经有9人回复
工科材料085601 279求调剂已经有9人回复
287求调剂已经有5人回复
22408 344分求调剂一志愿华电计算机技术已经有3人回复
289求调剂已经有3人回复
317求调剂已经有10人回复

letsplayagame

1楼 2015-11-04 22:55:55

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

wake_up

木虫 (著名写手)

应助: 124 (高中生)
金币: 4011.5
散金: 1654
红花: 24
帖子: 1113
在线: 424.2小时
虫号: 670119
注册: 2008-12-08
性别: GG
专业: 医学图像数据处理与分析

【答案】应助回帖

感谢参与，应助指数 +1

眼睛都看花了。你的自定义函数和@myfun名对应不起来。function一般也要化成等于0的形式。还有那些f()要起到什么作用？你有单步调试过没有？

发自小木虫Android客户端

赞一下

回复此楼

扫地扫地扫心地，心地不扫空扫地，人人都把心地扫，世上无处不净地．

2楼2015-11-05 01:03:57

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

qq642867921

银虫 (小有名气)

应助: 2 (幼儿园)
金币: 273.1
红花: 4
帖子: 146
在线: 30.5小时
虫号: 3860444
注册: 2015-05-10
性别: GG
专业: 钢铁冶金

引用回帖:

2楼: Originally posted by wake_up at 2015-11-05 01:03:57
眼睛都看花了。你的自定义函数和@myfun名对应不起来。function一般也要化成等于0的形式。还有那些f()要起到什么作用？你有单步调试过没有？

函数名称是弄错了，但是计算不出来，程序出错，以前我就计算过了的，可以算出结果。我的方程就是已知a1,a2,a3,b1,b2,b3然后计算f1-f29一共29个方程求出N1-N29，我另外一个很简单的程序都可以计算：
function f=myfun9(N)
global number;
number=number+1;
a1=0.45;a2=0.133333333;a3=0.35;b=0.066666667;%a1=SiO2,a2=Al2O3,a3=Cr2O3,b=MgO;
f(1)=N(1)+N(2)+N(3)+N(4)+N(5)+N(6)+N(7)+N(8)+N(9)-1;
f(2)=a3*(N(3)+N(7)+3*N(8))-a2*(N(4)+N(9));
f(3)=a3*(N(2)+N(5)+N(6)+2*N(8))-a1*(N(4)+N(9));
f(4)=a3*(0.2*N(1)+N(5)+2*N(6)+N(7)+N(9))-b*(N(4)+N(9));
f(5)=N(5)-6.680601877*N(1)*N(2);
f(6)=N(6)-44.56734705*N(1)*N(1)*N(2);
f(7)=N(7)-12.59153449*N(1)*N(3);
f(8)=N(8)-14.10214903*N(2)*N(2)*N(3)*N(3)*N(3);
f(9)=N(9)-6.684150735*N(1)*N(4);

N0=[0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1];
N=fsolve(@myfun9,N0)

赞一下

回复此楼

letsplayagame

3楼2015-11-05 09:31:56

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

wake_up

木虫 (著名写手)

应助: 124 (高中生)
金币: 4011.5
散金: 1654
红花: 24
帖子: 1113
在线: 424.2小时
虫号: 670119
注册: 2008-12-08
性别: GG
专业: 医学图像数据处理与分析

内容已删除

回复此楼

扫地扫地扫心地，心地不扫空扫地，人人都把心地扫，世上无处不净地．

4楼2015-11-05 10:09:14

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

dingd

铁杆木虫 (职业作家)

应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.7小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
qq642867921: 金币+10, ★★★很有帮助 2015-11-06 23:23:57

1stOpt试试，简单快捷的·多：

CODE:

Constant

a1=0.03,a2=0.12,a3=0.3,b1=0.35,b2=0.05,b3=0.15;//,%a1=Cr2O3,a2=Al2O3,a3=SiO2,b1=CaO,b2=MgO,b3=FeO;

Parameter N(29);

Function

N(1)+N(2)+N(3)+N(4)+N(5)+N(6)+N(7)+N(8)+N(9)+N(10)+N(11)+N(12)+N(13)+N(14)+N(15)+N(16)+N(17)+N(18)+N(19)+N(20)+N(21)+N(22)+N(23)+N(24)+N(25)+N(26)+N(27)+N(28)+N(29)-1;

a1*(0.5*N(1)+N(7)+2*N(8)+3*N(9)+N(10)+N(11)+3*N(12)+N(13)+12*N(14)+N(22)+2*N(24)+N(25)+3*N(26)+N(27)+2*N(28)+N(29))-b1*(N(4)+N(21)+N(22)+N(23));

a1*(0.5*N(2)+N(15)+2*N(16)+N(17)+N(23)+N(26)+N(27)+N(28)+N(29))-b2*(N(4)+N(21)+N(22)+N(23));

a1*(0.5*N(3)+N(19)+2*N(20)+N(21))-b3*(N(4)+N(21)+N(22)+N(23));

a1*(N(5)+N(10)+2*N(11)+N(12)+6*N(13)+7*N(14)+N(17)+3*N(18)+N(19)+N(24)+N(25))-a2*(N(4)+N(21)+N(22)+N(23));

a1*(N(6)+N(7)+N(8)+N(9)+N(15)+N(16)+2*N(18)+N(20)+N(24)+2*N(25)+2*N(26)+N(27)+2*N(28)+2*N(29))-a3*(N(4)+N(21)+N(22)+N(23));

N(7)-659.1868765*N(1)*N(6);

N(8)-48840.50342*N(1)^2*N(6);

N(9)-6410.811308*N(1)^3*N(6);

N(10)-30.59284903*N(1)*N(5);

N(11)-62.92544394*N(1)*N(5)^2;

N(12)-43.76471764*N(1)^3*N(5);

N(13)-262.9506479*N(1)*N(5)^6;

N(14)-55613987772*N(1)^12*N(5)^7;

N(15)-6.723640959*N(2)*N(6);

N(16)-44.56734705*N(2)^2*N(6);

N(17)-12.59153449*N(2)*N(5);

N(18)-14.10214903*N(5)^3*N(6)^2;

N(19)-0.0000927114*N(3)*N(5);

N(20)-4.209805425*N(3)^2*N(6);

N(21)-6.753781389*N(3)*N(4);

N(22)-5324.185271*N(1)*N(4);

N(23)-6.684150735*N(2)*N(4);

N(24)-170498.3476*N(1)^2*N(5)*N(6);

N(25)-59100.22934*N(1)*N(5)*N(6)^2;

N(26)-31859855.35*N(1)^3*N(2)*N(6)^2;

N(27)-1907.157919*N(1)*N(2)*N(6);

N(28)-0.053804824*N(1)^2*N(2)*N(6)^2;

N(29)-3593.935175*N(1)*N(2)*N(6)^2;

n1: 0.00740569433075613
n2: 0.0334451306656032
n3: 0.373025926308175
n4: 0.00107151740885337
n5: 0.0707709740557123
n6: 0.0359206838653582
n7: 0.175355314041228
n8: 0.0962179928931035
n9: 9.35308677531875E-5
n10: 0.0160339630815332
n11: 0.00233400841376471
n12: 1.25799584019101E-6
n13: 2.4466443104525E-7
n14: -5.57411935550417E-11
n15: 0.00807759375205905
n16: 0.00179071789055471
n17: 0.0298034629881549
n18: 6.44970151861671E-6
n19: 2.44752591688522E-6
n20: 0.0210418855773427
n21: 0.00269951190752645
n22: 0.0422491692366895
n23: 0.000239540173299834
n24: 0.023771222055233
n25: 0.0399667942132903
n26: 0.000558422925487669
n27: 0.0169679717120742
n28: 1.05593742617914E-10
n29: 0.00114857169843853

赞一下(1人)

回复此楼

5楼2015-11-05 11:14:28

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

相关版块跳转我要订阅楼主学员HBFXFr 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 工科材料085601 279求调剂 +7	困于星晨 2026-03-17	9/450	2026-03-20 17:38 by 无懈可击111
[考研] 0856调剂，是学校就去 +8	sllhht 2026-03-19	9/450	2026-03-20 14:25 by 无懈可击111
[考研] 求调剂 +4	Mqqqqqq 2026-03-19	4/200	2026-03-20 14:15 by 星空星月
[考研] 295复试调剂 +7	简木ChuFront 2026-03-19	7/350	2026-03-20 14:11 by 星空星月
[考研] 279分求调剂一志愿211 +9	chaojifeixia 2026-03-19	10/500	2026-03-20 12:29 by lature00
[考研] 招收调剂硕士 +4	lidianxing 2026-03-19	12/600	2026-03-20 12:25 by lidianxing
[考研] 085700资源与环境308求调剂 +10	墨墨漠 2026-03-18	11/550	2026-03-20 11:57 by 学员8dgXkO
[考研] 【考研调剂】化学专业 281分，一志愿四川大学，诚心求调剂 +6	吃吃吃才有意义 2026-03-19	6/300	2026-03-20 10:47 by 尽舜尧1
[考研] 材料学硕318求调剂 +5	February_Feb 2026-03-19	5/250	2026-03-19 23:51 by 23Postgrad
[考研] 一志愿苏州大学材料求调剂，总分315（英一） +3	sbdksD 2026-03-19	3/150	2026-03-19 23:21 by fmesaito
[考研] 266求调剂 +5	阳阳哇塞 2026-03-14	10/500	2026-03-19 15:08 by 阳阳哇塞
[考研] 化学求调剂 +3	临泽境llllll 2026-03-17	4/200	2026-03-19 13:59 by houyaoxu
[考研] 一志愿中海洋材料工程专硕330分求调剂 +7	小材化本科 2026-03-18	7/350	2026-03-19 10:46 by Linda Hu
[考研] 354求调剂 +4	Tyoumou 2026-03-18	7/350	2026-03-18 21:45 by Tyoumou
[考研] 268求调剂 +6	简单点0 2026-03-17	6/300	2026-03-18 09:04 by 无际的草原
[考研] 材料工程专硕274一志愿211求调剂 +6	薛云鹏 2026-03-15	6/300	2026-03-17 11:05 by 学员h26Tkc
[考研] 275求调剂 +4	太阳花天天开心 2026-03-16	4/200	2026-03-17 10:53 by 功夫疯狂
[考研] 一志愿，福州大学材料专硕339分求调剂 +3	木子momo青争 2026-03-15	3/150	2026-03-17 07:52 by laoshidan
[考研] 0703 物理化学调剂 +3	我可以上岸的对� 2026-03-13	5/250	2026-03-16 10:50 by 我可以上岸的对�
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3	嘉年新程 2026-03-15	3/150	2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿

24小时热门版块排行榜

qq642867921

[求助] matlab解多元方程组程序，以前我算出来了结果，现在怎么有问题，计算不出来 已有2人参与

» 猜你喜欢

wake_up

【答案】应助回帖

qq642867921

wake_up

dingd

【答案】应助回帖

[求助] matlab解多元方程组程序，以前我算出来了结果，现在怎么有问题，计算不出来已有2人参与