求任意一组数据的概率密度函数

我们在统计数据处理时,经常计算一个样本的概率密度估计,也就是说给出一组统计数据,要求你绘制出它的概率分布曲线,matlab的统计工具箱中有直接的函数  就是:Ksdensity核心平滑密度估计

[f,xi] = ksdensity(x)

计算样本向量x的概率密度估计,返回在xi点的概率密度f,此时我们使用plot(xi,f)就可以绘制出概率密度曲线。该函数,首先统计样本x在各个区间的概率(与hist有些相似),再自动选择xi,计算对应的xi点的概率密度

f = ksdensity(x,xi)

与上面的相似,只是这时xi我们帮Matlab选定了,ksdesity直接计算对应点的概率密度

1.   %by dynamic

2.   %see also http://www.matlabsky.com

3.   %contact me matlabsky@gmail.com

4.   09.2.21

5.   %

6.   %给一个随机样本

7.   x=[randn(30,1); 5+randn(30,1)];

8.   %计算出各点的概率密度

9.   [f,xi]=ksdensity(x);

10.  %绘制图形

11.  subplot(211)

12.  plot(x)

13.  title('样本数据(Sample Data)')

14.  subplot(212)

15.  plot(xi,f)

16.  title('概率密度分布(PDF)')



作者:来自外星球的柠檬甜橙
链接:https://www.jianshu.com/p/757e60657fde
来源:简书
著作权归作者所有。商业转载请联系作者获得授权,非商业转载请注明出处。

  • 4
    点赞
  • 17
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
要计算概率密度函数(PDF),可以使用scipy.stats模块中的norm.pdf方法。第一种方法是直接使用该方法来计算一组数据概率密度函数,代码如下: ```python import scipy.stats as st data = \[一组数据\] pdf_values = st.norm.pdf(data) ``` 第二种方法是自定义一个函数来计算概率密度函数,代码如下: ```python from math import exp, sqrt, pi def pdf(x, mean, var): return exp(-(x - mean) ** 2 / (2 * var ** 2)) / sqrt(2 * pi) * var data = \[一组数据\] mean = 平均值 var = 方差 pdf_values = \[pdf(x, mean, var) for x in data\] ``` 以上是两种常用的计算概率密度函数的方法。你可以根据具体的需选择其中一种方法来计算。 #### 引用[.reference_title] - *1* *3* [Python计算任意数据的分布函数(累积概率函数CDF)与概率密度(PDF),并绘图](https://blog.csdn.net/weixin_35757704/article/details/119541386)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^control_2,239^v3^insert_chatgpt"}} ] [.reference_item] - *2* [Python计算一组数据的PDF(概率密度函数)方法](https://blog.csdn.net/Littonlei/article/details/115242395)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^control_2,239^v3^insert_chatgpt"}} ] [.reference_item] [ .reference_list ]
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值