Аппроксимация табулированных функций

Аппроксимация табулированных функций

Реферат

Цель- решить задачу аппроксимации табулированной функции.

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

Предмет исследования: процесс аппроксимации табулированной функции методом наименьших квадратов. Основными методами исследования были выбраны общенаучные методы анализа и синтеза.

В соответствии с заданной целью, объектом и предметом исследования поставлены следующие задачи: необходимо описать метод аппроксимации табулированных функций методом наименьших квадратов, разработать алгоритм нахождения коэффициентов аппроксимирующего полинома и решить задачу аппроксимации в среде Delphi и MathCad

Метод наименьших квадратов применяется также для приближённого представления заданной функции другими (более простыми) функциями и часто оказывается полезным при обработке наблюдений. Разница в методе наименьших квадратов служит сумма квадратов отклонений действительных (экспериментальных) значений от теоретических. Выбираются такие значения параметров модели, при которых сумма квадратов разностей будет наименьшей - отсюда название метода.

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

Введение

Все существующие методы аппроксимации и интерполяции значений функции, по сути, являются различными комбинациями или модификациями этих трех различных логических подходов:

1.       Разложение функции в ряд Тейлора в окрестности некоторой точки;

2.       Интерполяция заданных значений функции;

.        Аппроксимация заданной функции (или ее дискретных значений).

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

Цель данного курсового исследования: решить задачу аппроксимации табулированной функции.

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

Предмет исследования: процесс аппроксимации табулированной функции методом наименьших квадратов.

В соответствии с заданной целью, объектом и предметом исследования поставлены следующие задачи: необходимо описать метод аппроксимации табулированных функций методом наименьших квадратов, разработать алгоритм нахождения коэффициентов аппроксимирующего полинома и решить задачу аппроксимации в среде Delphi и MathCad [2].

В данной курсовой работе надо найти коэффициенты трех полиномов аппроксимирующей функции, используя метод наименьших квадратов.

Метод наименьших квадратов применяется также для приближённого представления заданной функции другими (более простыми) функциями и часто оказывается полезным при обработке наблюдений.

Разница в методе наименьших квадратов служит сумма квадратов отклонений действительных (экспериментальных) значений от теоретических. Выбираются такие значения параметров модели, при которых сумма квадратов разностей будет наименьшей - отсюда название метода.

1.Теоретическая часть

.1 Обзор программных средств

В настоящее время для научно-технических расчетов на компьютерах все чаще и чаще используются не традиционные языки программирования и не электронные таблицы, а специальные математические программы типа Mathcad.

Математические пакеты Mathcad - позволяют специалистам в конкретной научно-технической области очень быстро освоить работу на компьютере и реализовать на них математические модели, не вдаваясь в тонкости программирования на традиционных языках.

Далее перечислим преимущества работы в среде математической программы Mathcad:

-       математические выражения в среде Mathcad записываются в их общепринятой нотации: числитель находится сверху, а знаменатель - внизу; в интеграле пределы интегрирования также расположены на своих привычных местах. Казалось бы, это все мелочи, никак не влияющие на вычислительный процесс. Программа должна быть понятной не только для компьютера, но и для человека. Пользователь, читая распечатку принтера или глядя на дисплей, видит, что данная величина записана в числителе и ее рост приводит к возрастанию всего выражения. А это очень важно при анализе математических моделей, форма и содержание которых едины [2].

-       в среде Mathcad процесс создания «программы» идет параллельно с ее отладкой. Пользователь, введя новое выражение, может не только сразу подсчитать, чему оно равно при определенных значениях переменных, но и построить график или поверхность, беглый взгляд на которые может безошибочно показать, где кроется ошибка, если она была допущена при вводе формул или при создании самой математической модели. «Отладочные» фрагменты можно оставить в готовом документе для того, чтобы, например, еще раз убедить воображаемого или реального оппонента в правильности модели;

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

1.       Решение алгебраических уравнений и систем (линейных и нелинейных);

2.       Решение обыкновенных дифференциальных уравнений и систем (задача Коши и краевая задача);

.        Решение дифференциальных уравнений в частных производных;

.        Статистическая обработка данных (интерполяция, экстраполяция, аппроксимация и многое другое);

