Moluch 251 c indd

Matlab-Script: 
Преобразуем уравнение (1) для программирования в Matlab-Script: 






1
1
1
1
1
1
1
1
.
rк
r
э
r
r
s
sx
sx
mx
sx
к
sy
my
sy
s
m
s
s
s
s
r k
l
l k
T
s i
i
u
i
i
r l
r
r
r
r



 



    













Переходим к оригиналу (s = d/dt): 






1
1
1
1
1
1
1
1
1
.
sx
rк
r
э
r
r
sx
mx
sx
к
sy
my
sy
s
m
s
s
s
s
s
di
r k
l
l k
i
u
i
i
dt
r l
r
r
r
r
T



 







  



















Переходим к конечным разностям (метод Эйлера): 


 
 
 
 
   
 
 
   
1
1
1
1
1
1
1
1
1
1
.
sx
sx
rк
r
э
sx
mx
sx
к
sy
s
m
s
s
r
r
my
sy
s
s
s
i i
i i
r k
l
i i
i
u i
i i i
dt
r l
r
r
l k
i
i
i i i
r
r
T







 


 














 







Отсюда ток 
i
sx
в Matlab-Script определится следующим образом [1]: 


 
 
 
 
   
 
 
   
1
1
1
1
1
1
1
1
1
.
rк
r
э
sx
sx
sx
mx
sx
к
sy
s
m
s
s
r
r
my
sy
s
s
s
r k
l
i i
i i
i i
i
u i
i i i
r l
r
r
l k
dt
i
i
i i i
r
r
T









 
 














 







(1') 
Си: 
Исключив (i) и (i+1) из (1’), получим выражение 
i
sx
на языке Си: 
1
1
1
1
1
1
1
1
.
rк
r
э
r
r
sx
sx
sx
mx
sx
к
sy
my
sy
s
m
s
s
s
s
s
r k
l
l k
dt
i
i
i
u
i
i
r l
r
r
r
r
T



 








  





 
 

 







Видимо, если циклы будут составлять несколько сотен тысяч, то это скажется на быстродействии. 
2. Определение статорного тока i
sy 
 
Уравнение для определения тока 
i
sy
в Simulink, полученное в работе [1]: 


 



1
1
1
2
3
4
5
1
.
1
s
rк
r
sy
my
sy
э
к
sx
mx
r
r
sx
m
s
r
r k
i
u
l
i
l k
i
l
T
s



 









 

 
   



  






 


(2) 
Структурная схема реализации уравнения (2) приведена на рис. 2. 
Аналогично преобразуем выражение тока 
i
sy
для Matlab-Script [1]:
 


 
 
 


   
 
 
   
1
1
1
1
1
1
1
1
1
1
.
rк
r
э
sy
sy
sy
my
sy
к
sx
s
m
s
s
r
r
mx
sx
s
s
s
r k
l
i i
i i
i i
i
u i
i i i
r l
r
r
l k
dt
i
i
i i i
r
r
T









 
 




  








 







(2') 
Из уравнения (2’) получим выражение 
i
sy
для программирования на Си: 
1
1
1
1
1
1
1
1
.
rк
r
э
r
r
sy
sy
sy
my
sy
к
sx
mx
sx
s
m
s
s
s
s
s
r k
l
l k
dt
i
i
i
u
i
i
r l
r
r
r
r
T



 








  





 
 

 







 
3. Определение потокосцепления 
ψ
mx
 
 
В работе [1] была получена структурная схема для определения потокосцепления 
ψ
mx
в Simulink (рис. 3) по следу-
ющему уравнению: 






1
1
5
4
2
3
1
1
.
1
rк
э
r
mx
srк
sx
sx
к
my
my
r
sy
m
s
s
r
r
m
r
l
l
r
i
u
l
i
l
l
l k
k
T
s





 
 










 




  
  













 











 
(3) 

---

“Young Scientist”
. 
#13 (251)
. 
March 2019
7
Physics
Рис.
 
1.
 Структурная схема для определения тока 
i
sx
 в Simulink


 



1
1
2
4
3
5
1
1
,
1
s
rк
r
sx
mx
sx
э
к
sy
my
r
r
sy
m
s
r
r k
i
u
l
i
l k
i
l
T
s



 









 



   



  






  

(1) 
где 
1
1
s
s
б
T
T

 
— постоянная времени статорной обмотки. 
Matlab-Script: 
Преобразуем уравнение (1) для программирования в Matlab-Script: 






1
1
1
1
1
1
1
1
.
rк
r
э
r
r
s
sx
sx
mx
sx
к
sy
my
sy
s
m
s
s
s
s
r k
l
l k
T
s i
i
u
i
i
r l
r
r
r
r



 



    













Переходим к оригиналу (s = d/dt): 






1
1
1
1
1
1
1
1
1
.
sx
rк
r
э
r
r
sx
mx
sx
к
sy
my
sy
s
m
s
s
s
s
s
di
r k
l
l k
i
u
i
i
dt
r l
r
r
r
r
T



 







  



















Переходим к конечным разностям (метод Эйлера): 


 
 
 
 
   
 
 
   
1
1
1
1
1
1
1
1
1
1
.
sx
sx
rк
r
э
sx
mx
sx
к
sy
s
m
s
s
r
r
my
sy
s
s
s
i i
i i
r k
l
i i
i
u i
i i i
dt
r l
r
r
l k
i
i
i i i
r
r
T







 


 














 







Отсюда ток 
i
sx
в Matlab-Script определится следующим образом [1]: 


 
 
 
 
   
 
 
   
1
1
1
1
1
1
1
1
1
.
rк
r
э
sx
sx
sx
mx
sx
к
sy
s
m
s
s
r
r
my
sy
s
s
s
r k
l
i i
i i
i i
i
u i
i i i
r l
r
r
l k
dt
i
i
i i i
r
r
T









 
 














 







(1') 
Си: 
Исключив (i) и (i+1) из (1’), получим выражение 
i
sx
на языке Си: 
1
1
1
1
1
1
1
1
.
rк
r
э
r
r
sx
sx
sx
mx
sx
к
sy
my
sy
s
m
s
s
s
s
s
r k
l
l k
dt
i
i
i
u
i
i
r l
r
r
r
r
T



 








  





 
 

 







Видимо, если циклы будут составлять несколько сотен тысяч, то это скажется на быстродействии.