Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- for i = 1:7
- Y(i) = input(['Podaj wartosc Y ', num2str(i),': ']);
- end
- for i = 1:4
- E(i) = input(['Podaj wartosc E ', num2str(i),': ']);
- end
- I(1) = input('Podaj wartosc I1: ');
- %Y(1) = 1;Y(2) = 2;Y(3) = 3;Y(4) = 4;Y(5) = 5;Y(6) = 6;Y(7) = 7;E(1) = 10;E(2) = 20;E(3) = 30;E(4) = 40;I(1) = 50;
- Y
- E
- I
- A = [
- (Y(1) + Y(2) + Y(7)), (-Y(2)), (-Y(1));
- (-Y(2)), (Y(2) + Y(3) + Y(4)), (-Y(3));
- (-Y(1)), (-Y(3)), (Y(1) + Y(3) + Y(5) + Y(6));
- ]
- B = [
- ((-E(1) .* Y(1)) - (E(2) .* Y(7)));
- I(1);
- ((-E(4) .* Y(6)) + (E(3) .* Y(5)) + (E(1) .* Y(1)));
- ]
- Ainv = inv(A)
- V = Ainv * B
- A * V
- I1 = ((V(1) - V(3) + E(1)) * Y(1));
- I2 = ((V(2) - V(1)) * Y(2));
- I3 = ((V(2) - V(3)) * Y(3));
- I4 = V(2) * Y(4);
- I5 = ((-V(3) + E(3)) * Y(5));
- I6 = ((V(3) + E(4)) * Y(6));
- I7 = ((V(1) + E(2)) * Y(7));
- IV1 = [
- -I1;
- I2;
- - I7;
- ]
- sum(IV1)
- IV2 = [
- -I2;
- -I3;
- I(1);
- -I4;
- ]
- sum(IV2)
- IV3 = [
- -I6;
- I1;
- I3;
- I5;
- ]
- sum(IV3)
Advertisement