24小时热门版块排行榜

>论坛更新日志 (4434)
>导师招生 (560)
>虫友互识 (396)
>文献求助 (320)
>招聘信息布告栏 (209)
>休闲灌水 (202)
>论文投稿 (190)
>学术会议 (182)
>考博 (131)
>基金申请 (117)
>博后之家 (109)
>硕博家园 (67)
>教师之家 (40)
>找工作 (35)
>绿色求助(高悬赏) (29)
>考研 (29)

返回列表

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

qiusj

金虫 (正式写手)

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

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

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

回复此楼

» 猜你喜欢

面上没中，邀请各位路过的虫友分析一下分数已经有10人回复
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急已经有9人回复
广西大学-广州大学招聘博士后欢迎广大优秀人才！！！已经有6人回复
两块石头已经有5人回复
哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生已经有3人回复
找导师已经有8人回复
国社科又开始会评了，不知道这次命运如何已经有14人回复
科研人应该花精力去思考如何解决问题，而不是去凝练问题已经有14人回复
要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够已经有17人回复
学科评审组评审是指会评吗？已经有5人回复

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

查看全部散金贴

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

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

龙缘浪子

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 434.1
帖子: 81
在线: 66.6小时
虫号: 728239

引用回帖:

2楼: Originally posted by Nkxz at 2012-10-16 15:35:36
用下面的matlab程序即可。其中X是需要被分组的样品的近红外光谱样品集，k是需要选出的样品个数。

function =kenstone(X,k)
tic
=size(X);
if k>=m | k<=0
h=errordlg('Wrongly specified nu ...

将笑脸替换为‘：’和‘）’

赞一下

回复此楼

15楼2018-06-25 19:09:00

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

查看全部 15 个回答

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');

赞一下(1人)

回复此楼

2楼2012-10-16 15:35:36

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