MATLAB в инженерных и научных расчетах


Продолжение таблицы 2.3 - часть 22


– число циклов вычисления определителя çА*(w)ç; аm

– частота собственных колебаний; dam

– шаг изменения частоты w; m – счетчик цикла вычисления определителя; l1,  l2,  l3,  l4 – длины стержней неразрезной балки; – коэффициент l фундаментальных функций; Х – вектор значений частоты w; Y – вектор значений определителя çА*(w)ç.

Текст программы

l1 = 4.0;  l2 = 6.0;  l3 = 3.0;  l4 = 1.0;  n = 16;  n1 = 300;  am = 0,01;

dam = 0.01;  X = zeros (n1,1);  Y = zeros (n1,1); 

for  m = 1 : n1

la = sqrt (am);  a

= zeros (n,n); 

a(1,2) = (sinh(la*l1) + sin (la*l1))/(2*la);

a(1,4) = - (si h(la*l1) - sin (la*l1))/(2*la^3);

a(2,2)=(cosh(la*l1)+cos(la*l1))/2; a(2,4) = - (cos h(la*l1) - cos (la*l1))/(2*la^2);

a(2,6) = - 1;  a(3,2) = la^ 4*a(1,4); a(3,4) = a(1,2); a(3,7) = - 1; a(4,1) = - 1;

a(4,2) = la^ 4*a(2,4); a(4,4) = a(2,2);

a(5,6) = (sinh(la*l2) + sin (la*l2))/(2*la);

a(5,7)= - (cosh(la*l2) -cos(la*l2))/(2*la^2);

a(5,8) = - (sinh(la*l2) - sin (la*l2))/(2*la^ 3);

a(6,6)=(cosh(la*l2)+cos(la*l2))/2; a(6,7) = - a(5,6); a(6,8) = a(5,7);

a(6,10) = - 1;  a(7,6) = la^ 4*a(5,8); a(7,7) = a(6,6); a(7,8) = a(5,6); a(7,11) = - 1;

a(8,3) = - 1;  a(8,6) = la^ 4*a(5,7); a(8,7) = - a(7,6); a(8,8) = a(6,6);

a(9,10) = (sinh(la*l3) + sin (la*l3))/(2*la);

a(9,11)= - (cosh(la*l3) -cos(la*l3))/(2*la^2);

a(9,12) = - (sinh(la*l3) - sin (la*l3))/(2*la^ 3);

a(10,10)=(cosh(la*l3)+cos(la*l3))/2; a(10,11) = - a(9,10); a(10,12) = a(9,11);

a(10,14) = - 1;  a(11,10) = la^ 4*a(9,12); a(11,11) = a(10,10); a(11,12) = a(9,10);

a(11,15) = - 1; a(12,5) = - 1;   a(12,10) = la^ 4*a(9,11); a(12,11) = - a(11,10);

a(12,12) = a(10,10); a(13,9) = - 1; a(13,14) = (sin h(la*l4) + sin (la*l4))/(2*la);

a(13,15)= - (cosh(la*l4) -cos(la*l4))/(2*la^2);




Начало  Назад  Вперед



Книжный магазин