MATLAB 关于交流电机及其故障仿真程序.rar

%按接地、内部短路状况调整绕组联接图，数组map_1存储调整后的绕组联接图信息 if((isc==0|isc==1)|isc2>1) column=fix(N_stator/sector); for(row=1:1:m_stator) for(i=1:1:(sector-1)) for(j=((i-1)*column+1):1:(i*column)) map_1((sector*(row-1)+i),j)=map(row,j); end for(j=((sector-1)*column+1):1:N_stator) map_1(sector*row,j)=map(row,j); end end end end if(isc==2&isc2==0) if(n_coil==0) for(row=1:1:m_stator) for(i=1:1:(sector-1)) for(j=((i-1)*n_sector+1):1:(i*n_sector)) map_1((sector*(row-1)+i),j)=map(row,j); end end for(j=((sector-1)*n_sector+1):1:N_stator) map_1((sector*row),j)=map(row,j); end end else for(row=1:1:m_stator) for(i=1:1:s_sector) for(j=((i-1)*n_sector+1):1:(i*n_sector)) map_1((sector*(row-1)+i),j)=map(row,j); end end for(j=(s_sector*n_sector+1):1:(s_sector*n_sector+n_coil)) map_1((sector*(row-1)+s_sector+1),j)=map(row,j); end for(j=(s_sector*n_sector+n_coil+1):1:((s_sector+1)*n_sector)) map_1((sector*(row-1)+s_sector+2),j)=map(row,j); end for(i=(s_sector+2):1:(sector-1)) for(j=((i-1)*n_sector+1):1:(i*n_sector)) map_1((sector*(row-1)+i+1),j)=map(row,j); end end end end end if(isc==0&isc2==1) if(n_coil==0&n_coil1==0) for(row=1:1:m_stator) for(i=1:1:(sector-1)) for(j=((i-1)*n_sector+1):1:(i*n_sector)) map_1((sector*(row-1)+i),j)=map(row,j); end end for(j=((sector-1)*n_sector+1):1:N_stator) map_1((sector*row),j)=map(row,j); end end end if((n_coil~=0&n_coil1==0)|(n_coil==n_coil1&n_coil~=0)) for(row=1:1:m_stator) for(i=1:1:s_sector) for(j=((i-1)*n_sector+1):1:(i*n_sector)) map_1((sector*(row-1)+i),j)=map(row,j); end end for(j=(s_sector*n_sector+1):1:(s_sector*n_sector+n_coil)) map_1((sector*(row-1)+s_sector+1),j)=map(row,j); end for(j=(s_sector*n_sector+n_coil+1):1:((s_sector+1)*n_sector)) map_1((sector*(row-1)+s_sector+2),j)=map(row,j); end for(i=(s_sector+2):1:(sector-1)) for(j=((i-1)*n_sector+1):1:(i*n_sector)) map_1((sector*(row-1)+i+1),j)=map(row,j); end end end end if(n_coil==0&n_coil1~=0) for(row=1:1:m_stator) for(i=1:1:s_sector1) for(j=((i-1)*n_sector+1):1:(i*n_sector)) map_1((sector*(row-1)+i),j)=map(row,j); end end for(j=(s_sector1*n_sector+1):1:(s_sector1*n_sector+n_coil1)) map_1((sector*(row-1)+s_sector1+1),j)=map(row,j); end for(j=(s_sector1*n_sector+n_coil1+1):1:((s_sector1+1)*n_sector)) map_1((sector*(row-1)+s_sector1+2),j)=map(row,j); end for(i=(s_sector1+2):1:(sector-1)) for(j=((i-1)*n_sector+1):1:(i*n_sector)) map_1((sector*(row-1)+i+1),j)=map(row,j); end end end end if(n_coil~=0&n_coil1~=0&n_coil~=n_coil1) if(s_sector~=s_sector1) if(s_sector<s_sector1) s_sector_a=s_sector; s_sector_b=s_sector1; n_coil_a=n_coil; n_coil_b=n_coil1; end if(s_sector>s_sector1) s_sector_a=s_sector1; s_sector_b=s_sector; n_coil_a=n_coil1; n_coil_b=n_coil; end for(row=1:1:m_stator) for(i=1:1:s_sector_a) for(j=((i-1)*n_sector+1):1:(i*n_sector)) map_1((sector*(row-1)+i),j)=map(row,j); end end for(j=(s_sector_a*n_sector+1):1:(s_sector_a*n_sector+n_coil_a)) map_1((sector*(row-1)+s_sector_a+1),j)=map(row,j); end for(j=(s_sector_a*n_sector+n_coil_a+1):1:((s_sector_a+1)*n_sector)) map_1((sector*(row-1)+s_sector_a+2),j)=map(row,j); end for(i=(s_sector_a+2):1:s_sector_b) for(j=((i-1)*n_sector+1):1:(i*n_sector)) map_1((sector*(row-1)+i+1),j)=map(row,j); end end for(j=(s_sector_b*n_sector+1):1:(s_sector_b*n_sector+n_coil_b)) map_1((sector*(row-1)+s_sector_b+2),j)=map(row,j); end for(j=(s_sector_b*n_sector+n_coil_b+1):1:((s_sector_b+1)*n_sector)) map_1((sector*(row-1)+s_sector_b+3),j)=map(row,j); end for(i=(s_sector_b+4):1:sector) for(j=((i-3)*n_sector+1):1:((i-2)*n_sector)) map_1((sector*(row-1)+i),j)=map(row,j); end end end end if(s_sector==s_sector1) s_sector_a=s_sector1; s_sector_b=s_sector; if(n_coil<n_coil1) n_coil_a=n_coil; n_coil_b=n_coil1; else n_coil_a=n_coil1; n_coil_b=n_coil; end for(row=1:1:m_stator) for(i=1:1:s_sector_a) for(j=((i-1)*n_sector+1):1:(i*n_sector)) map_1((sector*(row-1)+i),j)=map(row,j); end end for(j=(s_sector_a*n_sector+1):1:(s_sector_a*n_sector+n_coil_a)) map_1((sector*(row-1)+s_sector_a+1),j)=map(row,j); end for(j=(s_sector_a*n_sector+n_coil_a+1):1:(s_sector_a*n_sector+n_coil_b)) map_1((sector*(row-1)+s_sector_a+2),j)=map(row,j); end for(j=(s_sector_a*n_sector+n_coil_b+1):1:((s_sector_a+1)*n_sector)) map_1((sector*(row-1)+s_sector_a+3),j)=map(row,j); end for(i=(s_sector_a+4):1:sector) for(j=((i-3)*n_sector+1):1:((i-2)*n_sector)) map_1((sector*(row-1)+i),j)=map(row,j); end end end end end end if(isc==3) if(coil==0|coil==N_stator) for(row=1:1:m_stator) for(i=1:1:N_stator) map_1(row,i)=map(row,i); end end else column=coil; for(row=1:1:m_stator) for(i=1:1:(N_stator-column)) map_1((2*row-1),i)=map(row,i); end for(i=(N_stator-column+1):1:N_stator) map_1((2*row),i)=map(row,i); end end end end if(sector==1) map_1=map; end m_stator2=m_stator/3;%故障前的每相支路数 m_stator3=m_stator1/m_stator;%定子绕组每一支路所划分的段数 cam=ca/(N_stator0*m_stator2);%每一线圈对地等效电容，ca，cb，cc分别为每相绕组对地等效电容 cbm=cb/(N_stator0*m_stator2); ccm=cc/(N_stator0*m_stator2); ck=zeros((m_stator3+2),1);%一条支路上各段绕组所包含的匝数[0,ck,0]，中性点和机端为0 cs=zeros(4,(m_stator3+1));%定子每相绕组一条支路的电容分布A-B-C，机端－中性点 for(i=2:1:(m_stator3+1)) f