数据必须为可转换为双精度值的数值、日期时间、持续时间或数组。

sol=dsolve('DP1=lambda1-(lambda1+Mu1)*P1-lambda1*P2','DP2=lambda2-(lambda2+Mu2)*P2-lambda2*P1','P1(0)=0','P2(0)=0','x');
P1=simplify(sol.P1);
P2=simplify(sol.P2);
P0=1-P1-P2;
lambda1=1;
lambda2=1;
Mu1=2;
Mu2=4;
x=0:0.01:10;
P0=1-P1-P2;

plot（x,P0)

plot这一步的时候会报错，P0类型为1*1sym

仅供参考

1. sol=dsolve('DP1=lambda1-(lambda1+Mu1)*P1-lambda1*P2','DP2=lambda2-(lambda2+Mu2)*P2-lambda2*P1','P1(0)=0','P2(0)=0','x');
   
2. P1=simplify(sol.P1);
   
3. P2=simplify(sol.P2);
   
4. P0=1-P1-P2;
   
5. % lambda1=1;
   
6. % lambda2=1;
   
7. % Mu1=2;
   
8. % Mu2=4;
   
9. P1 = subs(P0,{'lambda1','lambda2','Mu1','Mu2'},{1,1,2,4});
   
10. x1=0:0.1:10;
    
11. for i=1:length(x1)
    
12.     P2(i)=double(subs(P1,'x',x1(i)));
    
13. end
    
14. plot(x1,P2)

复制代码

仅供参考

1. sol=dsolve('DP1=lambda1-(lambda1+Mu1)*P1-lambda1*P2','DP2=lambda2-(lambda2+Mu2)*P2-lambda2*P1','P1(0)=0','P2(0)=0','x');
   
2. P1=simplify(sol.P1);
   
3. P2=simplify(sol.P2);
   
4. P0=1-P1-P2;
   
5. % lambda1=1;
   
6. % lambda2=1;
   
7. % Mu1=2;
   
8. % Mu2=4;
   
9. P1 = subs(P0,{'lambda1','lambda2','Mu1','Mu2'},{1,1,2,4});
   
10. x1=0:0.1:10;
    
11. for i=1:length(x1)
    
12.     P2(i)=double(subs(P1,'x',x1(i)));
    
13. end
    
14. plot(x1,P2)

复制代码

20141303 发表于 2022-5-12 08:35
仅供参考

感谢回答