590

Імітаційне моделювання. Функція генератор випадкових (псевдо) чисел

Отчет о прохождении практики

Математика и математический анализ

Написати функцію генератор випадкових(псевдо) чисел. Дослідити поведінку ЛК на зміну параметрів. Побудувати графік. Реалізована функція xn+1=(axn+c) mod m. З вхідними параметрами. За допомогою функції і оримали масив 1000 значень. ПСЧ отримались в проміжку від 1 до 29.

Украинкский

2012-11-28

66.5 KB

4 чел.

PAGE  2

МІНІСТЕРСТВО ОСВІТИ І НАУКИ МОЛОДІ ТА СПОРТУ

КРИВОРІЗЬКИЙ ПЕДАГОГІЧНИЙ ІНСТИТУТ ДВНЗ КНУ

Кафедра інформатики і

прикладної математики

ЗВІТ

Імітаційне моделювання

Виконав:

студент 4-го курсу

групи  І-08

Крепчук О.О.

Керівник:

Євтєєв В.М.




м. Кривий Ріг

2012 р.

Завдання: Написати  функцію «генератор випадкових(псевдо) чисел». Дослідити поведінку ЛК на зміну параметрів. Побудувати графік.

Додаток функції майн.

n=1000;

x = randint(n,1);

  x(1)=100;

m=29;

  y=zeros(m,1);  

  a=1;

  c=20;

  for i=2:n

     x=LK(m,a,c,x,i);

  end;

  y=podshet(x,y,m,n);

  showplot(y);

Параметри ЛК

М = 29, краще брати м в яких нод = 1 і самому м.

С - повинно бути менше за м.

А – повинно бути більше нуля і менше М.

Додаток функції генератора.

 function x=LK(m,a,c,x,i)

x(i)=mod((a*x(i-1)+c),m);

end

Все як в лекціях. Реалізована функція xn+1=(axn+c) mod m. З вхідними параметрами. 

За допомогою функції і оримали масив 1000 значень. ПСЧ отримались в проміжку від 1 до 29 тобто до М.

Рис1. Вивід масиву наших ЛК чисел з нашими параметрами.

Інші параметри М=23, кількість в масиві 10000. а=1,с=21,першочергове значення 90.

Рис.2. Інші параметри.

І так з рис.2. видно, що генерує по рівномірному закону розподілу.

Додаток Експ закон.

function k = Expon( x,n )

for i=1:n

   k(i)=(-log(1-x(i))/2);

end;

end

Передаємо наш масив функції експ і вона реалізує закон Експ. З тими самими параметрами.

Додаток Ерленга закон.

%распределение Эрланга +

for i=1:10000

   sum=0;

   for y=1:l-1

       n=rand(1);

       sum=sum+(((lam*n)^y)/prod(1:y))*exp(-lam*n);

   end;

   a(i)=1-sum;

end;

Рис.3.

Додаток Норм. закон.

%нормальное распределение +

for i=1:10000

   a(i)=from+(to-from)*randn(1);

end;

Рис.4.

Висновок: Параметр м впливає на проміжок генерації. М також повинно мати нод 1 і саме себе. Розглянуті 4 закони розподілу. Проаналізовані параметри які входять в ЛК. М,а,с, і початковий параметр х0.

Припускали, що m>2, 1<а<m ( а=0 або а=1 брати безглуздо, а всі інші по модулю m еквівалентні числам з діапазону 0а<m), 0c<m і 0x0<m. Нам буде зручно виділити три випадки:

 xn+1=(a xn+c) mod m,  c0

(2)

xn+1=a xn mod m,  m – составе число.

(3)

xn+1=a xn mod m,  mпросте число.

(4)

Генератори (3), (4) є окремим випадком (2) і визначаються трійкою параметрів {m, a, x0}, причому насправді у випадку (4) вибір x0 байдужий, при 1 x0<m виходить один і той же генератор.

М повинно бути простим, тому що ми потрапимо на період. Тобто буде повторюватися наші ПСЧ.


 

А также другие работы, которые могут Вас заинтересовать

16684. НЕКОТОРЫЕ ПРОБЛЕМЫ ПРАВА НА НЕЖИЛОЕ ПОМЕЩЕНИЕ 50.09 KB
  НЕКОТОРЫЕ ПРОБЛЕМЫ ПРАВА НА НЕЖИЛОЕ ПОМЕЩЕНИЕ К.И. СКЛОВСКИЙ К.И. Скловский доктор юридических наук. Известно что такой объект как нежилое помещение является одним из самых распространенных в гражданском обороте а в сфере аренды недвижимости решительно пре
