gadproj.m

这是几何代数的matlab工具包

% 	PROJECTION, REJECTION
GAfigure; clc; %/
% 	PROJECTION, REJECTION
%
global x a A p r; %/
clf; %/
% 	Take two vectors x and a:
%.
x = e1 + e2 + e3;
a = e2/4 + e3/2;
draw(x,'r'); %/
random=unit(pi*e1+pi/exp(1)*e2+exp(1)*e3); %/
GAtext(0.7*x+0.1*unit(grade((random^x)/x,1)),'x'); %/
draw(a,'g'); %/
GAtext(0.7*a-0.1*unit(grade((random^a)/a,1)),'a'); %/
va = [-0.05 1.1 0 1.25 -0.25 1.5]; %/
axis(va); %/
axis('vis3d'); %/
axis off; %/
GAprompt; %/
GAorbiter(60,3); %/
axis(va); %/
GAprompt; %/
% 	projection of x onto a: 
%.
p = inner(x,a)/a 
draw(p,'m'); %/
draw(a,'g'); %/
DrawPolyline({x,p},'k'); %/
GAtext(0.7*p-0.1*unit(grade((random^p)/p,1)),'(x \bullet a) / a'); %/ 
axis(va); %/
GAprompt; %/
% 	rejection of x by a:
%.
r = (x^a)/a  
draw(r,'m'); %/
DrawPolyline({x,r},'k'); %/
GAtext(0.5*r-0.1*unit(grade((random^r)/r,1)),'r = (x \wedge a) / a'); %/ 
axis(va); %/
GAprompt; %/
axis(va); %/
GAorbiter(360,10); %/ 
axis(va); %/
GAprompt; %/
% 	Explanation: 
axis(va); %/
GAprompt; %/
%	span x^a,
DrawBivector(x,a,'y'); %/
GAtext((x+a)/2 + 0.1*(x^a)/I3,'x \wedge a'); %/
axis(va); %/
% 	reshape orthogonally,
GAprompt; %/
DrawBivector(r,a,'y'); %/
GAtext((r+a)/2 + 0.1*(x^a)/I3,'x \wedge a'); %/
axis(va); %/
GAprompt; %/
title('x \wedge a = r \wedge a = r \bullet a + r \wedge a = r a','FontSize',12);%/
% 	division by a then gives the rejection.
axis(va); %/
GAprompt; %/
%.
%	Projection and rejection formulas also work for bivectors:	
GAprompt; %/
clf; %/
x = e1/3 - e2/2 + e3;
A = e1^(e2/3 + e3/2);
draw(x,'r'); %/
axis off; %/
axis([-0.5 0.5 -0.8 0.2 -0.4 1]); %/
random=unit(pi*e1+pi/exp(1)*e2+exp(1)*e3); %/
GAtext(0.7*x-0.1*unit(grade((random^x)/x,1)),'x'); %/
draw(A,'g'); %/
GAtext(-0.4*unit(grade(inner(random,A)/A,1))+0.1*unit(A/I3),'A'); %/ 
va = [-0.5 0.5 -0.81 0.31 -0.4 1]; %/
axis(va); %/
axis('vis3d'); %/
GAprompt; %/
GAorbiter(360,4); %/
axis off; %/
% 	projection of x onto A: 
%.
p = inner(x,A)/A  
draw(p,'m'); %/
DrawPolyline({x,p},'k'); %/
GAtext(0.9*p+0.1*unit(grade((random^p)/p,1)),'(x \bullet A) / A'); %/ 
axis(va); %/
GAprompt; %/
% 	rejection of x by a:
%.
r = (x^A)/A  
draw(r,'m'); %/
DrawPolyline({x,r},'k'); %/
GAtext(0.7*r-0.1*unit(grade((random^r)/r,1)),'(x \wedge A) / A'); %/ 
axis(va); %/
GAprompt; %/
GAorbiter(340,10); %/