| 查看: 1790 | 回复: 14 | |||
| 当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖 | |||
[交流]
近红外kennard-stone选样本
|
|||
| 各位大侠,我最近做近红外,图谱扫描了,选择校正集合验证集想采用kennard-stone选样本,该如何操作,具体怎么弄?是用matlab吗,又该怎么弄? |
» 猜你喜欢
垃圾破二本职称评审标准
已经有19人回复
职称评审没过,求安慰
已经有53人回复
毕业后当辅导员了,天天各种学生超烦
已经有5人回复
26申博自荐
已经有3人回复
A期刊撤稿
已经有4人回复
» 抢金币啦!回帖就可以得到:
散金祈福
+1/612
郑州大学招收2026年外科学博士研究生
+1/461
DIY科研工具交流
+1/211
双一流南京医科大学招计算机、AI、统计、生物信息等方向26年9月入学博士
+1/183
湖北大学食品安全研究团队诚招博士后
+1/182
《春节催婚季,90后男生在线“招募”战友,一起过个轻松年》
+1/155
复旦大学化学系杰青/优青团队诚聘博士后
+3/150
Analytical Science Advances(Wiley出版社)长期征稿中...
+1/86
澳門大學 土木工程碩士和建造项目管理碩士招聘
+1/80
天津大学化学系吴立朋课题组申请考核制博士招生/博后招聘
+1/80
重庆大学杰青团队诚招2026年博士研究生
+2/46
博士招生
+1/32
南京-栖霞区-尧化门附件有房子出租吗?
+1/29
博士后招聘
+1/27
求硫醇的合成资源
+1/19
南京林业大学”申请-考核”制学术学位博士研究生招生
+1/5
实验室研究员
+1/4
山东大学集成电路学院博士招生1名
+1/4
想读博士(理工/化学类/新能源材料)求大佬捞一下 ?或者师兄师姐推荐/建议博导招生信息
+1/3
南京大学FinTech大模型实验室招募斯坦福国际联培博士生(2026)
+1/1
8楼2012-10-21 16:16:52
★
小木虫: 金币+0.5, 给个红包,谢谢回帖
小木虫: 金币+0.5, 给个红包,谢谢回帖
|
用下面的matlab程序即可。其中X是需要被分组的样品的近红外光谱样品集,k是需要选出的样品个数。 function [model,test]=kenstone(X,k) tic [m,n]=size(X); if k>=m | k<=0 h=errordlg('Wrongly specified number of objects to be selected to model set.','Error'); model=[]; if nargout==2 test=[]; end waitfor(h) return end x=[[1:size(X,1)]' X]; n=size(x,2); [i1,ind1]=min(fastdist(mean(x(:,2:n)),x(:,2:n))); model(1)=x(ind1,1); x(ind1, =[];[i2,ind2]=max(fastdist(X(model(1), ,x(:,2:n)));model(2)=x(ind2,1); x(ind2, =[];%h=waitbar(0,'Please wait ...'); %h=waitbar(0/k,h); for d=3:k [ii,ww]=max(min(fastdist(x(:,2:n),X(model, )));model(d)=x(ww,1); x(ww, =[];% h=waitbar(d/k,h); end if nargout==2 test=1:size(X,1); test(model)=[]; end %close(h); toc function D=fastdist(x,y) % Calculated Euclideam distances between two sets of objetcs D=((sum(y'.^2))'*ones(1,size(x,1)))+(ones(size(y,1),1)*(sum(x'.^2)))-2*(y*x'); |
2楼2012-10-16 15:35:36
★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
leecius: 金币+1, 谢谢参与。 2012-10-18 09:05:10
小木虫: 金币+0.5, 给个红包,谢谢回帖
leecius: 金币+1, 谢谢参与。 2012-10-18 09:05:10
|
不知道为什么会有笑脸出来,这个可能是个bug吧,我传到附件上了,你自己看吧。 http://g.zhubajie.com/urllink.php?id=128733297robb6l7cg4keoy0 |
3楼2012-10-16 15:39:09
4楼2012-10-17 21:19:32













回复此楼
=[];