what is the meaning of this error, "??? Undefined function or method 'Restart' for input arguments of type 'char' when using Matlab?

Question

ELISHA ANEBI 2023-7-8

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1993343-what-is-the-meaning-of-this-error-undefined-function-or-method-restart-for-input-arguments-o

评论： ELISHA ANEBI 2023-7-11

6 个评论
显示 4更早的评论隐藏 4更早的评论

ELISHA ANEBI 2023-7-8

编辑：Cris LaPierre 2023-7-8

在 MATLAB Online 中打开

Edit - fixed some apparent typos in the code

% The Matlab Codes for the forward bifurcation diagram
Rev_value=0.018:0.01:4;
Root_array=zeros (length (Rev_value), 2);
qI=0.001923; qA=0.00000004013; etaA=0.1213; etaQ=0.003808;  w=0.5925;
Lambda=7898250;theta=0.022; 
delta=0.125; mu=0.01119; pi=0.464360344; deltaQ=6.847e-4;  beta=1.086e-1;
qE=1.8113e-4; rhoQ=0.0815; a=0.16255; k=0.15; v1=0.71;  v2=0.29;alpha=0.57e-1; deltaI=0.00000000223; rhoA=0.1; rhoI=0.0666666; 
c1=mu+v1; c2=1-w; c3=mu+alpha+v2; c4=1-k; c5=qE+delta+mu; c6=rhoA+mu; c7=delta*(1-a); c8=rhoI+qI+deltaI+mu; c9=rhoA+deltaQ+mu; 
k1=delta*a*c8*c9+c7*k*qI*rhoQ+c8*k*qE*rhoQ; k2=c7*qI+c8*qE; b1=1-theta;
p2=(b1*c1*c2)*(a*c8*c9*delta*rhoA*pi+c4*c6*c7*qI*rhoQ*pi+c4*c6*c8*qE*rhoQ*pi+c8*k*qE*rhoA*rhoQ*pi+c7*k*qI*rhoA*rhoQ*pi)+(c2*c8*theta*v1*pi)*(c4*c6*qE*rhoQ+a*c9*delta*rhoA+k*qE*rhoA*rhoQ)+c2*c6*c7*c9*rhoI*theta*v1*pi+c3*c4*c6*c7*qI*rhoQ*theta*pi+c3*c7*k*qI*rhoA*rhoQ*theta*pi+(c3*c8)*(k*qE*rhoA*rhoQ*theta*pi+c4*c6*qE*rhoQ*theta*pi+a*c9*delta*rhoA*theta*pi)+b1*c4*c6*c7*qI*rhoQ*v2*pi+b1*c4*c6*c8*qE*rhoQ*v2*pi+b1*c7*k*qI*rhoA*rhoQ*v2*pi+b1*c8*k*qE*rhoA*rhoQ*v2*pi+a*b1*c8*c9*delta*rhoA*v2*pi+b1*c6*c7*c9*rhoI*v2*pi+c3*c6*c7*c9*rhoI*theta*pi+alpha*b1*c5*c6*c8*c9*pi+c2*c7*k*qI*rhoA*rhoQ*theta*v1*pi+c2*c4*c6*c7*qI*rhoQ*theta*v1*pi;
G2=b1*c1*c2*c6*mu*k2*pi+c2*theta*v1*c6*mu*k2*pi+b1*v2*c6*mu*k2*pi+c3*theta*c6*mu*k2*pi+c2*theta*c5*pi+b1*c1*c2*k1*pi*mu+c2*theta*v1*k1*pi*mu+b1*v2*k1*pi*mu+c3*theta*k1*pi*mu+b1*c1*c2*c6*c9*mu*c7*pi+c2*theta*v1*c6*c9*mu*c7*pi+b1*v2*c6*c9*mu*c7*pi+c3*theta*c6*c9*mu*c7*pi+b1*c1*c2*pi+c2*theta*v1*pi+b1*v2*pi+c3*theta*pi+c5*pi*b1-b1*c2*beta*etaQ*pi*mu*c6*k2-c2*theta*beta*etaQ*pi*mu*c6*k2-b1*c2*beta*c6*c9*c7*pi*mu-c2*theta*beta*c6*c9*c7*pi*mu-b1*c2*beta*etaA*pi*mu*k1-c2*theta*beta*etaA*pi*mu*k1;
U=(c1*c3-v2*v1)*mu/(alpha*c1*pi);
H1=b1*c2*c4*c6*c7*qI*rhoQ*pi+b1*c2*c4*c6*c8*qE*rhoQ*pi+b1*c2*c7*k*qI*rhoA*rhoQ*pi+b1*c2*c8*k*qE*rhoA*rhoQ*pi+c2*c4*c6*c7*qI*rhoQ*theta*pi+c2*c4*c6*c8*qE*rhoQ*theta*pi+c2*c7*k*qI*rhoA*rhoQ*theta*pi+c2*c8*k*qE*rhoA*rhoQ*theta*pi+b1*c2*c6*c7*c9*rhoI*pi+c2*c6*c7*c9*rhoI*theta*pi+a*c2*c8*c9*delta*rhoA*theta*pi+a*b1*c2*c8*c9*delta*rhoA*pi+b1*c2*c6*mu*k2*pi+c2*c6*mu*k2*theta*pi+b1*c2*k1*pi*mu+c2*theta*k1*pi*mu+b1*c2*c6*c9*mu*c7*pi+c2*theta*c6*c9*mu*c7*pi+b1*c2*pi+c2*pi*theta;
H2=p2+G2;
H3=alpha*c5*c6*c8*c9*theta*v1*pi+(alpha*b1*c1*c5*c6*c8*c9*pi)*(1-U*Rev_value)+(c5*pi)*(b1*v2+c3*theta+b1*c1+v1*theta)-(c2*theta*v1*beta*etaA*pi*mu+c3*theta*beta*etaA*pi*mu)*(k1)-(c2*theta*v1*beta*pi*mu+c3*theta*beta*pi*mu)*(c6*c9*c7)-(k2)*(c2*theta*v1*pi*mu*c6+c3*theta*pi*mu*c6)*(beta*etaQ);
hold on
for i=1:1:length(Rev_value)
    Rev=Rev_value(i);
    % bifurcation parameter
    %Coefficients of quadratic equation H1, H2, H3
    
    p=[H1,H2,H3];
    
    r =roots(p);
    len=length(r);
    for t=1:1:len
        if (imag(r(t))~=0) || (real(r(t))<0)
            Root_array(i,t)=0;
        else
            Root_array(i,t)=r(t);
        end
    end
