1 files changed, 66 insertions, 0 deletions
diff --git a/887/CH6/EX6.2/6_2.sce b/887/CH6/EX6.2/6_2.sce
new file mode 100755
index 000000000..ea53cc59a
--- /dev/null
+++ b/887/CH6/EX6.2/6_2.sce
@@ -0,0 +1,66 @@
+clc
+//ex6.2
+//given V_in(t)=3+2*cos(2000*%pi*t)+cos(4000*%pi*t-A), A=70*%pi/180
+//the three parts of V_in(t) are V_in_1=3, V_in_2=2*cos(2000*%pi*t),V_in_3=cos(4000*%pi*t-A)
+
+//first component V_1
+V_in_1=3;
+f_1=0;      //as omega is zero
+//equation of straight line of H_magnitude vs f is x+1000*y-4000=0
+H_1_max=(4000-f_1)/1000;      //magnitude of H(traansfer function)
+//equation of straight line of H_phase angle vs f is 6000*y=%pi*x (phase angle in radians)
+H_1_phi=%pi*f_1/6000;      //phase angle of H
+H_1=H_1_max*complex(cos(H_1_phi),sin(H_1_phi));
+V_out_1=H_1*V_in_1;
+V_out_1_R=real(V_out_1);      //real part
+V_out_1_I=imag(V_out_1);      //imaginary part
+V_out_1_max=sqrt((V_out_1_R^2)+(V_out_1_I^2));      //peak value
+V_out_1_phi=atan(V_out_1_I/V_out_1_R);      //phase angle
+
+//second component V_in_2
+V_in_2=2*complex(cos(0),sin(0));      //V_in_2 phasor
+w=2000*%pi;      //omega
+f_2=w/(2*%pi);      //frequency
+//equation of straight line of H_magnitude vs f is x+1000*y-4000=0
+H_2_max=(4000-f_2)/1000;      //magnitude of H(traansfer function)
+//equation of straight line of H_phase angle vs f is 6000*y=%pi*x (phase angle in radians)
+H_2_phi=%pi*f_2/6000;      //phase angle of H
+H_2=H_2_max*complex(cos(H_2_phi),sin(H_2_phi));
+V_out_2=H_2*V_in_2;
+V_out_2_R=real(V_out_2);      //real part
+V_out_2_I=imag(V_out_2);      //imaginary part
+V_out_2_max=sqrt((V_out_2_R^2)+(V_out_2_I^2));      //peak value
+V_out_2_phi=atan(V_out_2_I/V_out_2_R);      //phase angle
+
+//third component
+A=-70*%pi/180;      //-70 degrees in radians
+V_in_3=complex(cos(A),sin(A));      //V_in_3 phasor
+w=4000*%pi;      //omega
+f_3=w/(2*%pi);      //frequency
+//equation of straight line of H_magnitude vs f is x+1000*y-4000=0
+H_3_max=(4000-f_3)/1000;      //magnitude of H(traansfer function)
+//equation of straight line of H_phase angle vs f is 6000*y=%pi*x (phase angle in radians)
+H_3_phi=%pi*f_3/6000;      //phase angle of H
+H_3=H_3_max*complex(cos(H_3_phi),sin(H_3_phi));
+V_out_3=H_3*V_in_3;
+V_out_3_R=real(V_out_3);      //real part
+V_out_3_I=imag(V_out_3);      //imaginary part
+V_out_3_max=sqrt((V_out_3_R^2)+(V_out_3_I^2));      //peak value
+V_out_3_phi=atan(V_out_3_I/V_out_3_R);      //phase angle
+
+disp('Output voltage is Vout1+Vout2+Vout3 where')
+disp('')
+disp('FOR Vout1:')
+disp(V_out_1_max,'peak value in volts')
+disp(V_out_1_phi*180/%pi,'phase angle in degrees')
+disp(f_1,'with frequency in hertz')
+disp('')
+disp('FOR Vout2:')
+disp(V_out_2_max,'peak value in volts')
+disp(V_out_2_phi*180/%pi,'phase angle in degrees')
+disp(f_2,'with frequency in hertz')
+disp('')
+disp('FOR Vout3:')
+disp(V_out_3_max,'peak value in volts')
+disp(V_out_3_phi*180/%pi,'phase angle in degrees')
+disp(f_3,'with frequency in hertz')