1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
| [numRadii, numCurves] = size(P);
gridsOptions.number = 10;
ticksize = 3;
radius = 100;
th = linspace(0, 2*pi, numRadii+1);
x = NaN(numRadii*2+(numRadii-1), gridsOptions.number+1);
y = NaN(numRadii*2+(numRadii-1), gridsOptions.number+1);
x(1,:) = ticksize/2;
x(2,:) = -ticksize/2;
y(1,:) = radius * linspace(1, 0, gridsOptions.number+1);
y(2,:) = radius * linspace(1, 0, gridsOptions.number+1);
for n = 2:numel(th)-1
tmp = [cos(th(n)) sin(th(n)) ; -sin(th(n)) cos(th(n))]*[x(1,:) ; y(1,:)];
x(3*(n-1)+1,:) = tmp(1,:);
y(3*(n-1)+1,:) = tmp(2,:);
tmp = [cos(th(n)) sin(th(n)) ; -sin(th(n)) cos(th(n))]*[x(2,:) ; y(2,:)];
x(3*(n-1)+2,:) = tmp(1,:);
y(3*(n-1)+2,:) = tmp(2,:);
end
line(x(:,1:end-1), y(:,1:end-1), 'color', 'k') |
Partager