搜索
人大经济论坛 附件下载

附件下载

所在主题:
文件名:  1PIN.txt
资料下载链接地址: https://bbs.pinggu.org/a-1025810.html
附件大小:
请问下各位大牛:
小弟,在用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{:});
肿么办啊!
谢谢了


    熟悉论坛请点击新手指南
下载说明
1、论坛支持迅雷和网际快车等p2p多线程软件下载,请在上面选择下载通道单击右健下载即可。
2、论坛会定期自动批量更新下载地址,所以请不要浪费时间盗链论坛资源,盗链地址会很快失效。
3、本站为非盈利性质的学术交流网站,鼓励和保护原创作品,拒绝未经版权人许可的上传行为。本站如接到版权人发出的合格侵权通知,将积极的采取必要措施;同时,本站也将在技术手段和能力范围内,履行版权保护的注意义务。
(如有侵权,欢迎举报)
二维码

扫码加我 拉你入群

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

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

GMT+8, 2026-1-12 19:56