%efficiency of a pa when we consider the loss of the PA
%2007-02-16
%wujianfeng
%Sweden
clear all
close all
vd=0:0.3:42;
id=0:0.1:14;
vdd=48;
k=pi/4;
rsr=0.1;
rsh=10;
 for m=1:1:length(id)
for n=1:1:length(vd)
eff(m,n)=((vd(m)-id(n)*rsr)*id(n))/(vdd*k*(id(n)+(vd(m)/rsh)));
if eff(m,n)<=0
    eff(m,n)=0;
end
end
    eff1(m)=max(eff(m,:));
end
figure(1)
mesh(vd,id,eff);
xlabel('Vd');
ylabel('Id');
zlabel('eff');
figure(2)
c=contour(vd,id,eff,15);
clabel(c);
xlabel('Vd');
ylabel('Id');
figure(3)
c1=contour3(vd,id,eff,15);
clabel(c1);
xlabel('Vd');
ylabel('Id');
zlabel('eff');