PDF与CDF计算并绘图

kongfeng0824

1. clc;
   
2. clear;
   
3. [typ54416, desc54416, fmt54416] = xlsfinfo('D:\data\降水数据\站点数据\Beijing\miyun54416.xls');
   
4. [ndata54416, text54416, alldata54416] = xlsread('D:\data\降水数据\站点数据\Beijing\miyun54416.xls');
   
5. filename54416 = 'D:\data\降水数据\站点数据\Beijing\miyun54416.xls';
   
6. sheet = 1;
   
7. station54416=xlsread(filename54416,'A:A');
   
8. year54416 = xlsread(filename54416,'B:B');
   
9. month54416 = xlsread(filename54416,'C:C');
   
10. day54416 = xlsread(filename54416,'D:D');
    
11. precip54416 = xlsread(filename54416,'E:E');
    
12. y1=precip54416/10;
    
13. ymin=min(y1);
    
14. ymax=max(y1);
    
15. x=linspace(round(ymin),round(ymax),20);  %将最大最小区间分成20个等分点(19等分),然后分别计算各个区间的个数
    
16. yy=hist(y1,x);  %计算各个区间的个数
    
17. yy=yy/length(y1); %计算各个区间的个数密度
    
18. 
    
19. x=yy;
    
20. y1=poisspdf(x,2.5);
    
21. y2=poisspdf(x,5);
    
22. y3=poisspdf(x,10);
    
23. hold on
    
24. plot(x,y1,':r*')
    
25. plot(x,y2,':b*')
    
26. plot(x,y3,':g*')
    
27. hold off

复制代码

完整代码如上，具体的问题如下所示：

最近在处理降水的概率累计曲线和密度曲线：
遇到如下问题：
我在网上看到的cdf和pdf的用法如下所示：

x=-3:0.2:3;
y=normcdf(x,0,1);
plot(x,y)
x=-3:0.2:3;
y=normpdf(x,0,1);
plot(x,y)
我用的是根据“爱的待者”发表的帖子里面的：
y1=precip54416/10; %降水的数据
ymin=min(y1);
ymax=max(y1);
x=linspace(round(ymin),round(ymax),20);  %将最大最小区间分成20个等分点(19等分),然后分别计算各个区间的个数
yy=hist(y1,x);  %计算各个区间的个数
yy=yy/length(y1); %计算各个区间的个数密度
如果我想用normpdf和mormcdf得如何转换数据呢？
我用了一下的命令：
y=normcdf(yy,0,1);
figure(1);
plot(x,y);
hold on;
bar(x,yy) %画出概率密度分布图
figure(2);
y2=normpdf(yy,0,1);
plot(x,y2);
但是画出来的图不是那样的，结果如下图所示，不知道该如何修改？
请大神们指点！

kongfeng0824

自己先顶一下！

七十三棵树

同问，求高手解答！