MATLAB代码实现nakagami-m

以下是MATLAB代码实现nakagami-m分布的概率密度函数和累积分布函数：

function [pdf, cdf] = nakagami_m_distribution(x, m, omega)
% x: 自变量
% m: 形状参数
% omega: 尺度参数

pdf = 2 * m^m / (gamma(m) * omega^m) .* x.^(2*m-1) .* exp(-m*x.^2/omega);
cdf = gammainc(m, m*x.^2/omega, 'upper');

end

其中，`gamma`和`gammainc`分别是Gamma函数和不完全Gamma函数，需要先在MATLAB中导入`symbolic`工具箱。

使用方法如下：

x = 0:0.1:10;
m = 2;
omega = 1;

[pdf, cdf] = nakagami_m_distribution(x, m, omega);

figure;
subplot(2,1,1);
plot(x, pdf);
xlabel('x');
ylabel('pdf');
title(['nakagami-m distribution (m=', num2str(m), ', omega=', num2str(omega), ')']);

subplot(2,1,2);
plot(x, cdf);
xlabel('x');
ylabel('cdf');
title(['nakagami-m distribution (m=', num2str(m), ', omega=', num2str(omega), ')']);

输出结果如下图所示：