blob: 5165ac650323336004f3986ec3aa6a28afd74108 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
|
clc
//ex 6.1
// given V_in(t)=2*cos(2000*%pi*t+A), A=40*%pi/180
w=2000*%pi; //omega
f=w/(2*%pi); //frequency
A=40*%pi/180; //40 degrees in radians
//equation of straight line of H_magnitude vs f is x+1000*y-4000=0
H_max=(4000-f)/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_phi=%pi*f/6000; //phase angle of H
H=H_max*complex(cos(H_phi),sin(H_phi));
V_in=2*complex(cos(A),sin(A)); //input voltage phasor
V_out=H*V_in; //output voltage phasor
V_out_R=real(V_out); //real part
V_out_I=imag(V_out); //imaginary part
V_out_max=sqrt((V_out_R^2)+(V_out_I^2)); //peak value
V_out_phi=atan(V_out_I/V_out_R);
disp(V_out_max,'peak value of Vout in volts')
disp(V_out_phi*180/%pi,'phase angle of Vout in degrees')
disp(f,'with frequency equal to')
|
