kalmnsce.mdl

Nonstationary Channel Estimation using a Kalman Tracking Filter 卡尔曼滤波算法的一个一个用

Model {
  Name			  "kalmnsce"
  Version		  4.00
  Description		  "Nonstationary channel estimation using a Kalman tra"
"cking filter\n"
  SampleTimeColors	  off
  LibraryLinkDisplay	  "none"
  WideLines		  on
  ShowLineDimensions	  off
  ShowPortDataTypes	  off
  RecordCoverage	  off
  CovPath		  "/"
  CovSaveName		  "covdata"
  CovMetricSettings	  "d"
  CovNameIncrementing	  off
  CovHtmlReporting	  on
  BlockNameDataTip	  off
  BlockParametersDataTip  on
  BlockDescriptionStringDataTip	off
  ToolBar		  on
  StatusBar		  on
  BrowserShowLibraryLinks on
  BrowserLookUnderMasks	  off
  Created		  "Mon Jun 15 23:17:55 1998"
  Creator		  "don"
  UpdateHistory		  "UpdateHistoryNever"
  ModifiedByFormat	  "%<Auto>"
  LastModifiedBy	  "dlinebar"
  ModifiedDateFormat	  "%<Auto>"
  LastModifiedDate	  "Tue Jan  9 13:49:34 2001"
  ModelVersionFormat	  "1.%<AutoIncrement:105>"
  ConfigurationManager	  "none"
  SimParamPage		  "Solver"
  StartTime		  "0.0"
  StopTime		  "999999"
  SolverMode		  "SingleTasking"
  Solver		  "FixedStepDiscrete"
  RelTol		  "1e-3"
  AbsTol		  "1e-6"
  Refine		  "1"
  MaxStep		  "auto"
  MinStep		  "auto"
  MaxNumMinSteps	  "-1"
  InitialStep		  "auto"
  FixedStep		  "auto"
  MaxOrder		  5
  OutputOption		  "RefineOutputTimes"
  OutputTimes		  "[]"
  LoadExternalInput	  off
  ExternalInput		  "[t, u]"
  SaveTime		  on
  TimeSaveName		  "tout"
  SaveState		  off
  StateSaveName		  "xout"
  SaveOutput		  on
  OutputSaveName	  "yout"
  LoadInitialState	  off
  InitialState		  "xInitial"
  SaveFinalState	  off
  FinalStateName	  "xFinal"
  SaveFormat		  "Array"
  LimitDataPoints	  off
  MaxDataPoints		  "1000"
  Decimation		  "1"
  AlgebraicLoopMsg	  "warning"
  MinStepSizeMsg	  "warning"
  UnconnectedInputMsg	  "warning"
  UnconnectedOutputMsg	  "warning"
  UnconnectedLineMsg	  "warning"
  InheritedTsInSrcMsg	  "warning"
  SingleTaskRateTransMsg  "none"
  MultiTaskRateTransMsg	  "error"
  IntegerOverflowMsg	  "none"
  CheckForMatrixSingularity "none"
  UnnecessaryDatatypeConvMsg "none"
  Int32ToFloatConvMsg	  "warning"
  SignalLabelMismatchMsg  "none"
  LinearizationMsg	  "none"
  VectorMatrixConversionMsg "none"
  SfunCompatibilityCheckMsg "none"
  BlockPriorityViolationMsg "warning"
  ArrayBoundsChecking	  "none"
  ConsistencyChecking	  "none"
  ZeroCross		  on
  Profile		  off
  SimulationMode	  "normal"
  RTWSystemTargetFile	  "grt.tlc"
  RTWOptions		  " -aRollThreshold=2"
  RTWInlineParameters	  off
  RTWRetainRTWFile	  off
  RTWTemplateMakefile	  "grt_default_tmf"
  RTWMakeCommand	  "make_rtw"
  RTWGenerateCodeOnly	  off
  TLCProfiler		  off
  TLCDebug		  off
  TLCCoverage		  off
  AccelSystemTargetFile	  "accel.tlc"
  AccelTemplateMakefile	  "accel_default_tmf"
  AccelMakeCommand	  "make_rtw"
  ExtModeMexFile	  "ext_comm"
  ExtModeBatchMode	  off
  ExtModeTrigType	  "manual"
  ExtModeTrigMode	  "oneshot"
  ExtModeTrigPort	  "1"
  ExtModeTrigElement	  "any"
  ExtModeTrigDuration	  1000
  ExtModeTrigHoldOff	  0
  ExtModeTrigDelay	  0
  ExtModeTrigDirection	  "rising"
  ExtModeTrigLevel	  0
  ExtModeArchiveMode	  "off"
  ExtModeAutoIncOneShot	  off
  ExtModeIncDirWhenArm	  off
  ExtModeAddSuffixToVar	  off
  ExtModeWriteAllDataToWs off
  ExtModeArmWhenConnect	  off
  ExtModeLogAll		  on
  OptimizeBlockIOStorage  on
  BufferReuse		  on
  ParameterPooling	  on
  BlockReductionOpt	  off
  RTWExpressionDepthLimit 5
  BooleanDataType	  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		    "kalmnsce"
    Location		    [14, 77, 485, 343]
    Open		    on
    ModelBrowserVisibility  off
    ModelBrowserWidth	    200
    ScreenColor		    "white"
    PaperOrientation	    "landscape"
    PaperPositionMode	    "auto"
    PaperType		    "usletter"
    PaperUnits		    "inches"
    ZoomFactor		    "100"
    AutoZoom		    on
    ReportName		    "simulink-default.rpt"
    LastButtonDownPoint	    [0, 0]
    Block {
      BlockType		      Reference
      Name		      "Channel\nNoise"
      Ports		      [0, 1]
      Position		      [275, 30, 305, 60]
      Orientation	      "down"
      SourceBlock	      "simulink3/Sources/Band-Limited\nWhite Noise"
      SourceType	      "Continuous White Noise."
      Cov		      ".3"
      Ts		      "1"
      seed		      "12345"
      VectorParams1D	      on
    }
    Block {
      BlockType		      SubSystem
      Name		      "Demo Notes"
      Ports		      []
      Position		      [55, 169, 114, 212]
      BackgroundColor	      "cyan"
      DropShadow	      on
      ShowName		      off
      ShowPortLabels	      on
      TreatAsAtomicUnit	      off
      RTWSystemCode	      "Auto"
      RTWFcnNameOpts	      "Auto"
      RTWFileNameOpts	      "Auto"
      MaskDisplay	      "disp('Info')"
      MaskIconFrame	      on
      MaskIconOpaque	      on
      MaskIconRotate	      "none"
      MaskIconUnits	      "autoscale"
      System {
	Name			"Demo Notes"
	Location		[58, 381, 478, 700]
	Open			off
	ModelBrowserVisibility	off
	ModelBrowserWidth	212
	ScreenColor		"white"
	PaperOrientation	"landscape"
	PaperPositionMode	"auto"
	PaperType		"usletter"
	PaperUnits		"inches"
	ZoomFactor		"100"
	AutoZoom		on
	LastButtonDownPoint	[0, 0]
	Annotation {
	  Position		  [212, 199]
	  Text			  "\nThis demonstrates the ability of a Kalman"
" filter to track\nthe time-varying weights of a nonstationary 5th order\nFIR "
"filter.\n\n When run, a plot is made of each tap weight over time,\nwith the "
"\"true\" filter taps drawn in yellow, and the\nestimates of those tap weights"
" in magenta.  Each of \nthe 5 tap weights is plotted in a separate axis.\n\nF"
"or more information on Kalman filtering, see\nS. Haykin, \"Adaptive Filter Th"
"eory\",\n3rd Ed., Prentice Hall, 1996.\n"
	  FontName		  "Arial"
	  FontSize		  14
	}
	Annotation {
	  Position		  [217, 49]
	  Text			  "\nNonstationary Channel Estimation\nusing a"
" Kalman Tracking Filter\n"
	  FontName		  "Arial"
	  FontSize		  14
	}
      }
    }
    Block {
      BlockType		      Reference
      Name		      "Flip"
      Ports		      [1, 1]
      Position		      [295, 191, 330, 219]
      SourceBlock	      "dspindex/Flip"
      SourceType	      "Flip"
      dim		      "Columns"
      Port {
	PortNumber		1
	Name			"Estimated\nfilter"
	TestPoint		off
	LinearAnalysisOutput	off
	LinearAnalysisInput	off
	RTWStorageClass		"Auto"
      }
    }
    Block {
      BlockType		      Reference
      Name		      "Kalman\nAdaptive Filter"
      Ports		      [2, 2]
      Position		      [185, 161, 270, 219]
      SourceBlock	      "dspadpt3/Kalman\nAdaptive Filter"
      SourceType	      "Kalman Adaptive Filter"
      Ntaps		      "5"
      Jmin		      "0.3"
      q			      "0.1"
      ic0		      "0"
      K0		      "0.5"
      Adapt		      off
    }
    Block {
      BlockType		      SubSystem
      Name		      "Nonstationary\nChannel"
      Ports		      [1, 2]
      Position		      [180, 79, 255, 126]
      ShowName		      off
      ShowPortLabels	      on
      TreatAsAtomicUnit	      off
      RTWSystemCode	      "Auto"
      RTWFcnNameOpts	      "Auto"
      RTWFileNameOpts	      "Auto"
      MaskDescription	      "Nonstationary Nth Order Noisy Channel"
      MaskHelp		      "Implementation of a nonstationary filter.  Nons"
"tationarity parameter controls the filter tap nonstationarity, and may take v"
"alues in the range 0 <= a < 1.  The first system output is the actual filter "
"taps over time, and the second output is a vector of instantaneous tap weight"
" estimates (\"unknowns\")"
      MaskPromptString	      "Filter Order:|Nonstationarity Parameter (0 to 1"
"):|Sample Time:"
      MaskStyleString	      "edit,edit,edit"
      MaskTunableValueString  "on,on,on"
      MaskCallbackString      "||"
      MaskEnableString	      "on,on,on"
      MaskVisibilityString    "on,on,on"
      MaskToolTipString	      "on,on,on"
      MaskVariables	      "N=@1;a=@2;Ts=@3;"
      MaskDisplay	      "disp('Nonstationary\\nSystem')"
      MaskIconFrame	      on
      MaskIconOpaque	      on
      MaskIconRotate	      "none"
      MaskIconUnits	      "autoscale"
      MaskValueString	      "5|0.96|1"
      MaskVariableAliases     ",,"
      Port {
	PortNumber		2
	Name			"Actual\nfilter"
	TestPoint		off
	LinearAnalysisOutput	off
	LinearAnalysisInput	off
	RTWStorageClass		"Auto"
      }
      System {
	Name			"Nonstationary\nChannel"
	Location		[101, 535, 478, 705]
	Open			off
	ModelBrowserVisibility	off
	ModelBrowserWidth	200
	ScreenColor		"white"
	PaperOrientation	"landscape"
	PaperPositionMode	"auto"
	PaperType		"usletter"
	PaperUnits		"inches"
	ZoomFactor		"100"
	AutoZoom		on
	LastButtonDownPoint	[0, 0]
	Block {
	  BlockType		  Inport
	  Name			  "Input\nSignal"
	  Position		  [45, 20, 65, 40]
	  Port			  "1"
	  Interpolate		  on
	}
	Block {
	  BlockType		  Reference
	  Name			  "Average\nCoefficients\nTranspose Filter"
	  Ports			  [1, 1]
	  Position		  [95, 83, 165, 127]
	  SourceBlock		  "dsparch3/Direct-Form II\nTranspose Filter"
	  SourceType		  "Direct-Form II Transpose Filter"
	  num			  "1"
	  den			  "[1 -a]"
	  ic			  "0"
	}
	Block {
	  BlockType		  Reference
	  Name			  "Random \nTaps"
	  Ports			  [0, 1]
	  Position		  [40, 90, 70, 120]
	  SourceBlock		  "simulink3/Sources/Band-Limited\nWhite Noise"
	  SourceType		  "Continuous White Noise."
	  Cov			  "ones(N,1)"
	  Ts			  "Ts"
	  seed			  "rand(N,1)*1000"
	  VectorParams1D	  "on"
	}
	Block {
	  BlockType		  Reference
	  Name			  "Time-Varying FIR Filter"
	  Ports			  [2, 1]
	  Position		  [210, 18, 300, 62]
	  SourceBlock		  "dsparch3/Time-Varying Direct-\nForm II Tran"
"spose Filter"
	  SourceType		  "Time-Varying Direct-Form II Transpose Filte"
"r"
	  NumDen		  "All-Zero (FIR)"
	  ic			  "0"
	  filtCheck		  "on"
	  FiltPerFrame		  "One Filter Per Sample Time"
	}
	Block {
	  BlockType		  Outport
	  Name			  "Filter\nOutput"
	  Position		  [330, 30, 350, 50]
	  Port			  "1"
	  OutputWhenDisabled	  "held"
	  InitialOutput		  "0"
	}
	Block {
	  BlockType		  Outport
	  Name			  "Tap\nWeights"
	  Position		  [330, 95, 350, 115]
	  Port			  "2"
	  OutputWhenDisabled	  "held"
	  InitialOutput		  "0"
	}
	Line {