Xlabel('Час (в періодах малих власних коливань)','FontSize',12);

ylabel('Кут у градусах','FontSize',12);

subplot(2,4,1:2); plot(y(:,1)*180/pi,y(:,2));grid;

title('Фазовий портрет','FontSize',14);

Xlabel('Кут у градусах','FontSize',12); ylabel('Швидкість','FontSize',12);

subplot(2,4,3:4); axis('off');

h1= text(0.2,1.1,'Рух фізичного маятника','FontSize',14,'FontWeight','Bold');

h1=text(0.2, 1,'у відповідності до рівняння','FontSize',12);

h1=text(0,0.9,'fi" + 2*dz*fi'' + [1+nmy*sin(nu*t+ey)]*sin(fi) =','FontSize',14);

h1=text(0.45,0.8,' = - nmx*sin(nu*t+ex)*cos(fi)','FontSize',14);

h1=text(0,0.7,'за таких значень параметрів:','FontSize',12);

h1=text(0.45,0.6,sprintf('dz = %g',KM1(1)),'FontSize',12);

h1=text(0,0.5,sprintf('nmy = %g',KM1(2)),'FontSize',12);

h1=text(0.7,0.5,sprintf('nmx = %g',KM1(3)),'FontSize',12);

h1=text(0,0.4,sprintf('ey = %g град.',KM1(5)*180/pi),'FontSize',12);

h1=text(0.7,0.4,sprintf('ex = %g град.',KM1(6)*180/pi),'FontSize',12);

h1=text(0.45,0.3,sprintf('nu = %g',KM1(4)),'FontSize',12);

h1=text(0,0.2,'і початкових умов:','FontSize',12);

h1=text(0,0.1,sprintf('fi(0) = %g град.',fi0*180/pi),'FontSize',12);

h1=text(0.7,0.1,sprintf('fi''(0) = %g',fit0),'FontSize',12);

h1=text(0,0.05,'-------------------------------------------------------------------------------');

h1=text(0,-0.2,'-------------------------------------------------------------------------------');

h1=text(-0.05,-0.05,['Програма ',sprogram]);

h1=text(0.5,-0.05,'Автор - Лазарєв Ю.Ф., каф. ПСОН');

h1=text(-0.05,-0.15,['Виконав ',sname]);

tm=fix(clock); Tv=tm(4:5);

h1=text(0.65,-0.15,[sprintf(' %g:',Tv),' ',date]);

% Кінець програми FizmayatnU_Yadro

Файл FM2.m

function z=FM2(t,y);

% Процедура обчислення правих частин ДР маятника

% Лазарєв Ю.Ф. 1-02-2001

global MPFUN

z(1)=y(2);

z(2)=-sin(y(1))+feval(MPFUN,t,y);

z=z';

Файл MomFM1.m

function m=MomFM1(t,y);

% Процедура обчислення поточного значення моменту сил,

% що діють на маятник

% Лазарєв Ю.Ф. 1-02-2001

global KM1

m=-2*KM1(1)*y(2)-KM1(3)*sin(KM1(4)*t+KM1(6))*cos(y(1))-KM1(2)*sin(KM1(4)*t+KM1(5))*sin(y(1));

1.1. Залежність частоти власних коливань маятника від їх амплітуди

1.1.1. Мета роботи

Вивчити залежність частоти власних коливань маятника від їхньої амплітуди. Дослідити цю залежність шляхом експериментування на програмній моделі і порівняти її з аналогічними залежностями, одержаними шляхом теоретичних досліджень.

1.1.2. Сценарій роботи

1. Вивчити теоретичну частину роботи

2. Увійти у програмне оточення лабораторної роботи.

Встановити початкове значення кута відхилення маятника від вертикалі 179,9 градуса, початкову швидкість рівну нулю, амплітуди віброперевантажень рівними нулю і відносний коефіцієнт загасання рівним нулю. Зафіксувати результати моделювання на ЕОМ і, перш за все, величину періоду і форму коливань.

Зменшуючи початкове відхилення маятника, повторити експеримент, фіксуючи кожного разу період коливань, для 20-30 значень амплітуди від 179 до 0 градусів.

3. Побудувати експериментальну скелетну криву. Порівняти її з теоретичними. Зробити висновки.

4. Підготувати відповіді на запитання поділу 1.1.4.

5. Показати викладачеві одержані результати, пояснити їх. Розказати про зроблені висновки. Відповісти на запитання викладача.