matlab例程

开发平台：

Matlab

feature_x1_x2.m：源码内容

							clear all;
close all;
x1w1 = [-5.01 -5.43 1.08 0.86 -2.67 4.94 -2.51 -2.25 5.56 1.03];
x1w2 = [-0.91 1.3 -7.75 -5.47 6.14 3.6 5.37 7.18 -7.39 -7.5];
x2w1 = [-8.12 -3.48 -5.52 -3.78 0.63 3.29 2.09 -2.13 2.86 -3.33];
x2w2 = [-0.18 -2.06 -4.54 0.5 5.72 1.26 -4.63 1.46 1.17 -6.32];
figure,plot(x1w1,x2w1,'o',x1w2,x2w2,'r*');
mean11 = mean(x1w1);
mean21 = mean(x2w1);
sigma11 = std(x1w1);
sigma21 = std(x2w1);
sigma_x1_x2_w1 = sum((x1w1-mean11).*(x2w1-mean21))/length(x2w1);
mean12 = mean(x1w2);
sigma12 = std(x1w2);
mean22 = mean(x2w2);
sigma22 = std(x2w2);
sigma_x1_x2_w2 = sum((x1w2-mean12).*(x2w2-mean22))/length(x2w2);
%p(x1,x2|w1) and p(x1,x2|w2)
x = -10:0.1:10;
y = -10:0.1:10;
px1x2w1 = g2(mean11,sigma11,mean21,sigma21,sigma_x1_x2_w1,x,y);
px1x2w2 = g2(mean12,sigma12,mean22,sigma22,sigma_x1_x2_w2,x,y);
%Decision boundary
pw1 = 0.5;
pw2 = 0.5;
d = bi_dec2(px1x2w1,px1x2w2,pw1,pw2,x,y);
figure,imagesc(px1x2w1); %if you have good graphics card, please use mesh(x,y,px1x2w1);
axis image;
figure,imagesc(px1x2w2); %if you have good graphics card, please use mesh(x,y,px1x2w2);
axis image;
figure,imagesc(d);
axis image;

						