%按接地、内部短路状况调整绕组联接图，数组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
基于matlab关于交流电机及其故障仿真，同步发电机故障诊断实验
版权申诉
5星 · 超过95%的资源 126 浏览量
2022-04-17
15:46:35
上传
评论
收藏 9KB RAR 举报 
关于交流电机及其故障仿真，同步发电机故障诊断实验.rar （9个子文件） 
2.关于交流电机及其故障仿真，同步发电机故障诊断实验 
lamda.m 3KB test_input_TsingHua.m 2KB motor_simulation_myself.m 3KB generator_simulation_myself.m 3KB l_damp.m 2KB GeneratorML_Teach.m 2KB mabz.m 1KB mabTeach.m 2KB map_subroutine.m 9KB资源评论
liu180226203192024-04-07这个资源值得下载，资源内容详细全面，与描述一致，受益匪浅。
