资源描述:
《距离判别源代码清单》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、距离判别源代码清单OptionBase1'********************************************************************************************'函数部分'''********************************************************************************************Subputout(A()AsDouble,m%,n%)'m为矩阵的行,n为矩阵的列'***************************
2、*************************************************************'功能:输出矩阵'参数:a为待输出的矩阵;m为a矩阵的行;n为a矩阵的列'返回值:无返回值'****************************************************************************************Fori=1TomForj=1TonPrintTab(j*20);A(i,j);'Format(a(i,j),"0.0000");NextjPrintNextiEndSubSu
3、bpingjun(A()AsDouble,m%,qd%,zd%,pj()AsDouble)'****************************************************************************************'功能:求m列矩阵从第qd列到第zd列的平均'参数:a为待求矩阵;m为矩阵的行数;qd为起始列下标;zd终止列下标'返回值:pj返回最终结果,是一m为列向量'**********************************************************************
4、******************ReDimpj(m)AsDoubleFori=1TomForj=qdTozdpj(i)=pj(i)+A(j,i)Nextjpj(i)=pj(i)/(zd-qd+1)NextiEndSubSublicha(A()AsDouble,m%,qd%,zd%,lc()AsDouble)'****************************************************************************************'功能:求m列矩阵从第qd列到第zd列的离差'参数:a为待求矩阵;m为矩阵的
5、行数;qd为起始列下标;zd终止列下标'返回值:lc返回最终结果,是一m为列向量'****************************************************************************************ReDimlc(m,m)AsDouble'原来的lc矩阵可能已存入非零数据,所以要重新声明使其数据均为零Dimpj1()AsDoubleCallpingjun(A(),m,qd,zd,pj1())Forr=1TomFort=1TomForj=qdTozdsum=sum+(A(j,t)-pj1(t))*(A(
6、j,r)-pj1(r))Nextjlc(t,r)=sumsum=0Nextt,rEndSubSubchengfa(n()AsDouble,p()AsDouble,m()AsDouble)'****************************************************************************************'功能:求两矩阵的乘积'参数:n,p为待求乘积的两矩阵(n的列数必须与p的行数相同,否则结果无效)'返回值:m返回乘积结果的矩阵,行数与n的相同,列数与p的相同'********************
7、******************************************************************A=UBound(n,1):B=UBound(n,2):B=UBound(p,1):c=UBound(p,2)ReDimm(A,c)Fori=1ToAForj=1TocFork=1ToBm(i,j)=m(i,j)+n(i,k)*p(k,j)NextNextNextEndSubSubzhuanzhi(A()AsDouble,B()AsDouble)'******************************************
8、*******************