matlab做贝塔分布的例子,关于正态分布和贝塔分布的案例介绍

本文通过案例介绍了正态分布和贝塔分布的概念。

正态分布

正态分布,是一种非常常见的连续概率分布,其也叫做常态分布(normaldistribution),或者根据其前期的研究贡献者之一高斯的名字来称呼,高斯分布(Gaussiandistribution)。正态分布是自然科学与行为科学中的定量现象的一个方便模型。

各种各样的心理学测试结果和物理现象的观测值,比如光子计数等都被发现近似地服从正态分布。甚至生活中很多现象的表征结果也符合正态分布的分布规律。尽管这些现象的根本原因经常是未知的,甚至被采样的样本的原始群体分布并不服从正态分布,但这个变量的采样分布均值仍会近似服从正态分布。

正态分布的概率密度函数呈左右对称的钟形,其具体表达式为:

689aa313bfc7107283756415442344f0.png

因为正态分布是如此的常见而这个式子是如此的奇怪,我们打算重温高斯当年的推导过程,但部分细节不会那么严谨的证明,只是带领大家看看高斯当年的思路是如何的。

首先,高斯事先假定了如下条件,才得到了正态分布的连续密度函数。

即:误差分布导出的极大似然估计=算术平均值

这里我们把全部过程用直白的语言复述一遍。

  • 0
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
以下是使用 MATLAB 实现上述任务的步骤: 1. 安装 MATLAB 工具箱:Statistics and Machine Learning Toolbox,这个工具箱包含了威布尔分布贝塔分布正态分布函数。 2. 创建一个 MATLAB 脚本,导入所需的工具箱: ```matlab % 导入工具箱 addpath('C:\Program Files\MATLAB\R2021a\toolbox\stats') ``` 注意:将路径更改为您的 MATLAB 安装位置。 3. 定义随机场景的数量 n。 ```matlab n = 10; % 随机场景的数量 ``` 4. 定义风速、光伏和负荷的分布参数。这些参数根据您的实际情况进行调整。 ```matlab % 风速分布参数 alpha = 5; % 威布尔分布的形状参数 beta = 15; % 威布尔分布的比例参数 % 光伏分布参数 a = 2; % 贝塔分布的形状参数 b = 5; % 贝塔分布的比例参数 % 负荷分布参数 mu = 100; % 正态分布的均值 sigma = 20; % 正态分布的标准差 ``` 5. 生成随机场景的参数。使用拉丁超立方采样 (LHS) 方法生成参数,以确保场景的均匀分布。 ```matlab % 生成随机场景的参数 rng('default'); % 设置随机数种子 wind_speed = wblinv(lhsdesign(n, 1), alpha, beta); % 风速 solar_power = betainv(lhsdesign(n, 1), a, b); % 光伏 load = norminv(lhsdesign(n, 1), mu, sigma); % 负荷 ``` 6. 输出随机场景的参数: ```matlab % 输出随机场景的参数 disp('随机场景的参数:') disp(['风速:', num2str(wind_speed')]) disp(['光伏:', num2str(solar_power')]) disp(['负荷:', num2str(load')]) ``` 7. 可选:绘制随机场景的概率密度函数 (PDF)。 ```matlab % 绘制随机场景的概率密度函数 x = linspace(0, 30, 100); % 生成 x 值 figure; subplot(1, 3, 1); y1 = wblpdf(x, alpha, beta); % 风速的 PDF plot(x, y1); xlabel('风速'); ylabel('概率密度函数'); title('威布尔分布'); subplot(1, 3, 2); y2 = betapdf(x, a, b); % 光伏的 PDF plot(x, y2); xlabel('光伏'); ylabel('概率密度函数'); title('贝塔分布'); subplot(1, 3, 3); y3 = normpdf(x, mu, sigma); % 负荷的 PDF plot(x, y3); xlabel('负荷'); ylabel('概率密度函数'); title('正态分布'); ``` 8. 运行脚本,生成随机场景的参数和概率密度函数

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值