- 在线时间
- 15 小时
- 最后登录
- 2015-5-7
- 注册时间
- 2014-8-11
- 听众数
- 13
- 收听数
- 0
- 能力
- 0 分
- 体力
- 99 点
- 威望
- 0 点
- 阅读权限
- 20
- 积分
- 64
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 67
- 主题
- 3
- 精华
- 0
- 分享
- 0
- 好友
- 5
升级   62.11% TA的每日心情 | 衰 2015-5-7 10:21 |
---|
签到天数: 19 天 [LV.4]偶尔看看III
- 自我介绍
- I love there
 群组: Matlab讨论组 |
5体力
clc;clear;
A=csvread('c:\wk.txt');
m=6;
n=floor(length(A)/m);
for i=1:n;
y(i)=max(A([(i-1)*m+1:(i*m)]));
end
%hist(y,15);
x=y;
mypdf=@(x,a,b,u)(exp(-(1+b.*(x-u)/a).^(-1/b))./b.*(-(1+b.*(x-u)/a).^(-1/b-1)));
mycdf=@(x,a,b,u)(exp(-(1+b.*(x-u)/a).^(-1/b)));
[paramhat,paramint]=mle(x,'pdf',mypdf,'cdf',mycdf,'start',[.1,.5,.1])
Error using mlecustom>llf_pdfcdf (line 440)
The PDF function returned negative or zero values.
Error in fminsearch (line 191)
fv(:,1) = funfcn(x,varargin{:});
Error in mlecustom (line 183)
[phat,nll,err,output] = ...
Error in mle (line 229)
[phat, pci] = mlecustom(data,varargin{:});
Error in Untitled2 (line 12)
[paramhat,paramint]=mle(x,'pdf',mypdf,'cdf',mycdf,'start',[.1,.5,.1])
|
-
-
wk.txt
1.02 KB, 下载次数: 1, 下载积分: 体力 -2 点
最佳答案
madio 查看完整内容
这个和你的概率密度函数的定义有关,错误提示说这个函数返回了负值或者零值,这是不允许的,你好好推敲一下它的定义
Error using mlecustom>llf_pdfcdf (line 440)
The PDF function returned negative or zero values.
zan
|