end
f=1;
while (Root_array(f,1)==0 && Root_array(f,2)==0)
    f=f+1;
end
Rev_value_Cr=f;
for j=Rev_value_Cr:1:length(Rev_value)
    Root_array(j,:) =sort(Root_array(j,:));
end
f1=Rev_value_Cr;
while (Root_array(f1,1)~=0)
    f1=f1+1;
end
Rev_value_Cr2=f1;
Zero_1st=Rev_value(1,1:Rev_value_Cr2-1);
y_zero=zeros(1,length(Zero_1st));
Unstable =Rev_value(1,Rev_value_Cr:length(Rev_value));
figure (1)
plot(Unstable, Root_array(Rev_value_Cr:length(Rev_value),2),'r','LineWidth',3)
figure (2)
plot(Unstable,Root_array(Rev_value_Cr:length(Rev_value),2),'b','LineWidth',3)
ylabel('Force of Infection,(\lambda)','FontSize',12)
xlabel('reproduction number, R_o','FontSize',12)
hold off
figure (3)
plot(Rev_value,Root_array(:,1),'r',Rev_value,Root_array(:,2),'r',Rev_value,Root_array(:,2),'b','LineWidth',4)
ylabel('Force of Infection, \lambda','FontSize',12)
xlabel('reproduction number, R_o','FontSize',12)

ELISHA ANEBI 2023-7-8

This code is not ploting any graph as expected. Please what do you think is wrong?

Image Analyst 2023-7-8

