24小时热门版块排行榜    

查看: 809  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生 +3 灯灯灯灯DDDD 2026-09-03 3/150 2026-09-03 21:22 by 科研皮皮猪
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 YLHlRHNwYkce 2026-09-02 8/400 2026-09-03 20:40 by jyNB2vrbdUpR
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 YLHlRHNwYkce 2026-09-02 8/400 2026-09-03 20:31 by jyNB2vrbdUpR
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 uFOAtb1CDngh 2026-09-02 9/450 2026-09-03 20:20 by jyNB2vrbdUpR
[考博] 找导师 +7 yuanjiabao 2026-08-29 8/400 2026-09-03 20:15 by yinnelly
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 ero8OE6tv9cu 2026-09-02 5/250 2026-09-03 20:09 by jyNB2vrbdUpR
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ero8OE6tv9cu 2026-09-02 9/450 2026-09-03 20:07 by jyNB2vrbdUpR
[基金申请] 面上没中,邀请各位路过的虫友分析一下分数 +3 jackleilei 2026-09-03 3/150 2026-09-03 19:13 by god_tian
[硕博家园] 售一区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一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 ero8OE6tv9cu 2026-09-02 10/500 2026-09-03 16:27 by T0rGB46095mJ
[基金申请] 国社科又开始会评了,不知道这次命运如何 +10 雨打竹帘 2026-08-30 14/700 2026-09-03 13:17 by qsd10086
[基金申请] 要骂人了,新模版改版就是要淡化问题凝练这种虚的东西,结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 17/850 2026-09-01 19:51 by 超级无敌华子
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 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
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
信息提示
请填处理意见