minor editing of current code of displaying 3 timestamps to displaying multiple timestamps with for and pause command

1 次查看(过去 30 天)
I did this code that displays the first 3 timestamps for temperature difference, but now i need to use for loop and pause command to change it to show multiple timestamps and how it changes each step with the pause command. I have edited it but its still wrong.
%%first 3 timestamps
n=5;
dt=1000;
k=15;
rho=8055;
Cp=480;
deltax=1.6;
x=linspace(-0.8,0.8,n);
deltax = x(2)-x(1);
T=zeros(1,5); T(1,1)=40; T(1,5)=40;
figure
plot(x,T)
title('Temperature distribution for first 3 timestamps')
xlabel('x')
ylabel('Temperature')
hold on;
T_old = T;
T(2:end-1) = T_old(2:end-1) + k/(rho*Cp) * dt * (T_old(3:end)-2*T_old(2:end-1)+T_old(1:end-2))/(deltax.^2);
plot(x,T)
hold on
T_old = T;
T(2:end-1) = T_old(2:end-1) + k/(rho*Cp) * dt * (T_old(3:end)-2*T_old(2:end-1)+T_old(1:end-2))/(deltax.^2);
plot(x,T)
%%multiple timestamps
n=5;
dt=1000;
k=15;
rho=8055;
Cp=480;
deltax=1.6;
x=linspace(-0.8,0.8,n);
deltax = x(2)-x(1);
T=zeros(1,5); T(1,1)=40; T(1,5)=40;
t = 0;
dt = 1000;
for c = 1:t
for r = 1:t
T(r,c) = 1/(r+c-1);
pause(1000);
end
end
figure
plot(x,T)
title('Temperature distribution for first multiple timestamps')
xlabel('x')
ylabel('Temperature')

回答(1 个)

Ramtej
Ramtej 2023-9-14
Hi,
I assume that you are trying to plot Temperature distribution at every time step with dt=1000.
Since, I donot know your implementation details, I'm suggesting a possible workaround on how to plot your data at each time step.
%% define a figure and assign title, xlabel and ylabel to it
fig = figure;
title('Temperature distribution for first multiple timestamps')
xlabel('x')
ylabel('Temperature');
hold on % make sure use turn on the hold to plot the subsequent plots into the same figure
for t = 0:1000:350000
% calculate or update the temperature distrubution at the time t
T = yourFunction(t);
plot(fig,x,T);
pause(1); % pause the plotting for a second to observer the distribution at the time t
end
hold off
Hope this resolves your query!

类别

Help CenterFile Exchange 中查找有关 Annotations 的更多信息

产品


版本

R2020b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by