24小时热门版块排行榜    

查看: 4586  |  回复: 5

jhwhust

金虫 (小有名气)

[求助] 求助一个超定方程组最小二乘解问题 已有3人参与

我是学金属材料的,对matlab不是基本也就了解一点点皮毛,现在需要用这个工具处理几个数据,所以求高手们指点一下:
0.020644a+0.25b+c=2.3987e-06
0.015014a+0.157b+c=3.4030e-06
0.010322a+c=3.80862e-06
0.008692a+0.274b+c=3.54008e-06
0.008256a+0.16b+c=3.4693e-06
三个未知数,a,b,c,五个方程,超定方程组。
求这个超定方程组的最小二乘解,要求a>0。
  

我自己找了本matlab的书看了一下,可以用广义逆法和矩阵除法求解这个问题,但不知道怎么在程序里面定义a>0这个条件。
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

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

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
1stOpt求解:
CODE:
Function
0.020644*a+0.25*b+c=2.3987e-06;
0.015014*a+0.157*b+c=3.4030e-06;
0.010322*a+c=3.80862e-06;
0.008692*a+0.274*b+c=3.54008e-06;
0.008256*a+0.16*b+c=3.4693e-06;
a>=0;

a: 1.98880311304289E-23
b: -2.94201316951683E-6
c: 3.81878661509101E-6
2楼2014-07-19 17:02:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็

【答案】应助回帖

感谢参与,应助指数 +1
一看问题就知道楼主的分析出问题了;
你要求a>0,但是a只有在边界上才取到最小二乘
如果同时要求b>0,b也只有在边界上才最小二乘
所以,楼主这不是计算问题
ส็็็็็็็็็็็็็็็็็็็็
3楼2014-07-19 17:25:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jhwhust

金虫 (小有名气)

引用回帖:
3楼: Originally posted by cooooldog at 2014-07-19 17:25:27
一看问题就知道楼主的分析出问题了;
你要求a>0,但是a只有在边界上才取到最小二乘
如果同时要求b>0,b也只有在边界上才最小二乘
所以,楼主这不是计算问题

呃,a的值从金属材料研究的角度来说只能是正数;
对数值分析我是菜鸟啦,
所以不是很了解您所说的“不是计算问题”是何意?
还有“要求a>0,a只有在边界上才取得最小二乘”是您是怎么看出来的?
能否指点一二
4楼2014-07-19 18:57:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็

引用回帖:
4楼: Originally posted by jhwhust at 2014-07-19 18:57:14
呃,a的值从金属材料研究的角度来说只能是正数;
对数值分析我是菜鸟啦,
所以不是很了解您所说的“不是计算问题”是何意?
还有“要求a>0,a只有在边界上才取得最小二乘”是您是怎么看出来的?
能否指点一二 ...

很可能是,你用不恰当的或错误的分析,得到或列出了一个没有什么意义的方程;
从而对这个问题的计算变得没有意义
ส็็็็็็็็็็็็็็็็็็็็
5楼2014-07-19 22:13:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kabaisun

专家顾问 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
超定方程的最小二乘解是  INVERSE(A'A)*A'*b,解是唯一的,如果结果给出的解第一个元素是负的呢?怎么可以定义它是正的?不知道 @dingd 的解的原理。
严实诚勇
6楼2014-07-20 09:05:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jhwhust 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学调剂0703 +8 啊我我的 2026-03-11 8/400 2026-03-16 17:23 by 我的船我的海
[基金申请] 今年的国基金是打分制吗? 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 311求调剂 +5 26研0 2026-03-15 5/250 2026-03-16 16:21 by a不易
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
[考研] 材料与化工专硕调剂 +3 heming3743 2026-03-16 3/150 2026-03-16 15:05 by peike
[考研] 304求调剂 +6 小熊joy 2026-03-14 6/300 2026-03-16 12:59 by Iveryant
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考博] 欢迎申博同学联系 +3 天道酬勤2026686 2026-03-10 7/350 2026-03-15 19:03 by 天道酬勤2026686
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 265求调剂 +9 小木虫085600 2026-03-09 12/600 2026-03-14 01:11 by JourneyLucky
[考研] 一志愿湖师大化学289求调剂 +6 XMCMM3.14159 2026-03-10 6/300 2026-03-14 00:28 by JourneyLucky
[考研] 318求调剂 +3 李新光 2026-03-10 3/150 2026-03-14 00:21 by JourneyLucky
[考研] 337一志愿华南理工0805材料求调剂 +7 mysdl 2026-03-11 9/450 2026-03-13 22:43 by JourneyLucky
[考研] 290求调剂 +9 ADT 2026-03-11 9/450 2026-03-13 21:55 by JourneyLucky
[考研] 301求调剂 +6 Liyouyumairs 2026-03-11 6/300 2026-03-13 20:11 by JourneyLucky
[考研] 293求调剂 +3 世界首富 2026-03-11 3/150 2026-03-13 16:27 by JourneyLucky
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 085600材料与化工 309分请求调剂 +7 dtdxzxx 2026-03-12 8/400 2026-03-13 14:43 by jxchenghu
[考研] 材料301分求调剂 +5 Liyouyumairs 2026-03-12 5/250 2026-03-13 14:42 by JourneyLucky
[考研] 290求调剂 +3 柯淮然 2026-03-10 8/400 2026-03-11 13:48 by 柯淮然
信息提示
请填处理意见