Imagenes En 3D

IMÁGENES 3D close all clear all clc %Se grafican las figuras geometricas % Se ingrean los valores de la ecuacion de perspectiva X0=input('Ingresa elvalor de X0 = '); Y0=input('Ingresa el valor de Y0 = '); Z0=input('Ingresa el valor de Z0 = '); r1=input('Ingresa el valor de r1 = '); r2=r1; r3=r1;theta=input('Ingresa el valor de theta = '); alfa=input('Ingresa el valor de alfa = '); theta=theta*pi/180; alfa=alfa*pi/180; landa=input('Ingresa elvalor de landa = '); A = [1 5 5 1 1 1 5 1 5 5 1 1 5 5 5 5 1 1; 0 0 1 1 0 0 0 0 0 1 1 0 0 0 1 1 1 1; 0 0 0 0 0 4 4 4 4 4 4 4 4 0 0 4 4 0];plot3(A(1,:),A(2,:),A(3,:)) hold on grid on t = 1.5:-0.01:0; [X,Y,Z] = cylinder(t); X1=2+X; Y1=3+Y; Z1=5*Z; mesh(X1,Y1,Z1) axis([0 5 0 5 0 5]) axis squareAt=[A(1,:)',A(2,:)',A(3,:)']; [m n]=size(At); for i=1:m for j=1:nden=(-(At(i,1)-X0)*sin(theta)*sin(alfa)+(At(i,2)Y0)*cos(theta)*sin(alfa)-(At(i,3)-Z0)*cos(alfa)+r3+landa); xR(1,i)= landa*(((At(i,1)-X0)*cos(theta)+(At(i,2)-Y0)*sin(theta)r1))/den; yR(1,i)=landa*((At(i,1)-X0)*sin(theta)*cos(alfa)+(At(i,2)Y0)*cos(theta)*cos(alfa)+(At(i,3)-Z0)*sin(alfa)-r2)/den; end end [o p]=size(X1); for k=1:o for l=1:pden=(-(X1(k,l)-X0)*sin(theta)*sin(alfa)+(Y1(k,l)Y0)*cos(theta)*sin(alfa)-(Z1(k,l)-Z0)*cos(alfa)+r3+landa); xC(k,l)= landa*(((X1(k,l)-X0)*cos(theta)+(Y1(k,l)-Y0)*sin(theta)r1))/den; yC(k,l)=landa*((X1(k,l)-X0)*sin(theta)*cos(alfa)+(Y1(k,l)Y0)*cos(theta)*cos(alfa)+(Z1(k,l)-Z0)*sin(alfa)-r2)/den; end end

figure(2) plot(xR,yR,xC,yC) grid on