Смоделировать непрерывный источник сообщений в соответствии с параметрами, Закон распределения случайной величины

[digitalis]

Не туда добавлено !

program Lab1;

const NMAS=1000;
type tMass=array[1..NMAS] of Extended;

Procedure Otcenki(const x: tMass;var Mx,Dx:Real);
Var
  i:integer;
Begin
  Mx:=0; Dx:=0;
  For i:=1 to NMAS do Mx:=Mx+x[i];
  Mx:=Mx/NMAS;
  For i:=1 to NMAS do Dx:=Dx+Sqr(x[i]-Mx);
  Dx:=Dx/(NMAS-1);
end;

Var
  Mx,Dx,a,b,r:real;
  x:tMass;
  i :word;
  MaxVal,MinVal,M_b : real;
  hist: array [0..14] of integer;
  n:integer;
Begin
  Write('a  b:   ');
  ReadLn(a,b);
  Randomize;
  For i:=1 to NMAS do
   begin
    r:=Random;
    x[i]:=a+r*(b-a);
   end;
for i :=0 to 14 do hist[i] := 0 ;
    M_b := 15.0/(b-a) ;
    for i:=1 to NMAS do 
     begin
       n := Trunc((x[i]-a)*M_b);
       Inc (hist[n])
     end;
   MaxVal := x[i];
   MinVal := x[i];
   for i:=2 to NMAS do
     begin
       if x[i] < MinVal  then  MinVal := x[i];
       if x[i] > MaxVal  then  MaxVal := x[i];
     end;
  Otcenki(x,Mx,Dx);
  Writeln ('Min element is: ', MinVal:8:6,' Max element is: ', MaxVal:8:6);
  WriteLn(Mx:8:6,Dx:8:6);
  for i :=0 to 14 do WriteLn(hist[i]:3) ;
  ReadLn;
End.

В вывод я добавил аттрибуты формата, чтобы вместо чисел-портянок получить внятный для понимания вид.
Теперь дело за малым: запустить прогу (без всяких самодеятельных телодвижений) из #39 в FreePascal, если нормально нарисует, будем ее прикручивать к нашей задаче.

Ну все, я пошел за сухоньким: сегодня моя любимая Феррари одержала двойную победу! Под гимн Италии дерябнем с котиками.

[Мария34]

Да работает код из #39 на FreePascal.

[digitalis]

Отлично. Код прикрутим, если мой комп не загнется окончательно, как он пытался утром.

[Мария34]

С чего начать?

[Мария34]

Вы куда-то пропали.

[Мария34]

Вот вроде все сделала.

uses crt,graph;
const nmas=1000;
      n=14;
var x0,y0:integer;
type tmass=array[1..nmas] of extended;
type mas=array[0..n] of integer;
procedure histogram(a:mas);
var i,w:integer;
    mx,h:real;
    s:string;
begin
w:=10;{ширина столбцов}
mx:=a[0];
for i:=0 to n do
if a[i]>mx then mx:=a[i];{находим максимум}
h:=(y0-40)/mx;{масштаб по вертикали}
for i:=0 to n do
 begin
  setcolor(i+1);{цвет}
  setfillstyle(1,i+1);
  {рисуем прямоугольники с заданной шириной, шагом
    и высотой соответственно значению}
  bar(x0+4*w*i,y0-round(a[i]*h),x0+w*(4*(i+1)-2),y0);
  str(i+1,s);
  outtextXY(x0+4*w*i+15,y0+20,s);{выводим по низу номер}
  str(a[i],s);{преобразуем значения в строку}
  outtextXY(x0+4*w*i,y0-round(a[i]*h)-20,s);{выводим значения выше столбика}
 end;
end;

Procedure Otcenki(const x: tMass;var Mx,Dx:Real);
Var
  i:integer;
Begin
  Mx:=0; Dx:=0;
  For i:=1 to nmas do Mx:=Mx+x[i];
  Mx:=Mx/nmas;
  For i:=1 to nmas do Dx:=Dx+Sqr(x[i]-Mx);
  Dx:=Dx/(nmas-1);
end;

var x:tmass;
    hist:mas;
    MaxVal,MinVal,Mx,Dx,r,a,b,m_b:real;
    m,i,j,k,p:integer;
begin
clrscr;
write('a b:');
readln(a,b);
randomize;
 begin
  for j:=1 to nmas do
   begin
    r:=random;
    x[j]:=a+r*(b-a);
   end;
  for k:=0 to n do hist[k]:=0;
  m_b:=15/(b-a);
  for k:=1 to nmas do
   begin
    m:=trunc((x[k]-a)*m_b);
    inc(hist[m]);
   end;
   MaxVal := x[j];
   MinVal := x[j];
   for j:=2 to NMAS do
     begin
       if x[j] < MinVal  then  MinVal := x[j];
       if x[j] > MaxVal  then  MaxVal := x[j];
     end;
  Otcenki(x,Mx,Dx);
  Writeln ('Min element is: ', MinVal:8:6, ' Max element is: ', MaxVal:8:6);
  WriteLn('Мат.ожидание=',Mx:8:6,'Дисперсия=',Dx:8:6,'‘Среднеквад. отклонение',Sqrt(Dx):8:6);
  for k:=0 to n do
  write(hist[k]:3);
  readln;
  k:=0;
  initgraph(k,p,'');
  x0:=20;{левый край диаграммы}
  y0:=getmaxY-80;{нижний край}
  histogram(hist);
   begin
    readln;
    restorecrtmode;
   end;
 end;
readln
end.

Всем спасибо за помощь.

[digitalis]

Это не я пропал, это мой комп сломался. Починил.
Запустил эту прогу. Вроде работает. Только вместо подписей - крякозябры. Это только у меня в ТурбоПаскале или в ФриПаскале тоже? Разноцветные елочные огоньки можно убрать, рисовать одним цветом.
Ну поздравляю - курсач готов !

[Мария34]

Кракозябры это скорее FreePascal плохо реагирует на кирилицу.

[Мария34]

Уважаемый digitalis, согласитесь помочь еще с одним задание? Оно идет как дополнительное к этому. Правда не знаю сможете помочь нет посмотрите.

Суть задания:
Выполнить программное моделирование алгоритма поиска кодовых комбинаций и алгоритма кодирования методом Хаффмана. Для дискретного источника сообщений, смоделированного в моем первом задании, выполнить кодирование последовательности продолжительностью K = 100 и K = 1000 реализаций. Вычислить значение оценки средней длительности кодовых комбинаций для каждого из двух случаев

[Мария34]

Сам принцеп вроде поняла не могу найти в интернете даже примера чего то подобного.