rayleigh_eff.asv

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

s=1;
M=1000;
par=4;
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;
for tent=3:0.5:12
s=0.708;
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(n)=interp1(z,pdcint,abs(z));

pavg=mean(z.*z);
effavg=pavg/(mean(pdcint))