Материал: Спецглавы высшей математики. численные методы. Пантелеев И.Н

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

>>рlоt(Т,Z); % визуализация решения уравнения

%Вольтерра (рис. X.6)

Рис. X.6. Численное решение интегрального уравнения

t

 

 

 

1

, полученное на пространственно-временной

(t 2

s2

1)x(s)ds t 2

 

 

t

 

1

 

 

 

 

 

сетке, представленной на рис X.5.

При использовании квадратурных формул замкнутого типа и совпадающей системы узлов возникает проблема вычисления значения x1=x(t1)=x(s1)=x(a), которая не возникала при решении уравнений Вольтерра второго рода. Действительно, при i=1 равенство (X.25) теряет всякий смысл.

Для нахождения начального значения x1 продифференцируем уравнение (X.5) по t.

t

 

 

 

 

(X.31)

Q(t, t)x(t) Q (t, s)x(s)ds

f (t).

a

 

 

Положив в (X.31) t= а, получим

 

 

 

 

(X.32)

Q(a, a)x(a) f (a).

 

Откуда следует, что

x(a) f (a) Q(a, a)

следовательно,

x1 fQ(a)

11

,

. (X.33)

При использовании квадратурной формулы трапеций далее получаем:

 

 

 

 

 

 

 

h

 

 

 

 

 

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f2

 

h

 

Q21 x1

 

 

 

 

 

 

 

 

 

 

 

Q

21

x

1

 

Q

22

x

2

f

2

x

2

 

 

2

 

 

,

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q22

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

h

 

 

 

 

 

 

h

 

 

 

 

 

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f3

h

 

Q31 x1

 

hQ32 x2

 

 

Q

31

x

1

 

Q

32

x

2

 

Q

33

x

3

f

3

x

3

 

 

2

 

 

,

2

 

 

 

 

 

 

h

 

 

 

 

 

2

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q33

 

 

т. е. в общем случае при любому j= 2, 3,

 

 

 

 

 

 

 

 

 

2

 

 

 

 

j 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f j

Q j1 x1 h Q jk xk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x(s j

) x j

 

 

 

 

 

2

 

 

 

 

 

k 2

 

 

 

 

 

.

(X.34)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h

Q jj

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример X.4. Найти в пакете МАТ1АВ решение интегрального уравнения

t

1

 

(t 2 s2 1)x(s)ds t 2

,

 

t

1

 

 

используя замкнутые квадратурные формулы.

1. Создайте файл F11_4.m (листинг X.11), содержащий описание функции, возвращающей значения функции f(t). Листинг X.11. Файл F11_4.m

function z=F11_4(t) z=t.^2-1./t;

248

247

2.Создайте файл dF11_4.m (листинг X.12), содержащий описание функции, возвращающей значения производной функции f(t).

Листинг X.12. Файл dF11_4.m function z=dF11_4(t) z=2*t+1/t.^2;

3.Создайте файл Q11_4.m (листинг X.13), содержащий описание функции, возвращающей значения функции Q(t,s). Листинг X.13. Файл Q11_4.m

function z=Q11_4(t) z=t.^2+s.^2+1;

4. Создайте

файл

Volterra11_2.m

(листинг X.14),

содержащий

описание

функции, возвращающей решения

интегрального уравнения Вольтерра на совпадающей сетке.

Листинг X.14. Файл Volterra11_2.m function [T,Z]=Volterra11(t1,t2,N)

%задание временной сетки h=(t2-t1)/(N-1);

i=1:N; t(i)=t1+h*(i-1); s=t;

%вычисление значений ядра интегрального уравнения

%в узлах временной сетки

fоr i=1:N fоr j=1:N

q(i,j)=Q11_4(t(i),s(j)); end;

end;

%вычисление значений решения интегрального

%уравнения в соответствие с (X.33), (X.34) F=F11_4(t);

F(1)=dF11_4(t(1)); х(1)=F(1)/(h*q(1,1)); for m=2:N

s=0;

fоr j=2:m-1 s=s+q(m,j)*m(j);

еnd;

х(m) = (F(m) -h*s)/( h*q(m,m));

end; Т=t; Z=х;

5. Выполните следующую последовательность команд:

