24小时热门版块排行榜    

查看: 642  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +4 Tide man 2026-08-11 4/200 2026-08-12 00:39 by liuqidi
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +7 布布和一二 2026-08-10 9/450 2026-08-11 22:12 by 医学老男孩
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +8 Tide man 2026-08-10 9/450 2026-08-11 20:39 by beefly
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 分享一下我之前已中青C的计划书的filecode +3 布布和一二 2026-08-11 4/200 2026-08-11 15:20 by aasahr
[基金申请] 听说今天filecode变了 +26 布布和一二 2026-08-06 49/2450 2026-08-11 13:18 by WH3796
[基金申请] FileCode能看出啥? +7 要乐观耀哥 2026-08-10 17/850 2026-08-11 10:50 by 要乐观耀哥
[基金申请] filecode +9 documentary 2026-08-10 9/450 2026-08-11 10:08 by wk7465
[基金申请] 关于Filecode分析方法 +3 majunge000 2026-08-10 3/150 2026-08-10 15:46 by lch2012
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
信息提示
请填处理意见