.        Работа с векторами и матрицами (линейная алгебра и др.);

.        Поиск минимумов и максимумов функциональных зависимостей [3].

-       пакет Mathcad дополнен справочником по основным математическим и физико-химическим формулам и константам, которые можно автоматически переносить в документ без опасения внести в них искажения, увы, нередкие при ручной работе;

-       к пакету Mathcad можно приобрести те или иные электронные учебники по различным дисциплинам: решение обыкновенных дифференциальных уравнений, статистика, термодинамика, теория управления, сопротивление материалов и т.д. Прежде чем решать возникшую проблему, пользователь может изучить электронный учебник и перенести из него в свой документ нужные фрагменты, отдельные формулы и константы;

-       решая поставленную задачу, пользователь может вводить не только числовые значения переменных, но и дополнить их размерностями. При этом пользователь вправе выбирать и систему единиц (СИ, кг-м-с, г-см-с, британская), и конкретные размерности (мм, дюймы, футы и т.д.): система Mathcad в них сама разберется и выдаст ответ с заданной пользователем размерностью;

-       в систему Mathcad интегрированы средства символьной математики, что позволяет решать поставленные задачи (этап задачи) не только численно, но и аналитически;

-       не выходя из среды Mathcad, возможно открывать новые документы на других серверах и пользоваться теми преимуществами информационных технологий, предоставляемых Internet.

Кроме того, не следует забывать, что пакет Mathcad - это полноценное Windows-приложение. Решая поставленную задачу, можно в статике (через буфер обмена Windows) или в динамике (OLE-технологии) передать данные в среду другой программы (в среду языка fortran, например) и там решить часть задачи.

В Белоруссии Borland Delphi появляется в конце 1993 г. и сразу же завоевывает широкую популярность. Новые версии выходят практически каждый год. В них реализуются все новые мастера, компоненты и технологии программирования.

Действительно, процесс разработки в Delphi предельно упрощен. В первую очередь это относится к созданию интерфейса, на который уходит 80% времени разработки программы. Мы просто помещаем нужные компоненты на поверхность Windows-окна (в Delphi оно называется формой) и настраиваем их свойства с помощью специального инструмента (Object Inspector). С его помощью можно связать события этих компонентов (нажатие на кнопку, выбор мышью элемента в списке и т.д.) с кодом его обработки - и вот простое приложение готово. Причем разработчик получает в свое распоряжение мощные средства отладки (вплоть до пошагового выполнения команд процессора), удобную контекстную справочную систему (в том числе и по Microsoft API), средства коллективной работы над проектом, всего просто не перечислить [1].

Таким образом, мы можем использовать Delphi для создания как самых простых приложений, на разработку которых требуется 2-3 часа, так и серьезных корпоративных проектов, предназначенных для работы десятков и сотен пользователей. Причем для этого можно использовать самые последние веяния в мире компьютерных технологий с минимальными затратами времени и сил.

1.2 Постановка задачи

Задание 2 вариант 25. В данной курсовой работе надо:

1.       Разработать алгоритм нахождения коэффициентов трёх аппроксимирующих полиномов (многочленов) вида  для табулированной функции lny=f(x) (таблица 1.1). Степени полиномов 3, 4, 7. Построить блок-схему алгоритма. Создать программу на языке Delphi, реализующую разработанный алгоритм [1].

Таблица 1.1 - Исходные данные

Таблица

.        Рассчитать среднеквадратичные отклонения для каждого из трех случаев по формуле:

3.       Построить графики 3-х полученных приближающих функций в одной системе координат. На графике должны содержаться и исходные точки (х_i,y_i).

.        Решить задачу средствами MathCad.

Результаты решения задачи с помощью созданной программы и в среде MathCad нужно представить в виде построенных с помощью найденных коэффициентов трёх полиномов; таблицы, содержащей полученные с помощью найденных полиномов значения функции в точках х_i и среднеквадратичных отклонений [2].

1.3 Математическая модель

Пусть величина y является функцией аргумента x. Это означает, что любому значению x из области определения поставлено в соответствии значение y. Вместе с тем на практике часто неизвестна явная связь между y и x, т.е. невозможно записать эту связь в виде y=f(x). В некоторых случаях даже при известной зависимости y=f(x) она настолько громоздка (например, содержит трудно вычисляемые выражения, сложные интегралы и т.п.), что ее использование в практических расчетах затруднительно.

