24小时热门版块排行榜    

查看: 1291  |  回复: 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的回帖

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的回帖
查看全部 5 个回答

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的回帖

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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 i7NFEVbjQMM5 2026-08-16 5/250 2026-08-17 16:38 by CBFiACZS7HAK
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 EWi2p09MOOv8 2026-08-16 3/150 2026-08-17 16:35 by CBFiACZS7HAK
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +34 医学老男孩 2026-08-13 81/4050 2026-08-17 16:20 by yudaoqian88
[基金申请] 感觉是下周放榜了 +4 angus9576 2026-08-17 7/350 2026-08-17 16:11 by Vivilian
[基金申请] 时间戳变了,能看出什么问题? +4 基诺咪客 2026-08-17 4/200 2026-08-17 16:10 by Vivilian
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 6GojgJvkDudM 2026-08-16 4/200 2026-08-17 13:49 by EYK67XJMj64E
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 i7NFEVbjQMM5 2026-08-16 4/200 2026-08-17 12:10 by 4GBAYCdVQoK3
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 6GojgJvkDudM 2026-08-16 4/200 2026-08-17 11:46 by 4GBAYCdVQoK3
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 6GojgJvkDudM 2026-08-16 3/150 2026-08-17 09:49 by xLVPIRuSvCUe
[论文投稿] 岩土工程学报什么时候才能终审完呐 25+3 yeager111 2026-08-10 4/200 2026-08-16 20:42 by tfang
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] 重要来源:本周末出结果 +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
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
信息提示
请填处理意见