erfinv

erfinv(0.25)

ans = 
0.2253

对于 [-1,1] 之外的输入，erfinv 返回 NaN。对于 -1 和 1，erfinv 分别返回 -Inf 和 Inf。

erfinv([-2 -1 1 2])

ans = 1×4

   NaN  -Inf   Inf   NaN

求矩阵元素的逆误差函数。

M = [0 -0.5; 0.9 -0.2];
erfinv(M)

ans = 2×2

         0   -0.4769
    1.1631   -0.1791

绘制 -1 < x < 1 时的逆误差函数图。

x = -1:0.01:1;
y = erfinv(x);
plot(x,y)
grid on
xlabel('x')
ylabel('erfinv(x)')
title('Inverse Error Function for -1 < x < 1')

使用均匀分布的随机数生成高斯分布的随机数。要将均匀分布的随机数 $$x$$ 变换为高斯分布的随机数 $$y$$，请使用下面的变换

请注意，由于 x 采用 -1 + 2*rand(1,10000) 的形式，您可以通过使用 erfcinv 代替 erfinv 提高准确性。有关详细信息，请参阅 提示。

在区间 [-1,1] 内生成 10,000 个均匀分布的随机数。将它们转换为高斯分布的随机数。使用直方图说明这些数字遵照高斯分布形式。

rng('default')
x = -1 + 2*rand(1,10000);
y = sqrt(2)*erfinv(x);
h = histogram(y);