А лгоритм метода обращения разностной схемы:

1. Определяем , при ;

2. Находим ;

3. Определяем , при ;

4. Находим ;

5. По формуле (21) определяем ;

6. По формуле (20) находим , при ;

7. По формуле (21) определяем ;

8. По формуле (20) находим , при ;

и т.д.

Блок-схема срс 3

t:=k*h;

u0[k]:=f(t)

S[0]:=u0[0]

S[1]:=u1[0]

u1[k]:=0.5(u0[k]+u0[k+1])

u2[k]:=0.5(S[i]+S[i-2])*(u1[k+1]+u1[k])/S[i-2]-

-u0[k+1]*S[i]/S[i-2]

u0[k]:=u1[k];

u1[k]:=u2[k];

u0[n-i+1]:=u1[n-i+1]

S[i]:=S[i-2]*(u1[1]+S[i-1])/

(2S[i-2]-S[i-1]-u1[1]+2u0[1])

Лабораторная работа № 5.

В лабораторной работе № 5 надо рассчитать значения кажущегося сопротивления и построить график в логарифмическом масштабе, используя 3-точечный фильтр для установки Шлюмберже. Благодаря двукратному применению фильтра значения определяются в 6-ти точках на декаду логарифмического масштаба. Формула для расчета в логарифмическом масштабе

.

где

Коэффициенты фильтра 0,0148; -0,0814; 0,4018; -1,5716; 1,972; 0,1854; 0,1064; -0,0499; 0,0225.

Значения трансформанты сопротивления рассчитываются по соотношениям Пекериса.

Константы и переменные: nt=17; n=18; {число отсчетов}

R[1..10]; {удельные сопротивления}

T[1..nt]; {трансформанты сопротивления}

HS[1..9]; {мощности слоев}

x1, y, S, f, B: real; {х1- абсцисса первой точки измерения, у- текущая абсцисса, S – значение кажущегося сопротивления, f – шаг, B – текущие значения трансформанты сопротивления }

i9, i8, i, j, m : integer; {i9- число слоев, i8=i9-1, i, j, m- счетчики}

Procedure TI (var B: real); {вычисляет текущие значения трансформанты сопротивления}

Блок-схема процедуры ti

a1:=exp(u); a2:=(a1-1/a1)/(a1+1/a1);

B:=(B+a2*R[i])/(1+a2*B/R[i])

да

нет

B:=R[i9]

B:=R[i]

Блок-схема основной программы

F:=exp(ln(10)/6); y=x1*0,0105; B:=0

TI(B); T[nt]:=B; y:=y*f;

; вывод (S)

T[j]:=T[j+1]