MATLAB使用歐拉Euler法求解微分方程組 源程序代碼

1. clear;clc
   
2. c=2/3;%設(shè)置c的值
   
3. x(1)=0.1;%設(shè)置x初值為0.1
   
4. y(1)=0.3;%設(shè)置y初值為0.3
   
5. h=0.05; %設(shè)置步長為0.05
   
6. 
   
7. %for循環(huán)：根據(jù)Euler法，求解微分方程組
   
8. for i=1:1000
   
9. x(i+1)=x(i)+h*(x(i)*(c-x(i)/y(i)));
   
10. y(i+1)=y(i)+h*(y(i)*(1-y(i))-x(i)*y(i));
    
11. end
    
12. 
    
13. t=0:h:1000*h; %計(jì)算時(shí)間
    
14. plot(t,x) %繪制x曲線
    
15. hold on
    
16. plot(t,y,'r') %繪制y曲線
    
17. xlabel('time')%設(shè)置x軸的標(biāo)簽
    
18. ylabel('value') %設(shè)置y軸的標(biāo)簽
    
19. legend({'x','y'})%設(shè)置圖例
    
20. title('time evolution plot')%設(shè)置圖形的標(biāo)題
    
21. figure
    
22. plot(x,y) %繪制phase圖
    
23. title('phase plane plot')%設(shè)置標(biāo)題
    
24. xlabel('x') %設(shè)置x軸的標(biāo)簽
    
25. ylabel('y') %設(shè)置y軸的標(biāo)簽