24小时热门版块排行榜

>论坛更新日志 (4957)
>导师招生 (505)
>文献求助 (424)
>虫友互识 (292)
>硕博家园 (224)
>考博 (212)
>学术会议 (198)
>论文投稿 (194)
>教师之家 (169)
>找工作 (140)
>论文道贺祈福 (121)
>基金申请 (85)
>博后之家 (78)
>招聘信息布告栏 (65)
>休闲灌水 (61)
>考研 (52)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

qiusj

金虫 (正式写手)

应助: 4 (幼儿园)
金币: 4445.9
帖子: 587
在线: 272.4小时
虫号: 974224

[交流] 近红外kennard-stone选样本

各位大侠，我最近做近红外，图谱扫描了，选择校正集合验证集想采用kennard-stone选样本，该如何操作，具体怎么弄？是用matlab吗，又该怎么弄？

回复此楼

» 猜你喜欢

面上没中，邀请各位路过的虫友分析一下分数已经有4人回复
哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生已经有3人回复
找导师已经有8人回复
国社科又开始会评了，不知道这次命运如何已经有14人回复
麻烦专家们看看评委们的意见（F口面上）已经有13人回复
科研人应该花精力去思考如何解决问题，而不是去凝练问题已经有14人回复
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够已经有17人回复
学科评审组评审是指会评吗？已经有5人回复
小白求助投论文要求的highlights应该如何写已经有4人回复
面上意见出来了已经有23人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2012-10-14 10:16:06

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

Nkxz

铁杆木虫 (著名写手)

应助: 14 (小学生)
金币: 9349.4
帖子: 1251
在线: 267.7小时
虫号: 463532

★
小木虫: 金币+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');