- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
                        查看更多
                        
                    
                基于matlab的字符识别研究推荐
                    
PAGE  
PAGE  27
?
?
function Bayes2
%算法视线见《模式识别》P33-P44(各类样本的协方差不相等)
%为了提高实验样本测试的精度,故采用多次模拟求平均值的方法
N=input(实验模拟次数 N(N最好为奇数)= );
Result(1:3,1:3)=0;      %判别矩阵的初始化
for k=1:N             %控制程序模拟次数N
    %生成二维正态分布的样本2 X N 维的矩阵 
    X1=mvnrnd([1 2],[4 0;0 6],300);   %2 X N
    X2=mvnrnd([5 3],[5 0;0 1],200);
    X3=mvnrnd([4 7],[2 0;0 9],500);   %样本程序
    %%
    %测试样本
    X10=mvnrnd([1 2],[4 0;0 6],100);   %2 X N
    X20=mvnrnd([5 3],[5 0;0 1],100);
    X30=mvnrnd([4 7],[2 0;0 9],100);  
    %先验概率
    P(1)=length(X1)/(length(X1)+length(X2)+length(X3));
    P(2)=length(X2)/(length(X1)+length(X2)+length(X3));
    P(3)=length(X3)/(length(X1)+length(X2)+length(X3));
    %计算相关量  cov(X):协方差矩阵 Ave:均值
    %%
    W1=-1/2*inv(cov(X1)); W2=-1/2*inv(cov(X2)); W3=-1/2*inv(cov(X3));%
    Ave1=(sum(X1)/length(X1));Ave2=(sum(X2)/length(X2));
    Ave3=(sum(X3)/length(X3));%计算平均值(2维列向量)
    w1=inv(cov(X1))*Ave1;w2=inv(cov(X2))*Ave2;w3=inv(cov(X3))*Ave3;%2
    w10=-1/2*Ave1*inv(cov(X1))*Ave1-1/2*log(det(cov(X1)))+log(P(1));
    w20=-1/2*Ave2*inv(cov(X2))*Ave2-1/2*log(det(cov(X2)))+log(P(2));
    w30=-1/2*Ave3*inv(cov(X3))*Ave3-1/2*log(det(cov(X3)))+log(P(3));
    %%
    for i=1:3                                     
        for j=1:100                               
            if i==1
                g1=X10(:,j)*W1*X10(:,j)+w1*X10(:,j)+w10;  
                g2=X10(:,j)*W2*X10(:,j)+w2*X10(:,j)+w20;
                g3=X10(:,j)*W3*X10(:,j)+w3*X10(:,j)+w30;
                if g1=g2g1=g3     
                    Result(1,1)=Result(1,1)+1;
                elseif g2=g1g2=g3
                    Result(1,2)=Result(1,2)+1;%记录误判情况
                else
                    Result(1,3)=Result(1,3)+1;%记录误判情况
                end
            elseif i==2
                g1=X20(:,j)*W1*X20(:,j)+w1*X20(:,j)+w10;
                g2=X20(:,j)*W2*X20(:,j)+w2*X20(:,j)+w20;
                g3=X20(:,j)*W3*X20(:,j)+w3*X20(:,j)+w30;
                if g2=g1g2=g3
                    Result(2,2)=Result(2,2)+1;
       
                您可能关注的文档
- 图形的变化_七巧板 课件 七年级数学上册 初一课件.ppt
- 图形图象制作毕业论文 车厢广告设计.doc
- 图形编辑(三)ppt课件.ppt
- 图片基本处理方法ppt课件初中信息技术苏科版七年级全一册_1.ppt
- 图片基本处理方法ppt课件初中信息技术苏科版七年级全一册_2.ppt
- 图片基本处理方法ppt课件初中信息技术苏科版七年级全一册_3.ppt
- 图片基本处理方法ppt课件初中信息技术苏科课标版七年级全一册课件.ppt
- 图片基本处理方法ppt课件初中信息技术苏科课标版七年级全一册课件_1.ppt
- 图书购物系统毕业设计精选.doc
- 图片基本处理方法ppt课件初中信息技术苏科课标版七年级全一册课件_2.ppt
 原创力文档
原创力文档 
                        

文档评论(0)