Представляйте Распределение Коши Используя t Шкалу Местоположения

Этот пример показывает, как использовать t объект распределения вероятностей шкалы местоположения работать с распределением Коши с нестандартными значениями параметров.

Шаг 1. Создайте объект распределения вероятностей.

Создайте a t распределение вероятностей шкалы местоположения возражает со степенями свободы nu = 1. Задайте mu = 3, чтобы установить параметр положения, равный 3, и sigma = 1 устанавливать масштабный коэффициент, равный 1.

pd = makedist('tLocationScale','mu',3,'sigma',1,'nu',1)
pd = 
  tLocationScaleDistribution

  t Location-Scale distribution
       mu = 3
    sigma = 1
       nu = 1

Шаг 2. Вычислите описательную статистику.

Используйте объектные функции, чтобы вычислить описательную статистику для распределения Коши.

med = median(pd)
med = 3
r = iqr(pd)
r = 2
m = mean(pd)
m = NaN
s = std(pd)
s = Inf

Медиана распределения Коши равна своему параметру положения, и межквартильный размах равен два раза его масштабному коэффициенту. Его среднее и стандартное отклонение не определено.

Шаг 3. Вычислите и постройте PDF.

Вычислите и постройте PDF распределения Коши.

x = -20:1:20;
y = pdf(pd,x);
plot(x,y,'LineWidth',2)

Пик PDF сосредоточен в параметре положения mu = 3.

Шаг 4. Сгенерируйте вектор случайных чисел Коши.

Сгенерируйте вектор-столбец, содержащий 10 случайных чисел от распределения Коши с помощью функции random для t объект распределения вероятностей шкалы местоположения.

rng('default');  % For reproducibility
r = random(pd,10,1)
r = 10×1

    3.2678
    4.6547
    2.0604
    4.7322
    3.1810
    1.6649
    1.8471
    4.2466
    5.4647
    8.8874

Шаг 5. Сгенерируйте матрицу случайных чисел Коши.

Сгенерируйте матрицу 5 на 5 случайных чисел Коши.

r = random(pd,5,5)
r = 5×5

    2.2867    2.9692   -1.7003    5.5949    1.9806
    2.7421    2.7180    3.2210    2.4233    3.1394
    3.5966    3.9806    1.0182    6.4180    5.1367
    5.4791   15.6472    0.7558    2.8908    5.9031
    1.6863    4.0985    2.9934   13.9506    4.8792

Смотрите также

Похожие темы

Для просмотра документации необходимо авторизоваться на сайте