Plot at t = 0, 0.001, 0.01, and 10 in one plot

Question

Jennifer Yang 2018-8-1

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/412996-plot-at-t-0-0-001-0-01-and-10-in-one-plot

评论： Jennifer Yang 2018-8-1

Hello,

I'm solving for a time dependent pde. I was wondering how do I plot multiple lines on the same plot at specific time points (ex. at t=0,0.001,0.01,10)?

function simple_time_dependent_pdepe
    clear all; close all; clc;
    %
    D_ij = 1*10^-6; %Diffusion coefficient D (3.0*10^-7 cm^2/s -> 30 um^2/s)
    L0 = 1; %c0 [nM]
    x_f =0.02; %Length of domain [um]
    maxt = 10; %Max simulation time [s]
    %
    m = 0; %Parameter corresponding to the symmetry of the problem
    x = linspace(0,x_f,100); %xmesh
    t = linspace(0,maxt,100); %tspan
    %
    sol = pdepe(m,@DiffusionPDEfun,@DiffusionICfun,@DiffusionBCfun,x,t,[]);
    u = sol;
    %
    % Plotting
    hold all
    for n = linspace(1,length(t),10)
        plot(x,sol(n,:),'LineWidth',2)
    end
    title('Time Dependent')
    xlabel('Distance (\mum)')
    ylabel('Concentration (nM)')
    axis([0 x_f 0 L0])
    %
    function [c,f,s] = DiffusionPDEfun(x,t,u,dudx)
        D = D_ij;
        %
        %Rate constants
        k_1 = 0.25;
        R_L = -k_1.*u;
        %
        % PDE
        c = 1;
        f = D_ij.*dudx;
        s = R_L;
    end
    %
    function u0 = DiffusionICfun(x)
        u0 = 0;
    end
    %
    function [pl,ql,pr,qr] = DiffusionBCfun(xl,ul,xr,ur,t)
        c0 = L0;
        pl = ul-c0;
        ql = 0;
        pr = 0;
        qr = 1;
    end
end

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

请先登录，再进行评论。

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

Answer 1

KSSV 2018-8-1

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/412996-plot-at-t-0-0-001-0-01-and-10-in-one-plot#answer_331147

在 MATLAB Online 中打开

YOu should involve, space also while plotting. How about this approach?

[t,x,sol] =simple_time_dependent_pdepe() ;
[X,T] = meshgrid(x,t) ;
figure
surf(X,Y,sol)
ti = [0,0.001,0.01,10] ;
 xi = x ;
[Xi,Ti] = meshgrid(x,ti) ;
soli = interp2(X,T,sol,Xi,Ti) ;
figure(2)
surf(Xi,Ti,soli)