florenz_yz.m

一个小小的c语言编程

function florenz_yz(h_menu1_sub2_text_sigma,h_menu1_sub2_str_sigma,sigma,h_menu1_sub2_text_b,h_menu1_sub2_str_b,b,h_menu1_sub2_text_r,h_menu1_sub2_str_r,r,X0,Y0,Z0,t0,tf)
global X1 X2 X3 t                                                        %解方程并画出y-z的黑色平面轨线图                
global h_menu1_sub2_str_sigma h_menu1_sub2_str_b h_menu1_sub2_str_r
global h_menu1_sub2_text_sigma h_menu1_sub2_text_b h_menu1_sub2_text_r
global h_menu1_sub2_text_bifucationx1 h_menu1_sub2_text_bifucationx2

cla,set(h_menu1_sub2_text_sigma,'string',[h_menu1_sub2_str_sigma,sprintf('%1.4g\',sigma)]);
cla,set(h_menu1_sub2_text_b,'string',[h_menu1_sub2_str_b,sprintf('%1.4g\',b)]);
cla,set(h_menu1_sub2_text_r,'string',[h_menu1_sub2_str_r,sprintf('%1.4g\',r)]);
cla,set(h_menu1_sub2_text_bifucationx1,'string',['x1=',sprintf('%1.4g\',(b-1)/2)]);
cla,set(h_menu1_sub2_text_bifucationx2,'string',['x2=',sprintf('%1.4g\',-(b-1)/2)]);

tspan=[t0,tf];f0=[X0;Y0;Z0];
[t,X]=ode45('DYDt2',tspan,f0,[],sigma,b,r);
X1=X(:,1);X2=X(:,2);X3=X(:,3);
x1=-20;x2=20;y1=-30;y2=30;
grid on ,            %设置网格线    
hold on             %保持网格线
axis auto,          %自动容纳图形用axis auto，plot使用
plot(X2,X3,'k-'),hold on,