24小时热门版块排行榜    

查看: 1278  |  回复: 4

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,多谢啦,在线等~
回复此楼

» 猜你喜欢

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

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

查看全部散金贴

已阅   回复此楼   关注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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +32 医学老男孩 2026-08-13 74/3700 2026-08-16 21:05 by 医学老男孩
[论文投稿] 岩土工程学报什么时候才能终审完呐 25+3 yeager111 2026-08-10 4/200 2026-08-16 20:42 by tfang
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +9 Ldrop2023 2026-08-13 9/450 2026-08-16 13:24 by wukun_007
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +11 家与远方 2026-08-10 16/800 2026-08-16 10:28 by ray43
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
信息提示
请填处理意见