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

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


 

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

2175. Прибыль и рентабельность ОАО Ухтанефтегазгеология 60.11 KB
  Общие сведения о прибыли и рентабельности предприятия. Прибыль и рентабельность на примере предприятия ОАО Ухтанефтегазгеология, основные сведения о предприятий.
2176. Расчет редуктора и его составных частей 60.11 KB
  Краткое описание редуктора и технология его сборки. Выбор электродвигателя. Кинематический и силовой расчет привода. Конструктивные размеры червяка, червячного колеса и корпуса. Подбор подшипников качения.
2177. Етапи створення та класифікація інформаційних систем 25.89 KB
  Основні етапи створення інформаційних систем (ІС). Класифікація ІС. Структура комп’ютерних ІС та класифікація їх задач.
2178. CALS-технологии. И интегрированные автоматизированные системы управления 160.34 KB
  Интегрированные автоматизированные системы управления КИП. Базовые технологии управления данным и информационные модели. Преимущества использования CALS, общие представления об интегрированной информационной среде (ИИС)
2179. Изменение требований к осуществлению надзора за соблюдением законодательства РФ в области образования, государственного контроля качества образования 64.21 KB
  Подготовка образовательного учреждения к выездной проверке в сфере образования. Прохождение образовательным учреждением плановой выездной проверки в сфере образования. Основания для проведения контроля (надзора) в области образования.
2180. Расчет редуктора и ведомого вала с его деталями 125.94 KB
  Целью работы является закрепление теоретических знаний, полученных в ходе изучения материала и приобретения навыков создания современных технических систем.
2181. Воспитательная система школы 953.79 KB
  Содействие развитию личности и формированию компетенций обучающихся через создание для этого благоприятных условий деятельности. Создание условий формирования образованной , свободной, культурной, конкурентоспособной личности обучающегося, способной к творчеству, созиданию, само регуляции , непрерывному образованию , к жизнедеятельности в условиях рыночных отношений.
2182. Отношенческо-коммуникативный компонент воспитательной системы школы 195.24 KB
  Воспитатель личностно значим для воспитанников. Воспитатель социально положительно ориентирован. Толерантность сторон; Компетентность воспитателя. Учет базовых ценностей взаимодействующих субъектов, социума. Взаимодействие не должно приносить вреда другим личностям, обществу; Воспитание на личном примере.
2183. Энергетические системы 806.76 KB
  Общие сведения об энергетических системах. Режимы и параметры системы и сети. Трансформатор с расщепленной обмоткой низшего напряжения. Совместный расчет режима сетей нескольких номинальных напряжений. Расчеты режима линий с двусторонним питанием при различающихся напряжениях источников питания (по концам).