Версия для печати темы

Нажмите сюда для просмотра этой темы в обычном формате

Форум «Всё о Паскале» _ Задачи _ Решение СЛАУ методом Монте-Карло

Автор: MagicPRO 6.06.2007 20:58

Народ у меня расчетно-графическая робота, решение системы линейных уравнени методом Монте-Карло!
Может кто знает или у кого-то есть код решения этой задачи, принимается код на С++!
Заранее благодарен!

Автор: klem4 6.06.2007 21:57

В поиске смотрел ?

http://forum.pascal.net.ru/index.php?act=Search&CODE=show&searchid=a134e092020fdcdfe92a625ea01457fe&search_in=posts&result_type=topics&highlite=%2B%CC%EE%ED%F2%E5-%CA%E0%F0%EB%EE

Автор: MagicPRO 10.06.2007 19:41

Помогите пределать код с Basica на Pascal!

10 REM ~~~МОНТЕ_КАРЛО~~~
20 REM ВЫЧИСЛЕНИЯ КОРНЯ СИСТЕМЫ ЛИНЕЙНЫХ
УРАВНЕНИЙ
30 DIM A(20,2l),B(20,20),W(20)
40 INPUT "ВВЕДИТЕ ПОРЯДОК РЕШАЕМОЙ СИСТЕМЫ";N
50 PRINT "ВВЕДИТЕ ПО СТРОКАМ КОЭФФИЦИЕНТЫ
СИСТЕМЫ"
60 FOR I=1 ТО N
70 INPUT A(I,I):B(I,I)=ABS(A(I,I)) 80 FOR J=2 TO N
90 INPUT A(I,J):B(I,J)=B(I,J-I)+ABS(A(I,J))
100 NEXT J
110 INPUT A(I,N+1)
120 NEXT I
130 FOR I=1 TO N
140 W(I)=A(I,N+1)/(1-B(I,N))
150 NEXT I
160 INPUT "ВВЕДИТЕ НОМЕР РЕШАЕМОГО УРАВНЕНИЯ -В И КОЛИЧЕСТВО РЕАЛИЗАЦИИ СЛУЧАЙНОГО ПРОЦЕССА - М ";В,М
170 T=1:Y=0
180 S=B:Y=1
190 IF T>M THEN 260
200 C=RND(1)
210 FOR J=N TO 1 STEP -1
220 IF C<=B(S,J) THEN 240
230 IF J=N THEN T=T+1:Y=Y+V*W(S):GOTO 180 ELSE V=V*SNG(A(S,J+1)):S=J+1:GOTO 200
240 NEXT J
250 V=V*SGN(A(S,l)):S=1:GOTO 200
260 X=Y/M
270 PRINT "КОРЕНЬ ";B;" - ГО УРАВНЕНИЯ X=";X
280 END

Автор: MagicPRO 11.06.2007 13:40

Ура, я код переделал, но теперь есть ошибки с вводом, помогите их решить!
я думаю не одному мне этот метод понадобится, в инете я все обискал, только одну прогу платную нашел!

uses crt;

function SGN(a: real): integer;
begin
if a > 0 then
SGN := 1;
if a = 0 then
SGN := 0;
if a < 0 then
SGN := -1;
end;

var
a: array[1..20, 1..21] of real;
b: array[1..20, 1..20] of real;
w: array[1..20] of real;
n, d, m, t, s: byte;
i, j: word;
x, c, v, y: real;
label
label1, label2, label3;
begin
v := 0;
writeln('vvedite poradok system');
readln(n);
writeln('vvedite coeffs');
for I := 1 to N do
begin
readln(A[I, I]);
B[I, I] := ABS(A[I, I]);
for J := 2 to N do
begin
readln(A[I, J]);
B[I, J] := B[I, J - I] + ABS(A[I, J]);
end;
readln(A[I, N + 1]);
end;
readln;
for I := 1 to N do
W[i] := A[I, N + 1] / (1 - B[I, N]);
writeln('vvedite nomer uravnenia');
readln(d, m);
T := 1;
Y := 0;
label2:
S := d;
Y := 1;
if T > M then
goto label1;
label3:
C := random(1);
for J := N downto 1 do
begin
if C <= B[S, J] then
continue;
if J = N then
begin
T := T + 1;
Y := Y + V * W[S];
goto label2;
end
else
begin
V := V * SGN(A[S, J + 1]);
S := J + 1;
goto label3;
end;
end;
V := V * SGN(A[S, 1]);
S := 1;
goto label2;
label1:
X := Y / M;
writeln('coren ', d, ' - go uravnenia X=', X);
readln;
end.