Adrià Vilanova Martínez | ce8597a | 2021-09-21 18:26:59 +0200 | [diff] [blame^] | 1 | % Dibuix d'una banda de Moebius |
| 2 | % (codi amb Matlab i apte també per a Octave) |
| 3 | % Geometria Diferencial, FME, UPC |
| 4 | % 2021/02/12 |
| 5 | |
| 6 | % Paràmetres de la banda de Moebius |
| 7 | r = 11; % radi de la banda |
| 8 | t = linspace(0, 2*pi, 101); |
| 9 | zor = linspace(-1, 1, 21); % ample de la banda |
| 10 | |
| 11 | % Taula de punts de la banda |
| 12 | for k = 1:length(t) |
| 13 | A = [ |
| 14 | cos(t(k)) -sin(t(k)) 0 |
| 15 | sin(t(k)) cos(t(k)) 0 |
| 16 | 0 0 1 |
| 17 | ]; |
| 18 | B = [ |
| 19 | 1 0 0 |
| 20 | 0 cos(t(k)/2) -sin(t(k)/2) |
| 21 | 0 sin(t(k)/2) cos(t(k)/2) |
| 22 | ]; |
| 23 | |
| 24 | for l = 1:length(zor) |
| 25 | phi = A*(B*[0; r; zor(l)] + [0; r; 0]); |
| 26 | x(k, l) = phi(1); |
| 27 | y(k, l) = phi(2); |
| 28 | z(k, l) = phi(3); |
| 29 | end |
| 30 | end |
| 31 | |
| 32 | % Dibuixem la superfície |
| 33 | %xlabel('x'); |
| 34 | %ylabel('y'); |
| 35 | %zlabel('z'); |
| 36 | surf(x, y, z); |
| 37 | print("01.png", "-dpng") |
| 38 | %axis equal |