Граница решения рисования двух многомерных гауссовских

Я пытаюсь построить что-то похожее на следующее: введите здесь описание изображения

Я использую Матлаб. Я добился рисования контурных участков. Однако мне не удалось нарисовать дискриминант. Может ли кто-нибудь показать пример кода Matlab или дать какую-нибудь идею, чтобы нарисовать дискриминант?


person ipman    schedule 24.10.2013    source источник


Ответы (2)


Если вы знаете функцию плотности вероятности каждого из гауссов для данной точки (x,y), скажем, ее pdf1(x,y) и pdf2(x,y), то вы можете просто построить контурную линию f(x,y) := pdf1(x,y) > pdf2(x,y). Таким образом, вы определяете функцию f как 1 тогда и только тогда, когда pdf1(x,y)>pdf2(x,y). Таким образом, единственный контур будет размещен вдоль кривой, где pdf1(x,y)==pdf2(x,y) является границей решения (дискриминантом). Если вы хотите определить «хорошую» функцию, вы можете сделать это, просто установив f(x,y) = sgn( pdf1(x,y) - pdf2(x,y) ), и построение ее контурного графика приведет к точно такому же дискриминанту.

person lejlot    schedule 24.10.2013

Вот как я решил бы эту задачу аналитически: вы приравниваете эти две дискриминантные функции

g1(x)=x' W1 x + w1' x + w10
g2(x)=x' W2 x + w2' x + w20
g1(x) = g2(x)

==> x' (W2 - W1) x + (w2-w1)'x + w20 - w10

тогда я считаю, что W2 - W1 должна быть этой матрицей

W2-W1 = [a b; c d]

который затем, расширяя вектор x=[x1 x2]', мы получаем:

a x1^2 + (b+c) x1 x2 + d x2^2 + (w21-w11) x1 + (w22-w12) x2 + w20-w10 = 0

это уравнение эллипса, поэтому вы можете упростить его до следующей формы:

(x1 - a0)^2/h + (x2-b0)^2/g = r^2

Или вы можете предположить, что вы знаете диапазон x1, например, x1=[-2:0.1:2], а затем решить параболу

пример параболической границы решения

person Vahid Mirjalili    schedule 21.04.2014