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


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


А12

 

 

-1

 

12

 

 

 

 

-1

 

 

 

 

-l4А13

 l4А14

А11

 

 

 

 

13

 

 

 

 

 

 

 

 

-1

 

 

 

 

А12

-А13

-А14

14

 

 

 

 

 

 

 

 

 

 

 

 

-1

А11

-А12

-А13

15

 

 

 

 

 

 

 

 

 

 

 

 

 

 -l4А14

А11

А12

16

 

 

 

 

 

 

 

 

 

 

 

 

 

-l4А13

l4А14

А11

Наиболее просто корни могут быть найдены методом перебора в сочетании с процедурой вычисления определителя матрица А*(w). Организуется цикл вычисления определителя çА*(w) ç, значения которого вместе с частотой w

выводятся на печать в формате long e. При просмотре таблицы их значений определяются точки, где определитель изменяет знак. Эти точки и будут являться частотами собственных колебаний. Иногда бывают случаи, когда график функции d = çА*(w) çкасается оси w. Для контроля этих точек полезно построить график у

= d(w), просмотр которого позволяет определить интервалы, где находятся корни и существуют ли точки касания графика оси w. Далее корни могут быть уточнены повторными прогонами программы (без построения графика у = d(w). Рекомендуется начинать поиск корней с начального значения частоты

, с шагом
. Число вычислений п = 300 – 500 позволяет надежно определить первую частоту, старшие частоты определяются при расширении интервалов вычисления определителя çА*(w)ç. При решении данной задачи принимается, что l=EI=m=1.  

Обозначения переменных, принятых в программе

а - матрица А*(w); п – порядок матрицы А*(w); п1




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



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