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

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


 

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

41659. РАБОТА В ПРОГРАММНОЙ СРЕДЕ MICROSOFT OUTLOOK 757.34 KB
  Программная среда Microsoft Outlook пришла на смену разнообразным видам бумажных носителей которые использовали руководители и секретари для организации своей работы. Сегодня для организации документов и отправки почты планирования задач встреч событий и собраний ведения списка контактов а также учета всех выполненных работ используется программа Microsoft Outlook. Информация в среде Outlook организована в виде папок аналогичных по назначению своим бумажным предшественникам.
41660. Поверка средств измерений 39.3 KB
  Поверка средств измерений Цели и задачи работы: Изучение правил организации и порядка проведения поверки средств измерения. Краткие сведения из теории: Поверкой средств измерений называют совокупность действий выполняемых для определения и оценки погрешностей средств измерений. Вид поверки определяют в зависимости от того какой метрологической службой проведена поверка от характера поверки инспекционная экспертная каков этап работы средства измерений первичная периодическая внеочередная. Организацию и поверку средств измерений...
41661. Косвенные измерения. Определение показателей точности косвенных измерений 587.13 KB
  Косвенные измерения. Определение показателей точности косвенных измерений Цели и задачи работы: изучение методов измерения при которых искомое значение физической величины находят путем согласованных наблюдений других величин определяемых опытным путем связанных с искомой физической величиной известной зависимостью; ознакомление с правилами оценивания погрешностей косвенных измерений. При выполнении работы необходимо практически ознакомиться с системой допусков и посадок требованиями к точности линейных и угловых параметров изделий...
41662. Вставка и редактирование формул в редакторе WORD 73.64 KB
  Вставка и редактирование формул. Вставка формул. Вставка формул в редакторе WORD осуществляется с помощью формульного редактора. Вызов формульного редактора Eqution Editor из Word можно осуществить следующей последовательностью действий: поместите курсор в то место где должна быть вставлена формула; в меню вставка выберите команду обьект ; выберите закладку создание ; В окне тип обьекта выберите Microsoft Eqution 3.
41663. Теория электрической связи 263.74 KB
  Получение характеристик частотного модулятора при воздействии на его вход моногармонического сигнала. Напряжение смещения Есм являющееся постоянной составляющей модулирующего сигнала позволяет установить несущую частоту модулированного сигнала а переменная составляющая т. сам модулирующий сигнал поданный на гнезда КТ1 обеспечивает девиацию частоты fmx зависящую от амплитуды модулирующего сигнала. В схеме модулятора имеется блок автоматической регулировки усиления поддерживающий постоянную амплитуду ЧМ сигнала на схеме не показан.
41664. Исследование зависимости выходного напряжения усилительного каскада от амплитуды и частоты входного сигнала 155.55 KB
  Цель: Научиться определять и анализировать зависимости выходного напряжения усилительного каскада от амплитуды и частоты входного сигнала. Лабораторная работа №6 Тема: Исследование зависимости выходного напряжения усилительного каскада от амплитуды и частоты входного сигнала. Лабораторная работа №6 Тема: Исследование зависимости выходного напряжения усилительного каскада от амплитуды и частоты входного сигнала.
41667. Поиск данных и фильтрации в среде программирования Borland C++Builder 288.5 KB
  Компонент Table позволяет не только отображать, редактировать и упорядочивать данные, но и отсортировывать записи по определенным критериям. Фильтрация может задаваться свойствами Filter,Filtered и FilterOptions компонента Table.