Материал: Объектно-ориентированное программирование на примере численных методов

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

Ch: = ReadKey - возвращает считанный символ.

И, наконец, последняя, стандартная функция, используемая только в модуле myiter1. Это функция Abs() - возвращающая абсолютное значение аргумента. Она входит в состав стандартного модуля System. Возвращаемый тип данных тот же, что и у аргумента (real, integer).

Abs(x) - возвращает абсолютное значение х.


4. Описание нестандартных функций

Описание не стандартных функций сводится к описанию модулей.

Модуль hord1 содержит следующие функции и методы:(x:real):real; - принимает одно входное значение x, после чего функция вычисляет и возвращает значение исходной математической функции от данного x;.init данный метод инициализирует новый объект типа hord;.shag данный метод выполняет один математический шах по методу хорд в объекте hord;.poisk:real данная функция реализует цикл поиска корня линейного уравнения по методу хорд, после окончания поиска функция возвращает значение корня.

Модуль zeid1 - содержит одну процедуру (zeid). Которая выполняет решение системы линейных уравнений методом Зейделя.

Модуль roonge1 - содержит одну процедуру (roonge). Которая выполняет решение дифференциального уравнения методом Рунге-Кутта.


5. Описание интерфейса

Основная программа GLAV (использующая методы структурного программирования) работает следующим образом. Используя способ запроса, определяет дальнейший ход развития. При получении любого результата отличного от 1,2,3 вновь возвращается на начало программы. Тем самым, исключая возможность ошибочного ввода. При получении ответа соответствующего цифрам 1,2,3 передает управление одной из процедур описанных в не стандартных модулях пользователя. При этом выполнение главной программы практически заканчивается за исключением оператора выхода, который выполняется при вводе цифры 4. Управление передается соответственно одному из модулей (hord1, zeid1, roonge1). Каждый из перечисленных модулей по сути своей представляет отдельную программу, являющуюся составной частью другой. Это позволяет в зависимости от выбора пользователя выполнить тот или иной самостоятельный модуль, входящий в главную программу.

При выборе 1 управление передается модулю myzend (процедура zend) выполнение которого приводит к выходу из модуля в главную программу.

Аналогично построены и два других модуля входящих в программу GLAV.

Окно главной программы




Результат процедуры hord


Результат процедуы zeid




6. Численные примеры

Все выше перечисленное в данной курсовой работе проиллюстрировано на конкретных примерах.

1       Задано алгебраическое уравнение, 2х-3sin(2x) - 1=0, имеющее явно отрицательное значение при , и положительное значение при . При заданной точности вычисления , при решении методом хорд, дает .


a

b

f(a)

f(b)

x

f(x)

1

2

-1

5,270407

1

-1,7278923

1,463504

2

1,288182

5,270407

1,4635038

1,2881816

1,289957

2

-0,01791

5,270407

1,2899565

-0,0179125

1,292362

2

-0,00087

5,270407

1,2923616

-0,0008707

1,292478

2

-4,2E-05

5,270407

1,2924785

-4,153E-05

1,292484

2

-2E-06

5,270407

1,292484

-1,979E-06

1,292484

2

-9,4E-08

5,270407

1,2924843

-9,427E-08

1,292484

2

-4,5E-09

5,270407

1,2924843

-4,491E-09

1,292484

2

-2,1E-10

5,270407

1,2924843

-2,14E-10

1,292484

2

-1E-11

5,270407

1,2924843

1,292484

2

-4,9E-13

5,270407

1,2924843

-4,863E-13

1,292484

2

-2,3E-14

5,270407

1,2924843

-2,309E-14

1,292484

2

-1,1E-15

5,270407

1,2924843

-1,11E-15



Видно из графика, что корень уравнения приблизительно равен 1,3

2       Решение заданной системы

линейных уравнений осуществлялся с помощью Maple:




Заключение

В заключение хотелось бы отметить, что методы структурного программирования позволяют более мобильно и качественно проводить работу по модернизации программного обеспечения. Так достаточно усовершенствовать отдельный модуль откомпилировать его и главный модуль, для достижения, желаемого результата. Более того, текст главной программы более читаем, и занимает меньше места. Имеющиеся наработки в отдельных модулях можно использовать в других программах, не задумываясь об их работоспособности. Четко вычерченная структура позволяет яснее представить цели и задачи стоящие перед программистом. Не нужно все время перемещаться по тексту взад и вперед. В любой промежуток времени работы над очень большой программной появляется возможность сконцентрироваться на отдельном модуле. Открывается возможность совместной работы нескольких человек над одной громоздкой программной. Удачно разработанный модуль одного программиста становиться, доступен другим.

Все выше перечисленное позволяет сделать вывод. Что структурное программирование открывает не объятые просторы в развитии современной техники, и будущее, а вернее сказать уже настоящее именно за ним. Вместе с объектно-ориентированным программированием оно уже составляет авангард современной науки.


