【调剂】北京石油化工学院2024年16个专业接受调剂
查看: 889  |  回复: 6
【悬赏金币】回答本帖问题,作者工程师1995将赠送您 50 个金币

工程师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
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wurongjun

专家顾问 (职业作家)

用非线性拟合命令,nlinfit就可以!

发自小木虫Android客户端
善恶到头终有报,人间正道是沧桑.
2楼2019-04-17 10:26:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

工程师1995

金虫 (小有名气)

引用回帖:
2楼: Originally posted by wurongjun at 2019-04-17 10:26:10
用非线性拟合命令,nlinfit就可以!

我只能测到x(3)和gama的数据,x(1)和x(2)的数据没法测,这样也可以用nlinfit吗?
3楼2019-04-17 10:34:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hzlhm

至尊木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +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
QQ:2120156492
4楼2019-04-17 22:48:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

工程师1995

金虫 (小有名气)

引用回帖:
4楼: Originally posted by hzlhm at 2019-04-17 22:48:16
可以用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- ...

方程里的x(1)和x(2)并不是常数,x(1)和x(2)也都是随gama变化的,按照您的方法可以拟合吗?

发自小木虫Android客户端
5楼2019-04-17 23:25:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

把已知的数据都给出来。
6楼2019-04-30 10:52:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
7楼2019-05-10 11:29:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 工程师1995 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 家乡二本高校/沿海传统私企,如何抉择? 10+6 化学巷 2024-04-15 15/750 2024-04-19 20:43 by jiayouqf7
[基金申请] 基金和生小孩 +33 Ausy 2024-04-15 35/1750 2024-04-19 16:21 by feng6531
[教师之家] 美国教授开除中国研究生:我就不该录取你! +7 苏东坡二世 2024-04-15 7/350 2024-04-19 16:20 by feng6531
[基金申请] 院士建议:对勇于提出解决卡脖子问题新解决方案的学者加大资助力度 +7 zju2000 2024-04-16 10/500 2024-04-19 11:10 by dwuab
[考博] 25申博记录贴 +3 我属驴核动力驴 2024-04-18 4/200 2024-04-19 10:53 by 安塔瓦拉多
[论文投稿] CCS Chemistry投稿求助 30+3 wfqtriumph 2024-04-17 4/200 2024-04-19 08:56 by Bletilla
[论文投稿] 求助 100+3 guhangyu 2024-04-15 3/150 2024-04-19 08:48 by bobvan
[论文投稿] 国家重点研发计划论文可以发MDPI里的期刊吗 +4 lin_2018 2024-04-17 4/200 2024-04-19 08:46 by TsolmonZ
[高分子] 聚氨酯丙烯酸酯预聚体合成 +3 15692385990 2024-04-14 6/300 2024-04-18 19:43 by zhang1991
[论文投稿] 投稿求助 5+3 我是洲洲啊 2024-04-17 5/250 2024-04-18 17:13 by topedit
[论文投稿] MDPI初审7天,是不是被拒 20+6 applegirl5762 2024-04-15 8/400 2024-04-17 08:49 by bobvan
[考研] 332求调剂 +3 木叶下1999 2024-04-16 5/250 2024-04-17 00:56 by angeliar
[基金申请] 迟国泰通过向学生发放劳务费再回收的方式套取科学基金重点项目 +6 babu2015 2024-04-13 7/350 2024-04-16 20:32 by sundiv
[考博] 25年申博求助,射频方向 +4 SherlockAH 2024-04-15 10/500 2024-04-16 18:06 by ljysudaee
[考研] 329求调剂 +6 Kaylawander 2024-04-13 7/350 2024-04-16 12:00 by 风来花开1
[考博] 博士申请 +5 啦啦啦3737 2024-04-13 9/450 2024-04-16 07:07 by 1147wzzzz
[考研] 347求调剂 +4 安谷先生 2024-04-15 4/200 2024-04-15 18:27 by huao连
[材料综合] 请问哪里可以测试低压吸氢的PCT曲线 +3 yunshengcd 2024-04-14 5/250 2024-04-15 13:59 by zqdsb
[考研] 278求调剂 +4 月亮就蒜 2024-04-13 5/250 2024-04-14 23:03 by 永字号
[考研] 300求调剂 +5 @Xqy668800 2024-04-14 6/300 2024-04-14 23:03 by lincunhui
信息提示
请填处理意见