Студопедия Главная Случайная страница Обратная связь

Разделы: Автомобили Астрономия Биология География Дом и сад Другие языки Другое Информатика История Культура Литература Логика Математика Медицина Металлургия Механика Образование Охрана труда Педагогика Политика Право Психология Религия Риторика Социология Спорт Строительство Технология Туризм Физика Философия Финансы Химия Черчение Экология Экономика Электроника

Часть работы.СЛАУ.





[1] Постулат – утверждение принимаемое без доказательства, как аксиома. О справедливости того или иного постулата можно судить путем сравнения с экспериментом результатов, полученных при использовании того или иного постулата.

Оружие победы 2. Легкие танки

БТ-2
БТ-5
БТ-7  
Танк М3А1 Стюарт (ленд-лиз)  
Танк МКVII Тетрарх (ленд-лиз)  
Танк МК III Валентайн (ленд-лиз)  
Танк Т-26
Танк Т-50    
Танк Т-60
Т-70
САУ СУ-76
 

Метод деления отрезка пополам.

Интервал изоляции [-7;-6].

i a b c F(a) F(c) |b-a|
  -7 -6 -6,5   10,75967  
  -6,5 -6 -6,25 10,759668 0,796755 0,5
  -6,25 -6 -6,125 0,7967554 -3,34813 0,25
  -6,25 -6,125 -6,1875 0,7967554 -1,34018 0,125
  -6,25 -6,1875 -6,21875 0,7967554 -0,28821 0,0625
  -6,25 -6,21875 -6,23438 0,7967554 0,250101 0,03125
  -6,23438 -6,21875 -6,22656 0,2501008 -0,02009 0,015625
  -6,23438 -6,22656 -6,23047 0,2501008 0,114745 0,007813
  -6,23047 -6,22656 -6,22852 0,1147445 0,047261 0,003906
  -6,22852 -6,22656 -6,22754 0,0472615 0,013569 0,001953
  -6,22754 -6,22656 -6,22705 0,0135686 -0,00327 0,000977
c=(-7-6)/2=-6,5

f(a)=0,5-7-1-(-7-2)2-7=39

f(c)=0,5-6,5-1-(-6,5-2)2-6,5=10,75967

f(a)* f(c)>0

Интервал изоляции [1;2].

i a b c F(a) F(c) |b-a|
      1,5 -0,5 0,603553  
    1,5 1,25 -0,5 0,107948 0,5
    1,25 1,125 -0,5 -0,18212 0,25
  1,125 1,25 1,1875 -0,182123 -0,03359 0,125
  1,1875 1,25 1,21875 -0,033593 0,038053 0,0625
  1,1875 1,21875 1,203125 -0,033593 0,002449 0,03125
  1,1875 1,203125 1,195313 -0,033593 -0,01552 0,015625
  1,195313 1,203125 1,199219 -0,015518 -0,00652 0,007813
  1,199219 1,203125 1,201172 -0,006521 -0,00203 0,003906
  1,201172 1,203125 1,202148 -0,002033 0,000209 0,001953
  1,201172 1,202148 1,20166 -0,002033 -0,00091 0,000977
c=(1+2)/2=1,5

f(a)=0,51-1-(2-2)2+1=-0,5

f(c)=0,51,5-1-(1,5-2)2+1,5=0,603553

f(a)* f(c)<0

Интервал изоляции [3;4].

i a b c F(a) F(c) |b-a|
      3,5 1,125 0,338388  
  3,5   3,75 0,338388 -0,23817 0,5
  3,5 3,75 3,625 0,338388 0,065427 0,25
  3,625 3,75 3,6875 0,065427 -0,08254 0,125
  3,625 3,6875 3,65625 0,065427 -0,0076 0,0625
  3,625 3,65625 3,640625 0,065427 0,029154 0,03125
  3,640625 3,65625 3,648438 0,029154 0,010838 0,015625
  3,648438 3,65625 3,652344 0,010838 0,001635 0,007813
  3,652344 3,65625 3,654297 0,001635 -0,00298 0,003906
  3,652344 3,654297 3,65332 0,001635 -0,00067 0,001953
  3,652344 3,65332 3,652832 0,001635 0,000482 0,000977
c=(3+4)/2=3,5

f(a)=0,53-1-(3-2)2+3=1,125

f(c)=0,53,5-1-(3,5-2)2+3,5=0,338388

f(a)* f(c)>0

 

2. Метод Ньютона.

Условие начального приближения

f(x)*f''(x)>0

Интервал изоляции [-7;-6].

k xi Xi+1 –Xi F(xi) f'(xi)
  -7 0,559358   -69,7228
  -6,44064 0,193172 8,176259 -42,3264
  -6,24747 0,020124 0,7077 -35,1676
  -6,22735 0,000202 0,006952 -34,4783
  -6,22715   6,91E-07 -34,4715
x f f''
-7   59,497986
-6 -7 28,748993
 

 

