• 1475497788@qq.com
    了解作者
  • matlab
    开发工具
  • 9KB
    文件大小
  • rar
    文件格式
  • 0
    收藏次数
  • 10 积分
    下载积分
  • 12
    下载次数
  • 2018-09-25 16:40
    上传日期
MATLAB 关于交流电机及其故障仿真程序,清华大学博士生做的同步发电机故障诊断实验.
MATLAB 关于交流电机及其故障仿真程序.rar
  • 交流电机matlab
  • test_input_TsingHua.m
    1.7KB
  • l_damp.m
    1.8KB
  • lamda.m
    3KB
  • mabz.m
    1.2KB
  • GeneratorML_Teach.m
    2.1KB
  • generator_simulation_myself.m
    3.1KB
  • mabTeach.m
    1.9KB
  • map_subroutine.m
    8.6KB
  • motor_simulation_myself.m
    2.9KB
内容介绍
%按接地、内部短路状况调整绕组联接图,数组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
评论
    相关推荐