function[c, r]=stabil(s)% This function calculates the source and load stability circles% Usage: [center(2),radius(2)]=stabil(S(2, 2))%   where S(2, 2) is the two port S-parameter matrix%     center(1) is the center of the source stability circle%     center(2) is the center of the load stability circle%     radius(1) is the radius of the source stability circle%     radius(2) is the radius of the load stability circlej=sqrt(-1);del=s(1, 1)*s(2, 2)-s(1, 2)*s(2, 1);k=(1-abs(s(1, 1)) ^ 2-abs(s(2, 2)) ^ 2+abs(del) ^ 2)/(2*abs(s(1, 2)*s(2, 1)))den=abs(s(1, 1)) ^ 2-abs(del) ^ 2;c(1)=(conj(s(1, 1))-conj(del)*s(2, 2))/den;r(1)=abs(s(1, 2)*s(2, 1)/den);den=abs(s(2, 2)) ^ 2-abs(del) ^ 2;c(2)=(conj(s(2, 2))-conj(del)*s(1, 1))/den;r(2)=abs(s(1, 2)*s(2, 1)/den);disp('The center of the source and load circle are:')cdisp('The radius of the source and load circle are:')r

function[g]=conjmat(s)% This function calculates the source and load conjugate match% Usage: [gamma(2)]=conjmat(S(2, 2))%   where S(2, 2) is the two port S-parameter matrix%     gamma(1) is the source conjugate match%     gamma(2) is the load conjugate matchj=sqrt(-1);<a name="466"></a><a name="IDX-169"></a>del=s(1, 1)*s(2, 2)-s(1, 2)*s(2, 1);c1=s(1, 1)-del*conj(s(2, 2));b1=1+abs(s(1, 1)) ^ 2-abs(s(2, 2)) ^ 2-abs(del) ^ 2;c2=s(2, 2)-del*conj(s(1, 1));b2=1+abs(s(2, 2)) ^ 2-abs(s(1, 1)) ^ 2-abs(del) ^ 2;if ( b1 < 0)  temp=b1+sqrt(b1 ^ 2-4*abs(c1) ^ 2);else  temp=b1-sqrt(b1 ^ 2-4*abs(c1) ^ 2);endg(1)=conj(c1)*temp/(2*abs(c1) ^ 2);  if ( b2 < 0 )  temp=b2+sqrt(b2 ^ 2-4*abs(c2) ^ 2);else  temp=b2-sqrt(b2 ^ 2-4*abs(c2) ^ 2);endg(2)=conj(c2)*temp/(2*abs(c2) ^ 2);disp('The source...