20696

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

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

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

0001; J = [diffy1'x1' diffy1'x2' diffy1'x3' ; diffy2'x1' diffy2'x2' diffy2'x3' ; diffy3'x1' diffy3'x2' diffy3'x3' ]; p=[2;2;2]; x1=p1; x2=p2; x3=p3; dp=[inf;inf;inf]; while maxabsdp1:3 eps dp=[0;0;0]; Fk=[0;0;0]; Jk=evalJ; for i=1:3 Fki=evalFi:; end dp=invJkFk; p=pdp; x1=p1; x2=p2; x3=p3; end p 2.

Украинкский

2013-07-31

18.5 KB

3 чел.

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

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

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

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

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

10 варіант.  

 

Хід роботи

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. Результати виконання програми.

10 вар.    

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

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

Міністерство освіти і науки, молоді та спорту України

Черкаський політехнічний технікум

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

З дисципліни

«Чисельні методи»

Виконав:

студент групи ПС4-1

Марікуца Р.В.

Перевірив:

Русін О. О.

Черкаси 2013


 

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

50738. Финансовый контроль в бюджетных организациях 706 KB
  Цель и задачи работы обосновать значимость финансового контроля в комплексе государственных мероприятий РФ; провести анализ процесса финансового контроля, выявить проблемы, присущие этим процессам и обозначить возможные направления их решения
50739. Знаходження значення інтеграла по формулам Ньютона-Котеса 33.5 KB
  Мета. Навчитися знаходити значення інтеграла по формулам Ньютона-Котеса. Скласти програму. Устаткування: папір формату А4, ПК, С++.
50740. Знаходження інтеграла за формулами прямокутників 33.5 KB
  Мета. Навчитися знаходити значення інтегралу за формулами прямокутників. Скласти програму. Устаткування. папір формату А4, ПК, С++
50741. Знаходження інтегралу за формулами трапецій 31 KB
  Мета. навчитися знаходити значення інтегралу за формулами трапецій. Скласти програму. Устаткування: папір А4, ручка, ПК, програмне забезпечення С++.
50742. Обчислення інтегралу по формулі Сімпсона. Складання алгоритму 54 KB
  Мета. Навчитися обчислювати інтеграл по формулі Сімпсона; склаcти алгоритм. Устаткування: папір формату А4, ПК, С++
50743. Знаходження коренів нелінійного рівняння методом хорд 117.5 KB
  Мета. навчитися відокремлювати корені рівняння графічно та уточнювати методом хорд. Обладнання: лист формату А4, ручка, олівець, лінійка, програмне забезпечення С ++.
50744. Знаходження коренів нелінійного рівняння методом дотичних 133.5 KB
  Встановити корні рівняння графічно та уточнити один з них методом хорд з точністю до 0,001.Відділити корні рівняння аналітично та уточнити один зних методом хорд з точністю до 0,001.
50745. Знаходження коренів нелінійного рівняння комбінованим методом хорд та дотичних 28 KB
  Мета. Навчитися уточнювати корені нелінійного рівняння комбінованим методом хорд та дотичних. Обладнання. Лист формату А4, ручка , олівець, лінійка, програмне забезпечення С ++.
50746. Виготовлення та градуювання штучних термопар 167 KB
  Обладнання: термоелектроди для штучної термопари; джерело постійного струму; тигель з графітом; тигель з піском з електропідігріванням; вольтметр; термометр з робочим діапазоном температур що вимірюється 0 – 300 ˚С штатив; демонстраційні термопари різного типу та вимірювальні прилади для термопар. Для вимірювання вищих температур застосовують термопари які виготовляють з’єднанням між собою кінців різнорідних провідників. Для вимірювання температур до 1100 ˚С використовують термопари з неблагородних металів мідькопель залізокопель...