Предыдущая страница реферата | 1  2  3  4  5  6  7  8  9 | Следующая страница реферата

NN=0

KR1=((P-1)/8)*3+1

KR2=((P-1)/8)*5+1

KR3=((M-1)/4)*3+1 do IIY=2,P-1 do IIX=2,M if (NN.eq.N)then print *,'ПЕРЕПОЛНЕНИЕ МАССИВА Y' stop endif

c----------------------------------------------------------------- c проверка границы трубы с жидкостью c----------------------------------------------------------------- if ((IIY.ge.KR1).and.(IIY.le.KR2).and.(IIX.ge.KR3)) then q(IIX,IIY)=200.

c----------------------------------------------------------------- c проверка симметрии c----------------------------------------------------------------- elseif (((IIY.lt.KR1).or.(IIY.gt.KR2)).and.(IIX.eq.M)) then q(IIX,IIY)=6

NN=NN+1 y(1,NN)=IIX y(2,NN)=IIY y(3,NN)=2*HX y(4,NN)=KY y(5,NN)=0 y(6,NN)=KY

c----------------------------------------------------------------- c составление уравнений во внутренних точках фигуры c----------------------------------------------------------------- else q(IIX,IIY)=5

NN=NN+1 y(1,NN)=IIX y(2,NN)=IIY y(3,NN)=HX y(4,NN)=KY y(5,NN)=HX y(6,NN)=KY endif enddo enddo

c----------------------------------------------------------------- c присваивание начальных значений на границе фигуры c------------------------------------------------------------------

KY=0

KR1=P/2+1 do IIY=1,P if (IIY.le.KR1)then q(1,IIY)=0 else q(1,IIY)=500*KY-100 endif

KY=KY+K enddo do IIX=1,M q(IIX,1)=0 q(IIX,P)=100 enddo

N=NN end

ТЕСТ

Для тестирования составлю разностную систему с шагом вдоли оси X и Y=0.05

Неизвестные значения в узлах матрицы находящихся внутри фигуры высчитываются по формуле:
[pic]

Неизвестные значения в узлах матрицы находящихся на оси симметрии высчитываются по формуле:

[pic]

МЕТОД ЗЕЙДЕЛЯ.

Метод Зейделя относится к числу итерационных методов, в которых принципиально отсутствует фактор накопления погрешностей. Поэтому он широко применяется для решения больших систем уравнений. Будем рассматривать корни решаемой системы как компоненты некоторого вектора у . Основная идея всех итерационных методов заключается в том, что берется приближенное значение вектора у и по формулам, составленным на основании решаемых уравнений, вычисляется новое приближенное значение вектора у . Назовем эти приближенные значения y(k) и y(k+1) соответственно. Поскольку исходное приближение выбиралось произвольно, то у(k+1)в свою очередь может послужить исходным для получения по тем же формулам нового приближения y(k+2) .
Очевидно, этот процесс можно продолжать сколь угодно долго. Говорят, что процесс итераций сходится, если получаемая при этом последовательность векторов у(k) (к=0,1,2,...} имеет своим пределом вектор y,являющийся точным решением системы:

[pic]

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

Пусть дана линейная система

[pic]

Рекомендуем скачать другие рефераты по теме: реферат техника, решебник 6 класс виленкин.





Предыдущая страница реферата | 1  2  3  4  5  6  7  8  9 | Следующая страница реферата