readsonar.m

实现了fisher线性分类器的功能,源码使用matlab编程,通俗易懂

%%%%%%%%%%%%%%%%%%%%%%%%%%%% 读Sonar样本数据 %%%%%%%%%%%%%%%%%%%%%%%%%%%%
function [TrainSet,TestSet]=ReadSonar(train,test)
% 1.读训练样本
fp=fopen(train,'rt');
if(fp==-1) 
    display('File not exist!'); 
    return; 
end
%计算训练集样本记录个数
TrainSet.n=0;%训练集样本的总数
c=fscanf(fp,'%c',1);
while ~feof(fp)  
    if c=='M' | c=='R'
        TrainSet.n=TrainSet.n+1;  
    end
    c=fscanf(fp,'%c',1);
end
fseek(fp,0,'bof');
X=[1:60];
TrainSet.X=[];
TrainSet.CN=[];
TrainSet.n1=0;
TrainSet.n2=0;
for j=1:TrainSet.n
    for i=1:60
        X(i)=fscanf(fp,'%f');
        fscanf(fp,'%c',1);
    end
    c=fscanf(fp,'%c',1);
    if c=='M'
        cn=1;
        TrainSet.n1=TrainSet.n1+1;%计算M类的样本记录个数
    else
        cn=-1;
        TrainSet.n2=TrainSet.n2+1;%计算M类的样本记录个数
    end
    TrainSet.X=[TrainSet.X;X];%样本记录矩阵
    TrainSet.CN=[TrainSet.CN;cn];%每个样本记录对应的类别矩阵
end
fclose(fp);
TrainSet.tn=60;%维数

% 2.读测试样本
fp=fopen(test,'rt');
if(fp==-1) 
    display('File not exist!'); 
    return; 
end
TestSet.n=0;
c=fscanf(fp,'%c',1);
while ~feof(fp)  
    if c=='M' | c=='R'
        TestSet.n=TestSet.n+1;  
    end
    c=fscanf(fp,'%c',1);
end
fseek(fp,0,'bof');
X=[1:60];
TestSet.X=[];
TestSet.CN=[];
TestSet.n1=0;
TestSet.n2=0;
for j=1:TestSet.n
    for i=1:60
        X(i)=fscanf(fp,'%f');
        fscanf(fp,'%c',1);
    end
    c=fscanf(fp,'%c',1);
    if c=='M'
        TestSet.n1=TestSet.n1+1;
        cn=1;
    else
        cn=-1;
        TestSet.n2=TestSet.n2+1;
    end
    TestSet.X=[TestSet.X;X];
    TestSet.CN=[TestSet.CN;cn];
end
fclose(fp);
TestSet.tn=60;