% L/matlab/lampo.m   29.11.04

% Visualisoi a) pinta, b) lämpötilaprofiilit aikasiivuina.
% u(x,t)=100*sin(pi*x/80)*exp(-0.0001785*t)
% Rakennetaan lämpötilamatriisi U.
% Esim. 1 KRE s. 603
close all; clear
x=linspace(0,80,50);
t=linspace(0,400,20);
[X,T]=meshgrid(x,t);
U=100*sin(pi*X./80).*exp(-0.001785*T);
surf(x,t,U); shg
colorbar
figure
plot(x,U)
title('80 cm pitkän sauvan lämpötilajakaumia, alkulämpötila sin(pi*x/80)')
grid on;shg

% Esim. 2

close all; clear
x=linspace(0,80,50);
t=linspace(0,50,20);
[X,T]=meshgrid(x,t);
U=100*sin(3*pi*X./80).*exp(-0.01607*T);
surf(x,t,U); shg
colorbar
figure
plot(x,U); shg
'Ajanhetket: ',t     

% Esim. 3

close all; clear
L=2; c=1; 
b=4*L/pi^2*[1,0,-1/3^2,0,1/5^2,0,-1/7^2,0,1/9^2,0,-1/11^2];
T=linspace(0,1,10);
[U,x,t]=lamposin(L,c,b,T,40);
%[U,x,t]=lamposin(L,c,b,T,100);  % Kuva ei silmämääräisesti muutu lainkaan.
surf(x,t,U); shg
colorbar
figure
plot(x,U); 
title('Lämpötilat 0.1 s. välein 0:sta 1:een')
X=[0 L/2 L]; Y=[0 L/2 0]
hold on
plot(X,Y,'r'); shg


