24小时热门版块排行榜    

查看: 2554  |  回复: 12
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

rita128

新虫 (初入文坛)

[求助] 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;
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sd235634

新虫 (初入文坛)

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

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;
13楼2020-02-19 21:40:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

江心野鹤

专家顾问 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
是代码写错了,把原公式整上来啊

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 455445454.jpg
  • 2013-04-01 09:43:59, 35.32 K
雁去湘江烦寄语:此处有龙潜。
2楼2013-04-01 09:44:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

江心野鹤

专家顾问 (著名写手)

引用回帖:
2楼: Originally posted by 江心野鹤 at 2013-04-01 09:44:00
是代码写错了,把原公式整上来啊

额,怎么贴不上图片???
雁去湘江烦寄语:此处有龙潜。
3楼2013-04-01 09:45:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

江心野鹤

专家顾问 (著名写手)

【答案】应助回帖

★ ★
rita128: 金币+2 2013-04-02 18:47:44
1、parameter
2、exp()
由于没看到你表达式,不好改
雁去湘江烦寄语:此处有龙潜。
4楼2013-04-01 09:49:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +4 2h7du0nuhk 2026-02-07 4/200 2026-02-07 22:24 by ioveiuleh5
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 2h7du0nuhk 2026-02-07 4/200 2026-02-07 22:11 by ioveiuleh5
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 2h7du0nuhk 2026-02-07 5/250 2026-02-07 22:06 by ioveiuleh5
[教师之家] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 2h7du0nuhk 2026-02-07 5/250 2026-02-07 22:04 by ioveiuleh5
[教师之家] 有院领导为了换新车,用横向课题经费买了俩车 +7 瞬息宇宙 2026-02-04 7/350 2026-02-07 21:47 by tfang
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 2h7du0nuhk 2026-02-07 5/250 2026-02-07 21:46 by ioveiuleh5
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 2h7du0nuhk 2026-02-07 6/300 2026-02-07 21:44 by ioveiuleh5
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 2h7du0nuhk 2026-02-07 6/300 2026-02-07 21:31 by ioveiuleh5
[教师之家] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 2h7du0nuhk 2026-02-07 6/300 2026-02-07 21:26 by ioveiuleh5
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 2h7du0nuhk 2026-02-07 6/300 2026-02-07 21:11 by ioveiuleh5
[有机交流] 酰胺脱乙酰基 10+5 chibby 2026-02-03 12/600 2026-02-07 19:29 by 江东闲人
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全 +3 2h7du0nuhk 2026-02-07 3/150 2026-02-07 18:31 by 6ntj4f1rhx
[基金申请] 同年申请2项不同项目,第1个项目里不写第2个项目的信息,可以吗 +4 hitsdu 2026-02-06 4/200 2026-02-07 13:07 by jurkat.1640
[基金申请] 有时候真觉得大城市人没有县城人甚至个体户幸福 +9 苏东坡二世 2026-02-04 10/500 2026-02-07 12:37 by 小毛球
[考博] 天津大学招2026.09的博士生,欢迎大家推荐交流(博导是本人) +4 a793625982 2026-02-05 5/250 2026-02-07 10:57 by a793625982
[公派出国] CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生|MSCA经费充足、排名优 +4 雨念 2026-02-01 6/300 2026-02-06 23:32 by MelissaPon
[基金申请] 面上项目申报 +3 Tide man 2026-02-01 3/150 2026-02-05 22:56 by god_tian
[硕博家园] 博士延得我,科研能力直往上蹿 +7 偏振片 2026-02-02 7/350 2026-02-04 17:36 by 陈氏帝国
[基金申请] 面上基金申报没有其他的参与者成吗 +5 默默的小虫子 2026-01-31 5/250 2026-02-04 11:26 by 学员fOzRO9
[教师之家] 遇见不省心的家人很难过 +18 otani 2026-02-03 22/1100 2026-02-04 11:06 by tangmnt
信息提示
请填处理意见