rayleigh_eff_new.asv

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

n=0;
for par=4:0.5:10
    vmax=(10^(par/20));
%for vo=0:0.1:vmax
 %   if z(n)>mean(z)*(10^(par/20));
  %      z(n)=mean(z)*(10^(par/20));
   % end
%end
vdd=vmax;
m=0;
for tent=3:0.5:12
m=m+1;
    vtent(m)=tent;
    s1=1/(10^(tent/20));
   pavg=quad(rayleigh(vo,1),0,vmax);
    pdc=quad(eff_rayleigh(vo,vdd,s1,1),0,vmax)
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