testrake.m

来自「这是一个wcdma下行功率控制仿真的源代码」· M 代码 · 共 47 行

clear;
clc;
x = rayleigh(30000,30000);
averagePower=[0,-3,-6,-9];%multi-path average power dB
averagePower=10.^(averagePower/10);
averagePower=averagePower/sum(averagePower);

ray = sqrt(averagePower(1))*x;
y = ray(1:20);
%fading channel
z = kron(y,ones(1,128));

i = sqrt(-1);
x=bingen(5120);
x(find(x==0))=-1;
Len = length(x);
muxi = x(1:2:Len);
muxq = x(2:2:Len);
Len1 = length(muxi);
% source
output = zeros(1,Len1);
for count=1:Len1
    if muxi(count)==-1 && muxq(count)==-1
          output(count) = j;
    elseif muxi(count)==-1 && muxq(count)==1
          output(count) = 1;
    elseif muxi(count)==1 && muxq(count)==-1
          output(count) = -1;
    elseif muxi(count)==1 && muxq(count)==1
          output(count) = -j;
    end
end
% source through fading channel
obj = output.*z;
% add awgn noise
noise = complex_awgn(Len1,0);
%t = 1:Len1;
%plot(t,10*log10(abs(noise)),'--b');
obj = obj + noise;
% compesation fading channel
m = conj(z);
extz = abs(z).^2;
output
n = (obj.*m)./extz