24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

rita128

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 96.9
帖子: 21
在线: 10.3小时
虫号: 1129918
注册: 2010-10-23
性别: MM
专业: 无机合成和制备化学

[求助] sos!!1stopt解非线性方程组，一直出现常数定义错误的问题！求大神帮助！

第一次用，想解一个非线性方程组，不知道公式输入哪里出了错误，或是其他问题，请教大家帮忙，急迫~~~代码如下，就想解出X1，X2，X3.
//Parameters x1,x2,x3;

//Function 0.0167*exp2*x1-0.0159*expx1+1.0515*exp2*x2-0.1355*expx2+5.7379*exp2*x3-0.4025*expx3=0;
0.1723*exp2*x1-0.0443*expx1+1.0217*exp2*x2-0.1317*expx2+78.124*exp2*x3-5.7252*expx3=0;
0.03849*exp2*x1-0.0099*expx1+0.8214*exp2*x2-0.1058*expx2+4.723*exp2*x3-0.3461*expx3=0;

回复此楼

» 猜你喜欢

请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有6人回复
孩子确诊有中度注意力缺陷已经有14人回复
三甲基碘化亚砜的氧化反应已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复

1楼2013-03-31 18:35:09

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

sd235634

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 58.8
帖子: 21
在线: 20.1小时
虫号: 1833524
注册: 2012-05-25
专业: 运筹与管理

我也遇到了类似错误，不过看了下跟楼主的问题不一样，希望能有大佬也帮忙看看，一直提示常数定义有误

Title "test1";
//Parameters x1,y1[0,40],z1,x2[50,70],y2,z2;
Parameters ca[0,2],fb,pbb,paa,pba,pab,h_aa[0,1],h_ab[0,1],h_ba[0,1],h_bb[0,1],beta_a,beta_b;
//Constant
Constant m_aa = 1067;
Constant m_ab = 203;
Constant m_bb = 2357;
Constant fa = 0.38;
//Variable ;
//ConstStr
ConstStr  A = (h_aa - h_ab)*(h_ba - h_bb)
ConstStr  B = ((2*h_bb - (1-fa) * h_ba)*(h_aa - h_ab) + (2*h_ab - fa*(2*h_aa - h_ab))*(h_ba - h_bb))
ConstStr  C = (2*h_bb*(2*h_ab - fa*(2*h_aa - h_ab)) - 2*fa*h_ab*(h_ba - h_bb) - 2*(1-fa)*h_ba * h_ab)
ConstStr  D = - 4*fa*h_ab*h_bb
ConstStr  K = beta_b/(beta_a+beta_b)
ConstStr  Z = 1 - beta_b/(beta_a+beta_b)
//Function ;
Function
A*ca^3+B*ca^2+C*ca+D=0;
(pbb*((fb*h_bb * K)+ (fa*(1-h_bb)*Z) )) - (fb^2 * h_bb * K )=0;
(paa*((fa*h_aa * Z)+ (fb*(1-h_aa)*K) )) - (fa^2 * h_aa * Z )=0;
fb=1-fa;
beta_a - fa*h_aa/(h_aa*ca + h_ab * (2 - ca)) - fb*h_ba/(h_ba*ca + h_bb*(2-ca))=0;
beta_b - fb*h_bb/(h_ba*ca + h_bb * (2 - ca)) - fa*h_ab/(h_aa*ca + h_ab*(2-ca))=0;
pbb = m_bb/(m_bb+m_ab + m_aa);
paa = m_aa/(m_aa + m_ab + m_bb);
pba = m_ab/(m_aa+m_bb+m_ab);
pab = pba;
h_ab = 1- h_aa;
h_ba = 1- h_bb;