Why is my plot not showing?

Question

William Clark 2020-12-14

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/693055-why-is-my-plot-not-showing

回答： Star Strider 2020-12-14

I am trying to have Matlab create this plot, and no errors pop up in my code, but for some reason the figure will not show up and I don't know why.

tempdata1 = readcell('2020_Temp.xlsx','Range','A2:A201'); %Gather St. George data
tempSG = [tempdata1'];
tempStGeorge = tempSG(1:200);
function myODE
    tspan = [1 200];
    
    S0 = [226513, 0, 29, 0, 0, 0, 0, 0];
    
    [t S] = ode45(@g,tspan,S0);
    
    figure
    plot(t,S(:,3),'-b')
    legend('Active Cases')
    axis([1 200 0 1120])
    title('St. George')
end
function dSdt = g(t,S)
    alpha=0.5;   
    n=0.5;        
    s=1/5.1;      
    g_I=1/7;     
    g_A=1/7;     
    g_H=1/14;     
    p=0.025;      
    d=0.015;      
    lambda=20;   
    r=0.0003;        
    omega=0.5;    
    epsilon=0.5;
    
    dSdt = [-f(t)*(S(3)+n*S(4))*S(1)-f(t)*r*S(8)*S(1)
        f(t)*(S(3)+n*S(4))*S(1)-s*S(2)
        s*alpha*S(2)-p*S(3)-g_I*S(3)
        s*(1-alpha)*S(2)-g_A*S(4)
        p*S(3)-d*S(5)-g_H*S(5)
        g_I*S(3)+g_A*S(4)+g_H*S(5)
        d*S(5)
        lambda*epsilon*(S(3)+n*S(4))-omega*S(8)];
end
function beta = f(t)
    if (t <= 200)
        beta = (78-tempStGeorge(t+1))*(0.006)+0.18;
    end
end