Наиболее распространенным и практически важным случаем, когда вид связи между параметрами x и y неизвестен, является задание этой связи в виде некоторой таблицы {x_i y_i}. Это означает, что дискретному множеству значений аргумента {x_i} поставлено в соответствие множество значений функции {y_i} (i=0,1…n). Эти значения - либо результаты расчетов, либо экспериментальные данные. На практике нам могут понадобиться значение величины y и в других точках, отличных от узлов x_i. Однако получить эти значения можно лишь путем очень сложных расчетов или провидением дорогостоящих экспериментов. Таким образом, с точки зрения экономии времени и средств мы приходим к необходимости использования имеющихся табличных данных для приближенного вычисления искомого параметра y при любом значении (из некоторой области) определяющего параметра x, поскольку точная связь y=f(x) неизвестна.

Этой цели и служит задача о приближение (аппроксимации) функций: данную функцию f(x) требуется приближенно заменить (аппроксимировать) некоторой функцией g(x) так, чтобы отклонение (в некотором смысле) g(x) от f(x) в заданной области было минимальным. Функция g(x) при этом называется аппроксимирующей. Для практики весьма важен случай аппроксимации функции многочленом:

g(x)=a₀+a₁x+a₂x²+…+a_mx^m

При этом коэффициенты a_j будут подбираться так, чтобы достичь наименьшего отклонения многочлена от данной функции.

Если приближение строиться на заданном множестве точек {xi}, то аппроксимация называется точечной. К ней относятся интерполирование, среднеквадратичное приближение и др. При построении приближения на непрерывном множестве точек (например, на отрезке [a,b] аппроксимация называется непрерывной или интегральной). Одним из основных типов точечной аппроксимации является интерполирование. Оно состоит в следующем: для данной функции y=f(x) строим многочлен (1.2), принимающий в заданных точках x_i те же значения y_i, что и функция f(x), т.е. g(x_i)=y_i, i=0,1,…n. При этом предполагается, что среди значений x_i нет одинаковых, т.е. x_i¹x_k при этом i¹k. Точки x_i называются узлами интерполяции, а многочлен g(x) - интерполяционным многочленом.

Рисунок 1.1 - График интерполяционного многочлена

Таким образом, близость интерполяционного многочлена к заданной функции состоит в том, что их значения совпадают на заданной схеме точек (рисунок 1.1, пунктирной линией).

Максимальная степень интерполяционного многочлена m=n; в этом случае говорят о глобальной интерполяции.

При большом количестве узлов интерполяции получается высокая степень многочлена (1.2) в случае глобальной интерполяции, т.е. когда нужно уметь один интерполяционный многочлен для всего интервала изменения аргумента. Кроме того, табличные данные могли быть получены путем измерений и содержать ошибки. Построение аппроксимируемого многочлена с условием обязательного прохождения его графика через эти экспериментальные точки означало бы тщательное повторение допущенных при измерениях ошибок. Выход из этого положения может быть найден выбором такого многочлена, график которого проходит близко от данных точек.

Одним из таких видов является среднеквадратичное приближение функции с помощью многочлена (1.2). При этом m £ n; случай m = n соответствует интерполяции. На практике стараются подобрать аппроксимирующий многочлен как можно меньшей степени (как правило, m=1, 2, 3) [3].

Мерой отклонения многочлена g(x) от заданной функции f(x) на множестве точек (x_i,y_i) (i=0,1,…,n) при среднеквадратичном приближении является величина S, равная сумме квадратов разности между значениями многочлена и функции в данных точках:

Для построения аппроксимирующего многочлена нужно подобрать коэффициенты a₀, a₁,…,a_m так, чтобы величина S была наименьшей. В этом состоит метод наименьших квадратов[2].

;

После преобразования система слегка упростится:

 (1.5)

Введём следующие обозначения:

С учётом этих обозначений система (1.5) перепишется следующим образом:

     (1.6)

Среднеквадратичное отклонение вычисляется по формуле[2]:

.4 Алгоритм решения задачи

