| 查看: 2165 | 回复: 14 | |||
| 当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖 | |||
[交流]
近红外kennard-stone选样本
|
|||
| 各位大侠,我最近做近红外,图谱扫描了,选择校正集合验证集想采用kennard-stone选样本,该如何操作,具体怎么弄?是用matlab吗,又该怎么弄? |
» 猜你喜欢
科研人应该花精力去思考如何解决问题,而不是去凝练问题
已经有11人回复
要骂人了,新模版改版就是要淡化问题凝练这种虚的东西,结果有个评委还在说凝练得不够
已经有17人回复
基金系统什么内容也没有
已经有10人回复
学科评审组评审是指会评吗?
已经有5人回复
梦想
已经有8人回复
面上合作单位盖章
已经有9人回复
面上函评意见出来了,像什么等级?
已经有17人回复
为什么资助数各大高校都创新高,自己申请怎么就这么难
已经有14人回复
小白求助 投论文要求的highlights应该如何写
已经有4人回复
麻烦专家们看看评委们的意见(F口面上)
已经有12人回复
» 抢金币啦!回帖就可以得到:
云南大学丁军桥团队2027年博士招生
+1/810
【征稿中!】2026年机械、制造技术与材料国际会议 (MEMTM 2026)
+1/774
[求助]请问美国有哪些PC品牌厂商?
+1/189
哈尔滨工业大学(校本部)电信学院2027级博士招生
+1/181
2027级博士研究生——质谱分析
+5/135
广西大学作物学博士生招生
+2/50
同济大学环境学院肖倩课题组招收博士研究生(资格审核制)
+1/43
香港城市大学生物物理/软物质材料课题组2026-2027年招博士生1名
+2/28
评职称需要,卖发明专利了:D有技术储备需求的设备厂商可联系
+1/19
上海大学微电子学院杨军教授团队招收微电子/集成电路/材料/机械/微纳制造博士一名
+1/11
上海大学微电子学院杨军教授团队招收微电子/集成电路/材料/机械/微纳制造博士一名
+1/10
孩子爱看的纪录片,远比刷短视频更有价值
+1/9
哈工大马樱招收2027级计算机类/集成电路类博士生
+1/9
上海大学微电子学院杨军教授团队招收微电子/集成电路/材料/机械/微纳制造博士一名
+1/5
【科研诚信】武汉轻工大学何江玲3篇论文:涉嫌抄袭、成组数据重复、制造虚假显著
+1/5
国防补偿计划 推免硕士生及直博生
+1/4
义翘神州推出全长p-Tau217蛋白,助力阿尔茨海默病下一代生物标志物检测开发
+1/3
吉林大学高压与超硬材料全国重点实验室李达教授课题组招收2027申请考核制博士
+1/2
科研助理招聘
+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










回复此楼
=[];