代码汇总.m

一个小小的c语言编程

    'unit','normalized','position',[0.79,0.13,0.15,0.04],...
    'string','Z功率谱分析','callback','gonglvpu_z(X3)');

set(h_menu1_sub1_slider_sigma,'callback',[...			
   'sigma=get(gcbo,''value'');',...		
   'florenz_xy(h_menu1_sub1_text_sigma,h_menu1_sub1_str_sigma,sigma,h_menu1_sub1_text_b,h_menu1_sub1_str_b,b,h_menu1_sub1_text_r,h_menu1_sub1_str_r,r,X0,Y0,Z0,t0,tf)']);
set(h_menu1_sub1_slider_b,'callback',[...			
   'b=get(gcbo,''value'');',...		
   'florenz_xy(h_menu1_sub1_text_sigma,h_menu1_sub1_str_sigma,sigma,h_menu1_sub1_text_b,h_menu1_sub1_str_b,b,h_menu1_sub1_text_r,h_menu1_sub1_str_r,r,X0,Y0,Z0,t0,tf)']);
set(h_menu1_sub1_slider_r,'callback',[...			
   'r=get(gcbo,''value'');',...		
   'florenz_xy(h_menu1_sub1_text_sigma,h_menu1_sub1_str_sigma,sigma,h_menu1_sub1_text_b,h_menu1_sub1_str_b,b,h_menu1_sub1_text_r,h_menu1_sub1_str_r,r,X0,Y0,Z0,t0,tf)']);


function my_menu1_sub2                      %y-z平面轨线图和动态轨线图
global X1 X2 X3 t X0 Y0 Z0 t0 tf r sigma b time            %全局变量必须这样定义
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  
figure(4)
clf reset,
set(gcf,'unit','normalized','position',[0.002,0.01,0.94,0.92]);
set(gcf,'defaultuicontrolunits','normalized');
set(gcf,'defaultuicontrolfontsize',12);
set(gcf,'defaultuicontrolfontname','隶书');
set(gcf,'defaultuicontrolhorizontal','left');
str='lorenz function Y-Z phase';
set(gcf,'name',str,'numbertitle','off');	
h_axes=axes('position',[0.03,0.15,0.7,0.8]);	 
set(h_axes,'xlim',[-20,20],'ylim',[-30,30],'zlim',[0,50]);
h_menu1_sub2_str_sigma='sigma=';h_menu1_sub2_str_b='b=';h_menu1_sub2_str_r='r=';

h_menu1_sub2_text_1=uicontrol(gcf,'style','text',...
   'unit','normalized','position',[0.77,0.89,0.21,0.03],...
   'horizontal','left','string',{'y,z平面对应的奇点值'});
h_menu1_sub2_text_strangex=uicontrol(gcf,'style','text',...	       
   'unit','normalized','position',[0.77,0.85,0.10,0.03],...
   'horizontal','left','string',{'r*b*x/(b+x^2)'});
h_menu1_sub2_text_strangey=uicontrol(gcf,'style','text',...	       
   'unit','normalized','position',[0.88,0.85,0.10,0.03],...
  	'horizontal','left','string',{'r*x^2/(b+x^2)'});

h_menu1_sub2_text_2=uicontrol(gcf,'style','text',...
   'unit','normalized','position',[0.77,0.81,0.21,0.03],...
   'horizontal','left','string',{'焦-结分岔点'});
