模式识别及MATLAB实现【电子工业出版社】学习笔记

时肆知还

于 2023-12-13 22:50:44 发布

阅读量1k

点赞数 20

分类专栏：模式识别文章标签： matlab 学习笔记

本文链接：https://blog.csdn.net/H551493l/article/details/134983154

版权

模式识别专栏收录该内容

2 篇文章 0 订阅

订阅专栏

一、实验目的

1. 学习、理解基于最小错误率贝叶斯决策的基本原理；

2. 掌握贝叶斯决策 MATLAB 算法。

二、实验仪器

安装有 MATLAB 软件的计算机一台。

三、实验内容

1. 学习、理解基于最小错误率贝叶斯决策的基本原理；

2. 设计贝叶斯决策 MATLAB 算法，对一定数据进行分析，验证算法的正确性；

3. 对实验结果进行总结分析

四。实验代码

clear
iris=xlsread('C:\Users\admin\Desktop\iris.csv');
global meanx2 w1 w2 w3 meanx1 meanx3 test;
N=40;
 
for i=1:N
    for j=1:4
        w1(i,j)=iris(i,j+1);
    end
end
sumx1=sum(w1,1);
for i=1:4
    meanx1(1,i)=sumx1(1,i)/N;
end
 
for i=1:N
    for j=1:4
        w2(i,j)=iris(i+50,j+1);
    end
end
sumx2=sum(w2,1);
for i=1:4
    meanx2(1,i)=sumx2(1,i)/N;
end
 
for i=1:N
    for j=1:4
        w3(i,j)=iris(i+100,j+1);
    end
end
sumx3=sum(w3,1);
for i=1:4
    meanx3(1,i)=sumx3(1,i)/N;
end
 
z1(4,4)=0;
var1(4,4)=0;
for i=1:4
    for j=1:4
        for k=1:N
            z1(i,j)=z1(i,j)+(w1(k,i)-meanx1(1,i))*(w1(k,j)-meanx1(1,j));
        end
        var1(i,j)=z1(i,j)/(N-1);
    end
end
 
z2(4,4)=0;
var2(4,4)=0;
for i=1:4
    for j=1:4
        for k=1:N
            z2(i,j)=z2(i,j)+(w2(k,i)-meanx2(1,i))*(w2(k,j)-meanx2(1,j));
        end
        var2(i,j)=z2(i,j)/(N-1);
    end
end
 
z3(4,4)=0;
var3(4,4)=0;
for i=1:4
    for j=1:4
        for k=1:N
            z3(i,j)=z3(i,j)+(w3(k,i)-meanx3(1,i))*(w3(k,j)-meanx3(1,j));
        end
        var3(i,j)=z3(i,j)/(N-1);
    end
end
  Var1_inv=[];
  Var1_det=[];
  Var2_inv=[];
  Var2_det=[];
  Var3_inv=[];
  Var3_det=[];
 
var1_inv=inv(var1);
var2_inv=inv(var2);
var3_inv=inv(var3);
 
var1_det=det(var1);
var2_det=det(var2);
var3_det=det(var3);
 
M=10;
for i=1:M
    for j=1:4
        test(i,j)=iris(i+50,j+1);
    end
end
t1=0;
t2=0;
t3=0;
for i=1:M 
    x=test(i,1);
    y=test(i,2);
    z=test(i,3);
    h=test(i,4);
 
    g1=(-0.5)*([x,y,z,h]-meanx1)*var1_inv*([x,y,z,h]'-meanx1')-0.5*log(abs(var1_det))+log(w1);
    g2=(-0.5)*([x,y,z,h]-meanx2)*var2_inv*([x,y,z,h]'-meanx2')-0.5*log(abs(var2_det))+log(w2);
 
    if g1>g2
        t1=t1+1;
    else
        t2=t2+1;
    end
end 
plot3(w1,w2,w3)