24小时热门版块排行榜    

查看: 1003  |  回复: 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: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的回帖
查看全部 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的回帖

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

dingd

铁杆木虫 (职业作家)


abingchem(金币+3): 试了,不过感觉比较难写代码,结果也不太好 2011-06-02 10:42:01
可以试试1stOpt,拟合优化非常强大好用。
5楼2011-03-12 10:59:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料专硕306英一数二 +10 z1z2z3879 2026-03-16 13/650 2026-03-18 14:20 by 007_lilei
[考研] 070300化学319求调剂 +6 锦鲤0909 2026-03-17 6/300 2026-03-18 13:22 by Iveryant
[考研] 工科材料085601 279求调剂 +6 困于星晨 2026-03-17 6/300 2026-03-18 10:21 by kkcoco25
[基金申请] 被我言中:新模板不强调格式了,假专家开始管格式了 +4 beefly 2026-03-14 4/200 2026-03-17 22:04 by 黄鸟于飞Chao
[考研] 268求调剂 +7 好运连绵不绝 2026-03-12 8/400 2026-03-17 20:28 by xilongliang
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
[考研] 11408 一志愿西电,277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[基金申请] 今年的国基金是打分制吗? 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 机械专硕调剂 +3 笨笨兔子 2026-03-12 3/150 2026-03-15 20:02 by 栗子粥?
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 265求调剂 +4 威化饼07 2026-03-12 4/200 2026-03-14 17:23 by userper
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[考研] 土木第一志愿276求调剂,科研和技能十分丰富,求新兴方向的导师收留 +3 土木小天才 2026-03-12 3/150 2026-03-13 15:01 by JourneyLucky
信息提示
请填处理意见