СЛАУ по Гауссу, чтение массива из файла

@xfqrf · Регистрация: 25.04.2016

Author24 — интернет-сервис помощи студентам

Помогите, пожалуйста, исправить. Выдаёт ошибку External:SIGFTP, никак не понимаю, что он просит и как исправить.

Pascal

var
 
a: array[1..4, 1..4] of real;
b: array[1..4] of real;
arr: array [1..4, 1..5] of real;
x: array[1..5] of real;
n, i, j, k: integer; r, g: real;
f1,f2:text;
 
begin
assignfile(f1,'input.txt');
reset(f1);
assignfile(f2,'data.txt');
rewrite(f2);
n:=4;
 i := 0;
 while not Eof(f1) do begin
   j := 0;
   while not Eoln(f1) do begin
     Read(f1, arr[i,j]);
     Inc(j);
   end;
   ReadLn(f1);
   Inc(i);
 end;
 {for i:=1 to n do begin
   for j:=1 to n+1 do begin
A[i,j] :=arr[i,j+1] ;
B[i]:=arr[i,1];
end;end;
 for i:=1 to n do
 x[i]:=0;}
A[1,1]:=arr[1,2];  A[1, 2] := arr[1,3]; A[1, 3] := arr[1,4];A[1, 4] := arr[1,5];
A[2,1]:= arr[2,2]; A[2, 2] := arr[4,3]; A[2, 3] := arr[2,4];A[2, 4] := arr[2,5];
A[3,1]:= arr[3,2]; A[3, 2] := arr[4,3]; A[3, 3] := arr[3,4];A[3, 4] := arr[3,5];
A[4,1]:= arr[4,2]; A[4, 2] := arr[4,3]; A[4, 3] := arr[4,4];A[4, 4] := arr[4,5];
b[1] := arr[1,1]; b[2] := arr[2,1]; b[3] := arr[3,1];b[4] := arr[4,1];
 
{writeln('Матрица A: ');
for i := 1 to n do
begin
for j := 1 to n do
begin
write(a[i,j]:5:2);
a1[i,j]:=a[i,j];
end;
writeln;
end;
writeln;
Writeln('Матрица B: ');
for i:=1 to n do
writeln(b[i]);}
for k := 1 to n do { прямой ход Гаусса }
begin
for j := k + 1 to n do
begin
r := a[j, k] / a[k, k];
for i := k to n do
begin
a[j, i] := a[j, i] - r * a[k, i];
end;
b[j] := b[j] - r * b[k];
end;
end;
for k := n downto 1 do { обратный ход Гаусса }
begin
r := 0;
for j := k + 1 to n do
begin
g := a[k, j] * x[j];
r := r + g;
end;
if a[k,k]<>0 then
x[k] := (b[k]-r)/a[k,k]
//else  writeln(f2,'0!');
end;
//writeln('Корни системы:');
for i := 1 to n do
write(f2,'w[', i, ']=', x[i]:0:2, ' ');
close(f1);close(f2);
end.

ФедосеевПавел · 10.12.2016, 12:59

В строке 20 при чтении самого первого числа из файла чему будут равны индексы i и j? А какой диапазон индексов у переменной arr?

@xfqrf · 11.12.2016, 14:27 **[ТС]**

Точно, спасибо!

Добавлено через 17 часов 18 минут
Простите ещё раз, но не могу заставить его правильно считать. При вводе матрицы
25 1 1 1 1
20 4 3 2 1
25 4 1 1 4
22 4 1 4 1
должны быть ответы -14 12 13 14, а получаются 11,-4.67, -3.67, - 2.67.
что не так, скажите, пожалуйста?

ФедосеевПавел · 11.12.2016, 15:02

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

@xfqrf · 11.12.2016, 18:02 **[ТС]**

Проверил сами ходы. Отлично считает, ответы правильные. Но только если вводить с клавиатуры. Что-то не так с чтением из файла, но не понимаю что, пробовал делать другой цикл - всё поломал.

ФедосеевПавел · 11.12.2016, 20:17

Приведите последний вариант, который считает правильно из констант, но неправильно и файла.

@xfqrf · 11.12.2016, 20:40 **[ТС]**

Pascal

var
 
a: array[1..4, 1..4] of real;
b: array[1..4] of real;
//arr: array [1..4, 1..5] of real;
x: array[1..4] of real;
n, i, j, k: integer; r, g: real;
//f1,f2:text;
 
