matlab计算PDF和CDF出现的问题

爱的侍者

clear
clc
data1=load('g:\04\spring\uall.txt');
X=[0 25];
Y=[1 1];
y1=data1(:,1);
ymin=min(y1);
ymax=max(y1);
x=linspace(ymin,ymax,20);  %将最大最小区间分成20个等分点(19等分),然后分别计算各个区间的个数
yy=hist(y1,x);  %计算各个区间的个数
yy=yy/length(y1); %计算各个区间的个数
figure(1)
subplot(211)
bar(x,yy) %画出概率密度分布图
subplot(212)
s=0;
for i=2:length(x)
    s=[s,trapz(x([1:i]),yy([1:i]))];
end
plot(x,s,x,s,'*')
hold on
line(X,Y)


为啥计算的积分大于1 呀？？
搞不懂啦！~ 求解释

流觞曲水

应该是你上图各纵坐标之和大于1了。画累积分布图，最好直接利用累积分布函数，图光滑且不会大于1

kongfeng0824

流觞曲水 发表于 2013-2-22 10:42

                               
登录/注册后可看大图

应该是你上图各纵坐标之和大于1了。画累积分布图，最好直接利用累积分布函数，图光滑且不会大于1

不知道楼主数的累积分布函数，是直接使用内部函数pdf吗？