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

登录注册

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

回复此楼

» 猜你喜欢

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 个回答

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

[考研] 求调剂，一志愿:南京航空航天大学大学，080500材料科学与工程学硕，总分289分 +4	@taotao 2026-03-19	4/200	2026-03-20 22:14 by JourneyLucky
[考研] 一志愿苏州大学材料求调剂，总分315（英一） +5	sbdksD 2026-03-19	5/250	2026-03-20 22:10 by luoyongfeng
[考研] 中南大学化学学硕337求调剂 +3	niko- 2026-03-19	6/300	2026-03-20 21:58 by luoyongfeng
[考研] 一志愿北京化工大学0703化学318分，有科研经历，求调剂 +4	一瓶苯甲酸 2026-03-14	4/200	2026-03-20 20:36 by fen_rao
[考研] 一志愿吉林大学材料学硕321求调剂 +11	Ymlll 2026-03-18	15/750	2026-03-20 19:40 by 丁丁*
[考研] 086500 325 求调剂 +3	领带小熊 2026-03-19	3/150	2026-03-20 18:38 by 尽舜尧1
[考研] 广西大学家禽遗传育种课题组2026年硕士招生（接收计算机专业调剂） +3	123阿标 2026-03-17	3/150	2026-03-20 15:58 by 飞行琦
[考研] 08工学调剂 +5	用户573181 2026-03-20	5/250	2026-03-20 15:47 by xia_2003
[考研] 304求调剂 +5	曼殊2266 2026-03-18	5/250	2026-03-20 09:00 by ZHANG0tao
[考研] 332求调剂 +3	ydfyh 2026-03-17	3/150	2026-03-19 10:14 by 功夫疯狂
[考研] 材料专硕306英一数二 +10	z1z2z3879 2026-03-16	13/650	2026-03-18 14:20 by 007_lilei
[考研] 304求调剂 +12	小熊joy 2026-03-14	13/650	2026-03-18 12:34 by Linda Hu
[考研] 生物学071000 329分求调剂 +3	我爱生物生物爱� 2026-03-17	3/150	2026-03-18 10:12 by macy2011
[考研] 293求调剂 +11	zjl的号 2026-03-16	16/800	2026-03-18 08:10 by zhukairuo
[考研] 277调剂 +5	自由煎饼果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[考研] 290求调剂 +3	p asserby. 2026-03-15	4/200	2026-03-17 16:35 by wangkm
[考研] 302求调剂 +4	小贾同学123 2026-03-15	8/400	2026-03-17 10:33 by 小贾同学123
[考研] 304求调剂 +5	素年祭语 2026-03-15	5/250	2026-03-16 17:00 by 我的船我的海
[考研] 0856专硕279求调剂 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015
[考研] 289求调剂 +4	这么名字咋样 2026-03-14	6/300	2026-03-14 18:58 by userper

24小时热门版块排行榜

qq642867921

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

» 猜你喜欢

dingd

【答案】应助回帖

wake_up

【答案】应助回帖

qq642867921

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