24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2677  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 339求调剂 +8 hanwudada 2026-04-11 9/450 2026-04-12 15:36 by laoshidan
[考研] 求调剂 +16 张番茄不炒蛋 2026-04-10 17/850 2026-04-12 13:58 by 熬夜成!
[考研] 335求调剂 +13 想上岸呀!! 2026-04-12 13/650 2026-04-12 11:37 by 柴小白
[考研] 一志愿华中农微生物,288分,三年实验经历 +11 代fish 2026-04-09 11/550 2026-04-12 10:21 by Hayaay
[考研] 求调剂,一志愿材料科学与工程985,365分, +8 材化李可 2026-04-11 10/500 2026-04-12 08:42 by 852137818
[考研] 316求调剂 +5 想读研究生( ?∵ 2026-04-07 5/250 2026-04-12 00:43 by 蓝云思雨
[考研] 0860004 求调剂 309分 +9 Yin DY 2026-04-08 9/450 2026-04-11 22:55 by dongdian1
[考研] 药学专硕调剂 +8 ? 一路生?花? 2026-04-10 10/500 2026-04-11 21:21 by zhouxiaoyu
[考研] 085410 273分调剂 +4 X1999 2026-04-09 4/200 2026-04-11 13:05 by pies112
[考研] 考研调剂 +26 硕星赴 2026-04-09 27/1350 2026-04-10 22:24 by 猪会飞
[考研] 263能源动力专硕求调剂 +4 加大号饭盒袋 2026-04-10 4/200 2026-04-10 20:52 by gong120082
[论文投稿] mdpi小修rvr时间四五天了 20+3 哈哈high 2026-04-08 5/250 2026-04-10 16:02 by 北京莱茵润色
[考研] 273求调剂 +51 麦小叮当 2026-04-06 58/2900 2026-04-10 15:54 by jiajinhpu
[考研] 298求调剂 +13 钉叮咚冬瓜 2026-04-09 13/650 2026-04-10 15:49 by jiajinhpu
[考研] 083200 初试305分 求调剂 暂不考虑跨专业 +15 Claireyyyy 2026-04-09 15/750 2026-04-09 16:11 by zhuimr
[考研] 285求调剂 +12 AZMK 2026-04-05 18/900 2026-04-08 20:43 by 逆水乘风
[考研] 266调剂 +8 daya sun 2026-04-07 9/450 2026-04-08 20:27 by yutian743
[考研] 0703化学调剂 348分 +14 唉我超真没招了 2026-04-06 15/750 2026-04-08 19:16 by 我减肥1
[考研] 344求调剂 +11 魏子per 2026-04-07 11/550 2026-04-07 23:01 by JourneyLucky
[考研] 316求调剂 +4 15318418673 2026-04-07 4/200 2026-04-07 22:12 by hemengdong
信息提示
请填处理意见