mutistage_impedance0126.m

来自「Doherty功放的相关特性的仿真代码」· M 代码 · 共 64 行

%This is the efficiency about a multi-stage Doherty PA;
%wujianfeng
%2007-01-26
%Sweden
clear all
close all
gama1=2;
gama2=4;
vdd=1;
n=0;
rl=15;
for vo=0:0.01:vdd
    n=n+1;
    vout(n)=vo;
    if ((vo>=0)&(vo<=(vdd/gama2)))
        eff(n)=(pi/4)*(gama2*vo)/(vdd);
        rc(n)=gama2*gama2*rl;
        rp1(n)=2000;
        rp2(n)=2000;
        ic(n)=(vo/(gama2*rl));
        ip1(n)=0;
        ip2(n)=0;
    else if ((vo>(vdd/gama2))&(vo<=(vdd/gama1)))
            a1=vo/(vdd/gama1);
      rc(n)=gama1*gama2*rl/(a1);
      rp1(n)=(a1*gama1*gama1*gama2)*rl/(a1*gama2-gama1);
      rp2(n)=2000;
      ic(n)=(vo/(gama2*rl));
        ip1(n)=(vo/(gama1*rl))-(vdd/(gama1*gama2*rl));
        ip2(n)=0;
      eff(n)=(pi/4)*(gama1*gama2*(vo/vdd)*(vo/vdd))/((gama1+gama2)*(vo/vdd)-1);
        else if ((vo>(vdd/gama1))&(vo<=(vdd/1)))
                a2=vo/vdd;
                rc(n)=gama1*gama2*rl;
                rp1(n)=gama1*gama2*rl/(a2*gama2-1);
                rp2(n)=(a2*gama2*rl)/(a2*gama1-1);
                ic(n)=(vdd/(gama1*gama2*rl));
        ip1(n)=(vo/(gama1*rl))-(vdd/(gama1*gama2*rl));
        ip2(n)=(vo/rl)-(vdd/(gama1*rl));
      eff(n)=(pi/4)*(gama1)*(vo/vdd)*(vo/vdd)/((gama1+1)*(vo/vdd)-1);
            end
        end
    end
end
figure(1)
plot(vout,rc,'r-*');
hold on
plot(vout,rp1,'g-x');
hold on
plot(vout,rp2,'b-+');
legend('Carrier PA load impedance','Peak PA1 load impedance','Peak PA2 load impedance');
figure(2)
plot(vout,eff,'r');
figure(3)
plot(vout,ic,'r-*');
hold on
plot(vout,ip1,'g-x');
hold on
plot(vout,ip2,'b-+');
legend('Carrier PA current','Peak PA1 current','Peak PA2 current');