11681

Розвязання систем нелінійних рівнянь. Метод Ньютона

Лабораторная работа

Информатика, кибернетика и программирование

Лабораторна робота №4 Тема: Розв’язання систем нелінійних рівнянь. Метод Ньютона. Мета роботи: познайомитися з методами розв’язання систем нелінійних алгебраїчних рівнянь реалізувати заданий за варіантом метод у середовищі МatLAB. Завдання для виконання лаборат

Украинкский

2013-04-10

44.19 KB

27 чел.

Лабораторна робота №4

Тема: Розвязання систем нелінійних рівнянь. Метод Ньютона.

Мета роботи: познайомитися з методами розвязання систем нелінійних алгебраїчних рівнянь, реалізувати заданий за варіантом метод у середовищі МatLAB.

Завдання для виконання лабораторної роботи:

Розв’язати наступні рівняння методом Ньютона в середовищі МatLAB. Провести тестування створеної програми на прикладі, вибраному за варіантом.

2 варіант.  

Теоретичні відомості

Для розвязання нелінійних та трансцендентних рівнянь можуть застосовуватися звичайний ітераційний метод. Але при знаходженні розв‘язків збіжність ітераційного методу до конкретного розв‘язку залежить від початкових значень змінних.

Метод Ньютона оснований на знаходженні послідовності {[x1k,x2k,…,xnk]}, що збігається до розв‘язку (x1, x2, …, xn). Цей метод називають ітерацією нерухомої точки. Величина похідної в нерухомій точці визначає, чи буде ітераційний процес збіжним. Коли це правило застосовується для функції декількох змінних – похідні повинні бути частинними. Узагальненням “похідної” для системи функцій є матриця Якобі (Якобіан). Наприклад, для функцій трьох незалежних змінних f1(x,y,z), f2(x,y,z), f3(x,y,z) матриця Якобі має вигляд:

   (4.1)

Для функцій декількох змінних диференціал використовується, щоб показати, як змінення незалежних змінних вплине на залежні змінні. Наприклад, задані функції:

  (4.2)

Допустимо, що значення цих функцій відомі в точці (х0, у0, z0) і необхідно визначити їх значення в точці (x, y, z) віддаленій на ().

 (4.3)

де – диференціали залежних змінних, – диференціали незалежних змінних. Якщо змінення функції позначити dF, а змінення змінних dX, використовуючи векторне позначення можемо записати:

  (4.4)

Збіжність поблизу нерухомої точки. Ітерацію нерухомої точки визначаємо наступним чином:

    (4.5)

Теорема. Припустимо, що функції (2) та їх перші частинні похідні неперервні в області, в якій знаходиться нерухома точка (x, y, z). Якщо (х0, у0z0) достатньо близько розташована до точки (x, y, z) і виконуються умови:

  (4.6)

то ітерація збігається до нерухомої точки (x, y, z).

Метод Ньютона виконується за наступними етапами:

1етап – для здійснення обчислень сформуємо функцію:

     (4.7)

2 етап – обчислимо Якобіан:

 (4.8)

3 етап – розв‘яжемо систему рівнянь:

4 етап – обчислимо координати наступної точки – наступне наближення до розв‘язку має вигляд:

   (4.9)

Хід роботи

1. Складаю функцію для знаходження коренів

function nNewton(y1,y2,y3);

F=char(y1,y2,y3);

eps=0.0001;

J = [(diff(y1,'x1')) (diff(y1,'x2')) (diff(y1,'x3')) ;

    (diff(y2,'x1')) (diff(y2,'x2')) (diff(y2,'x3')) ;

    (diff(y3,'x1')) (diff(y3,'x2')) (diff(y3,'x3')) ];

p=[2;2;2];

x1=p(1);

x2=p(2);

x3=p(3);

dp=[inf;inf;inf];

while (max(abs(dp(1:3)))>eps)

   dp=[0;0;0];

   Fk=[0;0;0];

Jk=eval(J);

for i=1:3

    Fk(i)=eval(F(i,:));

end

dp=inv(Jk)*Fk;

p=p-dp;

x1=p(1);

    x2=p(2);

    x3=p(3);

end

p

2. Результати виконання програми.

2 вар.    

>> nNewton('(x1^2+2*(x2^2)+3*(x3^2))','(3*x1+x2^3+x3*8)','(5*(x1^2)+8*x2+7*(x3^2))')

p =

 1.0e-004 *

  -0.8573

  -0.0000

   0.3215

  

>> nNewton('(x1+2*(x2^2)+3*(x3^3))','(3*(x1^3)+x2+2*(x3^2))','(x1^2+8*(x2^3)+x3)')

