24小时热门版块排行榜    

查看: 805  |  回复: 3

小木虫zb

木虫 (正式写手)

[求助] 求用1sopt算个非线性方程组,并用原始数据比较

n1=1/0.33824,n2=1/0.44799, k1=0.45967,k2=0.17356,c0=0.5,1,2,4,6,8,10,15,
20,25,30这几个点,分别求这几个点上的c1, q1,c2,q2,然后作图(横坐标为c,纵坐标为q)方程如下:
cp1=(q1/k1)^n1*qp1/(qp1+qp2), cp2=(q2/k2)^n2*qp2/(qp1+qp2),qp1/q1+qp2/q2=1,
q1*n1=q2*n2,qp2=(c0-cp2)*10/0.1,qp1=(c0-cp1)*10/0.1
然后与原始数据比较:
cp1         qp1           cp2         qp2
                 0.12351        0.03243
0.03915        0.09443        0.29164        0.05851
0.12293        0.18251        0.80687        0.08522
0.53648        0.32369        2.05617        0.10752
1.16782        0.43388        3.4304        0.11203
2.01686        0.51311        4.85381        0.10956
2.93049        0.58315        6.22738        0.11417
5.70359        0.68868        9.87906        0.09473
8.64728        0.76995        13.29876        0.10828
11.81856        0.81884        16.88317        0.09841
15.18404        0.84011        20.65151        0.06237
回复此楼

» 猜你喜欢

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

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

dingd

铁杆木虫 (职业作家)

【答案】应助回帖


感谢参与,应助指数 +1
csgt0: 金币+1, 多谢指导 2013-05-09 14:19:11
没大看懂:
1:给出的cp1, qp1, cp2,qp2数据长度不一样
2:一共6个方程却有8个未知数(cp1, qp1, cp2,qp2,c1, q1,c2,q2)?
2楼2013-05-07 15:51:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小木虫zb

木虫 (正式写手)

引用回帖:
2楼: Originally posted by dingd at 2013-05-07 15:51:05
没大看懂:
1:给出的cp1, qp1, cp2,qp2数据长度不一样
2:一共6个方程却有8个未知数(cp1, qp1, cp2,qp2,c1, q1,c2,q2)?

少了两个空格,0.12351        0.03243对应cp2         qp2,比cp1,qp1多一个数据点,只有6个未知数,没有c1和c2
3楼2013-05-07 16:21:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫zb: 金币+10, ★★★★★最佳答案 2013-05-08 10:24:16
csgt0: 金币+1, 谢谢 2013-05-09 14:19:17
CODE:
Constant n1=1/0.33824,n2=1/0.44799, k1=0.45967,k2=0.17356;
LoopConstant c0=[0.5,1,2,4,6,8,10,15,20,25,30];
Function
cp1=(q1/k1)^n1*qp1/(qp1+qp2);
cp2=(q2/k2)^n2*qp2/(qp1+qp2);
qp1/q1+qp2/q2=1;
q1*n1=q2*n2;
qp2=(c0-cp2)*10/0.1;
qp1=(c0-cp1)*10/0.1;

结果:
c0        cp1        q1        qp1        qp2        cp2        q2
0.5        0.496449612805491        0.36901639360812        0.355038719450498        0.0185130624620537        0.499814869375379        0.488752525333678
1        0.995525055505814        0.468379827669076        0.447494449417464        0.0276621351986313        0.999723378648013        0.620356785083684
2        1.99436571697017        0.594637048202531        0.563428302980303        0.0413351638058763        1.99958664836194        0.787581158916241
4        3.99291360286306        0.755267767931132        0.708639713688471        0.0617576341907416        3.99938242365809        1.00033232987776
6        5.99190060341292        0.868904018739573        0.809939658699862        0.0780967481250867        5.99921903251874        1.15084056045924
8        7.99109743883292        0.959900128065062        0.89025611669791        0.0922416671377059        7.99907758332861        1.27136251797585
10        9.99042144518101        1.03709398242097        0.957855481886111        0.104949317828084        9.9989505068217        1.37360375109587
15        14.989062940042        1.19387370690032        1.09370599578367        0.132669514224926        14.9986733048577        1.58125437906163
20        19.9879873064828        1.31954573159967        1.20126935169745        0.156653980574166        19.9984334601942        1.74770367513914
25        24.9870829749811        1.42624040235786        1.29170250186046        0.178191951448142        24.9982180804855        1.88901796395246
30        29.9862954511155        1.51991775407721        1.37045488840877        0.197959693780442        29.9980204030621        2.01309115528236
4楼2013-05-08 09:27:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小木虫zb 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 YLHlRHNwYkce 2026-09-02 7/350 2026-09-03 18:50 by jyNB2vrbdUpR
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 YLHlRHNwYkce 2026-09-02 8/400 2026-09-03 18:50 by jyNB2vrbdUpR
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 YLHlRHNwYkce 2026-09-02 11/550 2026-09-03 18:42 by jyNB2vrbdUpR
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 rFEsUBKXRll0 2026-09-02 6/300 2026-09-03 16:49 by T0rGB46095mJ
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ero8OE6tv9cu 2026-09-02 8/400 2026-09-03 16:38 by T0rGB46095mJ
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ero8OE6tv9cu 2026-09-02 5/250 2026-09-03 16:38 by T0rGB46095mJ
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 ero8OE6tv9cu 2026-09-02 10/500 2026-09-03 16:27 by T0rGB46095mJ
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 ero8OE6tv9cu 2026-09-02 8/400 2026-09-03 14:48 by T0rGB46095mJ
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 YLHlRHNwYkce 2026-09-02 8/400 2026-09-03 12:03 by T0rGB46095mJ
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +8 gdd2018 2026-08-28 13/650 2026-09-03 06:22 by 羽毛枫f
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ero8OE6tv9cu 2026-09-02 4/200 2026-09-03 03:14 by rM1TE0WVDIIY
[基金申请] 科研人应该花精力去思考如何解决问题,而不是去凝练问题 +7 瞬息宇宙 2026-09-01 14/700 2026-09-02 17:33 by ma0526
[基金申请] 要骂人了,新模版改版就是要淡化问题凝练这种虚的东西,结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 10/500 2026-09-02 11:35 by 大不刘6
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 17/850 2026-09-01 19:51 by 超级无敌华子
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
信息提示
请填处理意见