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

De WikiLICC
Ir para: navegação, pesquisa
m
m
 
Linha 1: Linha 1:
−
N =10;
 
−
x1=-2;  x2= 2;
 
−
y1=-2;  y2= 2;
 
  
−
  for i=1:N
+
==Vetor tridimensional==
−
  for j=1:N
+
  function vetor( v )
−
    r = (i-1)/(N-1);
+
  param3d( [0;v(1)], [0;v(2)], [0;v(3)] )
−
    s = (j-1)/(N-1);     
+
endfunction
−
    x(i,j)=x1 *r + x2*(1-r);
+
 
−
    y(i,j)=y1 *s + y2*(1-s);
+
==Produto vetorial==
−
+
function [w]=prodvet( u, v )
−
    xv(i) =x(i,j);
+
 
−
    yv(j) =y(i,j);
+
  w(1)=u(2)*v(3) - u(3)*v(2);
−
    z(i,j)=1-x(i,j)-y(i,j);
+
  w(2)=u(3)*v(1) - u(1)*v(3);
−
  end
+
   w(3)=u(1)*v(2) - u(2)*v(1);
−
end
+
 
−
+
endfunction
−
surf(x,y, z,'facecol','interp')   
+
 
−
+
== Reta ==
−
  fig = gcf();
+
function reta( p, v )
−
fig.color_map=coolcolormap(32);
+
 
−
  contour2d( xv,yv, z,20)
+
  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

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
Disponível em "http://hunter.mat.ufrgs.br/mediawiki/index.php?title=Si2:Aula2&oldid=2985"