Технический форум - Генерирование вещественных чисел

Технический форум (http://www.tehnari.ru/)

- Delphi, Kylix and Pascal (http://www.tehnari.ru/f43/)

- - Генерирование вещественных чисел (http://www.tehnari.ru/f43/t87361/)

Виктория 24

10.04.2013 18:49

Генерирование вещественных чисел

Вложений: 2

Здравствуйте! Пишу программу Генератор случайных чисел Методом Середины Квадратов на Delphi.
Мат.ожидание высчитывает правильно (примерно "0,5"). (Цикл мат.ожидания пока заключила в скобки)
Скажите, правильно ли я написала код для вычисления значения дисперсии D = 1/12 (~ 0,083) согласно формуле. (Скрин.) А то что-то не то высчитывает...

Код:

begin

        sum2:= sum + sqr(chisla[j]-0.5);

        inc(j,1);

end;

     Edit6.Text:=floattostr((1/N)*sum2); // Дисперсия

Vladimir_S

10.04.2013 20:29

Цитата:

Сообщение от Виктория 24 (Сообщение 894724)

А то что-то не то высчитывает...

Естественно, не то. Потому что должно быть так:
sum2:= sum2 + sqr(chisla[j]-0.5);

interacia

10.04.2013 20:58

В куске представленного кода нет цикла. Подозреваю что вы его оставили за кадром. Если вычисления проходят в цикле то формула на рисунке реализована программно правильно.
Но вот только
Дисперсия равна разности между математическим ожиданием квадрата случайной величины и квадратом ее математического ожидания.
D(X)=M(X^2)-[M(X)]^2.

Проще, в данной формуле откуда то взялось деление на N.

Возможно так получиться

Код:

sum2:=0;

for j:=1 to N do

begin

        sum2:= sum2 + sqr(chisla[j]-0.5);

end;

в sum2 это будет дисперсия.

Виктория 24

10.04.2013 21:58

Вложений: 1

Цитата:

Естественно, не то. Потому что должно быть так:
sum2:= sum2 + sqr(chisla[j]-0.5);

Да, действительно, заработало!

Цитата:

Проще, в данной формуле откуда то взялось деление на N.Даже не знаю

что вам сказать в свое оправдание...
Может скриншот нашей лекции спасетtehno042
Вот что у меня получилось:

Код:

j:=1;

      sum:=0;

      sum2:=0;

      sum3:=0;

      //for j:=1 to N do

      while j <= High(chisla) do



      begin

        {sum:= sum + chisla[j];

        inc(j,1);

        Edit5.Text:=floattostr((1/N)*sum); // Мат.ожидание}



        sum2:= sum2 + sqr(chisla[j]-0.5);

        inc(j,1);

        Edit6.Text:=floattostr((1/N)*sum2); // Дисперсия



        {sum3:= sum3 + sqr(chisla[j]);

        inc(j,1);

        Edit7.Text:=floattostr((1/N)*sum3); // Сумма квадратов}

      end;

Но как-то кривовато работает, если считает все аспекты сразу!
А если по отдельности, то правильно.

Vladimir_S

11.04.2013 09:49

Цитата:

Сообщение от interacia (Сообщение 894804)

в sum2 это будет дисперсия.

Извините, сударь, но НИ В КОЕМ СЛУЧАЕ! В Sum2 появится лишь сумма квадратов отклонений случайной величины от матожидания, каковая отнюдь дисперсией не является. По "физической" сути этот самый второй момент (дисперсия) есть мера средневзвешенного отклонения реализаций случайной величины от математического ожидания. В куче книжек показано, что в случае равномерного распределения дискретной случайной величины математическое ожидание есть среднее арифметическое всех реализаций, а дисперсия - среднее значение квадрата отклонения, т.е. сумма квадратов отклонений, деленная на количество проб.

interacia

11.04.2013 21:06

Здравствуйте уважаемый Vladimir_S.
С N это действительно я пере умничал или не до умничал.:tehnari_ru_325:
В любом случаи спасибо а то так и помёр бы не просвещенным, за то теперь можно пойти спокойно умереть со стыда.:tehnari_ru_102:

Vladimir_S

11.04.2013 21:57

Цитата:

Сообщение от interacia (Сообщение 895233)

Ой! Что-то последнее время людей стал ненароком обижать... Извините. Кстати, сам, бывает, ляпы допускаю. И всегда прошу меня без экивоков харей тыкать, если кто заметит.
А умирать - всегда успеете. Не надо. Давайте еще поживем. Во благо нерадивых студентов.

Часовой пояс GMT +4, время: 09:01.