dessinp1.m

基于MATLAB电机磁场程序

Epai1=5*entref; Epai2=Epai1;
th=pi/24;
x1=Rs*sin(th); 
y1=Rs*cos(th);	
x2=1.25e-3;		
y2=sqrt(Rs*Rs-x2*x2);	
xr1=Rr*sin(beta1); 
yr1=Rr*cos(beta1); 
xr2=0; 
yr2=Rr; 
xr3=-xr1; 
yr3=yr1; 
xr4=-Rr*sin(beta); 
yr4=yr3; 
xr5=xr4; 
yr5=yr3-Epai1; 
xr6=xr4; 
yr6=-yr5; 
xr7=xr6; 
yr7=-yr4; 
xr8=xr3; 
yr8=-yr3; 
xr9=0; 
yr9=-Rr; 
xr10=xr1; 
yr10=-yr1; 
xr11=-xr7; 
yr11=yr7; 
xr12=-xr6; 
yr12=yr6; 
xr13=-xr5; 
yr13=yr5; 
xr14=xr13; 
yr14=yr4; 

xbr1=xr13+Epai2; 
ybr1=yr13; 
xbr2=xbr1; 
ybr2=-ybr1; 
xbr3=-xbr1; 
ybr3=ybr1; 
xbr4=xbr3; 
ybr4=-ybr3; 

%yam=(xr4-xr3)/2-xr4; yam=sqrt(yam^2+yr3^2); ram=1.25*(Rr-yam)/2; 


cus=[1		1			1			1		;
     Rc	   0		  -Rc			0		;
     0	  -Rc			0		   Rc		;
     0		Rc			0		  -Rc		;
     Rc		0		  -Rc			0		;
     2		2			2			2		;
     0		0			0			0		;
     0		0			0			0		;
     0		0			0			0		;
     Rc		Rc			Rc			Rc ]	;


st1 =[         1            2            1            2            1            2 
  		         x1  1.2500e-003  1.2500e-003  2.5000e-003  3.9500e-003 -3.9500e-003
  		         x2  1.2500e-003  2.5000e-003  3.9500e-003 -3.9500e-003 -2.5000e-003  
  		         y1  			 y2  3.9965e-002  4.1215e-002  4.9400e-002  4.9400e-002
  		         y2  3.9965e-002  4.1215e-002  4.9400e-002  4.9400e-002  4.1215e-002
               1            1            4            4            4            4
               2            2            2            2            2            2
            	0            0   1.2500e-003           0            0            0
            	0            0   4.1215e-002           0   4.9400e-002           0
      3.9385e-002           0   1.2500e-003           0   3.9500e-003           0];
            
st2=[          1           2              2            1 
    -2.5000e-003 -1.2500e-003   -1.2500e-003 -1.2500e-003 
    -1.2500e-003  1.2500e-003   -1.2500e-003          -x1
     4.1215e-002  3.9965e-002    3.9965e-002           y2
     3.9965e-002  3.9965e-002             y2           y1
               4           4              1            1
               2           1              2            2  
    -1.2500e-003           0              0            0  
     4.1215e-002           0              0            0           
     1.2500e-003           0              0   3.9385e-002];

gg=[st1 st2];

gs=gg; 
   
  for i=1:23;
      aes=pi/12*i;
      gg(6,3:8)=4+i;
      gs= [gs groto(gg,aes)];
  end
   
gsf=[cus gs];   
gr1= [1 xr1  xr2  yr1  yr2  3 1  0 0 Rr]';
gr2= [1 xr2  xr3  yr2  yr3  3 1  0 0 Rr]';
gr3= [2 xr3  xr4  yr3  yr4  3 1  0 0 0]';
gr4= [2 xr4  xr5  yr4  yr5  3 1  0 0 0]';
gr5= [2 xr5  xr6  yr5  yr6  3 28 0 0 0]';
gr6= [2 xr6  xr7  yr6  yr7  3 1  0 0 0]';
gr7= [2 xr7  xr8  yr7  yr8  3 1  0 0 0]';
gr8= [1 xr8  xr9  yr8  yr9  3 1  0 0 Rr]';
gr9= [1 xr9  xr10 yr9  yr10 3 1  0 0 Rr]';
gr10=[2 xr10 xr11 yr10 yr11 3 1  0 0 0]';
gr11=[2 xr11 xr12 yr11 yr12 3 1  0 0 0]';
gr12=[2 xr12 xr13 yr12 yr13 3 29  0 0 0]';
gr13=[2 xr13 xr14 yr13 yr14 3 1  0 0 0]';
gr14=[2 xr14 xr1  yr14 yr1  3 1  0 0 0]';

gr=[gr1 gr2 gr3 gr4 gr5 gr6 gr7 gr8 gr9 gr10 gr11 gr12 gr13 gr14];

bobr1=[2 xbr4 xr6  ybr4 yr6  28 1 0 0 0]';
bobr3=[2 xr5  xbr3 yr5  ybr3 28 1 0 0 0]';
bobr4=[2 xbr3 xbr4 ybr3 ybr4 28 1 0 0 0]';

bobr5=[2 xbr1 xr13 ybr1 yr13 29 1 0 0 0]';

bobr7=[2 xr12 xbr2 yr12 ybr2 29 1 0 0 0]';
bobr8=[2 xbr2 xbr1 ybr2 ybr1 29 1 0 0 0]';

bob1rot=[bobr1 bobr3 bobr4];
bob2rot=[bobr5 bobr7 bobr8];

bobrot=[bob1rot bob2rot];
ram=0; yam=0;
gam1=[1 ram   0   yam  yam+ram 30 3  0 yam ram]';
gam2=[1  0 -ram  yam+ram  yam  30 3  0 yam ram]';
gam3=[1 -ram  0   yam  yam-ram 30 3  0 yam ram]';
gam4=[1  0  ram  yam-ram  yam  30 3  0 yam ram]';
am=[gam1 gam2 gam3 gam4];
gam=am;
for i=1:2;
      amo=pi/11*i;
      gam= [gam groto(am,amo)];
end

 for i=1:2;
      amo=-pi/10*i;
      gam= [gam groto(am,amo)];
 end 
   gam=[gam groto(gam,pi)];
   gr=[gr bobrot];