f(-7)=0,5-7-1-(-7-2)2-7=39

f''(x)=0,5x *ln0.5-5*x+5;

f''(x)=0.5x*(ln0.5)2-5;

f'(-7)= 0,5-7 *ln0.5-5*(-7)+5 =-69.7228

f(-7)* f''(-7)= 39*59.497986>0

x1=x0- f(x)/ f'(x)=-7+39/(-69.7228)=-6.44064

|X1-X0|=-0,55936>0,001

 

Интервал изоляции [1;2].

k xi Xi+1 –Xi F(xi) f'(xi)
    0,188436 -0,5 2,653426
  1,188436 0,01355 -0,03142 2,318991
  1,201986 7,16E-05 -0,00016 2,294734
  1,202057   -4,6E-09 2,294605
x f f''
  -0,5 -1,759773
  1,25 -1,879887
 

 

f(1)=-0.5

f''(x)=0,5x *ln0.5-5*x+5;

f''(x)=0.5x*(ln0.5)2-5;

f'(1) =-1.759773

f(1)* f''(1)= -0.5*(-1.759773)>0

x1=x0- f(x)/ f'(x)=1+0.5/2.653426=1.188436

|X1-X0|>0,001

 

 

Интервал изоляции [3;4].

k xi Xi+1 –Xi F(xi) f'(xi)
    -0,30805 -0,9375 -3,04332
  3,691948 -0,0383 -0,09336 -2,43753
  3,653646 -0,00061 -0,00144 -2,36237
  3,653036   -3,6E-07 -2,36117
x f f''
  1,125 -1,939943
  -0,9375 -1,969972
 

 

f(4)=-0.9375

f''(x)=0,5x *ln0.5-5*x+5;

f''(x)=0.5x*(ln0.5)2-5;

f'(4) =-3.04332

f(4)* f''(4)= -0.9375*(-1.939943)>0

x1=x0- f(x)/ f'(x)=4+0.9375/(-3.04332)=3.691948

|X1-X0|>0,001

 

Метод Чебышева.

xk+1 =xk – (f(xk))/f'(xk) - (f(xk)* f''(xk))/ 2(f'(xk))3 ;

Интервал изоляции [-7;-6].

k xk Xk+1 –Xk F(xk) f'(xk) f''(xk)
  -7 0,562781   -69,7228 59,49799
  -6,43722 0,192484 8,031607 -42,1906 39,63389
  -6,24474 0,017682 0,61162 -35,0733 34,43372
  -6,22705 -9,4E-05 -0,0032 -34,4683 33,98989
  -6,22715   4,59E-05 -34,4715 33,99224
x1 =-7 – (39)/(-69.7228) - 39* 59.49799/ 2(-69.7228)3=-6.43722

Интервал изоляции [1;2].

k xk Xk+1 –Xk F(xk) f'(xk) f''(xk)
    0,164886 -0,5 2,653426 -1,75977
  1,164886 0,030778 -0,08653 2,361084 -1,78572
  1,195665 0,005303 -0,01471 2,306053 -1,79024
  1,200968 0,000904 -0,0025 2,296556 -1,79101
  1,201872   -0,00043 2,294937 -1,79114
x1 =1 – (-0.5)/(2,653426) – (-0.5)* (-1,75977)/ 2(2,653426)3=1,164886

Интервал изоляции [3;4].

k xk Xk+1 –Xk F(xk) f'(xk) f''(xk)
    -0,27529 -0,9375 -3,04332 -1,96997
  3,724709 -0,05873 -0,17427 -2,50185 -1,96366
  3,665978 -0,01066 -0,03072 -2,38656 -1,96215
  3,655322 -0,00188 -0,0054 -2,36566 -1,96187
  3,653439   -0,00095 -2,36196 -1,96182
x1 =4 – (-0,9375)/(-3,04332) -(-0,9375) * (-1.96997)/ 2(-3,04332)3=3,724709

Сравнение результатов, полученных по разным методам решения[3;4].

метод значение количество
решения корня итераций
деление отрезка пополам 3,652832  
Метод Чебышева 3,653439  
Ньютона 3,653036  
Более выгоден метод Ньютона т.к. он потребовал меньшее количество итераций.

Часть работы.СЛАУ.

Постановка задачи. Требуется найти решение системы m линейных уравнений, которая в общем виде записывается в виде

36х1-5x2-11х3 -19х4=-9

х1+33х2-11х3-20х4=-8

12+26х3-19х4=-7

11х1+4х2-5х3+21х4=-6

В матричном виде эта система уравнений записывается так:

А х=f, где

36 -5 -11 -19

А = 1 33 -11 -20 - матрица системы,

5 -1 26 -19

11 4 -5 21

-9 х1

b= -8 - вектор правых частей х= х2 - вектор неизвестных

-7 х3

-6 х4

 

Таким образом, задача состоит в том, чтобы при известных коэффициентах матрицы А и элементах вектора f найти такие значения (х12… хm) T, что при подстановке их в систему уравнений они превращаются в тождества.

