接受拒绝采样matlab代码

接受拒绝采样是一种概率分布采样方法，常用于生成服从非常规分布的随机变量。这种方法的基本思想是：首先生成一个服从某种比较容易处理的分布的随机变量，然后按照特定的条件接受或拒绝这个随机变量，最终得到所需的随机变量。Matlab代码实现如下：

function x = acceptReject(N, p, f, f_max)
%N为采样数，p为目标分布，f为仿真分布，f_max为仿真分布最大值
x = zeros(N,1); %初始化结果向量
for i = 1:N
accept = false; %初始化接受标志
while ~accept
x_prop = rand*f_max; %生成随机变量
u = rand; %生成接受概率
if u <= p(x_prop)/(f(x_prop)*f_max) %按照条件计算接受概率
x(i) = x_prop; %接受样本
accept = true; %修改接受标志
end
end
end

其中，N表示采样数，p为目标分布，f为仿真分布，f_max为仿真分布的最大值。首先，初始化结果向量x，然后对每个采样数执行循环。在循环内部，生成一个在仿真分布中的随机变量x_prop和一个介于0和1之间的接受概率u。如果接受概率满足p(x_prop)/(f(x_prop)*f_max)，则接受该样本，并将接受标志设为真。最终，返回结果向量x。

可以使用该函数生成服从各种非常规分布的随机变量，如beta分布、伽马分布、指数分布等。