3.4 Результати виконання лабораторної роботи

clear();

r = -0.9;

Mx = 10.0;

My = 20.0;

sigmaX = 0.6;

sigmaY = 0.8;

R1 = -0.6;

R2 = -0.7;

X = [11.8 11.2 10.0 8.8];

Ygr = My - sigmaX;

A = 0;

%_______________3.3_________________

P_K2 = func_erf((Ygr-My)/sigmaY);

P_K1 = 1 - P_K2;

%_______________3.4_________________

M_x_K1 = Mx + ( ( sigmaX*r*exp((-(Ygr-My)^2)/(2*sigmaY^2)))/((P_K1)*sqrt(2*pi)) );

M_x_K2 = Mx - ( ( sigmaX*r*exp((-(Ygr-My)^2)/(2*sigmaY^2)))/((P_K2)*sqrt(2*pi)) );

%_______________3.5_________________

Xa = [Mx-3*sigmaX Mx-2.5*sigmaX Mx-2*sigmaX Mx-1.5*sigmaX Mx-sigmaX Mx-0.8*sigmaX Mx-0.6*sigmaX Mx-0.4*sigmaX Mx-0.2*sigmaX Mx Mx+0.2*sigmaX Mx+0.4*sigmaX Mx+0.6*sigmaX Mx+0.8*sigmaX Mx+sigmaX Mx+1.5*sigmaX Mx+2*sigmaX Mx+2.5*sigmaX Mx+3*sigmaX];

for i = 1:length(Xa)

W_x(i) = (1/(sigmaX*(sqrt(2*pi))))*exp((-((Xa(i)-Mx)^2)/(2*(sigmaX^2))));

end

%_______________3.6_________________

Xa2K1 = [M_x_K1-3*sigmaX M_x_K1-2.5*sigmaX M_x_K1-2*sigmaX M_x_K1-1.5*sigmaX M_x_K1-sigmaX M_x_K1-0.8*sigmaX M_x_K1-0.6*sigmaX M_x_K1-0.4*sigmaX M_x_K1-0.2*sigmaX M_x_K1 M_x_K1+0.2*sigmaX M_x_K1+0.4*sigmaX M_x_K1+0.6*sigmaX M_x_K1+0.8*sigmaX M_x_K1+sigmaX M_x_K1+1.5*sigmaX M_x_K1+2*sigmaX M_x_K1+2.5*sigmaX M_x_K1+3*sigmaX];

Xa2K2 = [M_x_K2-3*sigmaX M_x_K1-2.5*sigmaX M_x_K2-2*sigmaX M_x_K1-1.5*sigmaX M_x_K2-sigmaX M_x_K2-0.8*sigmaX M_x_K2-0.6*sigmaX M_x_K2-0.4*sigmaX M_x_K2-0.2*sigmaX M_x_K2 M_x_K2+0.2*sigmaX M_x_K2+0.4*sigmaX M_x_K2+0.6*sigmaX M_x_K2+0.8*sigmaX M_x_K2+sigmaX M_x_K2+1.5*sigmaX M_x_K2+2*sigmaX M_x_K2+2.5*sigmaX M_x_K2+3*sigmaX];

for i = 1:length(Xa2K1)

W_x_K1_a(i) = (exp((-((Xa2K1(i)-Mx)^2))/(2*sigmaX^2))*(1-func_erf((sigmaX*(Ygr-My)-r*sigmaY*(Xa2K1(i)- Mx)) /

(sigmaX*sigmaY*sqrt(1-r^2))))) / (P_K1*sigmaX*sqrt(2*pi));

W_x_K2_a(i) = (exp((-((Xa2K2(i)-Mx)^2))/(2*sigmaX^2))*(func_erf((sigmaX*(Ygr-My)-r*sigmaY*(Xa2K2(i)-

Mx))/(sigmaX*sigmaY*sqrt(1-r^2))))) / (P_K2*sigmaX*sqrt(2*pi));

end

%_______________3.7_________________

Xkl = Mx + ((sigmaX*(Ygr-My)-A*sigmaX*sigmaY*sqrt(1-r^2))/(r*sigmaY));

