'Униполярная машина
'Униполярная машина представляет собой прямолинейный проводник с 
током,
'помещенный в магнитное поле
'Полагаем магнитное поле - перпендикулярно экрану;ток - от центра 
проводника
'к краю

screen 12
window (-320,-240) - (320,240):'поле экрана; центр экрана - 0,0
r=20:'длина проводника
mp=1:'масса проводника
f=2:'сила трения на краю токосъемника
'fp=4
'f=4
H=1:'напряженность магнитного поля
i=32:'сила тока в проводнике
dt=0.001:'квант времени
pi=4*atn(1):'Число Пи
fi=0:'Текущий угол ротора
fio=0:'Прежний угол
ff=0:'Угловая скорость ротора
fff=0:'Угловое ускорение ротора
?"Униполярная машина. <Esc> -выход"
def fnf(ff)=-sgn(ff)*(f*(-ff)):'сила трения,зависящая от скорости
circle (0,0),r+2,9
t=0
'z=0.02
per=8
t1=0*per
t2=.05*per
def fnt(t)=t-(int(t/per))*per
def fni(t)=i*(-1+sgn(t1-fnt(t)))*sgn(1-sgn(fnt(t)-t2))
?"Импульсы тока со скважностью=";int((t2/per)*10000)/100;"%"

krutim:
 t=t+dt
 circle (0,0),.2*r,12:
 'Сотрем изображение ротора в прежней позиции
 'x=r*cos(fio):y=r*sin(fio)
 'line (0,0) - (x,y),0
 'Нарисуем ротор в новом положении
 x=r*cos(fi):y=r*sin(fi)
 line (0,0) - (x,y),12
 if inkey$=chr$(27) then stop:'Прервать программу по <Esc>
 'сила, действующая на кусочек проводника длиной dr
 'df=fni(t)*H*dr
 'момент, действующий на кусочек проводника, нащодящегося на 
расстоянии
 'r от ценра; dm=df*r;тогда суммарный момент,m=summa(0-r)dm
 'Закон Ньютона в полярной системе координат -
 'МоментСилы=(УгловоеУскорение)*(МоментИнерции)
 'Определим функцию общего момента сил, на проводник

 def fnm(i)=(fni(t)*H*r/2)+(fnf(ff)*r)
' if abs(fnm(i))<abs(fnf(ff)) then locate 3,3:?"мало току":stop
 'момент инерции ротора in=масса ротора*(r/2)^2
 in=mp*((r/2)^2)
 'Угловое ускорение ротора - fff=(МоментСилы)/(МоментИнерции)
 fff=fnm(i)/in
 'Угловая скорость ротора ff=ff+fff*dt;где fff-угловое ускорение
 ff=ff+fff*dt

 'Приращение угла dfi=ff*dt; где ff-угловая скорость ротора
 dfi=ff*dt
 'Кладем старый угол ротора - предыдущий угол ротора
 fio=fi
 x=r*cos(fio):y=r*sin(fio)
 z=.01
 'delay z
 line(0,0)-(x,y),0
 fi=fi+dfi:'Угол через момент времени dt
 nn=abs(ff)

 'locate 2,1:?"     "
 'locate 2,2:?"Мощность-";f*r*2*pi*ff
 'locate 2,1:?"n=";int(3600*nn/(2*pi));"об/мин"
goto krutim