2.txt

matlab仿真程序 协同通信下的不同信道仿真程序

syms  fayi;
M=2;%调制指数，此处为BPSK调制;
for rt=0:2.5:25;   
r0=1/2*10^(rt/10);
r1k=r0;
r2k=r0;
rok=r1k+r2k;
rpk=r1k*r2k;
gpsk=(sin(pi/M))^2;
lk=rok+(rpk*gpsk)/(sin(fayi)^2);
f=(1+r0*gpsk/sin(fayi)^2)^-1*(rok*lk-4*rpk+4*rpk*(lk-rok)*log((lk+sqrt(lk^2-4*rpk))/(2*sqrt(rpk)))/sqrt(lk^2-4*rpk))/(lk^2-4*rpk);
Ps=vpa(int(f,fayi,0,pi*(M-1)/M))
end
------------
M=2;%调制指数，此处为BPSK调制;
for rt=0:2.5:25;   
d=0.01;
fayi=eps:d:pi*(M-1)/M;%eps代替0，无穷小
r0=1/2*10.^(rt/10);
r1k=r0;
r2k=r0;
rok=r1k+r2k;
rpk=r1k*r2k;
gpsk=(sin(pi/M))^2;
lk=rok+(rpk*gpsk)./(sin(fayi).^2);
f=(1+r0*gpsk./sin(fayi).^2).^-1.*(rok*lk-4*rpk+4*rpk*(lk-rok).*log((lk+sqrt(lk.^2-4*rpk))/(2*sqrt(rpk)))./sqrt(lk.^2-4*rpk))./(lk.^2-4*rpk);
ps=d*trapz(f)
end
-------------
M=2;%调制指数
K=1;%中继个数
rt=0:2.5:25;
r0=1/2*10.^(rt./10);
rpk=r0.^2;
rok=r0*2;
gpsk=sin(pi/M)^2;
syms fayi;
f=sin(fayi)^(2*K+2);
f2=sin(fayi)^(2*(K+1)+2);
wkm=(int(f,fayi,0,pi*(M-1)/M))/pi;
wkm2=(int(f2,fayi,0,pi*(M-1)/M))/pi;
Pl=((1+r0*gpsk).*(1+rpk*gpsk./rok)).^-1*wkm %k=1
Pl2=(((1+r0*gpsk).*(1+rpk*gpsk./rok)).*(1+rpk.*gpsk./rok)).^-1*wkm2 %k=2   
pb=(1-sqrt(r0./(1+r0)))/2