rayleigh_eff.m

Doherty功放的相关特性的仿真代码

s=1;
M=1000;
for par=4:0.5:10
x=randn(M,1)*s;
y=randn(M,1)*s;
z=abs(x+j*y);
z=z./max(z);
for n=1:1:length(z)
    if z(n)>mean(z)*(10^(par/20));
        z(n)=mean(z)*(10^(par/20));
    end
end
vdd=1;
m=0;
for tent=3:0.5:12
m=m+1;
    vtent(m)=tent;
    s=1/(10^(tent/20));
for n=1:1:length(z)
    vo=z(n);
    if vo<(vdd*s)
        effdh(n)=(pi/4)*vo/(s*vdd);
    elseif (vo==vdd)
        effdh(n)=pi/4;
    else
        effdh(n)=(pi*vo*vo)/(4*vdd*(vo+s*(vo-vdd)));
      end
   % end
pdcint(n)=vo*vo/effdh(n);
end
pdcinterp=interp1(z,pdcint,(abs(z)));

pavg=mean(z.*z);
effavg(m)=pavg/(mean(pdcint));
end
plot(vtent,effavg);
hold on
end