Список используемой литературы

1. Воробьев Г.Н., Бахвалов Н.С. «Численные методы». М.: Наука, 2006. 231 с.

. Ефимов А.В., Демидович Б.П. «Линейная алгебра и основы математического анализа». М.: Наука, 2011. 386 с.

. Бараненков Г.С., Демидович Б.П. «Задачи и упражнения по математическому анализу для ВТУЗОВ». М.: Наука, 2010. 184 с.

. Абрамов С.А., Зима Е.В. «Начало программирования на языке Паскаль». М.: Наука, 2007. 8 с.

. Епанешников А.Е., Красильников Ю.И. «Программирование в среде турбо Паскаль». М.: Центр МИФИ СП Диалог, 2008. 3-6 с.



Приложение

Метод хорд:myhord;

interfacecrt;hord=object, b, c, cp, eps:real;init;shag;poisk:real;;funct (x:real):real;funct (x:real):real;:=2*x-3*sin (2*x) - 1;hord.init;;(' Metod HORD');('Vvedite pogreshnost eps = ');(eps);('Vvedite nachalo otrezka: a = ');(a);('Vvedite konets otrezka: b = ');(b);:=b;:=b - (funct(b)*(a-b))/(funct(a) - funct(b));;hord.shag;(funct(a)*funct(c)>0)a:=cb:=c;:=c;:=b - (funct(b)*(a-b))/(funct(a) - funct(b));;hord.poisk:real;(abs (c-cp)>eps) do(funct(c)=0) then;;;;:=c;;.

Метод Зейделя:Dec_sys;

crt;=array [1.. 10,1..10] of real;=array [1..10] of real;_sys=object:integer;Zeidel (a:matrix; b:vector; x:vector; e:real);Dec_Zeidel;Init;Done;;

TDec_sys. Zeidel (a:matrix; b:vector; x:vector; e:real);i, j, flag:integer;, s2, s, v, m:real;

:=1;i:=1 to n do:= 0;j:=1 to n doj<>i then:=s+abs (a[i] [j]);s>=abs (a[i] [i]) then:=0;;

if (flag=1) then:=0;i:=1 to n do:=0;:=0;

j:=1 to i do:=s1+(a[i] [j]*x[j]);j:=i+1 to n do:=s2+a[i] [j] * x[j];

:=x[i];[i]:=x[i] - ((1 / a[i] [i])*(s1 + s2 - b[i]));

(abs (v-x[i]) >m) then:=abs (v - x[i]);;; until (m>=e);

('ђҐиҐ­ЁҐ бЁб⥬л:');i:=1 to n do('x', i, '=', x[i]:5:5);

writeln ('‘Ёб⥬ ­Ґ б室Ёвбп!');;TDec_sys. Dec_Zeidel;eps:real;:matrix;, x:vector;, j:integer;;('ђҐиҐ­ЁҐ бЁб⥬л га ў­Ґ­Ё©:');('a11ъx1+a12ъx2+a13ъx3=b1');('a21ъx1+a22ъx2+a23ъx3=b1');('a31ъx1+a32ъx2+a33ъx3=b1');

;('‚ўҐ¤ЁвҐ в®з­®бвм: '); readln(eps);

;('‚ўҐ¤ЁвҐ н «Ґ¬Ґ­вл а биЁаҐ­­®© ¬ ваЁжл бЁб⥬л:');i:=1 to n doj:=1 to n do('a', i, j, '=');(a[i] [j]);;;

;('‚ўҐ¤ЁвҐ н «Ґ¬Ґ­вл бв®«Ўж бў®Ў®¤­ле з «Ґ­®ў:');i:=1 to n do('b', i, '=');(b[i]);;

i:=1 to n do[i]:=0;;;

(a, b, x, eps);

;

;

TDec_sys. Init;:=3;;

TDec_sys. Done;;;

.

Метод Рунге-Кутта:rk;crt;TEiler=object, y, h, b, dx1, dx2, dx3, dx4, ddx:real;init;shag;poisk;;proizv (f, g:real):real;proizv (f, g:real):real;:=g*g-f*f;;TEiler.init;;:=0;:=0.5;:=0.1;:=1;;TEiler.shag;:=proizv (x, y);:=proizv (x+h/2, y+h*dx1/2);:=proizv (x+h/2, y+h*dx2/2);:=proizv (x+h, y+h*dx3);:=(dx1+2*dx2+2*dx3+dx4)/6;:=y+ddx*h;:=x+h;;TEiler.poisk;f:text;(f, '1.txt');(f);(x<=b) do('x=', x:3:3,' ', 'y=', y:3:3);(f, x, ' ', y);;;(f);;

end.