summaryrefslogtreecommitdiff
path: root/25/CH11/EX11.5/11_5.sce
diff options
context:
space:
mode:
Diffstat (limited to '25/CH11/EX11.5/11_5.sce')
-rwxr-xr-x25/CH11/EX11.5/11_5.sce20
1 files changed, 20 insertions, 0 deletions
diff --git a/25/CH11/EX11.5/11_5.sce b/25/CH11/EX11.5/11_5.sce
new file mode 100755
index 000000000..792860837
--- /dev/null
+++ b/25/CH11/EX11.5/11_5.sce
@@ -0,0 +1,20 @@
+//example:-11.5,page no.-617.
+//program to determine the stability of the transistor by calculating k and |delta|.
+s11=0.894*expm(%i*(-60.6)*%pi/180);
+s21=3.122*expm(%i*(123.6)*%pi/180);
+s12=0.02*expm(%i*(62.4)*%pi/180);
+s22=0.781*expm(%i*(-27.6)*%pi/180);
+delta=(s11*s22)-(s12*s21);
+[mag_delta,theta_delta]=polar(delta);
+k=(1+(abs(delta)^2)-(abs(s11)^2)-(abs(s22)^2))/(2*abs(s12*s21));
+Cl=conj(s22-delta*conj(s11))/(abs(s22)^2-abs(delta)^2);
+[mag_Cl,theta_Cl]=polar(Cl);
+Rl=abs(s12*s21)/(abs(s22)^2-abs(delta)^2);
+Cs=conj(s11-delta*conj(s22))/(abs(s11)^2-abs(delta)^2);
+[mag_Cs,theta_Cs]=polar(Cs);
+Rs=abs(s12*s21)/(abs(s11)^2-abs(delta)^2);
+disp([mag_Cl,theta_Cl])
+disp([mag_Cs,theta_Cs])
+disp(Rl)
+disp(Rs)
+disp("NOTE:-theta is in radian") \ No newline at end of file