Построение графика интеграла вероятности

@Шеда · Регистрация: 13.12.2013

Студворк — интернет-сервис помощи студентам

Здравствуйте!
Помогите, пожалуйста, с написанием программы (в MatLab) построения графика зависимости P(q).
Есть вот такая функция:

P=1-Ф[Ф^(-1)(1-t)-q]
Ф(z)=(1/sqrt(2*Pi))*int(exp(-z^2/2),’z’,-inf,z);
t=10^(-4);
q=0:0.1:16;

Ф^(-1) - функция, обратная к интегралу.

Заранее спасибо))

Зосима · 13.12.2013, 18:21

Сообщение от Шеда

Ф(z)=(1/sqrt(2*Pi))*int(exp(-z^2/2),’z’,-inf,z);

эта часть весьма напоминает нормальное распределение, с сигма=1, поэтому если вспомнить правило 3х сигм, то пожалуй можно записать:
F = @(z) 1/sqrt(2*pi))*quad('exp(-x.^2/2)',-100,z)
В ообщем на костылях получается так:

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
clear, clc
 
t=10^(-4);
q=0:0.1:16;
 
F = @(z)  1/sqrt(2*pi)*quad('exp(-x.^2/2)',-100,z);
z = 0:0.001:1; 
for i = 1:length(z) % считаем интергал для всех значений
    f(i) = F(z(i));
end
invF = @(x) spline(f,z, x);
invFt = invF(1-t);
for i = 1:length(q)
    P(i) = 1 - F( invFt - q(i) );
end
plot(q,P)
grid on
ylim([0 1.1])

*увы! интеграл не переваривает сразу массив пределов, поэтому пришлось делать циклы.

@Шеда · 15.12.2013, 12:34 **[ТС]**

Большое спасибо)))

@Шеда · 18.12.2013, 20:15 **[ТС]**

Здравствуйте! Подскажите, пожалуйста, по еще одной программке!
Что тут не правильно написано?

Matlab M
1
2
3
4
5
6
7
8
9
10
clear, clc
 
t=10^(-10);
q=0:0.1:16;
h=q*sqrt(2*log(1/t));
 
F = @(h) quad('(y/q^2*exp(-(y^2+q^2)/2*q^2)*besselj(0,y/q)',1000,h);
 
plot(q,F)
grid on

Зосима · 19.12.2013, 11:06

Первым делом повторюсь, что интеграл не умеет брать от переменных пределов (h - это массив), поэтому нужно организовать цикл. А во вторых внутри подинтегрального выражения неясно по какой переменной идет интегрирование и значения q он при такой записи не подставляет!
вот, немного причесал, но проверь еще раз формулу:

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
clear, clc
 
t=10^(-10);
q=0:0.1:16;
h=q*sqrt(2*log(1/t));
 
for i = 1:length(h)
    f = @(y) y./q(i).^2.*exp(-(y.^2+q(i)^2)./2.*q(i)^2).*besselj(0,y/q(i));
    F(i) = quad(f,1000,h(i));
end
plot(q,F)
grid on

Новые блоги и статьи Все статьи Все блоги /
Новый ноутбук volvo 07.12.2025 Всем привет. По скидке в "черную пятницу" взял себе новый ноутбук Lenovo ThinkBook 16 G7 на Амазоне: Ryzen 5 7533HS 64 Gb DDR5 1Tb NVMe 16" Full HD Display Win11 Pro	Музыка, написанная Искусственным Интеллектом volvo 04.12.2025 Всем привет. Некоторое время назад меня заинтересовало, что уже умеет ИИ в плане написания музыки для песен, и, собственно, исполнения этих самых песен. Стихов у нас много, уже вышли 4 книги, еще 3. . .	От async/await к виртуальным потокам в Python IndentationError 23.11.2025 Армин Ронахер поставил под сомнение async/ await. Создатель Flask заявляет: цветные функции - провал, виртуальные потоки - решение. Не threading-динозавры, а новое поколение лёгких потоков. Откат?. . .	Поиск "дружественных имён" СОМ портов Argus19 22.11.2025 Поиск "дружественных имён" СОМ портов На странице: https:/ / norseev. ru/ 2018/ 01/ 04/ comportlist_windows/ нашёл схожую тему. Там приведён код на С++, который показывает только имена СОМ портов, типа,. . .	Сколько Государство потратило денег на меня, обеспечивая инсулином. Programma_Boinc 20.11.2025 Сколько Государство потратило денег на меня, обеспечивая инсулином. Вот решила сделать интересный приблизительный подсчет, сколько государство потратило на меня денег на покупку инсулинов. . . .
Ломающие изменения в C#.NStar Alpha Etyuhibosecyu 20.11.2025 Уже можно не только тестировать, но и пользоваться C#. NStar - писать оконные приложения, содержащие надписи, кнопки, текстовые поля и даже изображения, например, моя игра "Три в ряд" написана на этом. . .	Мысли в слух kumehtar 18.11.2025 Кстати, совсем недавно имел разговор на тему медитаций с людьми. И обнаружил, что они вообще не понимают что такое медитация и зачем она нужна. Самые базовые вещи. Для них это - когда просто люди. . .	Создание Single Page Application на фреймах krapotkin 16.11.2025 Статья исключительно для начинающих. Подходы оригинальностью не блещут. В век Веб все очень привыкли к дизайну Single-Page-Application . Быстренько разберем подход "на фреймах". Мы делаем одну. . .	Фото: Daniel Greenwood kumehtar 13.11.2025	Расскажи мне о Мире, бродяга kumehtar 12.11.2025 — Расскажи мне о Мире, бродяга, Ты же видел моря и метели. Как сменялись короны и стяги, Как эпохи стрелою летели. - Этот мир — это крылья и горы, Снег и пламя, любовь и тревоги, И бескрайние. . .

@Шеда 0 / 0 / 0 Регистрация: 13.12.2013 Сообщений: 3

	Построение графика интеграла вероятности 13.12.2013, 17:29. Показов 6380. Ответов 4 Метки нет (Все метки) Здравствуйте! Помогите, пожалуйста, с написанием программы (в MatLab) построения графика зависимости P(q). Есть вот такая функция: P=1-Ф[Ф^(-1)(1-t)-q] Ф(z)=(1/sqrt(2Pi))int(exp(-z^2/2),’z’,-inf,z); t=10^(-4); q=0:0.1:16; Ф^(-1) - функция, обратная к интегралу. Заранее спасибо)) 0

@Шеда 0 / 0 / 0 Регистрация: 13.12.2013 Сообщений: 3
	15.12.2013, 12:34 [ТС]
	Большое спасибо))) 0