%_______________3.8_________________

P_reshK1 = func_erf( (Xkl-Mx)/sigmaX );

P_reshK2 = 1-P_reshK1;

%_______________3.9__________________

M_y_reshK1 = My - (((sigmaY*r)/(P_reshK1*sqrt(2*pi))) * exp(-((Xkl-Mx)^2/(2*sigmaX^2))));

M_y_reshK2 = My + (((sigmaY*r)/(P_reshK2*sqrt(2*pi))) * exp(-((Xkl-Mx)^2/(2*sigmaX^2))));

%_______________3.10_________________

Ya = [My-3*sigmaY My-2.5*sigmaY My-2*sigmaY My-1.5*sigmaY My-sigmaY My-0.8*sigmaY My-0.6*sigmaY My-0.4*sigmaY My-0.2*sigmaY My My+0.2*sigmaY My+0.4*sigmaY My+0.6*sigmaY My+0.8*sigmaY My+sigmaY My+1.5*sigmaY My+2*sigmaY My+2.5*sigmaY My+3*sigmaY];

for i = 1:length(Ya)

W_y(i) = (1/(sigmaY*(sqrt(2*pi))))*exp((-((Ya(i)-My)^2)/(2*(sigmaY^2))));

end

%_______________3.11_________________

Ya2K1 = [M_y_reshK1-3*sigmaX M_y_reshK1-2.5*sigmaX M_y_reshK1-2*sigmaX M_y_reshK1-1.5*sigmaX M_y_reshK1-sigmaX M_y_reshK1-0.8*sigmaX M_y_reshK1-0.6*sigmaX M_y_reshK1-0.4*sigmaX M_y_reshK1-0.2*sigmaX M_y_reshK1 M_y_reshK1+0.2*sigmaX M_y_reshK1+0.4*sigmaX M_y_reshK1+0.6*sigmaX M_y_reshK1+0.8*sigmaX M_y_reshK1+sigmaX M_y_reshK1+1.5*sigmaX M_y_reshK1+2*sigmaX M_y_reshK1+2.5*sigmaX M_y_reshK1+3*sigmaX];

Ya2K2 = [M_y_reshK2-3*sigmaX M_y_reshK2-2.5*sigmaX M_y_reshK2-2*sigmaX M_y_reshK2-1.5*sigmaX M_y_reshK2-sigmaX M_y_reshK2-0.8*sigmaX M_y_reshK2-0.6*sigmaX M_y_reshK2-0.4*sigmaX M_y_reshK2-0.2*sigmaX M_y_reshK2 M_y_reshK2+0.2*sigmaX M_y_reshK2+0.4*sigmaX M_y_reshK2+0.6*sigmaX M_y_reshK2+0.8*sigmaX M_y_reshK2+sigmaX M_y_reshK2+1.5*sigmaX M_y_reshK2+2*sigmaX M_y_reshK2+2.5*sigmaX M_y_reshK2+3*sigmaX];

for i = 1:length(Ya2K1)

W_y_reshK1(i) = (exp(-(Ya2K1(i)-My)^2/(2*sigmaY^2))*(func_erf((sigmaY*(Xkl-Mx)-r*sigmaX*(Ya2K1(i)-

My))/(sigmaX*sigmaY*sqrt(1-r^2)))))/(P_reshK1*sigmaY*sqrt(2*pi));

W_y_reshK2(i) = (exp(-(Ya2K2(i)-My)^2/(2*sigmaY^2))*(1-func_erf((sigmaY*(Xkl-Mx)-r*sigmaX*(Ya2K2(i)-

My))/(sigmaX*sigmaY*sqrt(1-r^2)))))/(P_reshK2*sigmaY*sqrt(2*pi));

end

%_______________3.12_________________

for i = 1:length(X)

W_x_K1(i) = (exp((-(X(i)-Mx)^2)/(2*sigmaX^2))*(1-func_erf((sigmaX*(Ygr-My)-r*sigmaY*(X(i)-Mx))/

(sigmaX*sigmaY*sqrt(1-r^2))))) / (P_K1*sigmaX*sqrt(2*pi));