1.       Начало алгоритма

2.       Для i=0,…,14

.        Обнуляем массивT_i=0;

.        Для i=0,…,7

.        Обнуляем массив C_i=0, A_i=0

.        Для i=0,…,7, для j=0,…,8

.        Обнуляем массив B_ij=0

.        Для i=1,…,m

.        Вычисляем коэффициенты массива T_j=T_j+Xi^j

.        Для j=0,…,n

.        Вычисляем коэффициенты массива C_ij= C_j+Yi*Xi^j

13.     T₀=m

Формируем расширенную матрицу

14.     Для i=0,…,n, для j=0,…,n

15.     B_ij=T_j+i

.        Для i=0,…,n

17.     B_in+1=C_i

Приводим её к треугольному виду:

18.     Для k=0,…,n-1, для i=k,…,n

19.     Bik = B_ik

.        Для j=k,…,n+1

.        Если i=k,

22.     то

23.     иначе

.        Для i=n,…,0

Вычисляем коэффициенты полинома

25.    

26.     Для i=n,…,0

27.     Вывод A_i

Вычисляем среднеквадратичное отклонение

28.     Delta=0

29.     Для i=1,..,m

30.     Yn=0, Для j=0,…,n

31.     Yn= Yn+A_i*Xi^j

32.     Вычисляем сумму квадратов разностей экспериментального значения функции и рассчитанного с помощью полинома:

.        Окончательное вычисление среднеквадратичного отклонения:

.        Конец алгоритма

1.5 Блок-схемы алгоритмов решения поставленной задачи

Рассмотрим блок-схему основного алгоритма решения задачи аппроксимации табулированной функции методом наименьших квадратов (рисунок 1.2).

Рис.

Рис.

Рисунок 1.2 -Блок-схемы алгоритма подпрограммы решения задачи

2. Практическая часть

.1 Решение поставленной задачи

Для решения поставленной задачи была создана программа в среде Delphi (Рисунок 2.1)[Приложение А].

Рисунок 2.1 - Интерфейс программы

Данная форма включает в себя: меню, исходные данные, графики полиномов, поле с полученными результатами.

Для того, чтобы загрузить исходные данные, нажимаем клавишу Исходные данные (Рисунок 2.1). В результате поле Исходные параметры Х и Y (Рисунок 2.1) заполняться экспериментальными данными.

В окнах Степени полинома вводим значения степеней полинома.

После ввода всех исходных данных щелкаем по кнопке Вычислить, в результате чего появятся графики трех полиномов и начальных параметров, а так же все проведенные вычисления и промежуточные числа выведены в окно рядом с графиками.

2.2 Перечень использованных в программе идентифика-торов

Таблица 2.1 - Идентификаторы программы

2.3 Текст программы

.3.1 Фрагмент основных модулей в Delphi

procedure TForm1.Button3Click(Sender: TObject);Pl;,Bik,Delta,Yn:extended;,i,k,q,j: integer;:=1; Pl:q=1 then n:=StrToInt(Edit1.Text)beginq=2 then n:=StrToInt(Edit2.Text)n:=StrToInt(Edit3.Text) end;i:=0 to 14 do T[i]:=0;i:=0 to 7 do begin C[i]:=0; A[i]:=0 end;i:=0 to 7 doj:=0 to 8 do B[i,j]:=0;.Lines.Add('Коэффициенты полинома степени ='+inttostr(n));i:=1 to 2*n do begin:=0;j:=1 to 40 do begin[j]:=StrToFloat(stringgrid1.Cells[1,j]);:=St+exp(i*ln(X[j]));[i]:=St;;i:=0 to n do begin:=0;j:=1 to 40 do begin[j]:=StrToFloat(stringgrid1.Cells[2,j]);:=St+exp(i*ln(X[j]))*Y[j];[i]:=St;;[0]:=40;i:=0 to n doj:=0 to n do[i,j]:=T[j+i];i:=0 to n do B[i,n+1]:=C[i];k:=0 to n-1 doi:=k to n do begin:=B[i,k];j:=k to n+1 doi=k then B[i,j]:=B[i,j]/BikB[i,j]:=B[i,j]/Bik-B[k,j]; end;i:=n downto 0 do begin[i]:=(B[i,n+1]-B[i,1]*A[1]-B[i,2]*A[2]-B[i,3]*A[3]-B[i,4]*A[4]-B[i,5]*A[5]-B[i,6]*A[6]-B[i,7]*A[7]-B[i,8]*A[8])/B[i,i];.Lines.Add('A['+inttostr(i)+']='+floattostr(A[i]));;

