garnissagep1.m

基于MATLAB电机磁场程序

npt=size(p,2);	ntrg=size(t,2);
K=sparse(npt,npt);  b=sparse(npt,1);
Ke=zeros(3);Me=zeros(3);

% D閒inition de J
J=sparse(1,ntrg);
indice1=find((t(4,:)>=4) & (t(4,:)<=7));
indice11=find((t(4,:)>=16) & (t(4,:)<=19));
indice2=find((t(4,:)>=12) & (t(4,:)<=15));
indice21=find((t(4,:)>=24) & (t(4,:)<=27));
indice3=find((t(4,:)>=20) & (t(4,:)<=23));
indice31=find((t(4,:)>=8) & (t(4,:)<=11));

indicer1=find(t(4,:)==28);
indicer2=find(t(4,:)==29);

ij=1;

J(indice1)=Jex1*ij*ones(1,length(indice1));
J(indice11)=Jex11*ij*ones(1,length(indice11));
J(indice2)=Jex2*ij*ones(1,length(indice2));
J(indice21)=Jex21*ij*ones(1,length(indice21));
J(indice3)=Jex3*ij*ones(1,length(indice3));
J(indice31)=Jex31*ij*ones(1,length(indice31));

J(indicer1)=Jexr1*ij*ones(1,length(indicer1));
J(indicer2)=-Jexr1*ij*ones(1,length(indicer2));

% D閒inition de nu
nu=nuo*ones(1,ntrg);
%rotor mur=1000;
indice=find(t(4,:)==3);
nu(indice)=nuo/murr*ones(1,length(indice));
%stator mur =1000;
indice=find(t(4,:)==2);
nu(indice)=nuo/murs*ones(1,length(indice));
sigma=sparse(1,ntrg);
% indice=find((t(4,:)==28)|(t(4,:)==29));
% sigma(indice)=sigmaexc*ones(size(indice));