16685. НЕКОТОРЫЕ ПРОБЛЕМЫ РЕСТИТУЦИИ 67.14 KB
  НЕКОТОРЫЕ ПРОБЛЕМЫ РЕСТИТУЦИИ К.И. СКЛОВСКИЙ К.И. Скловский доктор юридических наук. Вопросу о котором пойдет речь посвящено достаточное количество исследований тем не менее он не утрачивает своей остроты. Более того кажется что за последние годы он стал ещ
16686. ПРОБЛЕМЫ ИЗЪЯТИЯ ЗАЛОГА 12.83 KB
  ПРОБЛЕМЫ ИЗЪЯТИЯ ЗАЛОГА К.И. СКЛОВСКИЙ Скловский Константин Ильич доктор юридических наук адвокат профессор кафедры гражданского права КГАУ имеет около 120 научных публикаций. Главное качество залога право следования за вещью в силу которого залогодержате...
16687. НЕКОТОРЫЕ ПРОБЛЕМЫ ОСПАРИВАНИЯ В СУДЕ РЕАЛИЗАЦИИ ИМУЩЕСТВА НА ТОРГАХ 52.63 KB
  НЕКОТОРЫЕ ПРОБЛЕМЫ ОСПАРИВАНИЯ В СУДЕ РЕАЛИЗАЦИИ ИМУЩЕСТВА НА ТОРГАХ К.И. СКЛОВСКИЙ Скловский Константин Ильич доктор юридических наук адвокат профессор кафедры гражданского права КГАУ имеет около 120 научных публикаций. Нарастающее число споров о признан...
16688. ВЛИЯНИЕ ПРОЦЕДУРЫ ГОСУДАРСТВЕННОЙ РЕГИСТРАЦИИ НА ЧАСТНЫЕ ОТНОШЕНИЯ 29.53 KB
  О ВЛИЯНИИ ПРОЦЕДУРЫ ГОСУДАРСТВЕННОЙ РЕГИСТРАЦИИ НА ЧАСТНЫЕ ОТНОШЕНИЯ К.И. СКЛОВСКИЙ Скловский Константин Ильич доктор юридических наук адвокат профессор кафедры гражданского права КГАУ имеет около 120 научных публикаций. Соотношение публичных и частнопра...
16689. ОТНОШЕНИЯ СОБСТВЕННИКА С НЕЗАКОННЫМ ВЛАДЕЛЬЦЕМ И ПРИОБРЕТАТЕЛЬНАЯ ДАВНОСТЬ 19.08 KB
  ОТНОШЕНИЯ СОБСТВЕННИКА С НЕЗАКОННЫМ ВЛАДЕЛЬЦЕМ И ПРИОБРЕТАТЕЛЬНАЯ ДАВНОСТЬ К.И. СКЛОВСКИЙ Скловский Константин Ильич доктор юридических наук адвокат профессор кафедры гражданского права КГАУ имеет около 120 научных публикаций. Во владении физических и ос
16690. НЕГАТОРНЫЙ ИСК: ОТРИЦАНИЕ НЕ ДОЛЖНО БЫТЬ ОГУЛЬНЫМ 20.77 KB
  НЕГАТОРНЫЙ ИСК: ОТРИЦАНИЕ НЕ ДОЛЖНО БЫТЬ ОГУЛЬНЫМ К.И. СКЛОВСКИЙ Я отрицаю все и в этом суть моя признание героя популярной немецкой трагедии могло бы стать девизом отечественной правоприменительной практики в связи с неоправданно широким распространением не
16691. ПРИМЕНЕНИЕ НОРМ О ЗЛОУПОТРЕБЛЕНИИ ПРАВОМ В СУДЕБНОЙ ПРАКТИКЕ 21.93 KB
  О ПРИМЕНЕНИИ НОРМ О ЗЛОУПОТРЕБЛЕНИИ ПРАВОМ В СУДЕБНОЙ ПРАКТИКЕ К.И. СКЛОВСКИЙ Скловский Константин Ильич доктор юридических наук адвокат профессор кафедры гражданского права КГАУ имеет около 120 научных публикаций. Проблеме злоупотребления правом хорошо
16692. ЗАЩИТА АРЕНДАТОРА ОТ ВТОРОЙ АРЕНДЫ 9.49 KB
  О ЗАЩИТЕ АРЕНДАТОРА ОТ ВТОРОЙ АРЕНДЫ К.И. СКЛОВСКИЙ Скловский Константин Ильич доктор юридических наук адвокат профессор кафедры гражданского права КГАУ имеет около 120 научных публикаций. Собственник вправе неоднократно распорядиться вещью причем все его...