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 виходить один і той же генератор.

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


 

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

51145. Операционный блок микропрограммируемого процессора 53.57 KB
  Цель работы: изучение операционного блока на уровне структурной схемы, ознакомление с составом микрокоманд и порядком их выполнения, составление и отладка микропрограмм.
51146. Изучение микроконтроллера MSP430F2013 и адаптера eZ430-F2013 587.63 KB
  Микроконтроллеры широко применяются в автоматических системах во всех сферах жизни человека, например, в промышленности, на транспорте, в быту. Микроконтроллеры являются ядром системы управления роботов, мехатронных систем, автоматических устройств. Большой интерес вызывает активно развивающееся семейство микроконтроллеров MSP430 с флэш-памятью и сверхмалым энергопотреблением, производимых корпорацией Texas Instruments.
51147. Изучение способов адресации микроконтроллеров MSP430 204.88 KB
  Задание на лабораторную работу Разработка схемы алгоритма решения задачи Разработка программы Результаты вычислительного эксперимента Выводы Список литературы Введение Производительность микроконтроллера во многом определяется его центральным процессорным устройством ЦПУ. Разработка алгоритма Схема алгоритма...
51148. Программирование циклических алгоритмов 237.19 KB
  Последовательно вводятся N целых чисел. Найти сумму всех отрицательных среди них. Введи количество чисел Сумма отрицательных элементов Результат работы Задание...
51149. Изучение представления графической информации в WINDOWS 129.3 KB
  Цель работы: Написать программу, реализующую просмотр графического файла (формат BMP). Программа должна: загружать и выводить на экран произвольный файл (с использованием файловых функций); осуществлять проверку на допустимый формат файла;
51150. Построение и расчет временных параметров моделей СПУ 129.71 KB
  Согласно номеру своего варианта получите следующие исходные данные: время нормальной длительности каждой работы сетевой модели и описание упорядочения этих работ. В соответствии с сетевыми методиками: рассчитайте и отобразите на сетевом графике временные параметры событий: ранний и поздний срок свершения события резерв события; рассчитайте и представьте в таблице временные параметры работ: время раннего и позднего начала работ; время раннего и позднего окончания работ; полный и свободный резервы...