Nakagami分布，Nakagami信道的matlab仿真

瑞利，莱斯，nakagami三种信道的理论讲解参见瑞利、莱斯与Nakagami-m信道衰落模型
本文只给出matlab中nakagami信道的仿真代码。

N=10000; % 随机数个数

% Nakagami分布的参数
mu=1;
omega=2;

% 生成服从nakagami分布的随机数
pd=makedist('nakagami',mu,omega);
r=random(pd,N,1);

% 绘制生成的随机数的直方图，以概率密度形式
histogram(r,'Normalization','pdf','DisplayStyle','stairs');
hold on;

% nakagami分布的概率密度的理论表达式
f=@(x) (2*mu.^mu.*x.^(2*mu-1).*exp(-mu*x.^2./omega)./gamma(mu)./omega^mu).*(x>0);
fplot(f);
hold off;

运行结果：

可以看到，所生成的随机序列的直方图包络与理论曲线基本一直，说明这个序列是服从Nakagami分布的。
瑞利，莱斯信道等只需要改动下面的语句：
pd=makedist(‘nakagami’,mu,omega);

将名字换成需要的信道分布并修改参数即可。
信道的名字参考matlab帮助文档。

得到服从指定分布的随机序列后，与原信号相乘即可。
rx = tx.*r;

注：此例中的序列为N行1列，也可改为1行N列，与原信号要保持维度一致。