реклама на сайте
подробности

 
 
 
Reply to this topicStart new topic
> Генерация случайных чисел
DDA
сообщение Oct 5 2006, 07:33
Сообщение #1


Участник
*

Группа: Новичок
Сообщений: 15
Регистрация: 11-08-06
Пользователь №: 19 479



Вопросец такой. Как можно на компе генерить случайные числа с наперёд заданной плотностью распределения (не нормальная, не равномерная)
Go to the top of the page
 
+Quote Post
jorikdima
сообщение Oct 5 2006, 08:01
Сообщение #2


тут может быть ваша реклама
*****

Группа: Свой
Сообщений: 1 164
Регистрация: 15-03-06
Из: Санкт-Петербург/CA
Пользователь №: 15 280



Использовать МАТЛАБ, Statistiic Toolbox
Go to the top of the page
 
+Quote Post
gab
сообщение Oct 5 2006, 16:12
Сообщение #3


Местный
***

Группа: Свой
Сообщений: 376
Регистрация: 30-06-04
Из: Moskow
Пользователь №: 218



Цитата(jorikdima @ Oct 5 2006, 12:01) *
Использовать МАТЛАБ, Statistiic Toolbox

...или справку от Matlab. Если скорость нужна побольше.


--------------------
serpents on the way to paradise -
dying for love, fighting for ages.

Go to the top of the page
 
+Quote Post
vladv
сообщение Oct 6 2006, 20:25
Сообщение #4


Частый гость
**

Группа: Участник
Сообщений: 128
Регистрация: 7-06-06
Пользователь №: 17 825



Цитата(DDA @ Oct 5 2006, 11:33) *
Вопросец такой. Как можно на компе генерить случайные числа с наперёд заданной плотностью распределения (не нормальная, не равномерная)


Если есть генератор с равномерным распределением в диапазоне [0,1], то генератор с заданной
плотностью p(x), можно сделать с помощью функции F-1(x), которая является обратной функцией
интегрального распределения для заданной плотности (F(X)=P(x<X)): rand_p() = F-1(rand()).

Например, для распределения с "линейной" плотностью на отрезке x=[0,1]:
p(x) = 2*x
То, интегральная функция распределения будет:
F(X) = P(x<X) = X^2
Обратная функция для F(X) будет:
F-1(x) = x^(1/2)
Генератор:
rand()^(1/2)

Правда, вывести формулу аналитически можно далеко не для всех функций smile.gif.
Go to the top of the page
 
+Quote Post

Reply to this topicStart new topic
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0

 


RSS Текстовая версия Сейчас: 18th July 2025 - 08:39
Рейтинг@Mail.ru


Страница сгенерированна за 0.01387 секунд с 7
ELECTRONIX ©2004-2016