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 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40
|
maille=30
retraitx = 50
retraity = 50
x = 2200
y = 1300
xr = 2200 - (2*retraitx);
yr = 1300 - (2*retraity);
long_fin = floor(xr/maille);
larg_fin = floor(yr/maille);
ex = x/long_fin;
ey = y/larg_fin;
long = 1:long_fin;
larg = 1:larg_fin;
axis([-100 2500 -100 1500])
hold on
rectangle('Position',[0 0 x y])
rectangle('Position',[retraitx retraity xr yr],'Curvature',0.1,'EdgeColor', 'm')
pas=zeros(1,long_fin)
A = zeros([long_fin larg_fin])
for i=larg(1):larg(larg_fin)
esp = retraity:maille:yr + retraity;
pas(i) = (long(i)*maille)+retraitx;
for j=long(1):long(long_fin-1)
A(j,i) = esp(i);
endfor
a = A'(i,1)
plot(pas,a,'r+');
end |
Partager