h_menu1_sub2_text_bifucationx1=uicontrol(gcf,'style','text',...	       
   'unit','normalized','position',[0.77,0.77,0.07,0.05],...
  	'horizontal','left','string',['x1= ',sprintf('%1.5g\',(b-1)/2)]);

h_menu1_sub2_text_3=uicontrol(gcf,'style','text',...
   'unit','normalized','position',[0.77,0.73,0.21,0.03],...
   'horizontal','left','string',{'结-焦分岔点'});
h_menu1_sub2_text_bifucationx2=uicontrol(gcf,'style','text',...	       
   'unit','normalized','position',[0.77,0.69,0.07,0.05],...
  	'horizontal','left','string',['x2= ',sprintf('%1.5g\',-(b-1)/2)]);

h_menu1_sub2_text_2=uicontrol(gcf,'style','text',...	
   'position',[0.77,0.64,0.21,0.02],...
   'string',{'可更改参数设置'});

h_menu1_sub2_text_sigma=uicontrol(gcf,'style','text',...	
   'position',[0.77,0.61,0.21,0.03],...
   'string',['sigma=',sprintf('%1.4g\',10)]);
h_menu1_sub2_slider_sigma=uicontrol(gcf,'style','slider',...
   'position',[0.77,0.57,0.21,0.03],...
   'max',20,'min',5,...		
   'sliderstep',[0.01,0.05],...
   'Value',10);	

h_menu1_sub2_text_b=uicontrol(gcf,'style','text',...	
   'position',[0.77,0.53,0.21,0.03],...
   'string',['b=',sprintf('%1.4g\',8/3)]);
h_menu1_sub2_slider_b=uicontrol(gcf,'style','slider',...
   'position',[0.77,0.49,0.21,0.03],...
   'max',20,'min',2,...		
   'sliderstep',[0.01,0.05],...
   'Value',8/3);

h_menu1_sub2_text_r=uicontrol(gcf,'style','text',...	
   'position',[0.77,0.45,0.21,0.03],...
   'string',['r=',sprintf('%1.4g\',28)]);
h_menu1_sub2_slider_r=uicontrol(gcf,'style','slider',...
   'position',[0.77,0.41,0.21,0.03],...
   'max',40,'min',10,...		
   'sliderstep',[0.01,0.05],...
   'Value',28);

h_h_menu1_sub2_push1=uicontrol(gcf,'style','push',...
    'unit','normalized','position',[0.79,0.32,0.15,0.04],...
    'string','Y-Z平面投影图','callback',['plot(X2,X3),','grid on,','hold on,','xlabel(''y''),','ylabel(''z'')']);
h_menu1_sub2_push2=uicontrol(gcf,'style','push',...
    'unit','normalized','position',[0.79,0.28,0.15,0.04],...
    'string','Y-Z平面投影彗星图','callback','my_comet_yz(X2,X3,0.02,time)');
h_menu1_sub2_push3=uicontrol(gcf,'style','push',...
    'unit','normalized','position',[0.79,0.21,0.15,0.04],...
    'string','X功率谱分析','callback','gonglvpu_x(X1)');
h_menu1_sub2_push4=uicontrol(gcf,'style','push',...
    'unit','normalized','position',[0.79,0.17,0.15,0.04],...
    'string','Y功率谱分析','callback','gonglvpu_y(X2)');
h_menu1_sub2_push5=uicontrol(gcf,'style','push',...
    'unit','normalized','position',[0.79,0.13,0.15,0.04],...
    'string','Z功率谱分析','callback','gonglvpu_z(X3)');

set(h_menu1_sub2_slider_sigma,'callback',[...			
   'sigma=get(gcbo,''value'');',...		
   '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)']);
set(h_menu1_sub2_slider_b,'callback',[...			
   'b=get(gcbo,''value'');',...		
   '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)']);
set(h_menu1_sub2_slider_r,'callback',[...			
   'r=get(gcbo,''value'');',...		
   '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)']);


function my_menu1_sub3                      %x-z平面轨线图和动态轨线图
global X1 X2 X3 t X0 Y0 Z0 t0 tf r sigma b time            %全局变量必须这样定义
global h_menu1_sub3_str_sigma h_menu1_sub3_str_b h_menu1_sub3_str_r
global h_menu1_sub3_text_sigma h_menu1_sub3_text_b h_menu1_sub3_text_r
figure(6)
clf reset,
set(gcf,'unit','normalized','position',[0.002,0.01,0.94,0.92]);
set(gcf,'defaultuicontrolunits','normalized');
set(gcf,'defaultuicontrolfontsize',12);
set(gcf,'defaultuicontrolfontname','隶书');
set(gcf,'defaultuicontrolhorizontal','left');
str='lorenz function X-Z phase';
set(gcf,'name',str,'numbertitle','off');	
h_axes=axes('position',[0.03,0.15,0.7,0.8]);	 
set(h_axes,'xlim',[-20,20],'ylim',[-30,30],'zlim',[0,50]);
h_menu1_sub3_str_sigma='sigma=';h_menu1_sub3_str_b='b=';h_menu1_sub3_str_r='r=';

h_menu1_sub3_text_1=uicontrol(gcf,'style','text',...
   'unit','normalized','position',[0.77,0.89,0.21,0.03],...
   'horizontal','left','string',{'x,z平面对应的奇点值'});
h_menu1_sub3_text_strangex=uicontrol(gcf,'style','text',...	       
   'unit','normalized','position',[0.77,0.85,0.03,0.03],...
   'horizontal','left','string',{'y'});
h_menu1_sub3_text_strangez=uicontrol(gcf,'style','text',...	       
   'unit','normalized','position',[0.88,0.85,0.10,0.03],...
  	'horizontal','left','string',{'y^2/b'});

h_menu1_sub3_text_2=uicontrol(gcf,'style','text',...	
   'position',[0.77,0.64,0.21,0.02],...
   'string',{'可更改参数设置'});

h_menu1_sub3_text_sigma=uicontrol(gcf,'style','text',...	
   'position',[0.77,0.61,0.21,0.03],...
   'string',['sigma=',sprintf('%1.4g\',10)]);
h_menu1_sub3_slider_sigma=uicontrol(gcf,'style','slider',...
   'position',[0.77,0.57,0.21,0.03],...
   'max',20,'min',5,...		
   'sliderstep',[0.01,0.05],...
   'Value',10);	

h_menu1_sub3_text_b=uicontrol(gcf,'style','text',...	
   'position',[0.77,0.53,0.21,0.03],...
   'string',['b=',sprintf('%1.4g\',8/3)]);
h_menu1_sub3_slider_b=uicontrol(gcf,'style','slider',...
   'position',[0.77,0.49,0.21,0.03],...
   'max',20,'min',2,...		
   'sliderstep',[0.01,0.05],...
   'Value',8/3);

h_menu1_sub3_text_r=uicontrol(gcf,'style','text',...	
   'position',[0.77,0.45,0.21,0.03],...
   'string',['r=',sprintf('%1.4g\',28)]);
h_menu1_sub3_slider_r=uicontrol(gcf,'style','slider',...
   'position',[0.77,0.41,0.21,0.03],...
   'max',40,'min',10,...		
   'sliderstep',[0.01,0.05],...
   'Value',28);

h_menu1_sub3_push1=uicontrol(gcf,'style','push',...
    'unit','normalized','position',[0.79,0.32,0.15,0.04],...
    'string','X-Z平面投影图','callback',['plot(X1,X3),','grid on,','hold on,','xlabel(''x''),','ylabel(''z'')']);
h_menu1_sub3_push2=uicontrol(gcf,'style','push',...
    'unit','normalized','position',[0.79,0.28,0.15,0.04],...
    'string','X-Z平面投影彗星图','callback','my_comet_xz(X1,X3,0.02,time)');
h_menu1_sub3_push3=uicontrol(gcf,'style','push',...
    'unit','normalized','position',[0.79,0.21,0.15,0.04],...
    'string','X功率谱分析','callback','gonglvpu_x(X1)');
h_menu1_sub3_push4=uicontrol(gcf,'style','push',...
    'unit','normalized','position',[0.79,0.17,0.15,0.04],...
    'string','Y功率谱分析','callback','gonglvpu_y(X2)');
h_menu1_sub3_push5=uicontrol(gcf,'style','push',...
    'unit','normalized','position',[0.79,0.13,0.15,0.04],...
    'string','Z功率谱分析','callback','gonglvpu_z(X3)');

set(h_menu1_sub3_slider_sigma,'callback',[...			
   'sigma=get(gcbo,''value'');',...		
   'florenz_xz(h_menu1_sub3_text_sigma,h_menu1_sub3_str_sigma,sigma,h_menu1_sub3_text_b,h_menu1_sub3_str_b,b,h_menu1_sub3_text_r,h_menu1_sub3_str_r,r,X0,Y0,Z0,t0,tf)']);
set(h_menu1_sub3_slider_b,'callback',[...			
   'b=get(gcbo,''value'');',...		
   'florenz_xz(h_menu1_sub3_text_sigma,h_menu1_sub3_str_sigma,sigma,h_menu1_sub3_text_b,h_menu1_sub3_str_b,b,h_menu1_sub3_text_r,h_menu1_sub3_str_r,r,X0,Y0,Z0,t0,tf)']);
set(h_menu1_sub3_slider_r,'callback',[...			
   'r=get(gcbo,''value'');',...		
   'florenz_xz(h_menu1_sub3_text_sigma,h_menu1_sub3_str_sigma,sigma,h_menu1_sub3_text_b,h_menu1_sub3_str_b,b,h_menu1_sub3_text_r,h_menu1_sub3_str_r,r,X0,Y0,Z0,t0,tf)']);


function my_menu2_sub1
global X1 X2 X3 t X0 Y0 Z0 t0 tf r sigma b time 
global h_menu2_sub1_text_bifucationz1 h_menu2_sub1_text_bifucationz2 
global h_menu2_sub1_radio_1 h_menu2_sub1_radio_2 h_menu2_sub1_radio_3
figure(8)
clf reset,
set(gcf,'unit','normalized','position',[0.002,0.01,0.94,0.92]);
set(gcf,'defaultuicontrolunits','normalized');
set(gcf,'defaultuicontrolfontsize',12);
set(gcf,'defaultuicontrolfontname','隶书');
set(gcf,'defaultuicontrolhorizontal','left');
str='lorenz function X-Y strange place phase';
set(gcf,'name',str,'numbertitle','off');	
h_axes=axes('position',[0.03,0.15,0.7,0.8]);	 
set(h_axes,'xlim',[-20,20],'ylim',[-30,30],'zlim',[0,50]);

h_menu2_sub1_text_1=uicontrol(gcf,'style','text',...
   'unit','normalized','position',[0.77,0.89,0.21,0.03],...
   'horizontal','left','string',{'x,y平面对应的奇点值'});
h_menu2_sub1_text_strangex=uicontrol(gcf,'style','text',...	       
   'unit','normalized','position',[0.77,0.85,0.07,0.03],...
   'horizontal','left','string',{'0'});
h_menu2_sub1_text_strangey=uicontrol(gcf,'style','text',...	       
   'unit','normalized','position',[0.84,0.85,0.07,0.03],...
  	'horizontal','left','string',{'0'});

h_menu2_sub1_text_2=uicontrol(gcf,'style','text',...
   'unit','normalized','position',[0.77,0.81,0.21,0.03],...
   'horizontal','left','string',{'鞍-结分岔点'});
h_menu2_sub1_text_bifucationz1=uicontrol(gcf,'style','text',...	       
   'unit','normalized','position',[0.77,0.77,0.07,0.05],...
  	'horizontal','left','string',['z1= ',sprintf('%1.5g\',r-1)]);
h_menu2_sub1_text_3=uicontrol(gcf,'style','text',...
   'unit','normalized','position',[0.77,0.73,0.21,0.03],...
   'horizontal','left','string',{'结-焦分岔点'});
h_menu2_sub1_text_bifucationz2=uicontrol(gcf,'style','text',...	       
   'unit','normalized','position',[0.77,0.69,0.07,0.05],...
  	'horizontal','left','string',['z2= ',sprintf('%1.5g\',r-1+(sigma+1)^2/(4*sigma))]);

h_menu2_sub1_radio_1=uicontrol(gcf,'style','radio',...	
   'position',[0.77,0.52,0.21,0.04],...
   'string','鞍点区',...
   'callback',[...
       'set(h_menu2_sub1_radio_1,''value'',get(h_menu2_sub1_radio_1,''max'')),',...
        'set(h_menu2_sub1_radio_2,''value'',get(h_menu2_sub1_radio_2,''min'')),',...
        'set(h_menu2_sub1_radio_3,''value'',get(h_menu2_sub1_radio_3,''min'')),',...
    'x_y_an_space']);
h_menu2_sub1_radio_2=uicontrol(gcf,'style','radio',...	
   'position',[0.77,0.42,0.21,0.04],...
   'string','稳定结点区',...
   'callback',[...
       'set(h_menu2_sub1_radio_1,''value'',get(h_menu2_sub1_radio_1,''Min'')),',...
        'set(h_menu2_sub1_radio_2,''value'',get(h_menu2_sub1_radio_2,''Max'')),',...
        'set(h_menu2_sub1_radio_3,''value'',get(h_menu2_sub1_radio_3,''Min'')),',...
    'x_y_jie_space']);

h_menu2_sub1_radio_3=uicontrol(gcf,'style','radio',...	
   'position',[0.77,0.32,0.21,0.04],...
   'string','稳定焦点区',...
   'callback',[...
       'set(h_menu2_sub1_radio_1,''value'',get(h_menu2_sub1_radio_1,''min'')),',...
        'set(h_menu2_sub1_radio_2,''value'',get(h_menu2_sub1_radio_2,''min'')),',...
        'set(h_menu2_sub1_radio_3,''value'',get(h_menu2_sub1_radio_3,''max'')),',...
    'x_y_jiao_space']);

h_menu2_sub1_push=uicontrol(gcf,'style','push',...
    'unit','normalized','position',[0.79,0.25,0.15,0.04],...
    'string','X-Y轨线动态绘图','callback','my_comet_strange_xy(X1,X2,0.02,time)');


function my_menu2_sub2
global X1 X2 X3 t X0 Y0 Z0 t0 tf r sigma b time 
global h_menu2_sub2_text_bifucationx1 h_menu2_sub2_text_bifucationx2 
global h_menu2_sub2_radio_1 h_menu2_sub2_radio_2 h_menu2_sub2_radio_3
figure(10)
clf reset,
set(gcf,'unit','normalized','position',[0.002,0.01,0.94,0.92]);
set(gcf,'defaultuicontrolunits','normalized');
set(gcf,'defaultuicontrolfontsize',12);
set(gcf,'defaultuicontrolfontname','隶书');
set(gcf,'defaultuicontrolhorizontal','left');
str='lorenz function Y-Z strange place phase';
set(gcf,'name',str,'numbertitle','off');	
h_axes=axes('position',[0.03,0.15,0.7,0.8]);	 
set(h_axes,'xlim',[-20,20],'ylim',[-30,30],'zlim',[0,50]);

h_menu2_sub2_text_1=uicontrol(gcf,'style','text',...
   'unit','normalized','position',[0.77,0.89,0.21,0.03],...
   'horizontal','left','string',{'y,z平面对应的奇点值'});
h_menu2_sub2_text_strangey=uicontrol(gcf,'style','text',...	       
   'unit','normalized','position',[0.77,0.85,0.07,0.03],...
   'horizontal','left','string',{'r*b*x^2/(b+x^2)'});
h_menu2_sub2_text_strangez=uicontrol(gcf,'style','text',...	       
   'unit','normalized','position',[0.84,0.85,0.07,0.03],...
  	'horizontal','left','string',{'r*x^2/(b+x^2)'});
h_menu2_sub2_text_2=uicontrol(gcf,'style','text',...
   'unit','normalized','position',[0.77,0.81,0.21,0.03],...
   'horizontal','left','string',{'焦-结分岔点'});
h_menu2_sub2_text_bifucationx1=uicontrol(gcf,'style','text',...	       
   'unit','normalized','position',[0.77,0.77,0.07,0.05],...
  	'horizontal','left','string',['x1= ',sprintf('%1.5g\',(b-1)/2)]);
h_menu2_sub2_text_3=uicontrol(gcf,'style','text',...
   'unit','normalized','position',[0.77,0.73,0.21,0.03],...
   'horizontal','left','string',{'结-焦分岔点'});
h_menu2_sub2_text_bifucationx2=uicontrol(gcf,'style','text',...	       
   'unit','normalized','position',[0.77,0.69,0.07,0.05],...
  	'ho