270
Practical MATLAB
® Applications for Engineers
ANALYTICAL Solution
X
j C
j
j
C ( )
( / )
ϭ
ϭ
ϭϪ
1
1
10 1 40
4
·
·
Ω
X = 3 − j4 (Ω)
Then, Z = X/ /R 2 .
Z
X
X
ϭ ϩ
ϭ
Ϫ Њ
* 10
10
3 67
36
. ∠
Ω
MATLAB Solution
% Script file: imp _ Z
W = 10;
R1 = 3;
R2 = 10;
C = 1/40;
XC = -j/(W*C);
X = R1+XC;
Z = X*R2/(X+R2);
Zmag = abs(Z);
Phase = angle(Z);
Phase _ deg = Phase*360/(2*pi);
disp(‘^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^’)
disp(‘Magnitude of Z (in Ohms) and phase (in degrees) is given by :’);
Zmag
Phase _ deg
disp(‘^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^’)
The script fi le imp_Z is executed and the results are shown as follows.
FIGURE 3.51
Network of Example 3.3.
Z
R 2 = 10 Ω
R 1 = 3 Ω
C = 1/40 F
w = 10 rad/s
CRC_47760_CH003.indd 270
CRC_47760_CH003.indd 270
7/23/2008 1:27:44 PM
7/23/2008 1:27:44 PM
Practical MATLAB
® Applications for Engineers
ANALYTICAL Solution
X
j C
j
j
C ( )
( / )
ϭ
ϭ
ϭϪ
1
1
10 1 40
4
·
·
Ω
X = 3 − j4 (Ω)
Then, Z = X/ /R 2 .
Z
X
X
ϭ ϩ
ϭ
Ϫ Њ
* 10
10
3 67
36
. ∠
Ω
MATLAB Solution
% Script file: imp _ Z
W = 10;
R1 = 3;
R2 = 10;
C = 1/40;
XC = -j/(W*C);
X = R1+XC;
Z = X*R2/(X+R2);
Zmag = abs(Z);
Phase = angle(Z);
Phase _ deg = Phase*360/(2*pi);
disp(‘^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^’)
disp(‘Magnitude of Z (in Ohms) and phase (in degrees) is given by :’);
Zmag
Phase _ deg
disp(‘^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^’)
The script fi le imp_Z is executed and the results are shown as follows.
FIGURE 3.51
Network of Example 3.3.
Z
R 2 = 10 Ω
R 1 = 3 Ω
C = 1/40 F
w = 10 rad/s
CRC_47760_CH003.indd 270
CRC_47760_CH003.indd 270
7/23/2008 1:27:44 PM
7/23/2008 1:27:44 PM