Memo1.Lines.Add('Среднеквадратичное отклонение');

Delta:=0;i:=1 to 40 do begin:=0;j:=0 to n do begin:=1;k:=1 to j do St:=St*X[i];:=Yn+A[j]*St end;.Series[q].AddXY(X[i],Yn);.Series[0].AddXY(X[i],Y[i]);:=Delta+sqr(Yn-Y[i]); end;:=sqrt(Delta/40);.Lines.Add(floattostr(Delta));.Lines.Add('');q<3 then begin q:=q+1; goto Pl;;

end;

2.3.2 Фрагмент основных модулей в MathCad

Вычисляем коэффициенты системы уравнений:

Определяем матрицу коэффициентов и матрицу свободных членов системы

Находим коэффициенты полиномов, решив системы

2.4 Распечатка результатов

Рисунок 2.2 - Окно программы с результатами

2.5 Графическое представление результатов

Ниже представлены графики трех полиномов и начальных значений в Delphi (Рисунок 2.4) и в MathCad (Рисунок 2.5).

Рисунок 2.3 - Графики функций(в среде Delphi)

Рисунок 2.4 - Графики функций(в среде MathCad)

2.6 Анализ полученных результатов

При сравнении результатов полученных в среде MathCad с результатами полученными в среде Delphi мы можем наблюдать малую погрешность связанную с округлениями чисел (Таблица 2.2).

Таблица 2.2 - Сравнение результатов

Из выше изложенного подтверждает верность вычисления коэффициентов аппроксимирующего полинома методом наименьших квадратов в программе, созданной в среде программирования. А от степени полинома зависит точность аппроксимации: при высокой степени полинома среднеквадратичное отклонение становится меньше.

Заключение

аппроксимация табулированный математический mathcad

В данной курсовой работе была решена задача аппроксимации табулированных функций. Аппроксимация применяется для определения аналитического вида функции заданной таблично. Задача аппроксимации сводится к определению свободного параметра функции заданного вида, который обеспечит наилучшее приближение функции заданной таблично модельной аналитической функцией.

Наиболее распространенным методом аппроксимации экспериментальных данных является метод наименьших квадратов. Метод позволяет использовать аппроксимирующие функции произвольного вида и относится к группе глобальных методов.

Критерием близости в методе наименьших квадратов является требование минимальности суммы квадратов отклонений от аппроксимирующей функции до экспериментальных точек. Таким образом, не требуется, чтобы аппроксимирующая функция проходила через все заданные точки, что особенно важно при аппроксимации данных, заведомо содержащих погрешности.

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

Для получения данных результатов в курсовой работе был разработан алгоритм решения поставленной задачи. По разработанному алгоритму была создана программа на языке Delphi. Программа апробирована, результаты тестирования показывают работоспособность программы. Также задача была решена в среде Mathcad. Результаты, полученные в MathCad, полностью совпадают с результатами программы. А это говорит о правильности и верности вычислений созданной программы на языке Delphi.

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

.Наранович, О.И. Информатика. Задания и методические указания по выполнению, оформлению и защите курсовых работ для студентов дневной и заочной форм обучения специальностей 1-36 01 01, 1-36 01 03, 1-53 01 01. Часть 4./ О.И. Наранович, С.Г. Скобля. - Барановичи: БарГУ, 2005.

.Наранович, О.И. Информатика: методические указания и задания к лабораторным работам для студентов 2-го курса дневной формы обучения специальностей 40 01 02, 36 01 03, 36 01 01. Часть 3./ О.И. Наранович, С.Г. Скобля. - Барановичи: БарГУ, 2005.

.Бронштейн, И.Н. Справочник по высшей математике для инженеров и учащихся втузов / И.Н.Бронштейн. - М.: Наука , 1981г.

Приложение А

Текст программы заглавной формы

