| 所在主题： | |
| 文件名: 代码调试.zip | |
| 资料下载链接地址: https://bbs.pinggu.org/a-1856792.html | |
本附件包括：
|
|
| 附件大小: | |
|
代码如下：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) 结果出错，详细如下： ??? Conversion to double from sym is not possible. Error in ==> fminsearch at 175 fv(:,1) = funfcn(x,varargin{:}); Error in ==> my_mle1 at 3 [para,fv]=fminsearch(fun,para0,[],2,varargin{:}); 有没有会的大牛帮我解答下，万分感谢。详细的代码和数据我都放在附件里 |
|
熟悉论坛请点击新手指南
|
|
| 下载说明 | |
|
1、论坛支持迅雷和网际快车等p2p多线程软件下载，请在上面选择下载通道单击右健下载即可。 2、论坛会定期自动批量更新下载地址,所以请不要浪费时间盗链论坛资源,盗链地址会很快失效。 3、本站为非盈利性质的学术交流网站,鼓励和保护原创作品，拒绝未经版权人许可的上传行为。本站如接到版权人发出的合格侵权通知，将积极的采取必要措施；同时，本站也将在技术手段和能力范围内，履行版权保护的注意义务。 (如有侵权，欢迎举报) |
|
京ICP备16021002号-2 京B2-20170662号
京公网安备 11010802022788号
论坛法律顾问：王进律师
知识产权保护声明
免责及隐私声明
