Commit a6559747 authored by Conor McCoid's avatar Conor McCoid

Assistant work, MI: updates to solutions for series 2

parent 1bf4087c
......@@ -101,7 +101,7 @@
sur un domaine carré de côtés de longueur un. Son en-tête doit ressembler à
la suivante.
\VerbatimInput[firstline=0,lastline=7]{Heat2D.m}
\VerbatimInput[firstline=0,lastline=7]{serie02newSolutions/Heat2D.m}
\end{enumerate}
......
%% Vecteurs propres en dimension une
%% Espace propre (ex4)
G = [5/2, -1; 1, 0]; % matrice G
v0= [1;2]; v1= G*v0; v2= G*v1; % comportement de v0
u0= [2;1]; u1= G*u0; u2= G*u1; % comportement de u0
x0= linspace(0,2*pi,1000); x0= [x0, x0(:,1)];
x0= sqrt(5)*[cos(x0);sin(x0)]; % anneau
x1= G*x0; x2= G*x1; % comportement de l'anneau
figure(1)
hold on
plot(x0(1,:),x0(2,:),'b--',... % anneau
x1(1,:),x1(2,:),'r--',... % anneau apres 1 appl. de G
x2(1,:),x2(2,:),'m--') % anneau apres 2 appl. de G
quiver(zeros(1,4),zeros(1,4),...
[v0(1),u0(1),-v0(1),-u0(1)],... % v0 et u0
[v0(2),u0(2),-v0(2),-u0(2)],'b')
quiver(zeros(1,4),zeros(1,4),...
[v1(1),u1(1),-v1(1),-u1(1)],... % v0 et u0 apres 1 appl. de G
[v1(2),u1(2),-v1(2),-u1(2)],'r')
quiver(zeros(1,4),zeros(1,4),...
[v2(1),u2(1),-v2(1),-u2(1)],... % v0 et u0 apres 2 appl. de G
[v2(2),u2(2),-v2(2),-u2(2)],'m')
hold off
axis equal
line([0 0], ylim,'color','k') % les axes
line(xlim, [0 0],'color','k')
%% Vecteurs propres en dimension une (ex5)
n = 100; % discretization
x = linspace(0,1,n+2)'; % interval unitaire
......@@ -22,7 +50,7 @@ xlabel('x')
ylabel('Vecteurs propres')
axis([0,1,-1,1])
%% Vecteurs propres en dimension deux
%% Vecteurs propres en dimension deux (ex5)
n = 20; % discretization
x = linspace(0,1,n+2)'; % carr'e unitaire
......@@ -55,4 +83,6 @@ figure(3) % les vecteurs propres discrets
for i = 1:k^2
subplot(k,k,i)
surf(x(2:end-1),x(2:end-1)',reshape(V(:,i),n,n))
end
\ No newline at end of file
end
%% Conditions au bord (ex6)
\ No newline at end of file
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment