P449
Ans
Program
Here, we use current division to find i1 and then caculate Vth
%10.58
%Zth CD
z1=j*10;
z2=8-j*6;
Zth=z1*z2/(z1+z2) %parallel resistance
% Vth
vs=5.*exp(i*45/180*pi)
i1=z2/(z1+z2)*vs %using current division
Vth=i1*z1
disp('Ammplitude'),abs(Vth)
disp('Angle'), angle(Vth)*180/pi
Durjoy used node analysis, but he did not consider the direction of current.
So it is complex using Node analysis
Just using MESH analysis.
we use KVL to find the answer
% set current on 1/8F as i2, using kvl
% io*(4-j*4)-12+j4*i2=0; % loop1
%(j*2-j*8-j*4)*i2-j2*3*io+j*4*io=0 ;%loop2
A=[4-j*4 j*4;-j*2 -j*10]
B=[12 0]'
I=inv(A)*B
Vth=(-j*8)*I(2)
abs(Vth)
angle(Vth)*180/pi
10.67
z1=13-j*5;
z2=10;
z3=12;
z4=8+j*6;
z12=z1*z2/(z1+z2)%parrall
z34=z3*z4/(z3+z4)
Zth=z12+z34
vs=60.*exp(i*45/180*pi)
i1=vs/(z1+z2)
i2=vs/(z3+z4)
Vth=-i1*z1+i2*z3 %KVL method
abs(Vth)
angle(Vth)*180/pi
% you can use voltage difference
Vth=i1*z2-i2*z4
abs(Vth)
angle(Vth)*180/pi
% NORTON
IN=Vth/Zth
abs(IN)
angle(IN)*180/pi
Result
10.72
w=10^4;
vs=4;
zc=-j*10^5;
R=50*10^3;
vo=vs*zc/(R+zc)
io=vo/(100*10^3)
abs(io)
angle(io)*180/pi