| 查看: 1905 | 回复: 14 | |||
| 当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖 | |||
[交流]
近红外kennard-stone选样本
|
|||
| 各位大侠,我最近做近红外,图谱扫描了,选择校正集合验证集想采用kennard-stone选样本,该如何操作,具体怎么弄?是用matlab吗,又该怎么弄? |
» 猜你喜欢
309求调剂
已经有8人回复
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研经历
已经有8人回复
求调剂
已经有5人回复
求化学调剂
已经有5人回复
085600,材料与化工321分求调剂
已经有10人回复
0703 化学 求调剂,一志愿山东大学 342 分
已经有5人回复
一志愿南开大学0710生物学359求调剂
已经有3人回复
085600,专业课化工原理,320分求调剂
已经有4人回复
材料与化工272求调剂
已经有16人回复
290求调剂
已经有3人回复
» 抢金币啦!回帖就可以得到:
南方科技大学纳米-飞秒-超快动力学课题组---诚聘博士后/助理研究员
+1/189
【招聘】温州医科大学/附属医院/瓯江实验室李校堃院士团队诚聘博士后和科研人员
+1/187
南方科技大学纳米-飞秒-超快动力学课题组---诚聘博士后/助理研究员
+1/187
南昌航空大学材料学院教授/博导团队招收调剂研究生
+1/86
电力全国重点实验室双一流A类长江学者团队招2026年全日制博士1名(4月3日截止)
+2/68
西京学院土木水利 2026 级研究生招生
+2/40
欢迎化学、生物与医药、药学,生物、食品等相关专业的同学
+1/38
福建农林大学材料院招收调剂硕士生
+1/37
浙江省人民医院博士后岗位招聘
+1/36
浙江大学光电学院、极端光学实验室 阿秒方向 招聘科研助理2名
+1/12
西华大学材料学院表面科学与工程技术科研团队2026年招收研究生
+1/11
食品调剂-总分291-有中核论文导师一作、本人二作-一志愿南昌大学学硕-介绍详细
+1/8
【博士招生】南京林业大学车辆工程课题组招收26年9月入学“申请-考核制”博士研究生
+1/7
武汉工程大学化学与环境工程学院招收调剂生
+1/6
长江师范学院电子信息硕士和学科物理硕士调剂招生
+1/6
上海交通大学陈倩栎课题组招收2026年申请-考核博士生1名(能源材料方向)
+1/4
内蒙古科技大学铁电新能源材料与器件团队简介
+1/3
香港中文大学 生物医学学院/组织工程与再生医学研究所 招聘博士后/研究助理
+1/1
武汉理工大学船海与能源动力工程学院刘鹏飞院士团队博士后招聘启事
+1/1
南开大学环境科学与工程学院王鑫教授课题组招收申请考核制博士生
+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














回复此楼
=[];