W_x_K2(i) = (exp((-(X(i)-Mx)^2)/(2*sigmaX^2))*(1-func_erf((sigmaX*(Ygr-My)-r*sigmaY*(X(i)-Mx))/

(sigmaX*sigmaY*sqrt(1-r^2))))) / (P_K2*sigmaX*sqrt(2*pi));

Lambda(i) = W_x_K1(i)/W_x_K2(i);

end

%_______________3.13_________________

Porog = P_K2/P_K1;

%_______________3.15_________________

ra = [r 0.9*r 0.8*r 0.7*r 0.6*r 0.5*r 0.4*r 0.3*r 0.25*r 0.2*r 0.15*r ];

for i = 1:length(ra)

MR(i) = (sigmaY*ra(i)*exp(-((Xkl-Mx)^2/(ra(i)*sigmaX^2))))/(sqrt(2*pi)*P_reshK1*P_reshK2);

end

plot(MR,'o-');

%_______________3.16_________________

figure(2);

plot(Ya,W_y,'*-',Ya,W_y_reshK1,'o-',Ya,W_y_reshK2,'x-');

figure(3);

plot(Xa,W_x,'*-',Xa,W_x_K1_a,'o-',Xa,W_x_K2_a,'x-');

Рисунок 3.1 – Зв'язок модуля різниці умовних математичних сподівань зі ступенем тісноти

Рисунок 3.2 – Графік умовних та безумовної густин прогнозованого параметра

Рисунок 3.3 – Графік умовних та безумовної густин ознаки

W(x) =

0,00738641401989665 0,0292138341559476 0,0899849441886470 0,215862659443153

0,403284540865239 0,482819254602471 0,555374338153000 0,613783567172206

0,65173783292427 0,664903800669055 0,651737823292427 0,613783567172206

0,555374338153000 0,482819254602471 0,403284540865239 0,215862659443153

0,0899849441886470 0,0292138341559476 0,00738641401989665

W(x | K1) =

0,00313898543181332 0,0147925187389076 0,0542901867547479 0,155176844437108

0,345427949359645 0,443569762952807 0,547219300886291 0,648405284513300

0,737305822330868 0,802671795842526 0,832121788091824 0,813352416341321

0,738351845605966 0,610761753356202 0,451032256068835 0,114073831293116

0,0101378320018837 0,000280450602414319 2,27220807323020e-06

W(x | K2) =

1,48770856902627e-08 7,11758335330000e-16 0,000767153626010186 7,94909859916593e-09

0,270731239904521 0,497014320780731 0,766423241581132 1,00287728579928

1,12824998259098 1,10927820786343 0,971563326723257 0,773732089732803

0,571204174961792 0,397048947805725 0,262604051379422 0,0774855632369529

0,0177644828728965 0,00317164876958062 0,000441005243588190

W(y) =

0,00553981051492254 0,0219103756169607 0,0674887081414849 0,161896994582365

0,302463405648929 0,362114440951853 0,416530753614749 0,460337675379154

0,488803367469320 0,498677850501791 0,488803367469320 0,460337675379154

0,416530753614749 0,362114440951853 0,302463405648929 0,161896994582365

0,0674887081414849 0,0219103756169607 0,00553981051492254

W(y | решК1) =

0,00182758922455436 0,0179190984194965 0,0928269689313536 0,269534261623088

0,481863432485296 0,545650066257817 0,588413295588844 0,609037761224878

0,609034462788249 0,591293659550182 0,559178811726507 0,516095886515644

0,465360807092263 0,410147155049695 0,353401766217780 0,220695298490064

0,119739347157041 0,0564427048527591 0,0231156486738770.

W(y | решК2) =

0,00531641893884056 0,0184393406827766 0,0555605357550579 0,145277064316184

0,326592316601379 0,429705567095256 0,545887083685189 0,665103076031656

0,771022144583932 0,843118825916291 0,862205831494812 0,817956093136538

0,714709450355509 0,571661550359519 0,416416912781678 0,122066603969479

0,0185711619151053 0,00142597349646310 5,43294917960852e-05