| 查看: 1680 | 回复: 6 | |||
工程师1995金虫 (小有名气)
|
[求助]
求助MATLAB最小二乘 已有2人参与
|
|
求助MATLAB最小二乘法 (alfa*(x(1)-2*x(2))+1)*x(3)-lamda*gama*(1-x(2))=0; (alfa(*x(1)-2*x(2))+1)*x(1)-2*lamda*gama*x(3)=0; (1-alfa*x(2))*x(2)-alfa*(x(3))^2=0; 上面是一个方程组，从实验中测得了x(3)和gama之间的关系，有50组数据点。 我想要通过最小二乘的方法，回归出alfa，lamda的数值。 问题的难点在于这个方程组无法求解出x(3)关于gama的解析解。只能把三个方程放在一起求解数值解。 所以请大神帮帮我。已经折磨我两个星期了。@laosam280 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有9人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
国自然面上复盘~欢迎讨论
已经有15人回复
-
中青基了要发朋友圈吗？
已经有3人回复
-
国自然评审意见
已经有19人回复
wurongjun
专家顾问 (职业作家)
-
专家经验: +831 - 数学EPI: 9
- 应助: 791 (博后)
- 贵宾: 0.308
- 金币: 24609
- 散金: 310
- 红花: 75
- 帖子: 3004
- 在线: 881.4小时
- 虫号: 1368482
- 注册: 2011-08-14
- 性别: GG
- 专业: 计算数学与科学工程计算
- 管辖: 数学
2楼2019-04-17 10:26:10
工程师1995
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 682.9
- 散金: 3
- 帖子: 154
- 在线: 46小时
- 虫号: 2528602
- 注册: 2013-07-01
- 性别: GG
- 专业: 流体力学
3楼2019-04-17 10:34:00
hzlhm
至尊木虫 (著名写手)
- 应助: 387 (硕士)
- 金币: 17856.5
- 红花: 53
- 帖子: 2880
- 在线: 607.4小时
- 虫号: 1517335
- 注册: 2011-11-30
- 性别: GG
- 专业: 常微分方程与动力系统
【答案】应助回帖
感谢参与，应助指数 +1
|
可以用nlinfit函数来拟合alfa，lamda，x1，x2。其自定义拟合函数可以这样来考虑 function x3=func(x) %求助MATLAB最小二乘法 xalfa=x(1);lamda=x(2);x1=x(3);x2=x(4); %(alfa*(x1-2*x2)+1)*x3-lamda*gama*(1-x2)=0; %(alfa*(x1-2*x2)+1)*x1-2*lamda*gama*x3=0; %(1-alfa*x2)*x2-alfa*(x3)^2=0; x31=lamda*gama*(1-x2)/(alfa*(x1-2*x2)+1); x32=(alfa*(x1-2*x2)+1)*x1/(2*lamda*gama); x33=sqrt((1-alfa*x2)*x2/alfa); x3=(x31+x32+x33)/3; end |
4楼2019-04-17 22:48:16
工程师1995
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 682.9
- 散金: 3
- 帖子: 154
- 在线: 46小时
- 虫号: 2528602
- 注册: 2013-07-01
- 性别: GG
- 专业: 流体力学
5楼2019-04-17 23:25:07
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
6楼2019-04-30 10:52:30
7楼2019-05-10 11:29:14
