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