begin
{assign(f1,'input.txt');
reset(f1);
assign(f2,'data.txt');
rewrite(f2);}
n:=4;
 {i := 1;
 while not Eof(f1) do begin
   j := 1;
   while not Eoln(f1) do begin
     Read(f1, arr[i,j]);
     Inc(j);
   end;
   ReadLn(f1);
   Inc(i);
 end;
 for i:=1 to n do
 for j:=1 to n+1 do
 writeln(f2,arr[i,j]);
 for i:=1 to n do begin
   for j:=1 to n+1 do begin
A[i,j] :=arr[i,j+1] ;
B[i]:=arr[i,1];
end;end;
 }
A[1,1]:=1;//arr[1,2];
A[1, 2] := 1;//arr[1,3];
A[1, 3] := 1;//arr[1,4];
A[1, 4] := 1;//arr[1,5];
A[2,1]:= 4;//arr[2,2];
A[2, 2] := 3;//arr[4,3];
A[2, 3] := 2;//arr[2,4];
A[2, 4] := 1;//arr[2,5];
A[3,1]:= 4;//arr[3,2];
A[3, 2] := 1;//arr[4,3];
A[3, 3] := 1;//arr[3,4];
A[3, 4] := 4;//arr[3,5];
A[4,1]:= 4;//arr[4,2];
A[4, 2] := 1;//arr[4,3];
A[4, 3] := 4;//arr[4,4];
A[4, 4] := 1;//arr[4,5];
b[1] := 25;//arr[1,1];
b[2] := 20;//arr[2,1];
b[3] := 25;//arr[3,1];
b[4] := 22;//arr[4,1];
 
 
for k := 1 to n do { прямой ход Гаусса }
begin
for j := k + 1 to n do
begin
r := a[j, k] / a[k, k];
for i := k to n do
begin
a[j, i] := a[j, i] - r * a[k, i];
end;
b[j] := b[j] - r * b[k];
end;
end;
for k := n downto 1 do { обратный ход Гаусса }
begin
r := 0;
for j := k + 1 to n do
begin
g := a[k, j] * x[j];
r := r + g;
end;
x[k] := (b[k]-r)/a[k,k]
 
end;
 
for i := 1 to n do
write(x[i]:0:2,' ');
readln;readln;
//write(f2,'w[', i, ']=', x[i]:0:2, ' ');
//close(f1);close(f2);}
end.

ФедосеевПавел · 11.12.2016, 21:06

Сделайте формат файла следующим:

Код

4
1.00 1.00 1.00 1.00 25.00
4.00 3.00 2.00 1.00 20.00
4.00 1.00 1.00 4.00 25.00
4.00 1.00 4.00 1.00 22.00

Первой строкой - количество строк в системе
Далее - расширенная матрица - последний столбец - b.

Pascal

{.$define NoFile}
program SLAU;
 
type
  TMatrix = array[1..4, 1..4] of real;
  TArray  = array[1..4] of real;
var
{$ifdef NoFile}
  a: TMatrix = ((1, 1, 1, 1), (4, 3, 2, 1), (4, 1, 1, 4), (4, 1, 4, 1));
  b: TArray = (25, 20, 25, 22);
  n: integer = 4;
{$else}
  a: TMatrix;
  b: TArray;
  n: integer;
{$endif}
  x: TArray;
  i, j, k: integer;
  r, g: real;
  f1, f2: Text;
 
begin
{$ifndef NoFile}
  Assign(f1, 'input.txt');
  reset(f1);
  readln(f1, n);
  for i := 1 to n do
  begin
    for j := 1 to n do
      Read(f1, a[i, j]);
    readln(f1, b[i]);
  end;
  Close(f1);
{$endif}
 
  Assign(f2, 'data.txt');
  rewrite(f2);
 
  for i := 1 to n do
  begin
    for j := 1 to n do
      Write(f2, a[i, j]: 0: 2, ' ');
    writeln(f2, b[i]: 0: 2);
  end;
 
  for k := 1 to n do { прямой ход Гаусса }
  begin
    for j := k + 1 to n do
    begin
      r := a[j, k] / a[k, k];
      for i := k to n do
      begin
        a[j, i] := a[j, i] - r * a[k, i];
      end;
      b[j] := b[j] - r * b[k];
    end;
  end;
  for k := n downto 1 do { обратный ход Гаусса }
  begin
    r := 0;
    for j := k + 1 to n do
    begin
      g := a[k, j] * x[j];
      r := r + g;
    end;
    x[k] := (b[k] - r) / a[k, k];
 
  end;
 
  for i := 1 to n do
    Write(f2, x[i]: 0: 2, ' ');
  readln;
  //write(f2,'w[', i, ']=', x[i]:0:2, ' ');
  Close(f2);