It plots what you told it to plot. If you're plotting the wrong things, then probably your equations are wrong.

请先登录，再进行评论。

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

Answer 1

Sandeep Mishra 2023-7-10

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1993343-what-is-the-meaning-of-this-error-undefined-function-or-method-restart-for-input-arguments-o#answer_1270183

在 MATLAB Online 中打开

Hello ELISHA,

I tried running your commented code and found some different issue in your while loop.

Your code is running out of Root_array range, which is giving me error, you can stop that by checking width of Root_array as below

while (f<width(Root_array) && Root_array(f,1)==0 && Root_array(f,2)==0)
    f=f+1;
end
while (f1<width(Root_array) && Root_array(f1,1)~=0)
    f1=f1+1;
end

I also found that you are getting no plot in figure 1 and figure 2 because the plot inputs are empty (0×1 or 1×0 size)

Unstable =Rev_value(1,Rev_value_Cr:length(Rev_value))

unstable = 1×0 empty double row vector

Root_array(Rev_value_Cr:length(Rev_value),2)

Root_array = 0×1 empty double column vector

You can debug your code to fix these issues.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

ELISHA ANEBI 2023-7-11

Sandeep Mishra,

Thank you for your input. I copied the code you sent and replaced mine but it didn't still work. Please find below my equations and necessary parameters. I am ploting bifurcation graph- "Force of infection" against the "basic reproduction no Rev".

Please help me plot this graph.

% The Matlab Codes for the forward bifurcation diagram

Rev_value=0.018:0.01:4;

Root_array=zeros (length (Rev_value), 2);

qI=0.001923; qA=0.00000004013; etaA=0.1213; etaQ=0.003808; w=0.5925;

Lambda=0.1598643e-7; theta=0.022;

delta=0.125; mu=0.01119; pi=0.464360344; deltaQ=6.847e-4; beta=1.086e-1;

qE=1.8113e-4; rhoQ=0.0815; a=0.16255; k=0.15; v1=0.71; v2=0.29;alpha=0.57e-1; deltaI=0.00000000223; rhoA=0.1; rhoI=0.0666666;

c1=mu+v1; c2=1-w; c3=mu+alpha+v2;c4=1-k; c5=qE+delta+mu; c6=rhoA+mu; c7=delta*(1-a); c8=rhoI+qI+deltaI+mu; c9=rhoA+deltaQ+mu;

b1=1-theta;

B=delta*a*c8*c9+c7*k*qI*rhoQ+c8*k*qE*rhoQ;

G=qI*c7*c6+qE*c8*c6;

H1=c2*c5*c6*c8*c9;

H2=c5*c6*c8*c9*(c3*+c1*c2)-(b1*c2+c2*theta)*(c6*c9*c7*pi*beta+pi*G*beta*etaQ+pi*B*beta*etaA);

H3=c5*c6*c8*c9*(c3*c1-v1*v2)*(1-b1*Rev_value)-(c3*theta+c2*v1*theta)*(c6*c9*c7*pi*beta+pi*G*beta*etaQ+pi*B*beta*etaA);

hold on

for i=1:1:length(Rev_value);

Rev=Rev_value(i);

% bifurcation parameter

%Coefficients of quadratic equation H1, H2, H3

p=[H1,H2,H3];

r =roots(p);

len=length(r);

for t=1:1:len

请先登录，再进行评论。

what is the meaning of this error, "??? Undefined function or method 'Restart' for input arguments of type 'char' when using Matlab?

6 个评论
显示 4更早的评论隐藏 4更早的评论

回答（1 个）

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

另请参阅

类别

标签

Community Treasure Hunt

what is the meaning of this error, "??? Undefined function or method 'Restart' for input arguments of type 'char' when using Matlab?

6 个评论 显示 4更早的评论隐藏 4更早的评论

回答（1 个）

1 个评论 显示 -1更早的评论隐藏 -1更早的评论

另请参阅

类别

标签

Community Treasure Hunt

6 个评论
显示 4更早的评论隐藏 4更早的评论

1 个评论
显示 -1更早的评论隐藏 -1更早的评论