Revision 1105, 0.6 kB
(checked in by zimamiro, 14 years ago)
|
|
Rev | Line | |
---|
[1105] | 1 | function [box zero range eigvec]=make_box2(points)
|
---|
| 2 | eps=10^-6;
|
---|
| 3 | %OK
|
---|
| 4 |
|
---|
| 5 | N=size(points,1);
|
---|
| 6 | dim=size(points,2);
|
---|
| 7 | num_of_points=round(N^(1/dim));
|
---|
| 8 | eigvec=eye(2);
|
---|
| 9 |
|
---|
| 10 | range=(max(points)-min(points))';
|
---|
| 11 |
|
---|
| 12 |
|
---|
| 13 | box=zeros(N,dim);
|
---|
| 14 | pointer=min(points);
|
---|
| 15 | step=range/(num_of_points-1);
|
---|
| 16 | jump=[1 num_of_points ];
|
---|
| 17 | for i=1:N
|
---|
| 18 | for k=1:dim
|
---|
| 19 | box(i,k)=pointer(k)+step(k)*mod(floor((i-1)/jump(k)),num_of_points);
|
---|
| 20 | end
|
---|
| 21 | end
|
---|
| 22 |
|
---|
| 23 | range=ones(dim,1)./(2*range);
|
---|
| 24 | zero=box(1,:);
|
---|
| 25 |
|
---|
| 26 | %plot(box(:,1), box(:,2),'b.' ); hold on; plot(points(:,1), points(:,2), 'r.' ); hold off
|
---|
| 27 | end |
---|