p =

 1.0e-011 *

   0.0021

   0.1446

   0.0000

Висновок: виконавши лабораторну роботу, я розглянув ітераційні методи розв’язання систем нелінійних рівнянь, а саме метод Ньютона. Розв’язана запропоновану систему.


 

А также другие работы, которые могут Вас заинтересовать

40706. Направления развития сферы услуг в современных условиях 35 KB
  Направления развития сферы услуг в современных условиях. Рынок услуг сфера обращения или совокупность актов куплипродажи различных видов деятельности удовлетворяющие потребности отдельного человека группы людей организаций. Рынок нематериальных услуг как совокупность социальноэкономических институтов – это набор субъектов хозяйствования связей отношений норм социальноэкономической жизни которые имеют устойчивый характер и сохраняются в течение длительного времени. На рынке нематериальных услуг проявляются экономические отношения...
40707. Регион как территориальный социально-экономический комплекс. Основные характеристики 28.5 KB
  Она исследует экономические явления и процессы связанные с рыночным развитием хозяйства отдельных регионов и их включением в единое экономическое пространство России. Основные черты : узкая специализация регионов с доминированием какоголибо одного комплекса: рыбопромышленного мясомолочного военнопромышленного и т. Типы проблемных регионов: слаборазвитые хронически отстают от среднероссийского уровня экон.
40708. Устойчивое развитие регионов как фактор стабильности национальной экономики 76 KB
  И в городе и в районе и в субъекте Федерации необходимо осуществлять разработку стратегии развития оценивать положительные и отрицательные стороны этого развития выявлять факторы конкурентоспособности формировать общую и функциональные стратегии вести стратегический контроль проводить мероприятия городского и регионального маркетинга. Эти программы позволяют властным структурам координировать работу всех органов и лиц...
40709. Налоговое бремя экономики и его измерение 31.5 KB
  Проблема тяжести налогового бремени волнует различных агентов хозяйственной деятельности: государство как субъект управления хозяйственной деятельностью на своей территории и перераспределения доходов от нее в виде налогов в пользу прочих элементов государственной и социальной жизнедеятельности; организации и предприятия как объекты управляющего воздействия государства и субъекты собственно предпринимательской деятельности обеспечивающие формирование источника предпринимательского дохода и соответственно налогооблагаемой базы; ...
40710. Налоговая политика: сущность, содержание и механизмы реализации 26.5 KB
  Налоговая политика: сущность содержание и механизмы реализации. Налоговая политика комплекс правовых действий органов власти и управления определяющий целенаправленное применение налоговых законов. Налоговая политика является частью финансовой политики. Экономическая обоснованная налоговая политика преследует цель оптимизировать централизацию средств через налоговую систему.
40711. Инвестиции - источники формирования и объекты вложений (инвестирования) 32 KB
  Инвестиции источники формирования и объекты вложений инвестирования. Инвестиции от лат. Выделяют также производственные инвестиции направляемые на новое строительство реконструкцию расширение и техническое перевооружение действующих предприятий и интеллектуальные вкладываемые в создание интеллектуального духовного продукта; контролирующие прямые инвестиции обеспечивающие владение более чем 50 голосующих акций другой компании и неконтролирующие обеспечивающие владение менее чем 50 голосующих акций другой компании. В зависимости...
40712. Финансовый рынок: структура и механизм функционирования 28.5 KB
  Финансовый рынок – это сфера проявления экономических отношений между продавцами и покупателями финансовых денежных ресурсов и инвестиционных ценностей то есть инструментов образования финансовых ресурсов между их стоимостью и потребительной стоимости. Как и любой рынок финансовый рынок предназначен для установления непосредственных контактов между покупателями и продавцами финансовых ресурсов. Финансовая система Российской Федерации включает следующие звенья финансовых отношений: государственную бюджетную систему; внебюджетные...
40713. Рентабельность. Виды и методика расчета ее уровня 31.5 KB
  Рентабельность. РЕНТАБЕЛЬНОСТЬ от нем. Количественно рентабельность исчисляется как частное от деления прибыли на затраты расход ресурсов обеспечивающих получение прибыли. Рентабельность продукции определяется как отношение прибыли от ее реализации к себестоимости.
40714. Бизнес-планирование 37 KB
  Бизнеспланирование. В рыночной экономике бизнесплан является рабочим инструментом используемым во всех сферах предпринимательства. Основной целью разработки бизнесплана является планирование хозяйственной деятельности фирмы на ближайший и отдаленные периоды в соответствии с потребностями рынка и возможностями получения необходимых ресурсов. Другие цели разработки плана бизнеса могут быть различными например: уяснить степень реальности достижения намеченных результатов; доказать определенному кругу лиц целесообразность...