Необходимым и достаточным условием существования единственного решения СЛАУ является условие det А ≠0. В случае равенства нулю определителя матрица называется вырожденной и при этом СЛАУ либо не имеет решения, либо имеет их бесчисленное множество.

Все методы решения линейных алгебраических задач можно разбить на два класса: прямые(методы Крамера, обратной матрицы и Гаусса) и итерационные(методы простой итерации, Якоби и Гаусса – Зейделя).

Метод обратной матрицы. Если det A ≠0, то существует матрица А-1, обратная к данной. Умножим исходную систему уравнений на обратную, получим: А-1А х= А-1 f.

Известно, что произведение обратной матрицы на исходную дает единичную матрицу Е, следовательно, получаем Е х= А- 1f, х= А -1f

Решение СЛАУ свелось к умножению известной обратной матрицы на вектор правых частей. Т.о., задача решения СЛАУ и задача нахождения обратной матрицы связаны между собой, поэтому часто решение СЛАУ называют задачей обращения матрицы.

Метод простой итерации. Преобразуем исходную систему линейных уравнений А х=f к эквивалентной системе вида:

х=хα+β, где

х – искомый вектор, а α и β – некоторые новые матрица и вектор соответственно. В качестве первого приближения можно взять хi(0)=0. Следующие приближения находим по рекуррентным формулам

х(к+1)= α(к)+β, к=0,1,2,…

Такой итерационный процесс будем называть методом простых итераций(МПИ). Достаточным условием сходимости МПИ к решению системы при любом начальном векторе х(0) является требование || α;||<1, где || α;|| | - норма матрицы α;.

Существует несколько способов построения порождающей матрицы α;, для которой выполняется достаточное условие сходимости.

 

Метод Гаусса- Зейделя

В отличии от метода Якоби, в котором вычисление всех компонент вектора (k+1)-го приближения проводилось однообразно, в методе Гаусса- Зейделя для расчета i-й компоненты следующего приближения используется уже вычисленное на том, т.е. (k+1)-м шаге, новые значения первых i-1 компонент:

X1(k+1)=1/a11(f1-a11(k)-a13x3(k)-…-a1mxm(k))

X2(k+1)=1/a22(f2-a21(k+1)-a23x3(k)-…-a2mxm(k))

Xi(k+1)=1/aij(fi-ai1(k+1)-ai2x2(k+1)-…-ai,i-1xi-1(k+1)- ai,i+1xi+1(k)-…-a2mxm(k))

Xm(k+1)=1/amm(fm-am1(k+1)-am2x2(k+1)-…-am,m-1xm-1(k+1))

Или, в компактном виде:

Xi(k+1)=1/aij(fi- xj(k+1)- xj(k)), i=1,2…,m.

Достаточным условием сходимости этого метода, как и для метода Якоби, является условие диагонального преобладания:

/ aij/> , i=1,2…,m.







Дата добавления: 2015-10-15; просмотров: 511. Нарушение авторских прав; Мы поможем в написании вашей работы!




Расчетные и графические задания Равновесный объем - это объем, определяемый равенством спроса и предложения...


Кардиналистский и ординалистский подходы Кардиналистский (количественный подход) к анализу полезности основан на представлении о возможности измерения различных благ в условных единицах полезности...


Обзор компонентов Multisim Компоненты – это основа любой схемы, это все элементы, из которых она состоит. Multisim оперирует с двумя категориями...


Композиция из абстрактных геометрических фигур Данная композиция состоит из линий, штриховки, абстрактных геометрических форм...

Краткая психологическая характеристика возрастных периодов.Первый критический период развития ребенка — период новорожденности Психоаналитики говорят, что это первая травма, которую переживает ребенок, и она настолько сильна, что вся последую­щая жизнь проходит под знаком этой травмы...

РЕВМАТИЧЕСКИЕ БОЛЕЗНИ Ревматические болезни(или диффузные болезни соединительно ткани(ДБСТ))— это группа заболеваний, характеризующихся первичным системным поражением соединительной ткани в связи с нарушением иммунного гомеостаза...

Решение Постоянные издержки (FC) не зависят от изменения объёма производства, существуют постоянно...

Концептуальные модели труда учителя В отечественной литературе существует несколько подходов к пониманию профессиональной деятельности учителя, которые, дополняя друг друга, расширяют психологическое представление об эффективности профессионального труда учителя...

Конституционно-правовые нормы, их особенности и виды Характеристика отрасли права немыслима без уяснения особенностей составляющих ее норм...

Толкование Конституции Российской Федерации: виды, способы, юридическое значение Толкование права – это специальный вид юридической деятельности по раскрытию смыслового содержания правовых норм, необходимый в процессе как законотворчества, так и реализации права...

Studopedia.info - Студопедия - 2014-2025 год . (0.014 сек.) русская версия | украинская версия