24小时热门版块排行榜    

查看: 961  |  回复: 4

无双啊呜

新虫 (小有名气)

[求助] 求大神帮忙编写一个非线性拟合matlab程序

方程如下:y1=λ1*x1+a1*x1/(1+b1*x1+b2*x2)
                 y2=λ2*x2+a2*x2/(1+b1*x1+b2*x2)
程序最好能把图做出来,分别是x1~y1,x2~y2作图,跪谢!
x1                                   x2
0.10036985        0.097194278
0.255816948        0.24947162
0.469506086        0.462258883
0.751893727        0.743209506
1.008132022        0.994501615
1.991118914        1.970911398
2.934007491        2.914803876
4.2667603        4.175588371
5.348852996        5.231472081
y1                              y2
0.230062437        0.29335438
0.60943266        0.773713297
1.157600824        1.450597019
1.814596592        2.270642722
2.336650427        2.941861196
4.653100545        5.87118063
6.870687018        8.600178767
9.623398867        12.05369909
11.82017961        14.86530222
回复此楼

» 猜你喜欢

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

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

dingd

铁杆木虫 (职业作家)

1stOpt求解:
CODE:
Variable x1,x2,y1,y2;
SharedModel;
Function y1=r1*x1+a1*x1/(1+b1*x1+b2*x2);
         y2=r2*x2+a2*x2/(1+b1*x1+b2*x2);
Data;
x1        x2        y1        y2
0.10036985        0.097194278        0.230062437        0.29335438
0.255816948        0.24947162        0.60943266        0.773713297
0.469506086        0.462258883        1.157600824        1.450597019
0.751893727        0.743209506        1.814596592        2.270642722
1.008132022        0.994501615        2.336650427        2.941861196
1.991118914        1.970911398        4.653100545        5.87118063
2.934007491        2.914803876        6.870687018        8.600178767
4.2667603        4.175588371        9.623398867        12.05369909
5.348852996        5.231472081        11.82017961        14.86530222

均方差(RMSE):0.02807198483188
残差平方和(SSE):0.0141846539832234
相关系数(R): 0.999981218538354
相关系数之平方(R^2): 0.999962437429452
确定系数(DC): 0.999959595304795
F统计(F-Statistic): 12694.0543622117

参数                  最佳估算
--------------------        -------------
r1        3.48786912540209
a1        -1.08888004680341
b1        -0.524943279146923
b2        0.508491943305727
r2        4.06356282482197
a2        -1.04138783809056

====== 结果输出 ======

文件: 数据文件-1
No        目标 y1        计算 y1        目标 y2        计算 y2
1        0.230062437        0.240428064061444        0.29335438        0.29340646107397
2        0.60943266        0.611615479762708        0.773713297        0.752000811628464
3        1.157600824        1.12043989347115        1.450597019        1.39147159046462
4        1.814596592        1.78980755367499        2.270642722        2.23289591538678
5        2.336650427        2.3920616276863        2.941861196        2.98061668762716
6        4.653100545        4.67918017359755        5.87118063        5.86414524309084
7        6.870687018        6.84182596315963        8.600178767        8.62203816646627
8        9.623398867        9.62296019907896        12.05369909        12.0456680846732
9        11.82017961        11.8226564031366        14.86530222        14.8664374001913
2楼2013-12-02 19:46:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

水木青枫

银虫 (小有名气)

【答案】应助回帖

就是用光滑的曲线作图是吧?爪机无力,先mark下,下次补上

[ 发自小木虫客户端 ]
生活,就是心怀最大的善意在荆棘中穿行。即使被刺伤,亦不改初衷。
3楼2013-12-03 00:50:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

SavannahDawn

银虫 (初入文坛)

借楼主的贴问一下,如果不用1stOpt,单纯用Matlab的拟合或者优化功能(像solver)那样的,可以解出来么?
It's_never_too_late_-_in_fiction_or_in_life_-_to_revise.
4楼2013-12-03 09:41:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

Matlab本身也是一个编程软件,理论上可实现任何数值计算问题。依据个人能力了。
5楼2013-12-03 13:00:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 无双啊呜 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 oeyPlfmMUpOC 2026-09-26 3/150 2026-09-27 09:13 by 4IW0sJvtEqX8
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 GxESXJmptOtQ 2026-09-25 4/200 2026-09-27 06:44 by XgCC8uTcwILl
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 GxESXJmptOtQ 2026-09-25 4/200 2026-09-27 06:42 by XgCC8uTcwILl
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 GxESXJmptOtQ 2026-09-25 4/200 2026-09-27 06:39 by XgCC8uTcwILl
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 06:08 by XgCC8uTcwILl
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 06:05 by XgCC8uTcwILl
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 05:56 by XgCC8uTcwILl
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 05:56 by XgCC8uTcwILl
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 vQZoDrrm7mUF 2026-09-25 5/250 2026-09-27 05:53 by XgCC8uTcwILl
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 05:50 by XgCC8uTcwILl
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 vQZoDrrm7mUF 2026-09-25 5/250 2026-09-27 05:44 by XgCC8uTcwILl
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 vQZoDrrm7mUF 2026-09-25 3/150 2026-09-27 05:44 by XgCC8uTcwILl
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 vQZoDrrm7mUF 2026-09-25 4/200 2026-09-27 05:44 by XgCC8uTcwILl
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 vQZoDrrm7mUF 2026-09-25 6/300 2026-09-27 05:42 by XgCC8uTcwILl
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 PGTSwSC3F6nQ 2026-09-24 5/250 2026-09-25 11:19 by yhq9807
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 PGTSwSC3F6nQ 2026-09-24 7/350 2026-09-25 11:18 by yhq9807
[硕博家园] 上海工程技术大学 激光智能制造课题组 2027级博士研究生招生 +4 oxidpl 2026-09-20 6/300 2026-09-24 20:40 by oxidpl
[教师之家] 课题组招2027级博士 上海工程技术大学 激光智能制造方向 +4 oxidpl 2026-09-20 6/300 2026-09-24 20:37 by oxidpl
[基金申请] 师弟论文见刊大半年才想起来申请专利,还能抢救一下吗? +4 13108017953 2026-09-21 4/200 2026-09-24 10:06 by kudofaye
[有机交流] 有机合成以后会不会被AI改变?做科研的虫友怎么看 +5 zapen 2026-09-20 8/400 2026-09-23 17:21 by tianxiaoxian
信息提示
请填处理意见