Plot approach to steady state solution

Question

Jonathan Melia 2020-5-17

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/526125-plot-approach-to-steady-state-solution

评论： Jonathan Melia 2020-5-20

figure(1)
title('Exact Steady State Temperature Profile')
hold on
plot(r_profile,Up_exact,'r-');
xlabel('r')
ylabel('T')
for i = 1:ceil(tau/dt)
    Mplot(:,i) = U;
    U = Atinv*(D*U+dt*B);  
    [Uedge, redge] = GetBHboundary(theMesh,3,U);
    figure(2)
    plot(redge, Uedge,'k')
    axis([1 10 0 0.4])
    title('FE Temperature Profile')
    drawnow
    pause(0.05)
end

I have the following code which plots my steady state temperature profile in figure 1 and an animated transient temperature profile in figure 2. Is there a way i can have them both in one plot showing how the transient solution approaches the steady state?

Thanks :)

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

请先登录，再进行评论。

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

Answer 1

Asvin Kumar 2020-5-20

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/526125-plot-approach-to-steady-state-solution#answer_434132

在 MATLAB Online 中打开

Can you give this a try?

figure
title('Exact Steady State Temperature Profile')
hold on
plot(r_profile,Up_exact,'r-');
xlabel('r')
ylabel('T')
ax = gca;
for i = 1:ceil(tau/dt)
    Mplot(:,i) = U;
    U = Atinv*(D*U+dt*B);  
    [Uedge, redge] = GetBHboundary(theMesh,3,U);
    % figure(2)
    plot(ax, redge, Uedge,'k')
    axis([1 10 0 0.4])
    title('FE Temperature Profile')
    drawnow
    pause(0.05)
end