На главную

Построение векторной диаграммы параллельной цепи (RL + C)

>   restart:

>   with(plottools):

>   Ubx:=70:#Входное напряжение цепи
R:=200:#Активное сопротивление цепи
L:=2:C:=8*10^(-6):#Индуктивность и емкость (yL=yC)Равенство реактивных проводимостей.
#L:=1.237:C:=4*10^(-6):#Индуктивность и емкость (XL<XC)
#L:=1.237:C:=12*10^(-6):#Индуктивность и емкость (XL>XC)
 w:=314:#Круговая частота

>   xL:=evalf(w*L,4);xC:=evalf(1/(w*C),4);

xL := 628.

xC := 398.1

>   Zk:=R+I*w*L;#Cопротивление катушки индуктивности

>   Zc:=evalf(1/(I*w*C));#Сопротивление емкости
Yk:=1/Zk:#Проводимость катушки индуктивности
Yc:=1/Zc:#Проводимость емкости

>   Y:=Yk+Yc:#Проводимость всей цепи
k:=300:#к - масштабный коэффициент тока на векторной диаграмме
i:=Ubx*Y;#Ток в неразветвленной части цепи

     ik:=Ubx*Yk;#Ток в цепи катушки индуктивности
  ic:=Ubx*Yc;#Ток в цепи емкости

     Ur:=ik*R;#Напряжение на активном сопротивлении катушки индуктивности
  Ul:=ik*(I*w*L);#Напряжение на индуктивности

Zk := 200+628*I

Zc := -398.0891720*I

i := .3222954805e-1+.7463921913e-1*I

ik := 875/27149-5495/54298*I

ic := .1758400000*I

Ur := 175000/27149-549500/27149*I

Ul := 1725430/27149+549500/27149*I

#Построение векторной диаграммы

>   a:=0:
b:=abs(Ubx):
sm:=b/10:#sm- смещение надписей
wb:=b/60:wh:=b/15:hh:=0.1:#Параметры для вычерчивания вектора
l1 := arrow([a,a], [Re(Ubx),Im(Ubx)], wb,wh,hh, color=green):#Вектор входного напряжения
l2 := arrow([a,a], [Re(i)*k,Im(i)*k], wb,wh,hh, color=blue):#Вектор тока в неразветвленной части цепи
l3 := arrow([a,a], [Re(ik)*k,Im(ik)*k], wb,wh,hh, color=blue):#Вектор ik
l4 := arrow([a,a], [Re(ic)*k,Im(ic)*k], wb,wh,hh, color=blue):#Вектор ic
l5 := arrow([a,a], [Re(Ur),Im(Ur)], wb,wh,hh, color=red):#Вектор Ur
l6 := arrow([a,a], [Re(Ul),Im(Ul)], wb,wh,hh, color=red):
#Вектор Ul
l7 := arrow([Re(Ul),Im(Ul)], [Re(Ubx),Im(Ubx)], 0,0,0, color=blue):
l8 := arrow([Re(Ur),Im(Ur)], [Re(Ubx),Im(Ubx)], 0,0,0, color=blue):
l9 := arrow([Re(ik)*k,Im(ik)*k], [Re(i)*k,Im(i)*k], 0,0,0, color=blue):
l10 := arrow([Re(ic)*k,Im(ic)*k], [Re(i)*k,Im(i)*k], 0,0,0, color=blue):
t1:=textplot([Re(ik)*k+sm,Im(ik)*k,`Ik`],font=[COURIER, BOLD,14]):
t2:=textplot([Re(ic*k)+sm,Im(ic)*k,`Ic`],font=[COURIER, BOLD,14]):
t3:=textplot([Re(i)*k+sm,Im(i)*k+sm,`I `], font=[COURIER,BOLD,14]):
t4:=textplot([Re(Ur)-sm,Im(Ur),`Ur `], font=[COURIER,BOLD,14]):
t5:=textplot([Re(Ubx)+sm,Im(Ubx)-sm,`Uвх `], font=[COURIER,BOLD,14]):
t6:=textplot([Re(Ul)-sm,Im(Ul)+sm,`Ul `], font=[COURIER,BOLD,14]):

>   with(plots):plots[display](l1,l2,l3,l4,l5,l6,l7,l8,l9,l10,t1,t2,t3,t4,t5,t6,axes=none);
y:=abs(1/(R+I*w*L)+I*w*c):
z:=1/y:

>   i:=abs(Ubx)*y*10000:

>   

[Maple Plot]

>   

Hosted by uCoz