end.

@xfqrf · 11.12.2016, 21:43 **[ТС]**

Спасибо огромное! всё работает! и гораздо лучше, чем часы!

ФедосеевПавел · 11.12.2016, 21:48

Смотрите - там используется условная компиляция - когда куски кода или попадают или исключаются из компиляции.
Это, если будет преподаватель спрашивать.
В зависимости от того, определен символ NoFile (строка 1) параметры для вычислений будут взяты из констант или из файла.

@xfqrf · 11.12.2016, 22:19 **[ТС]**

Да, я понимаю. Спасибо. А без них можно? и тогда я удалю лишнее из var.

ФедосеевПавел · 11.12.2016, 22:24

Конечно, можно. Это же только для отладки.

ФедосеевПавел Модератор 8476 / 4335 / 1642 Регистрация: 01.02.2015 Сообщений: 13,461 Записей в блоге: 8
	10.12.2016, 12:59	2
	В строке 20 при чтении самого первого числа из файла чему будут равны индексы i и j? А какой диапазон индексов у переменной arr? 1

@xfqrf 0 / 0 / 0 Регистрация: 25.04.2016 Сообщений: 8
	11.12.2016, 14:27 [ТС]	3
	Точно, спасибо! Добавлено через 17 часов 18 минут Простите ещё раз, но не могу заставить его правильно считать. При вводе матрицы 25 1 1 1 1 20 4 3 2 1 25 4 1 1 4 22 4 1 4 1 должны быть ответы -14 12 13 14, а получаются 11,-4.67, -3.67, - 2.67. что не так, скажите, пожалуйста? 0

ФедосеевПавел Модератор 8476 / 4335 / 1642 Регистрация: 01.02.2015 Сообщений: 13,461 Записей в блоге: 8
	11.12.2016, 15:02	4
	Вечером вернусь и, если не забуду, посмотрю. Но вы не теряйте время - проверьте прямой и обратный ходы. 0

@xfqrf 0 / 0 / 0 Регистрация: 25.04.2016 Сообщений: 8
	11.12.2016, 18:02 [ТС]	5
	Проверил сами ходы. Отлично считает, ответы правильные. Но только если вводить с клавиатуры. Что-то не так с чтением из файла, но не понимаю что, пробовал делать другой цикл - всё поломал. 0

ФедосеевПавел Модератор 8476 / 4335 / 1642 Регистрация: 01.02.2015 Сообщений: 13,461 Записей в блоге: 8
	11.12.2016, 20:17	6
	Приведите последний вариант, который считает правильно из констант, но неправильно и файла. 0

@xfqrf 0 / 0 / 0 Регистрация: 25.04.2016 Сообщений: 8
	11.12.2016, 21:43 [ТС]	9
	Спасибо огромное! всё работает! и гораздо лучше, чем часы! 0

ФедосеевПавел Модератор 8476 / 4335 / 1642 Регистрация: 01.02.2015 Сообщений: 13,461 Записей в блоге: 8
	11.12.2016, 21:48	10
	Смотрите - там используется условная компиляция - когда куски кода или попадают или исключаются из компиляции. Это, если будет преподаватель спрашивать. В зависимости от того, определен символ NoFile (строка 1) параметры для вычислений будут взяты из констант или из файла. 0

@xfqrf 0 / 0 / 0 Регистрация: 25.04.2016 Сообщений: 8
	11.12.2016, 22:19 [ТС]	11
	Да, я понимаю. Спасибо. А без них можно? и тогда я удалю лишнее из var. 0

ФедосеевПавел Модератор 8476 / 4335 / 1642 Регистрация: 01.02.2015 Сообщений: 13,461 Записей в блоге: 8
	11.12.2016, 22:24	12
	Конечно, можно. Это же только для отладки. 1

СЛАУ по Гауссу, чтение массива из файла

Решение