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


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


a = zeros (16,16);  b =

zeros (16,1); X = zeros

(16,1); la = sqrt (0.1*0,315365);

l1 = 4.0;  l2 =

6.0;  l3 = 3.0;  l4 = 1.0; m = 20.0; q = 10.0; f1 = - 40.0;

f2 = 60.0; am = 2.0; af1 = 2.0; af2 = 1.0;

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

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

a(2,2) = (cosh(la*l1)+cos(la*l1))/2;

a(2,4) = - (cosh(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);

b(1,1) = - m* (cosh(la*(l1-am))- cos (la*(l1-am))) / (2*la^2);

b(2,1) = - m* (sinh(la*(l1-am)) + sin (la*(l1-am))) / (2*la);

b(3,1) =  m* (cosh(la*(l1-am)) + cos (la*(l1-am))) / 2;

b(4,1) =  m* la*(sinh(la*(l1-am)) - sin (la*(l1-am))) / 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);

b(5,1) = - q* (cosh(la*l2) + cos (la*l2) - 2) / (2*la^4);

b(6,1) = - q* (sinh(la*l2) - sin (la*l2)) / (2*la^3);

b(7,1) =  q* (cosh(la*l2) - cos (la*l2)) / (2*la^2);

b(8,1) =  q* (sinh(la*l2) + sin (la*l2)) / (2*la);

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);




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



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