- •1.Введение
- •2.Численное интегрирование
- •2.1 Метод трапеций
- •3.2 Схема алгоритма
- •3.3 Текст программы
- •3.4 Результат работы программы
- •3.Решение нелинейного уравнения
- •3.1 Решение нелинейного уравнения методом хорд-касательных.
- •3.2 График заданной функции
- •3.2 Текст программы.
- •3.3 Результат работы программы.
- •3.4 Алгоритм программы
- •4. Решение системы линейных алгебраических уравнений.
- •4.2. Описание метода
- •4.2. Схема алгоритма
- •4.3. Текст программы
- •4.4. Результаты работы программы
- •5. Проверка в программе wxMaxima
- •5.1Численное интегрирование
- •5.2 Нелинейное уравнение, слау
- •6.Заключение
- •7.Список использованной литературы
- •7.1 Реферативное описание прикладной программы Opera.
3.2 Схема алгоритма
Рис. 2.2
3.3 Текст программы
CLS
Const pi = 3.14
A = pi/3
B = pi/2
PRINT "e=";
INPUT e
n1 = 0
n2 = 1
DO
s1 = 0
s2 = 0
n1 = n1 + 1
h = (b - a) / n1
FOR x = a TO (b - h) STEP h
y1 = 1 / (sin(2*x)^2)
y2 = 1 / ((sin(2*x))^2+h)
s1 = s1 + ((y2 + y1) * h) / 2
NEXT x
n2 = n2 + 1
h = (b - a) / n2
FOR x = a TO (b - h) STEP h
y1 = 1 / (sin(2*x))^2
y2 = 1 / ((sin(2*x))^2 + h)
s2 = s2 + ((y2 + y1) * h) / 2
NEXT x
LOOP UNTIL ABS(s2 - s1) < e
PRINT "s="; s2
Проверка предоставлена в пункте 5
3.4 Результат работы программы
Рис. 2.3
3.Решение нелинейного уравнения
Общие сведения.
Уравнение типа F(x)=0 или x=f(x) называется нелинейным. Решить
уравнение это значит найти такое x, при котором уравнение превращается в
тождество. В общем случае уравнение может иметь 0; 1; 2;...∞ корней.
Рассмотренные ниже численные методы решения нелинейных уравнений
позволяют находить один корень на заданном интервале [a,b]. При этом на
интервале должен существовать только один корень.
Решение уравнения складывается из двух этапов:
1. отделения корня, т.е. грубой оценки значения корня, которое используется в
качестве начального приближения при его последующем уточнении;
2. уточнения корня, т.е. доведения его численного значения до заданной степени
точности.
Для отделения корня (определения начального приближения) следует
пользоваться графическим методом. Уточнение корня может быть произведено
одним из следующих методов.
3.1 Решение нелинейного уравнения методом хорд-касательных.
Описание метода:
Если в методе касательных производную функцию F’(Xi) заменить отношением конечных приращений, то получаем расчетную формулу для метода
X(i+1)=x(i)-
Исходные данные: x=tan(x)
3.2 График заданной функции
Рис 3.1
3.2 Текст программы.
REM
INPUT e
INPUT x2
IF ((tan(x2)-x2)*(1/cos(x2))^2-1))<0 GOTO 15
X3=x2-((tan(x2)-x2)/(1/cos(x2))^2-1)
IF ABS(x2-x3)<e THEN
GOTO 15
END IF
X1=X2
X2=X3
X3=X2-((tan(x2)-x2)*(x2-x3))-((tan(x2)-x2)-(tan(x1)-x1))
LOOP UNTIL ABS(x2-x3)<e
PRINT “x3=”;x3
Проверка предоставлена в пункте 5.
3.3 Результат работы программы.
Рис 3.2
3.4 Алгоритм программы
Рис. 3.3
4. Решение системы линейных алгебраических уравнений.
Реализовать программу на языке Q-Basic, которая производит вычисление корней системы линейных алгебраических уравнений (СЛАУ) методом Гаусса.
Исходная система:
-
21x1+3x2+4x3=20
x1+20x2+2x3=22
18x1+4x2+24x3=30
4.2. Описание метода
Метод Гаусса является одним из самых распространенных методов решения систем линейных уравнений. Этот метод еще называют методом последовательного исключения неизвестных.
Рассмотрим базовый вариант метода Гаусса, называемый схемой единственного деления.
ПРЯМОЙ ХОД состоит из n - 1 шагов исключения.
1-й шаг. Целью этого шага является исключение неизвестного x1 из уравнений с номерами i = 2, 3, …, n. Предположим, что коэффициент a11 ( 0. Будем называть его главным элементом 1-го шага.
Найдем величины qi1 = ai1/a11 (i = 2, 3, …, n), называемые множителями 1-го шага. Вычтем последовательно из второго, третьего, …, n-го уравнений системы первое уравнение, умноженное соответственно на q21, q31, …, qn1. Это позволит обратить в нуль коэффициенты при x1 во всех уравнениях, кроме первого. В результате получим эквивалентную систему
a11x1 + a12x2 + a13x3 + … + a1nxn = b1 ,
a22(1)x2 + a23(1)x3 + … + a2n(1)xn = b2(1) ,
a32(1)x2 + a33(1)x3 + … + a3n(1)xn = b3(1) ,
. . . . . . . .
an2(1)x2 + an3(1)x3 + … + ann(1)xn = bn(1) .
в которой aij(1) и bij(1) вычисляются по формулам
aij(1) = aij - qi1a1j , bi(1) = bi - qi1b1.
2-й шаг. Целью этого шага является исключение неизвестного x2 из уравнений с номерами i = 3, 4, …, n. Пусть a22(1) 0, где a22(1) – коэффициент, называемый главным (или ведущим) элементом 2-го шага. Вычислим множители 2-го шага
qi2 = ai2(1) / a22(1) (i = 3, 4, …, n)
и вычтем последовательно из третьего, четвертого, …, n-го уравнения системы второе уравнение, умноженное соответственно на q32, q42, …, qn2. В результате получим систему
a11x1 + a12x2 + a13x3 + … + a1nxn = b1 ,
a22(1)x2 + a23(1)x3 + … + a2n(1) = b2(1) ,
a33(2)x3 + … + a3n(2)xn = b3(2) ,
. . . . . . . . . . . . . . . . . .
. an3(2)x3 + … + ann(2)xn = bn(2) .
Здесь коэффициенты aij(2) и bij(2) вычисляются по формулам
aij(2) = aij(1) – qi2a2j(1) , bi(2) = bi(1) – qi2b2(1).
Аналогично проводятся остальные шаги. Опишем очередной r-й шаг.
r-й шаг. В предположении, что главный (ведущий) элемент r-го шага arr(r–1) отличен от нуля, вычислим множители r-го шага
qir = air(r–1) / arr(r–1) (i = r + 1, …, n)
и вычтем последовательно из (r + 1)-го, …, n-го уравнений полученной на предыдущем шаге системы r-e уравнение, умноженное соответственно на qr+1,r, qr+2,r, …, qnr.
После (n - 1)-го шага исключения получим систему уравнений
a11x1 + a12x2 + a13x3 + … + a1nxn = b1 ,
a22(1)x2 + a23(1)x3 + … + a2n(1)xn = b2(1) ,
a33(2)x3 + … + a3n(2)xn = b3(2) ,
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
ann(n–1)xn = bn(n–1) .
матрица A(n-1) которой является верхней треугольной. На этом вычисления прямого хода заканчиваются.
ОБРАТНЫЙ ХОД. Из последнего уравнения системы находим xn. Подставляя найденное значение xn в предпоследнее уравнение, получим xn–1. Осуществляя обратную подстановку, далее последовательно находим xn–1, xn–2, …, x1.
Вычисления неизвестных здесь проводятся по формулам
xn = bn(n–1) / ann(n–1),
xr = (bn(r–1) – ar,r+1(r–1)xr+1 – … – arn(r–1)xn) / ar-1r-1, (r = n –1, …, 1).
Вычисление множителей, а также обратная подстановка требуют деления на главные элементы arr(r–1). Поэтому если один из главных элементов оказывается равным нулю (или близким к нулю), то алгоритм не может быть реализован, необходимо принимать специальные меры.
i = 1,n,1
j = 1,n,1
i = 1,n,1