Well the last problem is part of this same assignment. Yes this is why I get frustrated it does take me this long. The problem with, the last problem is that I never got the mesh and nodal analysis to agree but I resolved to just send them in once I have everything completed and let them pick the bones out of that. Given the appalling standards of education these days no doubt I will get a pass just for making an effort!
Yes the explanation from Orolo is pretty good although I'm still a little confused about where his last 2 lines come from and have been trying to write down my own version of the calculations to see that I have it clear in my head. As much as I want to pass this course just to shut everybody up I would also like to learn something from it as I'm not your average student.
Next to the last line, Orolo is substituting in the value for e
1 in the numerator and for I'
2 in the denominator. Then there is some really ugly algebra in the last line. It takes a minute to see how the substitutions work out but they do.
I seem to get both mesh and nodal to agree but I'm only looking at Va and Vb.
wxMaxima code for Mesh Analysis
ratprint : false$
fpprintprec : 4$
eq1 : 0 = -V1 + Z1*I1 + Z4*(I1-I2) ;
eq2 : 0 = Z4*(I2-I1) + Z2*(I2-I4) + Z5*(I2-I3) ;
eq3 : 0 = Z5*(I3-I2) + Z3*I3 + V2 ;
eq4 : 0 = Z2*(I4-I2) + V3 ;
eq5 : VA = V1 - I1*Z1 ;
eq6 : VB = VA - V3 ;
eq7 : Z1 = 2 ;
eq8 : Z2 = -5*%i ;
eq9 : Z3 = 4 ;
eq10 : Z4 = -5*%i ;
eq11 : Z5 = 4*%i ;
eq12 : V1 = 120 ;
eq13 : V2 = 120*%i ;
eq14 : V3 = 14.14*%i + 14.14 ;
res : solve([eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8,eq9,eq10,eq11,eq12,eq13,eq14])$
results : expand(float(res))$
lngth : length(results[1])$
sorted : sort(results[1])$
print("")$
for i:1 thru lngth do
print(sorted[i])$
wxMaxima code for Nodal Analysis
ratprint : false$
fpprintprec : 4$
eq1 : 0 = -(V1-VA)/Z1 +(VA/Z4) + (VB / Z5) - (V2 - VB)/Z3;
eq2 : VA = VB + V3;
eq3 : V1 = 120;
eq4 : V2 = 120*%i;
eq5 : V3 = 14.14 + 14.14 * %i;
eq6 : Z1 = 2;
eq7 : Z3 = 4;
eq8 : Z4 = 0 -5 * %i;
eq9 : Z5 = 0 + 4 * %i;
res : solve([eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8,eq9])$
results : expand(float(res))$
lngth : length(results[1])$
sorted : sort(results[1])$
print("")$
for i:1 thru lngth do
print(sorted[i])$
For some strange reason, I can't do a copy and paste of the output. Nevertheless, for Va I get 86.38+j45.76 and for Vb 72.24+j31.62
One of the things I spent time on was getting the variables printed in sorted order with limited precision.
If you run out of things to do, you might try these in wxMaxima.