lorenz3d.mdl

运用非线性系统产生混沌吸引子的MATLAB源代码。

	  Points		  [50, 0]
	  DstBlock		  "-Bz+xy"
	  DstPort		  2
	}
	Line {
	  SrcBlock		  "Mux1"
	  SrcPort		  1
	  DstBlock		  "xyz"
	  DstPort		  1
	}
	Line {
	  SrcBlock		  "Rho * x"
	  SrcPort		  1
	  Points		  [45, 0]
	  DstBlock		  "Rx-y-xz"
	  DstPort		  3
	}
	Annotation {
	  Position		  [192, 42]
	  Text			  "Lorenz Attractor"
	}
      }
    }
    Block {
      BlockType		      SubSystem
      Name		      "Lorenz2"
      Ports		      [0, 1, 0, 0, 0]
      Position		      [40, 130, 80, 190]
      ShowPortLabels	      on
      MaskType		      "Lorenz Attractor System"
      MaskDescription	      "This autonomous system generates the lorentz at"
"tractor:\nxdot=sigma*(x-y); ydot=rho*x-y-x*z; zdot=-beta*z+x*y;\nParameters a"
"nd initial state can be selected from the mask."
      MaskHelp		      "This autonomous system generates the lorentz at"
"tractor.\nParameters and initial state can be selected from the mask.\n<p>\nT"
"he system equations are: \nxdot = sigma*(x-y),\nydot = rho*x-y-x*z,\nzdot = -"
"beta*z+x*y,"
      MaskPromptString	      "sigma|rho|beta|Initial State (x0, y0, y0)"
      MaskStyleString	      "edit,edit,edit,edit"
      MaskTunableValueString  "on,on,on,on"
      MaskCallbackString      "|||"
      MaskEnableString	      "on,on,on,on"
      MaskVisibilityString    "on,on,on,on"
      MaskVariables	      "s=@1;r=@2;b=@3;xyz0=@4;"
      MaskIconFrame	      on
      MaskIconOpaque	      on
      MaskIconRotate	      "none"
      MaskIconUnits	      "autoscale"
      MaskValueString	      "10|28|2+2/3|[2 2 2]"
      System {
	Name			"Lorenz2"
	Location		[95, 138, 546, 569]
	Open			off
	ModelBrowserVisibility	off
	ModelBrowserWidth	200
	ScreenColor		"white"
	PaperOrientation	"landscape"
	PaperPositionMode	"auto"
	PaperType		"usletter"
	PaperUnits		"inches"
	ZoomFactor		"100"
	AutoZoom		on
	Block {
	  BlockType		  Sum
	  Name			  "-Bz+xy"
	  Ports			  [2, 1, 0, 0, 0]
	  Position		  [150, 335, 170, 355]
	  NamePlacement		  "alternate"
	  IconShape		  "round"
	  Inputs		  "|-+"
	  SaturateOnIntegerOverflow on
	  Port {
	    PortNumber		    1
	    Name		    "zdot"
	    TestPoint		    off
	    RTWStorageClass	    "Auto"
	  }
	}
	Block {
	  BlockType		  Gain
	  Name			  "Beta * z"
	  Position		  [75, 325, 130, 365]
	  Gain			  "b"
	  SaturateOnIntegerOverflow on
	}
	Block {
	  BlockType		  Integrator
	  Name			  "Integrator"
	  Ports			  [1, 1, 0, 0, 0]
	  Position		  [245, 85, 275, 105]
	  ShowName		  off
	  ExternalReset		  "none"
	  InitialConditionSource  "internal"
	  InitialCondition	  "xyz0(1)"
	  LimitOutput		  off
	  UpperSaturationLimit	  "inf"
	  LowerSaturationLimit	  "-inf"
	  ShowSaturationPort	  off
	  ShowStatePort		  off
	  AbsoluteTolerance	  "auto"
	  Port {
	    PortNumber		    1
	    Name		    "x"
	    TestPoint		    off
	    RTWStorageClass	    "Auto"
	  }
	}
	Block {
	  BlockType		  Integrator
	  Name			  "Integrator1"
	  Ports			  [1, 1, 0, 0, 0]
	  Position		  [245, 210, 275, 230]
	  ShowName		  off
	  ExternalReset		  "none"
	  InitialConditionSource  "internal"
	  InitialCondition	  "xyz0(2)"
	  LimitOutput		  off
	  UpperSaturationLimit	  "inf"
	  LowerSaturationLimit	  "-inf"
	  ShowSaturationPort	  off
	  ShowStatePort		  off
	  AbsoluteTolerance	  "auto"
	  Port {
	    PortNumber		    1
	    Name		    "y"
	    TestPoint		    off
	    RTWStorageClass	    "Auto"
	  }
	}
	Block {
	  BlockType		  Integrator
	  Name			  "Integrator2"
	  Ports			  [1, 1, 0, 0, 0]
	  Position		  [245, 335, 275, 355]
	  ShowName		  off
	  ExternalReset		  "none"
	  InitialConditionSource  "internal"
	  InitialCondition	  "xyz0(3)"
	  LimitOutput		  off
	  UpperSaturationLimit	  "inf"
	  LowerSaturationLimit	  "-inf"
	  ShowSaturationPort	  off
	  ShowStatePort		  off
	  AbsoluteTolerance	  "auto"
	  Port {
	    PortNumber		    1
	    Name		    "z"
	    TestPoint		    off
	    RTWStorageClass	    "Auto"
	  }
	}
	Block {
	  BlockType		  Mux
	  Name			  "Mux1"
	  Ports			  [3, 1, 0, 0, 0]
	  Position		  [345, 33, 350, 407]
	  ShowName		  off
	  Inputs		  "3"
	  DisplayOption		  "bar"
	}
	Block {
	  BlockType		  Gain
	  Name			  "Rho * x"
	  Position		  [75, 254, 110, 276]
	  Gain			  "r"
	  SaturateOnIntegerOverflow on
	}
	Block {
	  BlockType		  Sum
	  Name			  "Rx-y-xz"
	  Ports			  [3, 1, 0, 0, 0]
	  Position		  [150, 210, 170, 230]
	  IconShape		  "round"
	  Inputs		  "--+"
	  SaturateOnIntegerOverflow on
	  Port {
	    PortNumber		    1
	    Name		    "ydot"
	    TestPoint		    off
	    RTWStorageClass	    "Auto"
	  }
	}
	Block {
	  BlockType		  Gain
	  Name			  "Sigma(y-x)"
	  Position		  [165, 81, 200, 109]
	  Gain			  "s"
	  SaturateOnIntegerOverflow on
	  Port {
	    PortNumber		    1
	    Name		    "xdot"
	    TestPoint		    off
	    RTWStorageClass	    "Auto"
	  }
	}
	Block {
	  BlockType		  Product
	  Name			  "x * y"
	  Ports			  [2, 1, 0, 0, 0]
	  Position		  [75, 380, 105, 400]
	  Inputs		  "2"
	  SaturateOnIntegerOverflow on
	}
	Block {
	  BlockType		  Product
	  Name			  "x * z"
	  Ports			  [2, 1, 0, 0, 0]
	  Position		  [75, 210, 105, 230]
	  Inputs		  "2"
	  SaturateOnIntegerOverflow on
	}
	Block {
	  BlockType		  Sum
	  Name			  "y - x"
	  Ports			  [2, 1, 0, 0, 0]
	  Position		  [80, 85, 100, 105]
	  IconShape		  "round"
	  Inputs		  "-+|"
	  SaturateOnIntegerOverflow on
	}
	Block {
	  BlockType		  Outport
	  Name			  "xyz"
	  Position		  [375, 213, 405, 227]
	  Port			  "1"
	  OutputWhenDisabled	  "held"
	  InitialOutput		  "[]"
	}
	Line {
	  SrcBlock		  "Mux1"
	  SrcPort		  1
	  DstBlock		  "xyz"
	  DstPort		  1
	}
	Line {
	  SrcBlock		  "x * y"
	  SrcPort		  1
	  Points		  [50, 0]
	  DstBlock		  "-Bz+xy"
	  DstPort		  2
	}
	Line {
	  SrcBlock		  "Beta * z"
	  SrcPort		  1
	  DstBlock		  "-Bz+xy"
	  DstPort		  1
	}
	Line {
	  Name			  "zdot"
	  Labels		  [-1, 1]
	  SrcBlock		  "-Bz+xy"
	  SrcPort		  1
	  DstBlock		  "Integrator2"
	  DstPort		  1
	}
	Line {
	  Name			  "x"
	  Labels		  [0, 1]
	  SrcBlock		  "Integrator"
	  SrcPort		  1
	  Points		  [0, 0; 35, 0]
	  Branch {
	    DstBlock		    "Mux1"
	    DstPort		    1
	  }
	  Branch {
	    Points		    [0, -30; -225, 0]
	    Branch {
	      Points		      [-40, 0; 0, 150]
	      Branch {
		Points			[0, 50]
		Branch {
		  Points		  [0, 120]
		  DstBlock		  "x * y"
		  DstPort		  1
		}
		Branch {
		  DstBlock		  "Rho * x"
		  DstPort		  1
		}
	      }
	      Branch {
		DstBlock		"x * z"
		DstPort			1
	      }
	    }
	    Branch {
	      DstBlock		      "y - x"
	      DstPort		      1
	    }
	  }
	}
	Line {
	  Name			  "y"
	  Labels		  [0, 1]
	  SrcBlock		  "Integrator1"
	  SrcPort		  1
	  Points		  [0, 0; 20, 0]
	  Branch {
	    DstBlock		    "Mux1"
	    DstPort		    2
	  }
	  Branch {
	    Points		    [0, -45; -140, 0]
	    Branch {
	      Points		      [-125, 0]
	      Branch {
		Points			[0, 220]
		DstBlock		"x * y"
		DstPort			2
	      }
	      Branch {
		Points			[0, -80]
		DstBlock		"y - x"
		DstPort			2
	      }
	    }
	    Branch {
	      DstBlock		      "Rx-y-xz"
	      DstPort		      1
	    }
	  }
	}
	Line {
	  Name			  "z"
	  Labels		  [0, 1]
	  SrcBlock		  "Integrator2"
	  SrcPort		  1
	  Points		  [0, 0; 25, 0]
	  Branch {
	    DstBlock		    "Mux1"
	    DstPort		    3
	  }
	  Branch {
	    Points		    [0, -45; -280, 0]
	    Branch {
	      Points		      [0, 45]
	      DstBlock		      "Beta * z"
	      DstPort		      1
	    }
	    Branch {
	      Points		      [0, -75]
	      DstBlock		      "x * z"
	      DstPort		      2
	    }
	  }
	}
	Line {
	  SrcBlock		  "x * z"
	  SrcPort		  1
	  DstBlock		  "Rx-y-xz"
	  DstPort		  2
	}
	Line {
	  Name			  "ydot"
	  Labels		  [-1, 1]
	  SrcBlock		  "Rx-y-xz"
	  SrcPort		  1
	  DstBlock		  "Integrator1"
	  DstPort		  1
	}
	Line {
	  SrcBlock		  "y - x"
	  SrcPort		  1
	  DstBlock		  "Sigma(y-x)"
	  DstPort		  1
	}
	Line {
	  Name			  "xdot"
	  Labels		  [1, 0]
	  SrcBlock		  "Sigma(y-x)"
	  SrcPort		  1
	  DstBlock		  "Integrator"
	  DstPort		  1
	}
	Line {
	  SrcBlock		  "Rho * x"
	  SrcPort		  1
	  Points		  [45, 0]
	  DstBlock		  "Rx-y-xz"
	  DstPort		  3
	}
	Annotation {
	  Position		  [192, 42]
	  Text			  "Lorenz Attractor"
	}
      }
    }
    Block {
      BlockType		      SubSystem
      Name		      "Lorenz3"
      Ports		      [0, 1, 0, 0, 0]
      Position		      [40, 220, 80, 280]
      ShowPortLabels	      on
      MaskType		      "Lorenz Attractor System"
      MaskDescription	      "This autonomous system generates the lorentz at"
"tractor:\nxdot=sigma*(x-y); ydot=rho*x-y-x*z; zdot=-beta*z+x*y;\nParameters a"
"nd initial state can be selected from the mask."
      MaskHelp		      "This autonomous system generates the lorentz at"
"tractor.\nParameters and initial state can be selected from the mask.\n<p>\nT"
"he system equations are: \nxdot = sigma*(x-y),\nydot = rho*x-y-x*z,\nzdot = -"
"beta*z+x*y,"
      MaskPromptString	      "sigma|rho|beta|Initial State (x0, y0, y0)"
      MaskStyleString	      "edit,edit,edit,edit"
      MaskTunableValueString  "on,on,on,on"
      MaskCallbackString      "|||"
      MaskEnableString	      "on,on,on,on"
      MaskVisibilityString    "on,on,on,on"
      MaskVariables	      "s=@1;r=@2;b=@3;xyz0=@4;"
      MaskIconFrame	      on
      MaskIconOpaque	      on
      MaskIconRotate	      "none"
      MaskIconUnits	      "autoscale"
      MaskValueString	      "10|28|2+2/3|[3 3 3]"
      System {
	Name			"Lorenz3"
	Location		[95, 138, 546, 569]
	Open			off
	ModelBrowserVisibility	off
	ModelBrowserWidth	200
	ScreenColor		"white"
	PaperOrientation	"landscape"
	PaperPositionMode	"auto"
	PaperType		"usletter"
	PaperUnits		"inches"
	ZoomFactor		"100"
	AutoZoom		on
	Block {
	  BlockType		  Sum
	  Name			  "-Bz+xy"
	  Ports			  [2, 1, 0, 0, 0]
	  Position		  [150, 335, 170, 355]
	  NamePlacement		  "alternate"
	  IconShape		  "round"
	  Inputs		  "|-+"
	  SaturateOnIntegerOverflow on
	  Port {
	    PortNumber		    1
	    Name		    "zdot"
	    TestPoint		    off
	    RTWStorageClass	    "Auto"
	  }
	}
	Block {
	  BlockType		  Gain
	  Name			  "Beta * z"
	  Position		  [75, 325, 130, 365]
	  Gain			  "b"
	  SaturateOnIntegerOverflow on
	}
	Block {
	  BlockType		  Integrator
	  Name			  "Integrator"
	  Ports			  [1, 1, 0, 0, 0]
	  Position		  [245, 85, 275, 105]
	  ShowName		  off
	  ExternalReset		  "none"
	  InitialConditionSource  "internal"
	  InitialCondition	  "xyz0(1)"
	  LimitOutput		  off
	  UpperSaturationLimit	  "inf"
	  LowerSaturationLimit	  "-inf"
	  ShowSaturationPort	  off
	  ShowStatePort		  off
	  AbsoluteTolerance	  "auto"
	  Port {
	    PortNumber		    1
	    Name		    "x"
	    TestPoint		    off
	    RTWStorageClass	    "Auto"
	  }
	}
	Block {
	  BlockType		  Integrator
	  Name			  "Integrator1"
	  Ports			  [1, 1, 0, 0, 0]
	  Position		  [245, 210, 275, 230]
	  ShowName		  off
	  ExternalReset		  "none"
	  InitialConditionSource  "internal"
	  InitialCondition	  "xyz0(2)"
	  LimitOutput		  off
	  UpperSaturationLimit	  "inf"
	  LowerSaturationLimit	  "-inf"
	  ShowSaturationPort	  off
	  ShowStatePort		  off
	  AbsoluteTolerance	  "auto"
	  Port {
	    PortNumber		    1
	    Name		    "y"
	    TestPoint		    off
	    RTWStorageClass	    "Auto"
	  }
	}
	Block {
	  BlockType		  Integrator
	  Name			  "Integrator2"
	  Ports			  [1, 1, 0, 0, 0]
	  Position		  [245, 335, 275, 355]
	  ShowName		  off