function [ pd ] =CAS_nak_m2(SNR,N,m,NN)

%N为衰落程度系数，mnakagami-m的衰落系数，SNR为信噪比，

pf=0.01;

SNR_integer=10^(SNR/10);%信噪比十进制形式

SNR_integer=1/SNR_integer*m^N;

u=1;

lamda=2*gammaincinv(1-pf,u);%判决门限根据虚警概率确定，见文中公式(2)

pp=(factorial(m-1))^N;

pp=1/pp;

s=zeros(1,length(NN)+1);

for n=0:NN; %为了计算方便，我们可以假设的

w=zeros(1,length(n+u-1)+1);

for k=0:n+u-1

w(k+1)=exp(-lamda/2)/(factorial(k))*(lamda/2)^k;

end

alpha=1-n;

ss=['[[',num2str(alpha),'],[]]'];

rr=['[[',num2str(m),',',num2str(m),'],[]]'];

s(n+1)=sum(w)*1/(factorial(n))*arrayfun(@(SNR_integer) double(feval(symengine,'meijerG',ss, rr,num2str(SNR_integer))), SNR_integer);

end

pd=pp*sum(s);

end

主程序;

clc

clear all

close all

m=2;

N=2;

NN=30;

SNR=[5:0.5:20];

pm=zeros(1,length(SNR));

for i=1:length(SNR)

pm(i)=1-CAS_nak_m2(SNR(i),N,m,NN);

end

figure;

semilogy(SNR,pm);

2016-10-5 16:38 上传

m=2,N=2,U=1,λ程序中给出计算方法

2016-10-5 16:43 上传

程序应该得到的结果