v
aus
komplexen Zahlen erzeugt wird, welcher beim Aufruf
plot(v); axis equal
ein reguläres Achteck zeichnet,
das im Kreis um mit dem
Radius 5 einbeschrieben ist.
Geben Sie auch die MATLAB-Befehle an zum Zeichnen des zugehörigen Umkreises.
k = 0:8; w = 2*pi/8*k z8 = 5*exp(i*w) plot(z8) hold on axis equal t = (0:0.005:1)*2*pi; kr = 5*exp(i*t); plot(kr,'k') hold off
A = [0 0 0]'; B = [0 -12 0]'; C = [5 -12 0]'; D = [5 0 0]' ; E = [0 0 5]' ; F = [5 0 5]'; Kl = [A B C D A E B E F C F D]; plot3(Kl(1,:),Kl(2,:),Kl(3,:),'m') axis equal eng = [ 0 0 1]', dkritg = 0 % Tests ob A,B,C,D in Grundflaeche dtA = eng'*A-dkritg, dtB = eng'*B-dkritg, dtC = eng'*C-dkritg, dtD = eng'*D-dkritg, % 2 Vektoren in Deckflaeche BCFE u = B-E, v = C-E N = cross(u,v) % = [ 0 -25 60] nN = norm(N) % = 65 endf = N/nN % = [0 -0.3846 0.9231] dkritdf = endf'*E % = 4.6154 % Tests ob B,C,F auch in Grundflaeche dtdB = endf'*B-dkritdf, dtdC = endf'*C-dkritdf, dtdF = endf'*F-dkritdf % Winkel direkt acos S-prod zwischen Einheitsvektoren wneig = acos(endf'*eng) % 0.3948 rad wneigd = wneig*180/pi % = 22.6154 Grad
Qo = [0 -6 -12 -6 0; 0 6 0 -6 0; 1 1 1 1 1]; % Mittelpunkt quadrat ist A+B+C+D/4 =[ -6 0]' Tz = [ 1 0 6; 0 1 0; 0 0 1] R = [0 1 0 ; -1 0 0; 0 0 1] Tb = [ 1 0 -6; 0 1 0; 0 0 1] Ttot = Tb*R*Tz Qt = Ttot*Qo plot(Qo(1,:),Qo(2,:),'g') hold on plot(Qt(1,:),Qt(2,:),'ro') axis equal ; hold off
t = (0:0.05:1)*pi % eine halbe Drehung % Vorschub ist 2 pro halbe Drehung, also ist g = 4 g = 4 z = t*g/(2*pi); % gleiches z fuer beide % Winkel 0 ist bei 5 0, also R = 5 R = 5 % Vorschub nach oben, also Gegenuhrzeiger in x-y Ebene ist rechtsdrehend xr = R*cos(t); yr = R*sin(t); % Uhrzeigersinn in x-y Ebene ist linksdrehend xl = R*cos(t); yl = -R*sin(t); plot3(xr,yr,z) ; hold on; plot3(xl,yl,z,'r') axis equal