全部版块 我的主页
论坛 提问 悬赏 求职 新闻 读书 功能一区 经管百科 爱问频道
7656 14
2012-01-04
请问下各位大牛:
     小弟,在用matlab编程计算知情交易者概率时遇到一个问题,麻烦各位解答一下
程序代码如下
function [para,fv]=my_mle1(fun,para0,varargin)
para0=para0(:);
[para,fv]=fminsearch(fun,para0,[],2,varargin{:});
fv=-fv;
end

function f=pin(para,num,B,S)
alpha=para(1);
mu=para(2);
eb=para(3);
es=para(4);
delta=para(5);
M=(min(B,S)+max(B,S))/2;
xs=es/(es+mu);
xb=eb/(eb+mu);
yy=(-eb-es + M*(log(xb) +log(xs)) + ...
    B*log(mu+eb)+S*log(mu+es) )+ ...
    log( alpha*(1-delta)*exp(-mu)*(xs.^(S-M)) .*(xb.^sym(-M)) + ...
    alpha* delta*exp(-mu) *(xb.^(B-M)).*(xs.^(-M)) + ...
    (1-alpha)*(xs.^(S-M)).*(xb.^(B-M)));
if num==1
    f=log(yy);
else
    f=-sum(log(yy));
end     
end

输入指令:B=x1PIN(:,3)
S=x1PIN(:,4)
[para,fv]=my_mle1('pin',[0.5;400;1000;1000;0.5],B,S)后
报错:
??? Error using ==> sym.sym>notimplemented at 2514
Function 'lt' is not implemented for MuPAD symbolic objects.
Error in ==> sym.sym>sym.lt at 803
            notimplemented('lt');
Error in ==> fminsearch at 323
    if fxr < fv(:,1)
Error in ==> my_mle1 at 3
[para,fv]=fminsearch(fun,para0,[],2,varargin{:});
肿么办啊!
谢谢了
附件列表

1PIN.txt

大小:6.14 KB

 马上下载

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2012-1-4 18:02:43
数据在附件里面
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2012-1-4 21:41:54
自己顶一下!!
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2012-1-5 19:38:20
用fminsearch函数的过程中出现如下问题:
Exiting: Maximum number of function evaluations has been exceeded
         - increase MaxFunEvals option.
         Current function value: -Inf
肿么办,怎么没高手肯出来帮忙呢!
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2012-1-7 12:31:01
怎么都是路人甲~路人乙~解决下问题吧
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2012-1-20 21:25:50
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群