>>t1=1; % левая граница интервала поиска решения

>>t2=1.3;%правая граница интервала поиска решения

>>N=300; % число интервалов разбиения отрезка

%[t1,t2]

>>[Т 2]=Vо1tеrrа11_2 (t1,t2,N) ; % численного

% решения интегрального уравнения Вольтерра

>> рlо1; (Т,Z); % визуализация решения уравнения % Вольтерра (рис. X.7)

Рис. X.7. Численное решение интегрального уравнения

t

 

 

 

1

, полученное при использовании замкнутых

(t 2

s2

1)x(s)ds t 2

 

 

t

 

1

 

 

 

 

 

квадратурных формул

249

250

 

ЗАКЛЮЧЕНИЕ

Система MatLab представляет собой уникальный сплав универсальных программных и алгоритмических средств с широкой гаммой специализированных приложений. Входной язык и среда программирования MatLab очень близки к современным системам визуального программирования на базе универсальных алгоритмических языков типа C++, Java, Objekt Pascal. По ряду аспектов MatLab уступает указанным системам – режим интерпретации, небольшой запас визуальных компонентов. Однако с его библиотекой численных методов ни по объему, ни по качеству не может сравниться ни одна из систем программирования. На базе MatLab созданы многочисленные расширения, обеспечивающие моделирование и анализ систем в разнообразных сферах человеческой деятельности.

Многие учебные заведения у нас и за рубежом используют MatLab при подготовке специалистов широкого профиля. Для современного инженера и научно-технического работника MatLab является незаменимым инструментом моделирования и исследования различных прикладных систем, прежде всего за счет использования готовых решений. Но не менее важно научиться создавать новые приложения, используя программные и алгоритмические средства MatLab, а также возможность объединения модулей, разработанных в разных системах программирования. Использование системы MatLab в учебном процессе позволяет сблизить дисциплины основанные на математическом моделировании и численных методах.

251

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1.Амосов А.А. Вычислительные методы для инженеров / А.А. Амосов, Ю.А. Дубинский, Н.А. Копчёнова. – М.: Высш.

шк., 1994.

2.Бахвалов Н.С. Численные методы / Н.С. Бахвалов, Н.П. Жидков, Г.М. Кобельков. – М.: Наука, 1987.

3.Боглаев Ю.П. Вычислительная математика и программирование / Ю.П. Боглаев. – М.: Высш. шк., 1990.

4. Васильев Ф.П. Численные методы решения экстремальных задач / Ф.П. Васильев. – М.: Наука, 1998.

5.Волков Е.А. Численные методы / Е.А. Волков. – М.:

Наука, 1982.

6.Демидович Б.П. Основы вычислительной математики / Б.П. Демидович, И.А. Марон. – М.: Наука, 1966.

7.Калиткин Н.Н. Численные методы / Н.Н. Калиткин. -

М.: Наука, 1974.

8. Кетков Ю.Л. MATLAB 6.x: программирование численных методов / Ю.Л. Кетков, А.Ю. Кетков, М. Шульц. CПб.: БХВ-Петербург, 2004.

9.Марчук Г.И. Методы вычислительной математики / Г.И. Марчук. - М.: Наука, 1977.

10.Пирумов У.Г. Численные методы / У.Г. Пирумов. – М.:

Дрофа, 2003.

11.Поршнев С.В. Вычислительная математика / С.В. Поршнев. - Курс лекций. – СПб.: БХВ-Петербург, 2004.

12.Самарский А.А. Введение в численные методы / А.А. Самарский. – М.: Наука, 1982.

13.Самарский А.А. Численные методы / А.А. Самарский, А.В. Гулин. – М.: Наука, 1989.

14.Тихонов А.Н. Вводные лекции по прикладной математике / А.Н. Тихонов, Д.П. Костромаров. М.: Наука, 1984.

15.Турчак Л.И. Основы численных методов / Л.И. Турчак, П.В. Плотников. – М.: ФИЗМАТЛИТ, 2005.

252

ОГЛАВЛЕНИЕ

 

4.4. Метод Ньютона ………………………………….

93

 

 

4.5. Решение систем нелинейных уравнений

103

