script
%              Расчет характеристик схем включения  двигателя постоянного  тола
%                      последовательного возбуждения ДПВ. 
%                             
 prog='redpv';                                   
%                                              ДВИГАТЕЛЬ
   
%Каталожные данные двигателя (при ПВ=40%):

tip='D32'
Un=220;               %Напряжение номинальное  ,В	             	           
Pn=9500;             %Мощность номинальная , Вт
nn=760;               %Скорость номинальная ь,об/мин
In=53;                   %Ток  номинальный,   А
Mmax=677;         %Момент максимальный , Нм
%;Jdv=10.25;      %Момент инерции якоря,кгм^2

%Обмоточные данные двигателя

roa=0.2;               %Сопротивление якоря невыключаемое ,Ом
rdp=0.08;             %Cопротивление дополнительных полюсов, Ом	
rov=0.0972;         %Cопротивление обмотки возбуждения,Ом
ra=roa+rdp+rov;  %= 0.2+0.08+0,0972=0.377; 
       
%Базовые значения 

wn=nn/9.55
Mn=Pn/wn

%Каталожные кривые в границах прямолинейных участков

B=[1  2  3  4  5  6  7  8  9  10];
I5=[15  20	32	40	47	53	60	80	100 150	] ;          
n5=[1850 1400 1000 880	810	760	720	630	550	400];
M5=[ 1.5  2.5   6   8	10	12	14.5  22	29  48];

% Пересчет каталожных данных в систему СИ

 M51(B)=[M5]*9.81; %Момент на валу, Нм
 %M51(B)
 w5(B)=[n5]/9.55; %Скорость, рад/с 
 %w5(B)
 kF(B)=(Un-I5*ra)./w5(B);
 %kF(B)
 M52(B)=kF(B).*I5; %Момент электромагнитный, нм
% M52(B)
 M53(B)=2*M52(B) - M51(B); % Момент на валу тормозной, Нм
 %M53(B)

% Перевод в о.е.
   I51(B)=I5(B)/In;             % Ток якоря
  w51(B)=w5(B)/wn;        %Скорость
  M511(B)=M51(B)/Mn;    %Момент на валу
  M521(B) =M52(B)/Mn;   %Момент электромагнитный
  M531(B)= M53(B)/Mn;   % Момент на валу тормозной
  
  %tab=[I51(5) M511(5)  w51(5) M521(5)  M531(5)]
 %tab1=[I51(7) M511(7)  w51(7) M521(7)  M531(7)] 
 
 figure (1)  % Каталожные и расчетные характеристики
 
 axis ([0   3   0   6])       %рамка
 
plot  (I51(B),M511(B),'*-k', I51(B),w51(B),'*-b',I51(B),M521(B),'*-r', I51(B),M531(B),'*-g')
 hold on
 grid on
     xlabel('I')
    ylabel('w,M') 
     title('w=f(I)-b, M=f(I)-r');  
     
   %Режим работы:в заданной точке 
   
  % ПРОВЕРИТЬ ЗНАК СКОРОСТИ wzad
  
 % rr=1; U1= 110       %- двигательный режим при R1=var  
  rr=2; R2= 2.8        % - двигательный режим при  U= Un и R2=var
 
 % УСТАНОВИТЬ ЗНАК СКОРОСТИ wzad
 
 % rr=3; R3=8           % - торможение противовключением при  U= Un и  R3=var
 % rr=4; R4=2.3       % - динамическое с самовозбуждением   при R4=var 
 % rr=5; R5= 2.4      % - динамическое с независимым   возбуждением при Ioв = Iн и R5=var    
            

 % rr=5  
  
 % Расчет параметров схем включения
  
 Mzad=0.8  
wzad=0.5   

if  Mzad*wzad>0   rg=1,           %двигательный режим 
else rg=2 ;         %тормозной режим
end

U=0; R=ra;
            
 
 if  rg==1    D=B;  D=1;
 
       while M511(D)<Mzad, D=D+1;        end 
 
     kM51= (I51(D)- I51(D-1)) / (M511(D) - M511(D-1));
     Izad=I51(D-1)+kM51*(Mzad-M511(D-1)) 
  
      kw51= (w51(D)- w51(D-1)) / (M511(D) - M511(D-1));
     wezad=w51(D-1)+kw51*(Mzad-M511(D-1)); 
   
   Dw=wzad/wezad
   Ee=Un-Izad*In*ra
   
