应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 概率论与统计 » 求这个模型参数估计与优化的Matlab代码
71/1返回列表
查看: 1201  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

ranyouhua

新虫 (初入文坛)

[求助] 求这个模型参数估计与优化的Matlab代码 已有1人参与

模型:
Y=a*T1*P1*W*E*P2/(1+P2/b)+c*T2+d

输入数据：
T1= [0.9902 0.9696 0.9851 0.9341 0.9787 0.9502 0.8941 0.9987 0.9987 0.9897 0.9946 0.7887]
T2= [18.261 22.945 22.074 24.285 22.474 23.742 25.379 20.631 19.377 21.735 18.711 11.493]
P1= [0.5961 0.5777 0.5989 0.6213 0.6258 0.6478 0.6416 0.6554 0.6514 0.6553 1.0000 1.0000]
P2= [19265.5 18877 18052.5 15476 20851.5 21657 22012 6481 21738 20884 17597.5 18969]
W= [0.8755 0.8492 0.8791 0.9119 0.9176 0.9511 0.9417 0.9625 0.9552 0.9615 0.9525 0.9180]
E= [0.2927 0.2817 0.3708 0.4006 0.3684 0.4377 0.4311 0.4939 0.4847 0.4762 0.4602 0.4043]

Y的观测数据：
Y= [-6.23 -7.24 2.21 6.1 15.18 16.38 13.63 2.53 10.51 4.08 0.73 3.51]

需要根据观测估计参数：a, b, c, d，参数的范围a[0.05, 0.25]，b[300, 20000]，c[0.01, 0.3]，d[-0.5, 2]，假设参数的先验分布是正态分布。
请问用matlab 怎么实现a, b, c, d的估计？非常感谢大家的帮助。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-07-13 16:58:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

蝶恋花lwp

铁虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
看你用什么方法了是ols还是mle！

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
2楼2015-07-13 19:09:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主
似乎公式有问题
回复此楼
3楼2015-07-13 20:08:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ranyouhua

新虫 (初入文坛)
引用回帖:
3楼: Originally posted by feixiaolin at 2015-07-13 20:08:01
似乎公式有问题

什么问题啊
回复此楼
4楼2015-07-13 22:12:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主
引用回帖:
4楼: Originally posted by ranyouhua at 2015-07-13 22:12:53
什么问题啊

Y与E、W有较强相关性，可得不错的拟合效果；但与T1、T2、P1、P2相关性很差。
我用1stop 软件试了一下，用你的公式则出不来结果。
赞一下(1人) 回复此楼
5楼2015-07-13 22:28:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ranyouhua

新虫 (初入文坛)
内容已删除
回复此楼
6楼2015-07-14 08:51:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主
内容已删除
回复此楼
7楼2015-07-14 10:16:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ranyouhua 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 cu9Nq1xK233Z 2026-08-03 14/700 2026-08-04 14:44 by jzdxtiUmJ6PL
[基金申请] 面上提前没消息，有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
[基金申请] 纯娱乐，不喜欢勿喷 +6 Tide man 2026-08-04 9/450 2026-08-04 14:22 by biomichael
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 Ocx3OKOzGU4q 2026-08-04 5/250 2026-08-04 13:31 by jzdxtiUmJ6PL
[公派出国] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 Ocx3OKOzGU4q 2026-08-04 6/300 2026-08-04 13:28 by jzdxtiUmJ6PL
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 Ocx3OKOzGU4q 2026-08-03 8/400 2026-08-04 13:15 by jzdxtiUmJ6PL
[基金申请] 娱乐 +4 Tide man 2026-08-03 4/200 2026-08-04 11:51 by wgch518
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 Ocx3OKOzGU4q 2026-08-03 6/300 2026-08-04 11:34 by MWRYsWbJw4DE
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 Ocx3OKOzGU4q 2026-08-03 7/350 2026-08-04 11:31 by MWRYsWbJw4DE
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 f1dmPx051StA 2026-08-03 7/350 2026-08-04 09:31 by v1TJC4qJauAB
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 09:17 by v1TJC4qJauAB
[有机交流] 一个有机合成实验室都需要哪些设备？ 50+3 kf2781974 2026-07-31 10/500 2026-08-04 08:05 by 88817753
[找工作] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 07:37 by k7OM8YghWbkC
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 cu9Nq1xK233Z 2026-08-03 6/300 2026-08-04 06:10 by k7OM8YghWbkC
[公派出国] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 04:46 by k7OM8YghWbkC
[基金申请] 什么时候能放榜呀？ +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[高分子] UV压敏胶开发 +3 ichall 2026-07-30 5/250 2026-08-03 14:40 by Sunrisepay
[基金申请] 2026年国自然面上资助率 +16 布布和一二 2026-07-30 22/1100 2026-08-03 09:20 by gy116024
[高分子] HXDI做水性聚氨酯乳液，是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
信息提示
关闭
请填处理意见
关闭