查看单个帖子
旧 2019-12-10, 20:48   #1
poster
高级会员
 
注册日期: 2019-11-21
帖子: 3,006
声望力: 66
poster 正向着好的方向发展
帖子 帮助理解确定的积分

我正在尝试使用MATLAB和C#将书中的函数转换为代码。

我首先尝试使函数在MATLAB中正常工作。

以下是说明:



变量是:

xt and m can be ignored. zMax = Maximum Sensor Range (100) zkt = Sensor Measurement (49) zkt* = What sensor measurement should have been (50) oHit = Std Deviation of my measurement (5) 我已经在MATLAB中编写了第一个公式N(zkt; zkt *,oHit):

hitProbabilty = (1/sqrt( 2*pi * (oHit^2) ))... * exp(-0.5 * (((zkt- zktStar) ^ 2) / (oHit^2)) ); 这给了我期望的高斯曲线。

我下面的定积分有问题,我不明白如何将其转换为实数,因为我从代码中得到了可怕的值,这是这样的:

func = @(x) hitProbabilty * zkt * x; normaliser = quad(func, 0, max) ^ -1; hitProbabilty = normaliser * hitProbabilty; 有人可以帮我这个积分吗?它应该使我的曲线归一化,但是会发疯....(我正在针对zkt 0:1:100进行此操作,其他所有条件都相同,并绘制出应该输出的概率。)



回答:

您应该使用误差函数 ERF (在基本MATLAB中可用)

编辑1:

@Jim Brissom所述,累积分布函数(CDF)通过以下方式与误差函数相关:

normcdf(X) = (1 + erf(X/sqrt(2)) / 2 , where X~N(0,1) 请注意, NORMCDF需要统计工具箱

编辑2:

我认为看到这些评论会有一个小混乱。.上面仅计算归一化因子,因此,如果要计算某个值范围内的最终概率,则应执行以下操作:

zMax = 100; %# Maximum Sensor Range zktStar = 50; %# What sensor measurement should have been oHit = 5; %# Std Deviation of my measurement %# p(0
poster 当前离线   回复时引用此帖