Augmenting Latin Hypercube Points in MATLAB -

-

i know commands create lhs design , augment more points later on if model not enough? example, first create 50 points lhs design, add more points (perhaps in batches of 20) incrementally until model accurate enough. example:

set1=lhsdesign(5,5); %5x5 matrix %use of set 1, determine more points needed  set2=%some command adds 20 points set1 make 25x5 matrix

the difficulty in running lhsdesign again 20 new points not take account of original points. there work done using original points , generating new set of points scratch wastes work new points not included in new set.

i able make function solves problem. not sure final matrix true latin hypercube, adds required number of points given points , moves them closest available open 'channel' (sub range no point exists) if need be.

usage follows:

x1=lhsdesign(200,17); xf=lhsaugment(x1,200);

which adds 200 points x1 set of points, resulting in xf being 400x17 matrix. function follows:

function xf = lhsaugment(x1,npoi) %function xf = lhsaugment(x1,npoi) %function augment given latin hypercube x1 number of points, %npoi. length changed, i.e. points added length. %the original points left unctouched , appear first in output %xf. size of xf [size(x1,1)+npoi size(x1,2)]. x2=lhsdesign(npoi,size(x1,2)); npoi=size(x2,1); opoi=size(x1,1); tpoi=npoi+opoi; fint=1/tpoi; i=1:tpoi     cbound(i,:)=[(i-1)*fint i*fint]; end xf=zeros(tpoi,size(x1,2)); bx1=zeros(size(x1)); bx2=zeros(size(x2)); bf=zeros(tpoi,size(x1,2)); if=zeros(1,size(x1,2)); imove=0; i=1:opoi     j=1:size(cbound,1)         l=1:size(x1,2)             if (x1(i,l)>cbound(j,1))&&(x1(i,l)<=cbound(j,2))&&(bf(j,l)==0)                 if(1,l)=if(1,l)+1;                 xf(if(1,l),l)=x1(i,l);                 bx1(i,l)=1;                 bf(j,l)=1;             elseif (x1(i,l)>cbound(j,1))&&(x1(i,l)<=cbound(j,2))&&(bf(j,l)~=0)                 imin=size(cbound,1);                 pmin=size(cbound,1);                 m=j:-1:1                     if (bf(m,l)==0)                         imin=m;                         pmin=j-m;                         break                     end                 end                 m=j:size(cbound,1)                     if (bf(m,l)==0)&&(m-j<pmin)                         imin=m;                         pmin=j+m;                         break                     end                 end                 if(1,l)=if(1,l)+1;                 xf(if(1,l),l)=x1(i,l);                 bx1(i,l)=1;                 bf(imin,l)=1;             end         end     end end i=1:npoi     j=1:size(cbound,1)         l=1:size(x2,2)             if (x2(i,l)>cbound(j,1))&&(x2(i,l)<=cbound(j,2))&&(bf(j,l)==0)                 if(1,l)=if(1,l)+1;                 xf(if(1,l),l)=x2(i,l);                 bx2(i,l)=1;                 bf(j,l)=1;             elseif (x2(i,l)>cbound(j,1))&&(x2(i,l)<=cbound(j,2))&&(bf(j,l)~=0)                 imin=size(cbound,1);                 pmin=size(cbound,1);                 m=j:-1:1                     if (bf(m,l)==0)                         imin=m;                         pmin=j-m;                         break                     end                 end                 m=j:size(cbound,1)                     if (bf(m,l)==0)&&(m-j<pmin)                         imin=m;                         pmin=j+m;                         break                     end                 end                 if(1,l)=if(1,l)+1;                 xf(if(1,l),l)=(x2(i,l)-(floor(x2(i,l)/fint)*fint))+((imin-1)*fint);                 bx2(i,l)=1;                 bf(imin,l)=1;                 if l==1                 imove=imove+1;                 end             end         end     end     end

Comments

Post a Comment

Popular posts from this blog

node.js - Mongoose: Cast to ObjectId failed for value on newly created object after setting the value -

-
this question has answer here: what's mongoose error cast objectid failed value xxx @ path “_id”? 7 answers even though has been asked many times because entity has no valid objectid field. i came across situation objectid set value can casted objectid still error. the schema: var serviceschema = parammongoose.schema({ servicename:{ type: parammongoose.schema.types.objectid , ref: 'servicenames',required:true } }); pseudo code: servicefromdb=new service({servicename:'some name'}); servicefromdb.servicename='000000000000000000000001'; servicefromdb.save(function(paramerror,paramdata){ if(paramerror){ console.log('but but...',servicefromdb,paramerror) } }); the output of code is: but but... { servicename: 000000000000000000000001, _id: 55079a90286f49280364f78b } { [casterror: cast objectid faile...
Read more

gradle error "Cannot convert the provided notation to a File or URI" -

-
i have defined 2 system properties in gradle.properties: systemprop.builddir=build systemprop.cfgdir=build\\cfg and have following task defined in build.gradle: task clean(group:'clean',description:'clean project') << { ant.sequential { delete(dir: system.properties.builddir) mkdir(dir: system.properties.builddir) delete(dir: system.properties.cfgdir) mkdir(dir: system.properties.cfgdir) } } this generates following error: execution failed task ':clean'. > cannot convert provided notation file or uri: {dir=build}. following types/formats supported: - string or charsequence path, e.g 'src/main/java' or '/usr/include' - string or charsequence uri, e.g 'file:/usr/include' - file instance. - uri or url instance. but equivalent block of code not generate error , works expected: task clean(group:'clean',description:'clean project') << { ...
Read more

python - NameError: name 'subprocess' is not defined -

-
i'm trying make rng website: http://jamesdotcom.com/?p=417 . when run following code: captureimage = subprocess.popen( [ "fswebcam", "-r", "356x292", "-d", "/dev/video0", "static.jpg", "--skip", "10" ], stdout=devnull, stderr=devnull) captureimage.communicate() i following error message: nameerror: name 'subprocess' not defined what should do? thanks have imported first? import subprocess captureimage = subprocess.popen(["fswebcam", ...) subprocess included in python standard library since python 2.4. should not necessary install it.
Read more