diff options
Diffstat (limited to '1430/CH11/EX11.1/exa11_1.sce')
-rw-r--r-- | 1430/CH11/EX11.1/exa11_1.sce | 25 |
1 files changed, 25 insertions, 0 deletions
diff --git a/1430/CH11/EX11.1/exa11_1.sce b/1430/CH11/EX11.1/exa11_1.sce new file mode 100644 index 000000000..a1b9ee28b --- /dev/null +++ b/1430/CH11/EX11.1/exa11_1.sce @@ -0,0 +1,25 @@ +// Example 11.1
+// A Frequency-Selective Network
+// v_in(t)=10*cos(20*t)+10*cos(300*t)
+R=8;
+L=0.2;
+s=%s;
+H_s= R/(s*L+R); // H(s)=V_out/V_in , applying KVL in figure 11.1
+// Selecting input frequency to
+omega1= 20;
+H_omega1= horner(H_s,%i*omega1);
+a_omega1=abs(H_omega1);// amplitude ratio
+theta1_r=atan(imag(H_omega1),real(H_omega1)); // Phase shift in radian
+theta1_d=atan(imag(H_omega1),real(H_omega1))*(180/%pi); // Phase shift in degree
+// Selecting input frequency to
+omega2=300;
+H_omega2= horner(H_s,%i*omega2);
+a_omega2=abs(H_omega2);// amplitude ratio
+theta2_d=atan(imag(H_omega2),real(H_omega2))*(180/%pi);// Phase shift in degree
+theta2_r=atan(imag(H_omega2),real(H_omega2));// Phase shift in radians
+t=0:0.001:5
+v_out=a_omega1*10*cos(omega1*t+theta1_r)+a_omega2*10*cos(omega2*t+theta2_r)
+plot(t,v_out);
+xlabel('t');
+ylabel('v_out(t)')
+title('Steady State Output Voltage Waveform')
|