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

Question

Ellie Matlab 2022-6-28

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1749760-minor-editing-of-current-code-of-displaying-3-timestamps-to-displaying-multiple-timestamps-with-for

回答： Ramtej 2023-9-14

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')

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Ramtej 2023-9-14

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1749760-minor-editing-of-current-code-of-displaying-3-timestamps-to-displaying-multiple-timestamps-with-for#answer_1309941

在 MATLAB Online 中打开

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