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


 

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

23040. Операційні підсилювачі (негативний зворотний зв`язок) 68.5 KB
  Вступ Операційний підсилювач – це диференційний підсилювач постійного струму який в ідеалі має нескінченний коефіцієнт підсилення за напругою і нульову вихідну напругу за відсутністю сигналу на вході великий вхідний опір і малий вихідний а також необмежену смугу частот сигналів що підсилюються. Мета роботи ознайомитись із властивостями операційних підсилювачів опанувати способи підсилення електричних сигналів в ОП охопленому негативним зворотним зв`язком та способи виконання математичних операцій за допомогою ОП. Операційні...
23041. Пасивні RC-фільтри 129.5 KB
  Пасивний чотириполюсник не містить у собі джерела енергії; потужність що виділяється в елементі кола підключеного до виходу чотириполюсника менше потужності що споживається від джерела сигналу підключеного до входу чотириполюсника; на виході такого чотириполюсника ніколи не буває гармонік яких би не було у поданому на його вхід сигналі якщо цей чотириполюсника створений на базі лінійних елементів. Функцію перетворення будьякого чотириполюсника можна подати кількома варіантами в залежності від способу впливу...
23042. Напівпровідникові діоди. Вольт-амперна характеристика (ВАХ) 83.5 KB
  Вольтамперна характеристика ВАХ – це залежність величини струму ІД крізь pn перехід діода від величини і полярності напруги UД прикладеної до діода. Виконання роботи передбачає використання осцилографа як характериографа з метою одержання на екрані двоканального осцилографа зображення ВАХ діода а також побудову ВАХ шляхом вимірювання деякої кількості величин струму ІД що відповідають певним величинам та полярності напруги UД і представленням результату у вигляді графіка. Залежність струму крізь діод від прикладеної до...
23043. Транзистори 88 KB
  Вихідна вольтамперна характеристика ВАХ біполярного транзистора – це залежність величини струму колектора ІК від напруги між колектором та емітером UКЕ при певному струмі бази ІБ або напруги між базою та емітером UБЕ . Вихідна вольтамперна характеристика ВАХ польового транзистора – це залежність величини струму стока ІС від напруги між стоком та витоком UСВ при певній напрузі між затвором та витоком UЗВ . Виконання роботи передбачає використання осцилографа як характериографа з метою одержання на екрані двоканального...
23044. ПІДСИЛЮВАЧІ НА ТРАНЗИСТОРАХ 103 KB
  Він є лише керувальним пристроєм а збільшення потужності сигналу відбувається за рахунок зовнішнього джерела напруги струмом в колі якого й керує транзистор. Характер зміни вхідного сигналу повинен передаватися на вихід без помітних спотворень. Кажуть що має місце інверсія фази сигналу. Як випливає з рівняння ЕберсаМола [1] імпеданс для малого сигналу з боку емітера при фіксованій напрузі на базі дорівнює rе = kT еIк 5 де k – стала Больцмана Т – абсолютна температура е – заряд електрона Iк – струм колектора.
23045. Дешифратори та мультиплексори 1.3 MB
  Це здійснюється аналогічно заданню параметрів елементів схеми за допомогою редактора пробних сигналів Stimulus Editor. Це робиться аналогічно заданню мітки вузла схеми причому в описі шини слід перерахувати через кому мітки усіх вузлів що входять у шину Альтернативна можливість полягає у використанні конструкцій типу BUS[1n] де BUS – ім’я шини BUS[1]BUS[n] – відповідні мітки вузлів. Пакет OrCAD дозволяє провести суто цифрове моделювання для даного вузла схеми якщо до цього вузла під’єднані лише цифрові входи та виходи. Зазначимо що...
23046. Тригери 1.45 MB
  1 зображено схему найпростішого RSтригера на елементах 2ІНЕ серії 74 із зворотнім зв’язком. Встановлення тригера в 10 відбувається при подачі нуля на NSNR при цьому протилежний вхід повинен бути встановлений в одиницю. Подача двох нулів є забороненою комбінацією при якій стан тригера буде невизначеним. Для даної схеми тригера доцільно у початковий момент встановити режим зберігання інформації потім у деякий момент подати імпульс встановлення 1 потім імпульс встановлення 0 після цього знову використати режим зберігання і нарешті...
23047. Регістри та лічильники 1.83 MB
  Виведіть графік залежності вхідних Reset Shift Info та вихідних Q0Q2 цифрових сигналів регістра від часу та поясніть ці залежності. Виведіть відповідні графіки для вхідних та вихідних сигналів та поясніть ці залежності. Джерела пробних сигналів підберіть таким чином щоб регістр послідовно виконав операції: а паралельного запису числа 0101; б перетворення цього числа на послідовний код; в послідовного запису числа 1010; Проведіть моделювання для цієї схеми. Виведіть відповідні графіки для вхідних та вихідних сигналів та поясніть ці...
23048. Імпульсні цифрові схеми 2.62 MB
  Формувачі імпульсів. 1 зображено схему формувача імпульсів на логічних елементах ІНЕ. 1 Недоліком цієї схеми є те що для формування імпульсів досить великої тривалості потрібно використати велику кількість логічних елементв. 2 Для формування імпульсів з синусоїдальної напруги часто застосовується тригер Шмітта рис.