sim12xx1.mdl

solve power flows by gauss seidel method ans newton Rapson

Model {
  Name			  "sim12ex11"
  Version		  2.09
  SimParamPage		  Solver
  SampleTimeColors	  off
  InvariantConstants	  off
  WideVectorLines	  off
  ShowLineWidths	  off
  PaperOrientation	  landscape
  PaperType		  usletter
  PaperUnits		  inches
  StartTime		  "0.0"
  StopTime		  "2"
  Solver		  ode45
  RelTol		  "1e-3"
  AbsTol		  "1e-3"
  Refine		  "1"
  MaxStep		  ".010"
  InitialStep		  "auto"
  FixedStep		  "auto"
  MaxOrder		  5
  OutputOption		  RefineOutputTimes
  OutputTimes		  "[]"
  LoadExternalInput	  off
  ExternalInput		  "[t, u]"
  SaveTime		  off
  TimeSaveName		  "tout"
  SaveState		  off
  StateSaveName		  "xout"
  SaveOutput		  off
  OutputSaveName	  "yout"
  LoadInitialState	  off
  InitialState		  "xInitial"
  SaveFinalState	  off
  FinalStateName	  "xFinal"
  LimitMaxRows		  off
  MaxRows		  "1000"
  Decimation		  "1"
  AlgebraicLoopMsg	  warning
  MinStepSizeMsg	  warning
  UnconnectedInputMsg	  warning
  UnconnectedOutputMsg	  warning
  UnconnectedLineMsg	  warning
  ConsistencyChecking	  off
  ZeroCross		  on
  SimulationMode	  normal
  RTWSystemTargetFile	  "grt.tlc"
  RTWInlineParameters	  off
  RTWRetainRTWFile	  off
  RTWTemplateMakefile	  "grt_vc.tmf"
  RTWMakeCommand	  "make_rtw"
  RTWGenerateCodeOnly	  off
  ExtModeMexFile	  "ext_comm"
  ExtModeBatchMode	  off
  BlockDefaults {
    Orientation		    right
    ForegroundColor	    black
    BackgroundColor	    white
    DropShadow		    off
    NamePlacement	    normal
    FontName		    "Helvetica"
    FontSize		    10
    FontWeight		    normal
    FontAngle		    normal
    ShowName		    on
  }
  AnnotationDefaults {
    HorizontalAlignment	    center
    VerticalAlignment	    middle
    ForegroundColor	    black
    BackgroundColor	    white
    DropShadow		    off
    FontName		    "Helvetica"
    FontSize		    10
    FontWeight		    normal
    FontAngle		    normal
  }
  LineDefaults {
    FontName		    "Helvetica"
    FontSize		    9
    FontWeight		    normal
    FontAngle		    normal
  }
  System {
    Name		    "sim12ex11"
    Location		    [147, 280, 837, 694]
    Open		    on
    ScreenColor		    white
    Block {
      BlockType		      Demux
      Name		      "Demux"
      Ports		      [1, 3, 0, 0, 0]
      Position		      [425, 96, 480, 164]
      Outputs		      "3"
    }
    Block {
      BlockType		      StateSpace
      Name		      "Feedback Gain \nusing LQR Design."
      Position		      [295, 225, 365, 285]
      Orientation	      left
      A			      "[]"
      B			      "[]"
      C			      "[]"
      D			      "k"
      X0		      "0"
      MaskType		      "Regulator"
      MaskDescription	      "Linear quadratic regulator.\n(Requires Control "
			      "Toolbox)"
      MaskHelp		      "Implements a linear state quadratic regulator. "
			      "Mask and unmask this block to see how it works."
      MaskPromptString	      "Plant A matrix:|Plant B matrix:|Plant C matrix:|"
			      "Plant D matrix:|Output weighting, y'Qy:|Input "
			      "weighting, u'Ru:"
      MaskStyleString	      "edit,edit,edit,edit,edit,edit"
      MaskInitialization      "[k]=lqry(@1,@2,@3,@4,@5,@6);"
      MaskDisplay	      "disp('K')"
      MaskIconFrame	      on
      MaskIconOpaque	      on
      MaskIconRotate	      none
      MaskIconUnits	      autoscale
      MaskValueString	      "[-5   0   -100;  2   -2    0;   0    0.1   "
			      "-0.08]|[0;   0;    -0.1]|eye(3)|zeros(3,1)|[20, "
			      "0  0;  0  10   0;  0   0   5]      |eye(1,1)*.15"
    }
    Block {
      BlockType		      Scope
      Name		      "Scope"
      Ports		      [1, 0, 0, 0, 0]
      Position		      [540, 126, 580, 174]
      Floating		      off
      Location		      [188, 365, 512, 604]
      Open		      off
      Grid		      on
      TickLabels	      on
      ZoomMode		      on
      TimeRange		      "auto"
      YMin		      "-5"
      YMax		      "5"
      SaveToWorkspace	      off
      SaveName		      "ScopeData"
      LimitMaxRows	      on
      MaxRows		      "5000"
      Decimation	      "1"
      SampleInput	      off
      SampleTime	      "0"
    }
    Block {
      BlockType		      StateSpace
      Name		      "State-Space"
      Position		      [260, 99, 360, 161]
      A			      "[-5  0  -100;  2  -2  0;  0  0.1  -0.08]"
      B			      "[0;   0 ;  -0.1]"
      C			      "eye(3)"
      D			      "zeros(3, 1)"
      X0		      "0"
    }
    Block {
      BlockType		      Step
      Name		      "Step Input"
      Position		      [55, 95, 95, 135]
      Time		      "10"
      Before		      "0.2"
      After		      ".2"
    }
    Block {
      BlockType		      Sum
      Name		      "Sum"
      Ports		      [2, 1, 0, 0, 0]
      Position		      [185, 104, 210, 151]
      Inputs		      "+-"
    }
    Line {
      SrcBlock		      "Sum"
      SrcPort		      1
      DstBlock		      "State-Space"
      DstPort		      1
    }
    Line {
      SrcBlock		      "Step Input"
      SrcPort		      1
      DstBlock		      "Sum"
      DstPort		      1
    }
    Line {
      SrcBlock		      "Demux"
      SrcPort		      3
      DstBlock		      "Scope"
      DstPort		      1
    }
    Line {
      SrcBlock		      "State-Space"
      SrcPort		      1
      Points		      [30, 0]
      Branch {
	DstBlock		"Demux"
	DstPort			1
      }
      Branch {
	Points			[0, 125]
	DstBlock		"Feedback Gain \nusing LQR Design."
	DstPort			1
      }
    }
    Line {
      SrcBlock		      "Feedback Gain \nusing LQR Design."
      SrcPort		      1
      Points		      [-135, 0; 0, -115]
      DstBlock		      "Sum"
      DstPort		      2
    }
    Annotation {
      Position		      [348, 335]
      FontSize		      12
      Text		      "Simulation block diagram for Example 12.11"
    }
    Annotation {
      Position		      [512, 132]
      VerticalAlignment	      top
      Text		      "Dw"
    }
    Annotation {
      Position		      [322, 47]
      VerticalAlignment	      top
      FontSize		      12
      Text		      "LQG Demo. Start simulation.\nTry varying "
			      "weighting parameters."
    }
  }
}