buduicheng.m

用MATLAB编的配电网自动化的一些程序

function q=buduicheng
syms zs_1 zs_2 zs_0 zf_1 zf_2 zf_0 u_a u_b u_c m_0 m_1 m_2 m_3 m_4 l mm tt
m_1=symadd(zf_1,zf_2)
m_0=symadd(m_1,zf_0)
m_3=symadd(symdiv(m_0,3),zs_1)
m_4=symdiv(zs_1,m_3)
u_a=symsub(1,m_4)%a xiangjiediduanlu
%u_sanxiang=

%u_a=symdiv(zs_1,symadd(zs_1,symdiv(symadd(symadd(zf_1,zf_2),zf_0)),3))))

u_a_0=subs(u_a,{zf_1,zf_2,zf_0},{symmul(0.278,l),symmul(0.278,l),symmul(1.478,l)})
u_a_1=subs(u_a_0,{zs_1,zs_2,zs_0},{0.663,0.663,0.663})
simple(u_a_1)
%u_a_2=numeric(u_a_1)
kk=solve(mm-u_a_1,l)
%solve(tt-kk,mm)

t=0:0.01:0.99
u_a_2=numeric(subs(kk,mm,t))
plot(t,u_a_2)

xlabel('临界电压(p.u.)')
ylabel('临界距离(km)')
title('不同a值的临界距离与临界电压的关系曲线')
axis([0,1,0,30])


syms u l_crit z zs a m n k_1 k_2
m=symmul(zs,u)
n=symmul(z,symsub(1,u))
k_2=symsub(1,symmul(u,symsub(1,cos(a))))
k_1=symmul(u,cos(a))

l_crit=symmul(symdiv(m,n),k_2)
l_tu1=subs(l_crit,{zs,z,a},{0.663,0.278,0})
l_tu2=subs(l_crit,{zs,z,a},{0.663,0.278,-30*pi/180})
l_tu3=subs(l_crit,{zs,z,a},{0.663,0.278,-60*pi/180})
l_tu4=subs(l_crit,{zs,z,a},{0.663,0.278,120*pi/180})
y=qujiandiv(0.663,0.278,0.05)

t=0:0.01:0.99

%y=sym([symmul(y(1,1),l_tu1);symmul(y(2,1),l_tu1);symmul(y(1,3),l_tu1)])
l_tu1_0=numeric(subs(l_tu1,u,t))
%plot(t,y)
plot(t,l_tu1_0,'r',t,y(1,1)*l_tu1_0,t,y(1,3)*l_tu1_0)
xlabel('临界电压(p.u.)')
ylabel('临界距离(km)')
title('不同a值的临界距离与临界电压的关系曲线')
axis([0,1,0,10])
l_tu2_0=numeric(subs(l_tu2,u,t))
l_tu3_0=numeric(subs(l_tu3,u,t))
l_tu4_0=numeric(subs(l_tu4,u,t))
plot(t,u_a_2,'*r',t,l_tu1_0,'b')
%plot(t,l_tu4_0,'r',t,l_tu1_0,t,l_tu3_0,t,l_tu2_0)
%plot(t,l_tu1,'r',t,l_tu2_0,':b',t,l_tu3_0,'-.k',t,l_tu4_0,'--m')
%grid
xlabel('临界电压(p.u.)')
ylabel('临界距离(km)')
title('a相短路接地与三相接地的临界距离与临界电压的关系曲线')
axis([0,1,0,30])