interface, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,, StdCtrls, ExtCtrls, TeeProcs, TeEngine, Chart, XPMan, Grids,;= class(TForm): TStringGrid;: TButton;: TXPManifest;: TChart;: TMemo;: TButton;: TEdit;: TEdit;: TEdit;: TLabel;: TLineSeries;: TLineSeries;: TLineSeries;: TLineSeries;: TButton;FormCreate(Sender: TObject);Button1Click(Sender: TObject);Button2Click(Sender: TObject);Button3Click(Sender: TObject);

{ Private declarations }

{ Public declarations };: TForm1;,C,A,X,Y:array [0..100] of extended;:array [0..40,0..40] of extended;

{$R *.dfm}TForm1.FormCreate(Sender: TObject);i: integer;.Clear;.Cells[1,0]:=' Xi';.Cells[2,0]:=' ln|Yi|';.Cells[3,0]:=' Yi';i:=1 to 40 do stringgrid1.Cells[0,i]:=IntToStr(i);.Text:='3';.Text:='4';.Text:='7';;TForm1.Button3Click(Sender: TObject);Pl;,Bik,Delta,Yn:extended;,i,k,q,j: integer;:=1; Pl:q=1 then n:=StrToInt(Edit1.Text)beginq=2 then n:=StrToInt(Edit2.Text)n:=StrToInt(Edit3.Text) end;i:=0 to 14 do T[i]:=0;i:=0 to 7 do begin C[i]:=0; A[i]:=0 end;i:=0 to 7 doj:=0 to 8 do B[i,j]:=0;.Lines.Add('Коэффициенты полинома степени ='+inttostr(n));i:=1 to 2*n do begin:=0;j:=1 to 40 do begin[j]:=StrToFloat(stringgrid1.Cells[1,j]);:=St+exp(i*ln(X[j]));[i]:=St;;i:=0 to n do begin:=0;j:=1 to 40 do begin[j]:=StrToFloat(stringgrid1.Cells[2,j]);:=St+exp(i*ln(X[j]))*Y[j];[i]:=St;;[0]:=40;i:=0 to n doj:=0 to n do[i,j]:=T[j+i];i:=0 to n do B[i,n+1]:=C[i];k:=0 to n-1 doi:=k to n do begin:=B[i,k];j:=k to n+1 doi=k then B[i,j]:=B[i,j]/BikB[i,j]:=B[i,j]/Bik-B[k,j]; end;i:=n downto 0 do begin[i]:=(B[i,n+1]-B[i,1]*A[1]-B[i,2]*A[2]-B[i,3]*A[3]-B[i,4]*A[4]-B[i,5]*A[5]-B[i,6]*A[6]-B[i,7]*A[7]-B[i,8]*A[8])/B[i,i];.Lines.Add('A['+inttostr(i)+']='+floattostr(A[i]));

end;.Lines.Add('Среднеквадратичное отклонение');

Delta:=0;i:=1 to 40 do begin:=0;j:=0 to n do begin:=1;k:=1 to j do St:=St*X[i];:=Yn+A[j]*St end;.Series[q].AddXY(X[i],Yn);.Series[0].AddXY(X[i],Y[i]);:=Delta+sqr(Yn-Y[i]); end;:=sqrt(Delta/40);.Lines.Add(floattostr(Delta));.Lines.Add('');q<3 then begin q:=q+1; goto Pl;;;TForm1.Button2Click(Sender: TObject);.Close;;TForm1.Button1Click(Sender: TObject);F:TextFile;: string;: integer;: array [0..40] of extended;:real;(F, 'Xi.txt');(F);i:=1 to 40 do begin Readln(F, S);.Cells[1,i]:=S;;(F);(F, 'Yi.txt');(F);i:=1 to 40 do begin Readln(F, S);[i]:=StrToFloat(S);:=ln(L[i]);.Cells[2,i]:=FloatToStr(n);.Cells[3,i]:=S;;(F);.Enabled:=True;

end;.

Приложение Б

Текст программы в MathCad

Рис.

Вычисляем коэффициенты системы уравнений:

Определяем матрицу коэффициентов и матрицу свободных членов системы

Находим коэффициенты полиномов, решив системы