24小时热门版块排行榜    

查看: 641  |  回复: 3

greetrunners

金虫 (著名写手)

[求助] 求助 用 非线性最小二乘法求解

>> a1=[0.0393
0.0649
0.0991
0.1292
0.1365
0.1762
0.1929
0.2043
0.2162
0.2438
0.3212
];
>> a2=[0.0286
0.0799
0.1665
0.1945
0.2224
0.2858
0.3002
0.3092
0.3123
0.3196
0.3298
];
>> a3=[0.9321
0.8552
0.7344
0.6763
0.6411
0.5380
0.5069
0.4865
0.4715
0.4366
0.3490
];
>> b1=[0.9126
0.8545
0.7341
0.6766
0.6403
0.5708
0.5442
0.5372
0.5176
0.5002
0.4300
];
>> b2=[0.0299
0.0808
0.1762
0.2094
0.2356
0.2796
0.2992
0.3088
0.3146
0.3230
0.3293
];
>> b3=[0.0575
0.0647
0.0897
0.1140
0.1241
0.1496
0.1566
0.1540
0.1678
0.1768
0.2407
];
目标函数  min{(a1.*k1-b1)^2+(a2.*k2-b2)^2+(a3.*k3-b3)^2}
利用 最小二乘法 lsqnonlin函数求 k1、 k2、 k3
回复此楼

» 猜你喜欢

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

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

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
greetrunners: 金币+20, ★★★很有帮助 2012-08-24 15:49:36
dbb627: 金币+1, 感谢应助 2012-08-25 21:02:24
用1stOpt试试:
CODE:
Constant a1=[0.0393,0.0649,0.0991,0.1292,0.1365,0.1762,0.1929,0.2043,0.2162,0.2438,0.3212],
         a2=[0.0286,0.0799,0.1665,0.1945,0.2224,0.2858,0.3002,0.3092,0.3123,0.3196,0.3298],
         a3=[0.9321,0.8552,0.7344,0.6763,0.6411,0.5380,0.5069,0.4865,0.4715,0.4366,0.3490],
         b1=[0.9126,0.8545,0.7341,0.6766,0.6403,0.5708,0.5442,0.5372,0.5176,0.5002,0.4300],
         b2=[0.0299,0.0808,0.1762,0.2094,0.2356,0.2796,0.2992,0.3088,0.3146,0.3230,0.3293],
         b3=[0.0575,0.0647,0.0897,0.1140,0.1241,0.1496,0.1566,0.1540,0.1678,0.1768,0.2407];
MinFunction Sum(a1,b1)((a1*k1-b1)^2)+Sum(a2,b2)((a2*k2-b2)^2)+Sum(a3,b3)((a3*k3-b3)^2);

k1: 2.77154928053109
k2: 1.01003065554773
k3: 0.186457075839164

» 本帖已获得的红花(最新10朵)

2楼2012-08-24 09:40:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

greetrunners

金虫 (著名写手)

送鲜花一朵
引用回帖:
2楼: Originally posted by dingd at 2012-08-24 09:40:46
用1stOpt试试:

Constant a1=,
         a2=,
         a3=,
         b1=,
         b2=,
         b3=;
MinFunction Sum(a1,b1)((a1*k1-b1)^2)+Sum(a2,b2)((a2*k2-b2)^2)+Sum(a3,b3)((a3*k3-b3)^2);

...

请问 1stOpt  从那儿 下载的 。
努力!努力!努力!
3楼2012-08-24 15:48:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

greetrunners

金虫 (著名写手)

引用回帖:
2楼: Originally posted by dingd at 2012-08-24 09:40:46
用1stOpt试试:

Constant a1=,
         a2=,
         a3=,
         b1=,
         b2=,
         b3=;
MinFunction Sum(a1,b1)((a1*k1-b1)^2)+Sum(a2,b2)((a2*k2-b2)^2)+Sum(a3,b3)((a3*k3-b3)^2);

...

你可以用 lsqnonlin  函数 解决吗 ?
努力!努力!努力!
4楼2012-08-24 15:48:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 greetrunners 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +17 Lanmanbaby 2026-08-09 29/1450 2026-08-12 00:24 by zhanghaozhu
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +8 Tide man 2026-08-10 9/450 2026-08-11 20:39 by beefly
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +3 majunge000 2026-08-11 4/200 2026-08-11 20:13 by lch2012
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +3 工厂打螺丝 2026-08-08 3/150 2026-08-11 13:42 by 会议编辑
[基金申请] 听说今天filecode变了 +26 布布和一二 2026-08-06 49/2450 2026-08-11 13:18 by WH3796
[基金申请] filecode +9 documentary 2026-08-10 9/450 2026-08-11 10:08 by wk7465
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 关于Filecode分析方法 +3 majunge000 2026-08-10 3/150 2026-08-10 15:46 by lch2012
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
信息提示
请填处理意见