function[Sww]=fisher2D（xy）
load?data.mat
%%%將樣本點和測試點畫出;
plot（Sample1（:1）Sample1（:2）‘g*‘）
hold?on;
plot（Sample2（:1）Sample2（:2）‘b.‘）
plot（xy‘r*‘）;
%先求樣本的平均值
u1=mean（Sample1）;
u2=mean（Sample2）;
%求兩樣本的內間方差s1s2；
s1=[00;00];
for?i=1:1000;
????s1=s1+（Sample1（i:）-u1）‘*（Sample1（i:）-u1）;
end

s2=[00;00];
for?i=1:1000;
????s2=s2+（Sample2（i:）-u2）‘*（Sample2（i:）-u2）;
end
%求兩樣本的內間方差;
Sw=s1+s2;
%求兩樣本的投影向量;
w=inv（Sw）*（u1-u2）‘;
%設置以固定閾值；
for?i=1:1000;
x1（i）=Sample1（i:）*w;
x2（i）=Sample2（i:）*w;
end;
w0=（mean（x1）+mean（x2））/2;
if（[xy]*w-w0>0）
????disp（‘該點屬于sample1類.‘）
end;
if（[xy]*w-w0<0）
???disp（‘該點屬于Sample2類.‘）
end;
%m=-10*w（1:）:10*w（1:）;
%n=-10*w（2:）:10*w（2:）;
line（[-1000*w（1）1000*w（1）][-1000*w（2）1000*w（2）]）;
hold?off

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件???????30712??2011-04-15?10:57??fisher\data.mat
?????文件?????????701??2011-04-20?22:28??fisher\fisher2D.asv
?????文件?????????851??2011-04-20?22:45??fisher\fisher2D.m
?????文件??????????71??2011-04-15?18:59??fisher\main.m
?????目錄???????????0??2011-05-04?20:58??fisher\