| 查看: 1779 | 回复: 11 | |||
| 本帖产生 1 个 计算强帖 ，点击这里进行查看 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
求Matlab 高手改进程序 可付费
|
|||
|
用MATLAB自己写了个MLE 估计， 由于数据量比较大 发现如果完全不出错跑完该程序需要60天， 求高手帮忙改进。 抱歉 目前金币紧缺 如果有高手能解决问题 可付费 先谢谢啦  [ Last edited by suyan200326 on 2012-7-13 at 16:20 ] |
回复此楼» 猜你喜欢
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有12人回复
-
物理学I论文润色/翻译怎么收费?
已经有68人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有13人回复
-
桂林理工大学物理学专业招收调剂
已经有18人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助高手将以下C语言翻译为Matlab语言，急急急！
已经有12人回复
- matlab 气泡边缘提取 求助高手！！
已经有5人回复
- 求帮忙啊~ matlab一多项式结果太长，怎么写到txt里？
已经有14人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- 大规模优化数值计算，梯度的计算方式讨论
已经有24人回复
- 采用matlab中nlinfit拟合，出现问题了，不知道是哪错了，该怎么解决，求高手指点一下
已经有9人回复
- 【求助】请各位高手告诉我：这种想法能不能用matlab实现？200金币求助
已经有16人回复
- 【求助】【求助】寻求matlab年轻高手
已经有14人回复
- 【求助】用matlab编程，需要将数据自检分类（高手进来指点下）
已经有5人回复
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 计算强帖: 12
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
4楼2012-07-19 21:24:22
onesupeng
金虫 (职业作家)
- 计算强帖: 13
- 应助: 256 (大学生)
- 贵宾: 1.36
- 金币: 2336.2
- 散金: 9228
- 红花: 92
- 帖子: 4583
- 在线: 1303.8小时
- 虫号: 394701
- 注册: 2007-06-07
- 专业: 流体力学
2楼2012-07-17 08:24:25
3楼2012-07-19 18:40:13
★
fegg7502: 金币+1, 多谢交流 2012-07-23 18:27:51
fegg7502: 金币+1, 多谢交流 2012-07-23 18:27:51
|
是在是太感激你了，最近被这个问题折磨得已经抑郁了。 机器配置很一般 学校的电脑 2.13GH,2.99GB 目前CODE 如下： /*likelihood function*/ function llf=likelihoodfunction(theta, ETS, TS,DS,V0) %global n ETS TS DS; mu=theta(1); sigma=theta(2); n=length(ETS); %V0=10; lf=0; for i=1:n Et=ETS(i); t=TS(i); D=DS(i); %Vt=invoptionprice(Et,D,sigma,t) Vt=fsolve(@(Vt) Vt*normcdf((log(Vt/D)+0.5*sigma^2*t)/(sigma*t^0.5))-D*normcdf((log(Vt/D)+0.5*sigma^2*t)/(sigma*t^0.5)-sigma*t^0.5)-Et, V0); dt=log(Vt/D)+0.5*sigma^2*t; if i>1 f=-0.5*log(2*pi)-0.5*log(sigma^2)-log(Vt)-log(normcdf(dt))-0.5/sigma^2*(log(Vt/Vt_1)-mu)^2; lf=lf+f; end Vt_1=Vt; end llf=-lf; /*MLE*/ Data=xlsread('data1.xls') DS=Data(:,1); ETS=Data(:,2); LENGTH=length(ETS); TS=Data(:,3); %TS=ones(LENGTH,1); theta0=[0,0.1]; V0=10000; options = optimset('Display','iter','MaxIter',(1e+003),'MaxFunEvals',1e+4,'TolFun',10^(-10),'TolX',10^(-10)); [theta,fval]=fminsearch(@(theta) likelihoodfunction(theta,ETS, TS,DS,V0) , theta0,options); n=length(ETS); mu=theta(1); sigma=theta(2); VS=[]; for i=1:n Et=ETS(i); t=TS(i); D=DS(i); Vt=fsolve(@(Vt) Vt*normcdf((log(Vt/D)+0.5*sigma^2*t)/(sigma*t^0.5))-D*normcdf((log(Vt/D)+0.5*sigma^2*t)/(sigma*t^0.5)-sigma*t^0.5)-Et, V0); VS=[VS,Vt]; end |
5楼2012-07-19 23:34:10
20