aa.rar_doos

% s11=0.72*exp(-i*116*pi/180); % s21=2.60*exp(i*76*pi/180); % s12=0.03*exp(i*57*pi/180); % s22=0.73*exp(-i*54*pi/180); % %input,take impedance Load 13.3648 +53.1389i % s11=0.444*exp(-i*155*pi/180); % s21=1.280*exp(i*64.915*pi/180); % s12=0.14*exp(i*9.224*pi/180); % s22=0.631*exp(-i*78.084*pi/180); %output,take impedance Source 2.9943 +12.1942i s11=0.769*exp(-i*155.8*pi/180); s21=1.157*exp(i*33.3*pi/180); s12=0.11*exp(i*-20.3*pi/180); s22=0.419*exp(i*122.289*pi/180); det=s11*s22-s12*s21 B1=1+abs(s11)^2-abs(s22)^2-abs(det)^2; B2=1+abs(s22)^2-abs(s11)^2-abs(det)^2; C1=s11-s22'*det; C2=s22-s11'*det; gammaS1=(B1+sqrt(B1^2-4*abs(C1)^2))/(2*C1) % MgammaS1=abs(gammaS1) % PmgammaS1=angle(gammaS1)*180/pi impedanceS1=(1+gammaS1)/(1-gammaS1)*50 gammaS2=(B1-sqrt(B1^2-4*abs(C1)^2))/(2*C1) % MgammaS2=abs(gammaS2) % PmgammaS2=angle(gammaS2)*180/pi impedanceS2=(1+gammaS2)/(1-gammaS2)*50 gammaL1=(B2+sqrt(B2^2-4*abs(C2)^2))/(2*C2) % MgammaL1=abs(gammaL1) % PmgammaL1=angle(gammaL1)*180/pi impedanceL1=(1+gammaL1)/(1-gammaL1)*50 gammaL2=(B2-sqrt(B2^2-4*abs(C2)^2))/(2*C2) % Mgammal2=abs(gammaS2) % PmgammaL2=angle(gammaL2)*180/pi impedanceL2=(1+gammaL2)/(1-gammaL2)*50 %---------- f=10^9*[15,35] z=[2.153, 0.727] c=1./(2*pi*(f.*z)*50) %---- length=186*10^-6; Lt=6.621106686656013e-7; Ct=1.1302751243809054e-10; c1=0.15*10^-12; wBrag=1/(length*sqrt(Lt*(Ct+c1)/length))