Mudanças entre as edições de "Si2:Aula2"

De WikiLICC
Ir para: navegação, pesquisa
m (Criou página com 'N =10; x1=-2; x2= 2; y1=-2; y2= 2; for i=1:N for j=1:N r = (i-1)/(N-1); s = (j-1)/(N-1); x(i,j)=x1 *r + x2*(1-r); y(i,j)=y1 *s + y2*(1-s); xv(i) =x(i…')
 
m
 
(Uma revisão intermediária pelo mesmo usuário não está sendo mostrada)
Linha 1: Linha 1:
N =10;
 
x1=-2;
 
x2= 2;
 
y1=-2;
 
y2= 2;
 
  
 +
==Vetor tridimensional==
 +
function vetor( v )
 +
  param3d( [0;v(1)], [0;v(2)], [0;v(3)] )
 +
endfunction
  
for i=1:N
+
==Produto vetorial==
  for j=1:N
+
function [w]=prodvet( u, v )
    r = (i-1)/(N-1);
+
 
    s = (j-1)/(N-1);     
+
  w(1)=u(2)*v(3) - u(3)*v(2);
    x(i,j)=x1 *r + x2*(1-r);
+
   w(2)=u(3)*v(1) - u(1)*v(3);
    y(i,j)=y1 *s + y2*(1-s);
+
  w(3)=u(1)*v(2) - u(2)*v(1);
 +
 
 +
endfunction
  
     xv(i) =x(i,j);
+
== Reta ==
     yv(j) =y(i,j);
+
function reta( p, v )
     z(i,j)=1-x(i,j)-y(i,j);
+
 
   end
+
  for i=1:30 
end
+
     t=(i-1)/(30-1)*10-5;
 +
    x(i)= p(1)+v(1)*t;
 +
     y(i)= p(2)+v(2)*t;
 +
     z(i)= p(3)+v(3)*t;
 +
   end  
 +
  param3d( x,y,z)
 +
   
 +
endfunction
  
surf(x,y, z,'facecol','interp')
+
==Projeção==
 
+
function [p]=projecao( v, b )
fig = gcf();
+
 
fig.color_map=coolcolormap(32);
+
  p=((v'*b) / norm(b)^2 )* b  ;
contour2d( xv,yv, z,20)
+
endfunction

Edição atual tal como às 18h51min de 31 de agosto de 2010

Vetor tridimensional

function vetor( v )
  param3d( [0;v(1)], [0;v(2)], [0;v(3)] )
endfunction

Produto vetorial

function [w]=prodvet( u, v )
  
  w(1)=u(2)*v(3) - u(3)*v(2);
  w(2)=u(3)*v(1) - u(1)*v(3);
  w(3)=u(1)*v(2) - u(2)*v(1);
 
endfunction

Reta

function reta( p, v )
 
 for i=1:30  
   t=(i-1)/(30-1)*10-5;
   x(i)= p(1)+v(1)*t;
   y(i)= p(2)+v(2)*t;
   z(i)= p(3)+v(3)*t;
 end   
 param3d( x,y,z)
   
endfunction

Projeção

function [p]=projecao( v, b )
  
 p=((v'*b) / norm(b)^2 )* b  ;
endfunction