参考基于matlab的字符识别研究.docVIP

下载本文档

6
0
约2.12万字
约 40页
2018-06-08 发布于贵州
举报
版权申诉

参考基于matlab的字符识别研究.doc

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的字符识别研究

? ? 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(:,