24小时热门版块排行榜

返回列表

amorphous

铁杆木虫 (著名写手)

小木虫之超级奶爸

应助: 7 (幼儿园)
贵宾: 0.15
金币: 8567.3
红花: 1
帖子: 1290
在线: 160.1小时
虫号: 431986
注册: 2007-08-12
专业: 催化化学

[求助] matlab解非线性方程组，求助！

帮忙看一下下列程序,如何求得实数解.

function q=zxxfun(n)
logkp=[11.3,6.2,14.8,12.4,19.5,17.1,22.5,2.4,9.8];
p=1;
p0=1;
nn10=1;
nn20=2;
q(1)=logkp(1)-(log10(n(4))+log10(n(5))+2*log10(sum(n))-3*log10(n(2))-log10(n(1))-2*log10(p/p0));
q(2)=logkp(2)-(2*log10(n(4))+log10(n(6))+3*log10(sum(n))-4*log10(n(2))-2*log10(n(1))-3*log10(p/p0));
q(3)=logkp(3)-(2*log10(n(4))+log10(n(7))+4*log10(sum(n))-5*log10(n(2))-2*log10(n(1))-4*log10(p/p0));
q(4)=logkp(4)-(3*log10(n(4))+log10(n(8))+5*log10(sum(n))-6*log10(n(2))-3*log10(n(1))-5*log10(p/p0));
q(5)=logkp(5)-(3*log10(n(4))+log10(n(9))+6*log10(sum(n))-7*log10(n(2))-3*log10(n(1))-6*log10(p/p0));
q(6)=logkp(6)-(4*log10(n(4))+log10(n(10))+7*log10(sum(n))-8*log10(n(2))-4*log10(n(1))-7*log10(p/p0));
q(7)=logkp(7)-(4*log10(n(4))+log10(n(11))+8*log10(sum(n))-9*log10(n(2))-4*log10(n(1))-8*log10(p/p0));
q(8)=logkp(8)-(log10(n(2))+log10(n(3))-log10(n(1))-log10(n(4)));
q(9)=logkp(9)-(log10(n(3))+log10(sum(n))-2*log10(n(1))-log10(p/p0));
q(10)=nn10-(n(1)+2*n(3)+n(4));
q(11)=2*nn20-(2*n(2)+2*n(4)+4*n(5)+4*n(6)+6*n(7)+6*n(8)+8*n(9)+8*n(10)+10*n(11));

clear

options=optimset('Display','iter');
n=fsolve(@zxxfun,[0.5 0.5 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1,0.1],options);

回复此楼

» 猜你喜欢

寻求一种能扛住强氧化性腐蚀性的容器密封件已经有7人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有8人回复
申请2026年博士已经有6人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有5人回复
2025冷门绝学什么时候出结果已经有7人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有6人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有7人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复

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

求助！matlab用fsolve函数求解非线性方程组的问题！已经有19人回复
matlab求解方程组“” 已经有4人回复
matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较已经有4人回复
Matlab求解偏微分方程组已经有22人回复
matlab求解非线性方程组已经有16人回复
matlab数值求解边界条件微分方程组已经有7人回复
高金求助matlab解微分方程组已经有12人回复
求matlab求解一个方程组的问题已经有6人回复
求高人指点用matlab求解非线性方程组，解决了追加100金币；已经有11人回复
matlab的fsove 命令求解非线性方程组已经有6人回复
【求助】MATLAB 解方程组【已解决】已经有7人回复
【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题已经有21人回复
【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！已经有18人回复
【求助】matlab怎么求解偏微分方程组啊，先谢谢了已经有13人回复
【求助】求教matlab解非线性方程组已经有9人回复

1楼2012-05-24 22:21:13

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

dingd

铁杆木虫 (职业作家)

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

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
amorphous: 金币+100, ★★★★★最佳答案, 谢谢。 2012-05-25 15:16:06

不知道理解对否，用1stOpt计算：

CODE:

Constant logkp=[11.3,6.2,14.8,12.4,19.5,17.1,22.5,2.4,9.8];

Constant p=1,p0=1,nn10=1,nn20=2;

ConstStr f=sum(i=1:11)(n[i]);

Function logkp[1]-(log(n4)+log(n5)+2*log(f)-3*log(n2)-log(n1)-2*log(p/p0));

         logkp[2]-(2*log(n4)+log(n6)+3*log(f)-4*log(n2)-2*log(n1)-3*log(p/p0));

         logkp[3]-(2*log(n4)+log(n7)+4*log(f)-5*log(n2)-2*log(n1)-4*log(p/p0));

         logkp[4]-(3*log(n4)+log(n8)+5*log(f)-6*log(n2)-3*log(n1)-5*log(p/p0));

         logkp[5]-(3*log(n4)+log(n9)+6*log(f)-7*log(n2)-3*log(n1)-6*log(p/p0));

         logkp[6]-(4*log(n4)+log(n10)+7*log(f)-8*log(n2)-4*log(n1)-7*log(p/p0));

         logkp[7]-(4*log(n4)+log(n11)+8*log(f)-9*log(n2)-4*log(n1)-8*log(p/p0));

         logkp[8]-(log(n2)+log(n3)-log(n1)-log(n4));

         logkp[9]-(log(n3)+log(f)-2*log(n1)-log(p/p0));

         nn10-(n1+2*n3+n4);

         2*nn20-(2*n2+2*n4+4*n5+4*n6+6*n7+6*n8+8*n9+8*n10+10*n11);

n4: 0.818320026230474
n5: 0.585572164827173
n1: 4.65521821183815E-6
n2: 0.0105340963919806
n3: 0.0908376592736609
n6: 1.85174177533351E-13
n7: 5.15897663173426E-7
n8: 8.17642694378187E-17
n9: 7.20355202419344E-12
n10: 1.14168605646827E-21
n11: 2.00691894813261E-18