if rr==1  
    
   % U1 = Dw*Ee+Izad*In*ra 
    
Ezad(B)=(U1-ra.*I5(B));
Eezad(B)=(Un-ra.*I5(B));
ke2(B)=Ezad(B)./Eezad(B);
w52(B)=w5(B).*ke2(B)/wn;

figure (2)
plot  (I51(B),w51(B),'*-k', I51(B),w52(B), '*-b', Izad,wzad, 'o r')
end

if rr==2   
    
    %R2= (Un-Dw*Ee)/(Izad*In)  

Ezad(B)=(Un-R2.*I5(B));
Eezad(B)=(Un-ra.*I5(B));
ke3(B)=Ezad(B)./Eezad(B);
w53(B)=w5(B).*ke3(B)/wn;

figure (2)
plot ( I51(B),w51(B),'*-k', I51(B),w52(B), '*-b', I51(B),w53(B), '*--b', Izad,wzad, 'o r')

end 

 else
 D1=B;
      D1=1;
       while M531(D1)<Mzad,  D1=D1+1;    end            
 
     kM511= (I51(D1)- I51(D1-1)) / (M531(D1) - M531(D1-1));         
     Izad1=I51(D1-1)+kM511*(Mzad-M531(D1-1)) 
  
     kw51= (w51(D1)- w51(D1-1)) / (M531(D1) - M531(D1-1));
     wezad1=w51(D1-1)+kw51*(Mzad-M531(D1-1)) 
     
   wzad1 =wzad;
   Dw1=wzad1/wezad1;
   Ee1=Un-Izad1*In*ra;
    
if rr==3
    
    %R3= (Un-Dw1*Ee1)/(Izad1*In)  

Ezad(B)=(Un-R3.*I5(B));
Eezad(B)=(Un-ra.*I5(B));
ke4(B)=Ezad(B)./Eezad(B);
w54(B)=w5(B).*ke4(B)/wn;

figure (2)
plot ( I51(B),w51(B),'*-r', I51(B),w52(B), '*-b', I51(B),w53(B), '*--b', I51(B),w54(B), '*-.m' , Izad1,wzad1, 'or')
end

if rr==4   
    %R4= (-Dw1*Ee1)/(Izad1*In)  

Ezad(B)=-(R4.*I5(B));
Eezad(B)=(Un-ra.*I5(B));
ke5(B)=Ezad(B)./Eezad(B);
w55(B)=w5(B).*ke5(B)/wn;

figure (2)
plot ( I51(B),w51(B),'*-r', I51(B),w52(B), '*-b', I51(B),w53(B), '*--b', I51(B),w54(B), '*-k', I51(B), w55(B),'*-.m' ,  Izad1,wzad1, 'or' )

end 

if rr==5   
    % Iv=In; wezad=wn; 
    %R5= -(Un-In*ra)*wzad1/(Izad1*In)  

Ezad(B)=(-R5.*I5(B));

Eezad(B)=(Un-ra.*In);
ke6(B)=Ezad(B)./Eezad(B);
w56(B)=ke6(B);

figure (2)
plot ( I51(B),w51(B),'*-r', I51(B),w52(B), '*-b', I51(B),w53(B), '*--b', I51(B),w54(B), '*-k', I51(B),w55(B), '*--m', I51(B),w56(B), '*-m', Izad,wzad2, 'o r', Izad1,wzad1, 'or' )

end 

 end
  axis ([ 0   3   -1.5   3])       %рамка

 %plot  (I51(B),M511(B),'*-k', I51(B),w51(B),'*-b',I51(B),M521(B),'*-r', I51(B),M531(B),'*-g',)
% figure (2)
%plot ( I51(B),w52(B),'-b', I51(B),w53(B),'-b', Izad,wzad, 'or') 
  
 hold on
 grid on
     xlabel('I')
    ylabel('w,M') 
     title('w=f(I)-b, M=f(I)-r');      
     
     
for i=-0.5:0.001:3 plot(i,0,'k'); end;  %оси
for i=-4.0:0.001:3 plot(0,i,'k'); end;

 %  y-Жё   m-Фи   c-Го   r-Кр   g-Зе   b-Си   w-Бе   k-Чё
 
