testjakes.m

这是一个wcdma下行功率控制仿真的源代码

clear;
clc;
sampletime = 1/(3.84 * 10^6);
totaltime = 2560* sampletime;
averagePower=[0,-3,-6,-9];%multi-path average power dB
averagePower=10.^(averagePower/10);
average = 1/(sum(sqrt(averagePower)).^2);
averagePower = averagePower*average
%averagePower=averagePower/sum(averagePower);% avergePower unitary
rand('state',sum(100*clock));
u = rand(1,5120)>0.5;
[u1,u2] = SerailToParallel(u);
[I_SymbolsTx,Q_SymbolsTx] = QPSKEncoder(u1,u2);
Symbols = complex(I_SymbolsTx,Q_SymbolsTx)
symbos_power = mean(abs(Symbols).^2)

fader1 = RayleighGen(0,60,totaltime,sampletime,averagePower(1));
fader1_power = mean(abs(fader1).^2)
sf1 = Symbols.*fader1;
sf1_power = mean(abs(sf1).^2)
%fakefader1 = sqrt(averagePower(1))*RayleighGen(0,60,totaltime,sampletime,1);
%m = mean(abs(fader1).^2)
fader2 = RayleighGen(0,60,totaltime,sampletime,averagePower(2));
fader2_power = mean(abs(fader2).^2)
sf2 = Symbols.*fader2;
sf2_power = mean(abs(sf2).^2)
%fakefader2 = sqrt(averagePower(2))*RayleighGen(0,60,totaltime,sampletime,1);
%m = mean(abs(fader2).^2)
fader3 = RayleighGen(0,60,totaltime,sampletime,averagePower(3));
fader3_power = mean(abs(fader3).^2)
sf3 = Symbols.*fader3;
sf3_power = mean(abs(sf3).^2)
%fakefader3 = sqrt(averagePower(3))*RayleighGen(0,60,totaltime,sampletime,1);
%m = mean(abs(fader3).^2)
fader4 = RayleighGen(0,60,totaltime,sampletime,averagePower(4));
fader4_power = mean(abs(fader4).^2)
sf4 = Symbols.*fader4;
sf4_power = mean(abs(sf4).^2)
%fakefader4 = sqrt(averagePower(4))*RayleighGen(0,60,totaltime,sampletime,1);
%m = mean(abs(fader4).^2)
sf = sf1+sf2+sf3+sf4;
sf_power = mean(abs(sf).^2)
%fader = fader1+fader2+fader3+fader4;
%fakefader = (sqrt(averagePower(1))+sqrt(averagePower(2))+sqrt(averagePower(3))+sqrt(averagePower(4)))*RayleighGen(0,60,totaltime,sampletime,1);
%m = mean(abs(fader).^2)
%n = mean(abs(fakefader).^2)

envlope = abs(fader1);
subplot(2,2,1);
cdfplot(envlope);
n = mean(envlope.^2) %均方根值


step = 0.1;
t=0:step:5.5;
h=hist(envlope,t);
pdfc=h/(step*sum(h));
subplot(2,2,2);
plot(t,pdfc,':b.');
xlabel('接收信号包络r');
ylabel('p(r)');

envlope=envlope(1:2560);
pws = 10*log10(envlope.^2);
sizes = 1:length(envlope);
t = sizes* 1000/(3.84 * 10^6);
subplot(2,2,3);
plot(t,pws,'--b');
xlabel('时间ms');
ylabel('干扰电平(dB,参考rms)');

%subplot(2,2,4);
%plot(sizes,envlope);