diff options
| author | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
|---|---|---|
| committer | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
| commit | 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch) | |
| tree | dbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /1430/CH9/EX9.11/exa9_11.sce | |
| parent | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff) | |
| download | Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.gz Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.bz2 Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.zip | |
initial commit / add all books
Diffstat (limited to '1430/CH9/EX9.11/exa9_11.sce')
| -rw-r--r-- | 1430/CH9/EX9.11/exa9_11.sce | 32 |
1 files changed, 32 insertions, 0 deletions
diff --git a/1430/CH9/EX9.11/exa9_11.sce b/1430/CH9/EX9.11/exa9_11.sce new file mode 100644 index 000000000..c3c60fcea --- /dev/null +++ b/1430/CH9/EX9.11/exa9_11.sce @@ -0,0 +1,32 @@ +//Example 9.11
+// Underdamped Zero-Input Response
+// form figure 9.25
+L=0.1;
+R=5;
+C=1/640;
+alpha=R/(2*L);
+omega_0=sqrt(1/(L*C));
+//Characteristic Values
+p1=-alpha+sqrt(alpha^2-omega_0^2);
+omega_d=sqrt(omega_0^2-alpha^2);
+p2=p1'; // Complex conjugate
+V_s1=30; // t<0
+V_s2=0;//t>0
+// using initial conditions we find
+i_L_aft=0;// i(0^+)=0
+i_L_aft_d=-30/L; // i'(0^+)=0
+I_ss= 0; // when capacitor becomes fully charge before t<0
+//Using complex matrix equation
+P=[1,1;p1,p2];
+I=[i_L_aft-I_ss;i_L_aft_d]
+A=P\I
+A_1=A(1);
+A_1_m=abs(A_1);
+phase_A_1=atan(imag(A_1),real(A_1))*(180/%pi);
+t=0:0.001:0.5
+i_L=2*A_1_m*exp(-alpha*t).*cos(omega_d*t+phase_A_1);
+plot(t,i_L)
+xlabel('t')
+ylabel('i_L(t)')
+title('Current Waveform')
+
|