Pascal. Програмирование итерационных процессов

LizarD · 13.06.2011, 23:52

Воть: U[i]=(1/x[i])+(1/R)*("знак суммы от бесконечности до n=0"(((n+1)*(a-b))/(n*a+b*(n+1)))*(c*x[i]/R)^n
При R=2.5 a=5.0 b=1.5 c=2.0
x1=0.01 x2=0.02 x3=0.05 x4=0.07 x5=0.1 x6=0.2 x7=0.5 x8=0.7 x9=1.0
Извеняюсь, что формула в корявом виде, если посоветуете как более читаемо выложить её, так и сделаю

Основная загвоздка в написании U[i]

Моего понимания хватило только на следующее:

делее предпологается использовать repeat...until...

Благодарю за прочтение, а тем более желание помочь

Vladimir_S · 14.06.2011, 10:16

Если я правильно понял, и формула выглядит так:

, и при этом (чего Вы не сообщили, а в условии должно быть) вычисление требуется произвести с заданной точностью (е), то я бы это дело решил так:

Код:

CONST
 R=2.5;
 a=5.0;
 b=1.5;
 c=2.0;
 x:Array[1..9] of Real=(0.01, 0.02, 0.05, 0.07, 0.1, 0.2, 0.5, 0.7, 1.0);
 e=0.000001;

VAR
 n:Array[1..9] of Word;
 U:Array[1..9] of Real;
 i,j:Integer;
 Sum,Q:Real;

BEGIN
 For i:=1 to 9 do
  Begin
   Sum:=(a-b)/b;
   j:=0;
   Repeat
    j:=j+1;
    Q:=(j+1)*(a-b)/(j*a+b*(j+1))*Exp(j*Ln(c*x[i]/R));
    Sum:=Sum+Q;
   Until Abs(Q)<e;
   n[i]:=j;
   U[i]:=(1.0/x[i])+Sum/R;
   Writeln('x[',i,']= ',x[i]:4:2,'    U[',i,']= ',U[i]:9:5,'   n= ',n[i]:3);
  End;
 Readln
END.

Выводятся x[i], U[i] и число членов суммы.

LizarD · 14.06.2011, 23:20

Владимир, именно так!
Про точность я и не подумал.

В знаминателе (j-1)

Благодарю, Владимир!

Vladimir_S · 15.06.2011, 09:20

Цитата:

Сообщение от LizarD

В знаминателе (j-1)

Ох! Виноват...

LizarD · 15.06.2011, 09:35

Эти ошибки учат бить внимательным не только вас!

Спасибо...

Vladimir_S · 15.06.2011, 11:14

Цитата:

Сообщение от LizarD

Эти ошибки учат бить внимательным не только вас!

Спасибо...

Ага. Только тогда - еще одно исправление:
строку
S:=(a-b)/b;
следует заменить на
S:=(b-a)/b;
Потому что при j=0 в знаменателе не b, а (-b).

13.06.2011, 23:52	#1 (permalink)
LizarD Новичок Регистрация: 11.06.2011 Сообщений: 6 Сказал(а) спасибо: 0 Поблагодарили 0 раз(а) в 0 сообщениях Репутация: 10	Pascal. Програмирование итерационных процессов Воть: U[i]=(1/x[i])+(1/R)("знак суммы от бесконечности до n=0"(((n+1)(a-b))/(na+b(n+1)))(cx[i]/R)^n При R=2.5 a=5.0 b=1.5 c=2.0 x1=0.01 x2=0.02 x3=0.05 x4=0.07 x5=0.1 x6=0.2 x7=0.5 x8=0.7 x9=1.0 Извеняюсь, что формула в корявом виде, если посоветуете как более читаемо выложить её, так и сделаю Основная загвоздка в написании U[i] Моего понимания хватило только на следующее: делее предпологается использовать repeat...until... Благодарю за прочтение, а тем более желание помочь

14.06.2011, 10:16	#2 (permalink)
Vladimir_S Специалист Регистрация: 27.08.2008 Адрес: Санкт-Петербург Сообщений: 27,807 Сказал(а) спасибо: 340 Поблагодарили 583 раз(а) в 208 сообщениях Репутация: 113184	Если я правильно понял, и формула выглядит так: , и при этом (чего Вы не сообщили, а в условии должно быть) вычисление требуется произвести с заданной точностью (е), то я бы это дело решил так: Код: CONST R=2.5; a=5.0; b=1.5; c=2.0; x:Array[1..9] of Real=(0.01, 0.02, 0.05, 0.07, 0.1, 0.2, 0.5, 0.7, 1.0); e=0.000001; VAR n:Array[1..9] of Word; U:Array[1..9] of Real; i,j:Integer; Sum,Q:Real; BEGIN For i:=1 to 9 do Begin Sum:=(a-b)/b; j:=0; Repeat j:=j+1; Q:=(j+1)(a-b)/(ja+b(j+1))Exp(jLn(cx[i]/R)); Sum:=Sum+Q; Until Abs(Q)<e; n[i]:=j; U[i]:=(1.0/x[i])+Sum/R; Writeln('x[',i,']= ',x[i]:4:2,' U[',i,']= ',U[i]:9:5,' n= ',n[i]:3); End; Readln END. Выводятся x[i], U[i] и число членов суммы.

13.06.2011, 23:52
Helpmaster Member Регистрация: 08.03.2016 Сообщений: 0	Ваш вопрос обязательно решится, а пока вы можете почитать схожие топики Слишком много процессов в диспетчере задач Самовольный запуск процессов opera.exe и др.

14.06.2011, 23:20	#3 (permalink)
LizarD Новичок Регистрация: 11.06.2011 Сообщений: 6 Сказал(а) спасибо: 0 Поблагодарили 0 раз(а) в 0 сообщениях Репутация: 10	Владимир, именно так! Про точность я и не подумал. В знаминателе (j-1) Благодарю, Владимир!

15.06.2011, 09:35	#5 (permalink)
LizarD Новичок Регистрация: 11.06.2011 Сообщений: 6 Сказал(а) спасибо: 0 Поблагодарили 0 раз(а) в 0 сообщениях Репутация: 10	Эти ошибки учат бить внимательным не только вас! Спасибо...

Ads
Яндекс Member Регистрация: 31.10.2006 Сообщений: 40200 Записей в дневнике: 0 Сказал(а) спасибо: 0 Поблагодарили 0 раз(а) в 0 сообщениях Репутация: 55070