Page 129 - 363_

130

окрему процедуру. Для цього спочатку дещо перетворимо початкове рівняння,
записуючи його у вигляді:
    sin( )  ( S  , , ),

де штрих – позначення похідної за безрозмірним часом  =  0t,

mgl
 0  ,
J

а через S(, ,  ') позначено деяку задану функцію безрозмірного часу, куту

повороту маятника та його безрозмірної швидкості:
d
  .
d
У поданому випадку ця функція набуває вигляду:

, t ( S  ,  )  2   n mx  sin(     x ) cos  n my  sin(     y )  sin ,

причому безрозмірні величини  і  визначаються виразами:

R 
  ;   .
2 mgl J   0

Така безрозмірна форма подання рівнянь є зручнішою, бо дозволяє
скоротити кількість параметрів (у нашому випадку замість трьох розмірних

параметрів J , R та mgl залишається один – ), а також подавати розв'язок

рівняння у більш узагальненій формі.
Винесення обчислення моментів зовнішніх дій в окрему обчислювальну

процедуру також дозволяє зробити процедуру правих частин рівняння

маятника більш загальною, якщо звернення до процедури обчислення
моментів здійснювати також через функцію feval.

Утворимо процедуру MomFml, яка обчислюватиме моменти сил, що
діють на маятник:

function m = MomFM1(t, y);
% Обчислення Моментів сил, що діють на Фізичний Маятник.
% Здійснює розрахунок моменту “m” сил за формулою:
% m =-2*dz*y(2) - (nmx*sin(nu*t+ex)*cos(y(l)) +•••
% + nmy*sin(nu*t+ey)*sin(y(l)),
% Коефіцієнти передаються у процедуру через глобальний вектор
% KM l=[dz, nmy, nmx, nu, ey, ex]
global KM1
m = -2*KM l(l)*y(2)- (KM l(3)*sin(KM l(4)*t+KM l(6))*cos(y(l)) +...
KM l(2)*sin(KMl(4)*t+KM l(5))*sin(y(l)));
% Кінець процедури MomFM1
Відтепер слід перебудувати процедуру правих частин. Назвемо цей