Не вижу дельта-пика в плотности распределении функции ksdensity

@SpaceQuester · Регистрация: 19.01.2016

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

Всем привет!

Есть такой работающий модельный код:

Matlab M

omega_cur = 0;
omega_total = [];
for i = 1:36
  omega_cur = 0.05;
  omega_total = [omega_total omega_cur'];
end
start_density_base = 0;
d_density_base     = 0.000001;
end_density_base   = 0.15;
density_base = (start_density_base : d_density_base : end_density_base)';
[f] = ksdensity(omega_total, density_base);
f = f / max(f);
plot(density_base, f, 'r', 'LineWidth', 2);
xlim([0 0.15]);
ylim([0 1.01]);

В теории здесь должен получиться дельта пик на частоте 0.05. Но его нет. Подскажите, в чем проблема.

@Krasme · 13.07.2018, 13:37

Matlab M

1	ylim([0.999 1.001]);

@SpaceQuester · 13.07.2018, 13:54 **[ТС]**

Это всё равно не дельта пик

Нужен чистый дельта-пик.

@Krasme · 13.07.2018, 14:01

Сообщение от SpaceQuester

чистый дельта-пик.

что это значит?

@SpaceQuester · 13.07.2018, 14:23 **[ТС]**

Вот как надо:

Matlab M

omega_cur = 0;
omega_total = [];
for i = 1:36
  omega_cur = 0.05;
  omega_total = [omega_total omega_cur'];
end
start_density_base = 0;
d_density_base     = 0.000001;
end_density_base   = 0.15;
density_base = (start_density_base : d_density_base : end_density_base)';
[f] = ksdensity(omega_total, density_base, 'Bandwidth', 0.0001);
f = f / max(f);
plot(density_base, f, 'r', 'LineWidth', 2);
xlim([0 0.15]);
ylim([0 1.01]);

Ответ заключался в том что бы задать принудительно полосу фильтрации:

Matlab M

1	'Bandwidth', 0.0001

Я думал функция автоматом сама поймёт что мне нужно. Но MATLAB оказался строг со мной

@SpaceQuester 2 / 2 / 0 Регистрация: 19.01.2016 Сообщений: 223
	13.07.2018, 13:54 [ТС]	3
	Это всё равно не дельта пик Нужен чистый дельта-пик. 0

@Krasme 6829 / 4889 / 2065 Регистрация: 02.02.2014 Сообщений: 13,046
	13.07.2018, 14:01	4
	Сообщение от SpaceQuester чистый дельта-пик. что это значит? 0

	13.07.2018, 14:23

Опции темы