На главную

>    [Maple OLE 2.0 Object]

>    restart;
with(plottools):

>    a:=i1+i2=i3;
b:=z1*i1+z3*i3=e1;
c:=z2*i2+z3*i3=e2;

a := i1+i2 = i3

b := z1*i1+z3*i3 = e1

c := z2*i2+z3*i3 = e2

>   

>    eqns:={a,b,c}:#Система  уравнений
soln:=solve(eqns,{i1,i2,i3});#решение системы уравнений

soln := {i2 = (-z3*e1+z1*e2+e2*z3)/(z2*z1+z2*z3+z1*z3), i3 = (z1*e2+z2*e1)/(z2*z1+z2*z3+z1*z3), i1 = (z3*e1-e2*z3+z2*e1)/(z2*z1+z2*z3+z1*z3)}

>    assign(soln);

>    i1;#Приравниваем переменным их значения

>    i2;

>    i3;

>   

(z3*e1-e2*z3+z2*e1)/(z2*z1+z2*z3+z1*z3)

(-z3*e1+z1*e2+e2*z3)/(z2*z1+z2*z3+z1*z3)

(z1*e2+z2*e1)/(z2*z1+z2*z3+z1*z3)

>    e1:=e1m*exp(I*Pi/3);#Комплексное значение напряжения e1

e1 := e1m*(1/2+1/2*I*3^(1/2))

>    # I - мнимая единица

>    e2:=e2m*exp(-I*Pi/6);#Комплексное значение напряжения e2

e2 := e2m*(1/2*3^(1/2)-1/2*I)

>    e1m:=280:#Численные значения величин

>    e2m:=170:

>    z1:=60:

>    z2:=40:

>    #z3:=30-I*35;

>    z3:=30+I*35;

z3 := 30+35*I

>   

>   

>   

>    abs(evalf(e1));#Проверка значения e1

>    evalf(e1,3);#Значение е1 - алгебраич. запись комплексного числа

280.0000000

140.+242.*I

>    abs(evalf(e2));#Проверка значения e2

>    evalf(e2,3);#Значение е2 - алгебраич. запись комплексного числа

170.0000000

147.-85.*I

>   

>    u3:=z3*i3:

>    evalf(u3,7);#Напряжение на участке ab

89.83246+60.87816*I

>   

>   

>   

>   

>    U3:=abs(evalf(u3));#Абсолютное значение U на участке ab

U3 := 108.5173732

>    evalf(i1,3);#Комплексное выражение тока i1

.844+3.02*I

>    I1:=abs(%);#Величина тока i1

I1 := 3.135719375

>    evalf(i2,3);#Комплексное выражение тока i2

1.42-3.64*I

>    I2:=abs(%);#Величина тока i2

I2 := 3.907172891

>    evalf(i3,3);#Комплексное выражение тока i3

2.26-.620*I

>    I3:=abs(%);#Величина тока i3

I3 := 2.343501654

>    tgu3:=evalf(Im(u3)/Re(u3));#Тангенс угла напряжения Uab

tgu3 := .6776853175

>    tgi3:=evalf(Im(i3)/Re(i3));#Тангенс угла тока i3

tgi3 := -.2730773893

>   

>   

>   

>    Fiu3:=evalf(arctan(tgu3));#Фазовый угол Uab

Fiu3 := .5955921569

>    Fii3:=evalf(arctan(tgi3));#Фазовый угол i3

Fii3 := -.2665778979

>    cosab:=cos(Fiu3-Fii3);#Косинус угла ab

cosab := .6507913734

>    FiabRad:=Fiu3-Fii3;#Угол ab в радианах

>   

FiabRad := .8621700548

>    FiabGrad:=evalf(FiabRad*180/Pi);#Угол ab в градусах

FiabGrad := 49.39870536

>    P:=evalf(abs(i3)*abs(u3)*cosab);#Мощность на участке ab

P := 166.2496982

>    a:=0:b:=220:wb:=b/40:wh:=b/10:hh:=0.1:
#wb - width of the body of the arrow, wh- width of the head of the #arrow
#hh - height of the head of the arrow as a ratio of the length of the #body
l1 := arrow([a,a], [Re(i1)*100,Im(i1)*100], wb,wh,hh, color=green):
t1:=textplot([Re(i1)*100+20,Im(i1)*100,`i 1`]):
l2 := arrow([a,a], [Re(i2)*100,Im(i2)*100], wb,wh,hh, color=green):
t2:=textplot([Re(i2)*100+20,Im(i2)*100,`i 2`]):
l3 := arrow([a,a], [Re(i3)*100,Im(i3)*100], wb,wh,hh, color=green):
t3:=textplot([Re(i3)*100+20,Im(i3)*100,`i 3`]):
E1:= arrow([a,a], [Re(e1),Im(e1)], wb,wh,hh, color=red):
t4:=textplot([Re(e1)+20,Im(e1),`e 1`]):
E2:= arrow([a,a], [Re(e2),Im(e2)], wb,wh,hh, color=red):
t5:=textplot([Re(e2)+20,Im(e2),`e 2`]):
Uab:= arrow([a,a], [Re(u3),Im(u3)], wb,wh,hh, color=red):
t6:=textplot([Re(u3)+20,Im(u3),`U ab`]):

>   

>   

>   

>    with(plots):#Только здесь и не ранее!
plots[display](l1,l2,l3,E1,E2,Uab,t1,t2,t3,t4,t5,t6,axes=none,
scaling=constrained);

Warning, the names arrow and changecoords have been redefined

[Maple Plot]

>   

>   

>   

>   

>   

>   

>   

>   

>   

>   

>   

>   

>   

>   

>   

>   

>   

>   

>