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

登录注册

应《网络安全法》要求，自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

回复此楼

» 猜你喜欢

308求调剂已经有3人回复
321求调剂已经有11人回复
一志愿中海洋材料工程专硕330分求调剂已经有8人回复
330求调剂已经有4人回复
304求调剂已经有7人回复
一志愿南京理工大学085701资源与环境302分求调剂已经有6人回复
南京大学化学376求调剂已经有3人回复
考研调剂求学校推荐已经有5人回复
323求调剂已经有3人回复
317求调剂已经有12人回复

letsplayagame

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

已阅回复此楼关注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的回帖

查看全部 5 个回答

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的回帖

查看全部 5 个回答

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

[考研] 321求调剂 +9	何润采123 2026-03-18	11/550	2026-03-20 23:19 by JourneyLucky
[考研] 一志愿武汉理工材料工程专硕调剂 +9	Doleres 2026-03-19	9/450	2026-03-20 22:36 by JourneyLucky
[考研] 一志愿中南化学（0703）总分337求调剂 +8	niko- 2026-03-19	9/450	2026-03-20 21:57 by luoyongfeng
[考研] 290求调剂 +7	^O^乜 2026-03-19	7/350	2026-03-20 21:43 by JourneyLucky
[考研] 求调剂一志愿南京航空航天大学289分 +3	@taotao 2026-03-19	3/150	2026-03-20 21:34 by JourneyLucky
[考研] 274求调剂 +9	S.H1 2026-03-18	9/450	2026-03-20 20:54 by zhukairuo
[考研] 一志愿华中农业071010，总分320求调剂 +3	困困困困坤坤 2026-03-20	3/150	2026-03-20 20:38 by 学员8dgXkO
[考研] 319求调剂 +3	小力气珂珂 2026-03-20	3/150	2026-03-20 19:47 by JourneyLucky
[考研] 085700资源与环境308求调剂 +11	墨墨漠 2026-03-18	12/600	2026-03-20 19:43 by 丁丁*
[考研] 298-一志愿中国农业大学-求调剂 +9	手机用户 2026-03-17	9/450	2026-03-20 14:24 by 无懈可击111
[考研] 085600材料与化工求调剂 +6	绪幸与子 2026-03-17	6/300	2026-03-19 13:27 by houyaoxu
[考研] 346求调剂[0856] +3	WayneLim327 2026-03-16	6/300	2026-03-19 11:21 by WayneLim327
[考研] 材料专业求调剂 +5	hanamiko 2026-03-18	5/250	2026-03-18 20:19 by 星空星月
[考研] 344求调剂 +6	knight344 2026-03-16	7/350	2026-03-18 20:13 by walc
[考研] 0854，计算机类招收调剂 +3	胡辣汤放糖 2026-03-15	6/300	2026-03-18 12:09 by 上岸上岸……..
[考研] 334求调剂 +3	志存高远意在机� 2026-03-16	3/150	2026-03-18 08:34 by lm4875102
[考研] 277调剂 +5	自由煎饼果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[考研] 318求调剂 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu
[考研] 080500，材料学硕302分求调剂学校 +4	初识可乐 2026-03-14	5/250	2026-03-14 21:08 by peike
[考研] 297一志愿上交085600求调剂 +5	指尖八千里 2026-03-14	5/250	2026-03-14 17:26 by a不易

24小时热门版块排行榜

qq642867921

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

» 猜你喜欢

dingd

【答案】应助回帖

wake_up

【答案】应助回帖

qq642867921

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