努力解决问题

23日视图(30天)
雅苒
雅苒 2023年5月18日在2:10
移动: Torsten 2023年5月18日12:20
我有解决这个问题:
实现定制的概率分布(50分)
经常在工程,你画一个定制的随机样本概率distribu -
优化选择函数p (x)。MATLAB提供的功能,比如兰德,randn画样本
统一的和正态分布,分别。但是,什么如果你的概率函数
其他的吗?
(一)这个问题,你将实现一个方法来画样本:
p (x) =
(2
25 x + 2
5如果x5][0,
0否则。(1)
这样做,首先显式计算累积概率密度函数:
P (x) =
Zx
−∞
杜p (u)。(2)
注意:你所做的不需要为这部分使用MATLAB。在你的报告显示结果。
(b)范围内的累积概率密度函数是[0,1],因为需求
R∞
−∞p (x) dx = 1。与这个事实,你可以样品p (x)首先画一个随机样本吗
y从均匀分布域[0,1](这是什么兰特),然后得到x
使用x = P−1 (y)。创建一个名为myRand的功能,它提供了一个随机样本x的
概率密度函数p (x)。
函数x = myRand ()
%输出:
% x:随机样本来自p (x)
(c)画One hundred.,000个样本从你myRand功能和创建一个结果——的柱状图
荷兰国际集团(ing)MATLAB的直方图函数。设置规范化属性的直方图
pdf。指//www.tianjin-qmedu.com/help/matlab/ref/matlab.graphics。
chart.primitive.histogram。html查找如何做到这一点。也,在的柱状图,
情节p (x)在一个不同的颜色。截图是不会被接受的。评论你的im -
实现接近概率密度函数。,例如,图1显示了Y数量作为一个
函数(X)数量。我们观察到结果上几句话。”
我解决了一部分然后写了下面的代码来解决:
% % b部分
函数x = myRand ()
y = rand ();%生成一个随机数在0和1之间
如果y < = 0%处理这个案子当y = 0
x = 0;
elseify > = 1% y = 1时处理这个案子
x = 5;
其他的
x = (5/2) + sqrt (25 * y + 25/4);%转化CDF实验组的方程
结束
% % c部分
%生成100000个样本
样品= 0 (100000);
i = 1∶
样品(i) = myRand ();
结束
%建立直方图
图;
直方图(样本,“归一化”,“pdf”);
持有;
%画出原始概率密度函数p (x)
x = linspace (0, 1000);
p = (-2/25) * x + 2/5;
情节(x, p,“r”,“线宽”2);
%添加标签和标题
包含(“x”);
ylabel (的概率密度);
标题(柱状图的随机抽样和p (x)的);
%显示情节
网格;
传奇(“样本”,“p (x)”);
结束
然而,而不是得到一个直方图,p (x)策划的总体趋势,我有一个柱状图,基本上是相反的方向(插入)yoiu能帮我。谢谢你!
1评论
沃尔特·罗伯森
沃尔特·罗伯森 2023年5月18日2时38分
编辑:沃尔特·罗伯森 曾于2023年5月18日
或许你感兴趣知道 没有一个 的浮点MATLAB可以生成随机数生成技术支持万博1manbetx 完全 0或 完全 1(除非你想算上一代正态分布)。所以当你使用rand()你知道值将被严格> 0 < 1——你 其他的 是唯一的条款,可以调用。
这也暗示你可以vectorize随机数生成。

登录置评。

答案(1)

Torsten
Torsten 2023年5月18日12:03
移动:Torsten 2023年5月18日12:20
你倒提供是错误的。
在myRand,集合x
x = 5 * (1-sqrt (1 y))
而不是
x = (5/2) + sqrt (25 * y + 25/4);
——沃尔特说——你不需要这样的if语句,因为y总是满足0 < y < 1。

类别

找到更多的在描述性统计荷兰国际集团(ing)ydF4y2Ba帮助中心文件交换

标签

s manbetx 845


释放

R2020b

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!