exm020602.m

来自「matlab电子书籍」· M 代码 · 共 20 行

syms t
A=sym([sqrt(3)/2,1/2;1/2,sqrt(3)/2])
G=[cos(t),-sin(t);sin(t),cos(t)];			
GA=G*A 

An=subs(GA,t,110/180*pi);						
Op=[0;0];													
Ad=double(A);											
v1=[Op,Ad(:,1)]';v2=[Op,Ad(:,2)]';		
u1=[Op,An(:,1)]';u2=[Op,An(:,2)]';		
plot(v1(:,1),v1(:,2),'--k',v2(:,1),v2(:,2),'b')
axis([-1,1,-1,1]),axis square
hold on
hu=plot(u1(:,1),u1(:,2),'--k',u2(:,1),u2(:,2),'b');
set(hu,'LineWidth',4)
title('Givens Rotation')
legend(' v1',' v2',' u1',' u2','Location','South')
hold off 
grid on