Введение …………………………………………………..

3

средствами пакета MATLAB …………………..

I. Теория погрешностей ………………………………….

5

V. Интерполирование функций …………………………

108

1.1.Общие сведения об источниках погрешностей,

 

5.1. Постановка задачи ………………………………

108

их классификация……………………………………

5

5.2. Интерполяционный полином Лагранжа ……….

113

1.2. Виды погрешностей ……………………………..

5

5.3. Интерполяционный полином Ньютона

116

1.3. Абсолютнаяи относительная погрешности,

 

для равноотстоящих узлов ………………………

формы записи данных ……………………………

7

5.3.1 . Конечныеразности……………………………..

116

1.4. Вычислительная погрешность …………………..

9

5.3.2. ПерваяинтерполяционнаяформулаНьютона

116

1.5. Понятия о погрешности машинных вычислений

11

5.3.3. ВтораяинтерполяционнаяформулаНьютона

121

II. Решение уравнений с одной переменной …………..

16

5.4. Погрешность интерполяции …………………….

124

2.1. Общие сведенья и основные определения ……

16

5.5. Сплайн-интерполяция …………………………… 125

2.2. Отделение корней ……………………………….

17

VI. Численное дифференцирование и интегрирование … 134

2.3. Метод половинного деления ……………………

17

6.1. Численное дифференцирование функций,

134

2.4. Метод простой итерации ………………………..

22

заданных аналитически …………………………

2.5. Преобразование уравнения к

 

6.2. Особенности задачи численного

138

итерационному виду ……………………………

26

дифференцирования функций, заданных таблицей

III. Методы решения систем линейных алгебраических

 

6.3. Интегрирование функций, заданных

138

уравнений ……………………………………………..

37

аналитически …………………………………

3.1. Общие сведения и основные определения …………

6.4. Погрешность численного интегрирования……….

147

3.2. МетодГауссаиегореализациявпакетеMATLAB

37

6.5. Вычисление интегралов методом Монте-Карло.. 150

38

3.3. Вычисление определителей …………………….

42

VII. Методы обработки экспериментальных данных ….

154

3.4. Решение систем линейных уравнений методом

 

7.1. Метод наименьших квадратов ………………….

154

простой итерации ……………………………….

46

7.2. Нахождение приближающей функции в виде

158

3.5. Метод Зейделя …………………………………

53

линейной функции и квадратичного трехчлена..

3.6. Решение систем линейных уравнений

 

7.3. Аппроксимация функцией произвольного вида..

167

средствами пакета MATLAB …………………..

58

VIII. Численные методы решения

171

IV. Методы решения систем нелинейных уравнений ….

60

обыкновенных дифференциальных уравнений ….

4.1. Векторная запись нелинейных систем. Метод

 

8.1. Общие сведения и определения ………………..

171

простых итераций ……………………………….

60

8.2. Метод Пикара ……………………………………

174

4.2. Метод Ньютона решения систем нелинейных

 

8.3. Метод Эйлера ……………………………………

176

уравнений ………………………………………..

64

8.4. Метод Рунге-Кутта ………………………………

183

4.3. Решение нелинейных систем методами спуска … 69

 

 

253

254

 

8.5.Средства пакета MATLAB для решения обыкновенных дифференциальных уравнений … 190

IX. Численные методы решения дифференциальных уравнений в частных производных …………………. 194

9.1.Общие сведения и классификация уравнений в частных производных …………………………… 194

9.2.Численные методы решения эллиптических уравнений ………………………………………… 197

9.3.Явные разностные схемы для уравнений параболического и эллиптического типов ……… 206

9.4.Неявная разностная схема для уравнения параболического типа …………………………… 213

9.5.Решение уравнений с частными производными

методом Монте-Карло …………………………... 219 X.Численные методы решения интегральных уравнений.. 227

10.1.Общие сведения об интегральных уравнениях .. 227

10.2.Квадратурный метод решения интегральных уравнений Фредгольма ………………………… 233

10.3.Квадратурный метод решения интегральных уравнений Вольтера ……………………………. 240

Заключение ……………………………………………… 251

Библиографический список …………………………. 252

255

Источник: https://studfile.net/preview/16565850/