应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 计算数学 » 【求助】请教高手如何拟合这样的两组数据
51/1返回列表
查看: 998  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

abingchem

木虫 (著名写手)

[交流] 【求助】请教高手如何拟合这样的两组数据

现在已知一组数据，分别是y1,y2,x1,x2,x3,欲求解下面方程中的A1,A2,A3，请教高手给个思路！
y1=A1*x2^2 + A2*x3^2 + (A1 + A2 - A3)*x2*x3
y2=A1*x1^2 + A3*x3^2 + (A1 + A3 - A2)*x1*x3
目前知道的用线性拟合，在matlab回归可以得到两组A1,A2,A3，但两组结果不一样。因为这两个方程用的是一组A1,A2,A3数据，所以请教高手怎么样才可以将两个方程统一起来求解得到一组A1,A2,A3，多谢啦，在线等~
回复此楼

» 猜你喜欢

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

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2011-03-01 20:27:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leedobb

金虫 (正式写手)
abingchem(金币+3): 恩，回头自己编程序试试 2011-03-02 13:26:36
知道最小二乘法的原理吗？知道的话，拟合起来不难啊
构造函数

F= sum(A1*X2...-y1(k))^2 +sum(A1*X1(k）..-y2(k)) ^2
然后由dF/dA1=0,dF/dA2=0,dF=dA3=0得到三个方程，

解一下就得到答案了。

有一点烦，自己编吧
赞一下 回复此楼
2楼2011-03-02 10:45:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abingchem

木虫 (著名写手)
引用回帖:
Originally posted by leedobb at 2011-03-02 10:45:29:
知道最小二乘法的原理吗？知道的话，拟合起来不难啊
构造函数

F= sum(A1*X2...-y1(k))^2 +sum(A1*X1(k）..-y2(k)) ^2
然后由dF/dA1=0,dF/dA2=0,dF=dA3=0得到三个方程，

解一下就得到答案了。

有一点烦 ...

恩，是个思路，之前想到过matlab这样的代码，不知道结果是不是一样的：

XX1=[X2.^2+X2.*X3   X3.^2+X2.*X3   (-1)*X2.*X3];
XX2 = [X1.^2+X1.*X3   (-1)*X1.*X3    X3.^2+X1.*X3];
AA=[XX1;XX2]\[lnf1;lnf2];

这样得出来的结果还可以，但还是不理想，高手继续啊
赞一下 回复此楼
3楼2011-03-02 11:35:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leedobb

金虫 (正式写手)
abingchem(金币+3): 我这个都是实验数据，好像没有权重的差别~ 2011-03-02 13:36:35
引用回帖:
Originally posted by abingchem at 2011-03-02 11:35:30:
恩，是个思路，之前想到过matlab这样的代码，不知道结果是不是一样的：

XX1=[X2.^2+X2.*X3   X3.^2+X2.*X3   (-1)*X2.*X3];
XX2 = [X1.^2+X1.*X3   (-1)*X1.*X3    X3.^2+X1.*X3];
AA=[XX1;XX2]\[lnf1;lnf ...

改进不理想的方式：
你自己判断一下曲线的哪个区域的形状是最希望被保持的，然后在这段加大权重

那么F函数可变成

F = sum ( (yk-fk) ^2 * f(x) )  其中f(x)就是权重函数，权重函数自己随意定啦。

另外有人希望y数值高的地方误大点也没事，那么上式可进一步改成：

F = sum ( （(yk-fk)/yk^a） ^2 * f(x) ) :   a=0~1调节a看看情况了。
赞一下 回复此楼
4楼2011-03-02 12:12:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)
abingchem(金币+3): 试了，不过感觉比较难写代码，结果也不太好 2011-06-02 10:42:01
可以试试1stOpt，拟合优化非常强大好用。
赞一下 回复此楼
5楼2011-03-12 10:59:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 abingchem 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[基金申请] 国自科面上基金字体 +4 iwuli 2026-03-12 5/250 2026-03-15 17:07 by 风云无泪
[考研] 材料专硕326求调剂 +4 墨煜姒莘 2026-03-15 4/200 2026-03-15 11:02 by dyw
[考研] 309求调剂 +4 花与叶@ 2026-03-10 4/200 2026-03-14 21:26 by a不易
[考研] 265求调剂 +4 威化饼07 2026-03-12 4/200 2026-03-14 17:23 by userper
[基金申请] 现在如何回避去年的某一个专家，不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 328求调剂 +3 5201314Lsy！ 2026-03-13 6/300 2026-03-14 15:31 by hyswxzs
[考研] 求调剂 +5 鹤遨予卿 2026-03-09 5/250 2026-03-14 02:44 by JourneyLucky
[考研] 求调剂！ +4 朔朔话 2026-03-09 4/200 2026-03-14 01:38 by JourneyLucky
[考研] 材料工程专硕，一志愿中国矿业大学，总分314，求调剂 +5 无懈可击的巨人 2026-03-10 5/250 2026-03-14 00:37 by JourneyLucky
[考研] 085600调剂 +5 漾漾123sun 2026-03-12 5/250 2026-03-13 22:06 by 星空星月
[考研] 290求调剂 +9 ADT 2026-03-11 9/450 2026-03-13 21:55 by JourneyLucky
[考研] 281求调剂 +9 Koxui 2026-03-12 11/550 2026-03-13 20:50 by Koxui
[考研] 【考研调剂求收留】 +3 Ceciilia 2026-03-11 3/150 2026-03-13 20:18 by JourneyLucky
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 工科调剂 +4 Jiang191123！ 2026-03-11 4/200 2026-03-13 15:15 by Miko19
[考研] 285求调剂 +4 ytter 2026-03-12 4/200 2026-03-13 14:48 by jxchenghu
[考研] 0856化学工程280分求调剂 +4 shenzxsn 2026-03-11 4/200 2026-03-13 11:55 by ymwdoctor
[考研] 大连大学化学专业研究生调剂 +3 琪久. 2026-03-10 8/400 2026-03-11 10:02 by 琪久.
信息提示
关闭
请填处理意见
关闭