Shrikant: I dont recall how we adjusted the model for different
amplitude step stimuli. I know that we adjusted it for age. But if
I use a 2D step I get a good response but if I use a 1D step I get a
double humped response for the 25 year old. Then if I decrease the
phasic vel gain from 4 to 3 and increase the tonic position signal
from 4 to 6 the response looks OK, but I dont recall doing any of this before.
I looked in several of our model papers but could not find what I
was looking for. I have to get matlab licence updated for our web
model to work again.
Any ideas on how to change the model for different amplitudes? Ive
attached the model for a 25 year old.
Thanks,
Cliff
--
EYE-MOVEMENT mailing list ([log in to unmask])
N.B. Replies are sent to the list, not the sender
To unsubscribe, etc. see http://www.jiscmail.ac.uk/files/eye-movement/introduction.html
Other queries to list owner at [log in to unmask]
Model {
Name "simplified25"
Version 5.0
SaveDefaultBlockParams on
SampleTimeColors off
LibraryLinkDisplay "none"
WideLines off
ShowLineDimensions off
ShowPortDataTypes off
ShowLoopsOnError on
IgnoreBidirectionalLines off
ShowStorageClass off
ExecutionOrder off
RecordCoverage off
CovPath "/"
CovSaveName "covdata"
CovMetricSettings "dw"
CovNameIncrementing off
CovHtmlReporting on
covSaveCumulativeToWorkspaceVar on
CovSaveSingleToWorkspaceVar on
CovCumulativeVarName "covCumulativeData"
CovCumulativeReport off
DataTypeOverride "UseLocalSettings"
MinMaxOverflowLogging "UseLocalSettings"
MinMaxOverflowArchiveMode "Overwrite"
BlockNameDataTip off
BlockParametersDataTip off
BlockDescriptionStringDataTip off
ToolBar on
StatusBar on
BrowserShowLibraryLinks off
BrowserLookUnderMasks off
Created "Wed Oct 13 10:36:02 2004"
UpdateHistory "UpdateHistoryNever"
ModifiedByFormat "%<Auto>"
LastModifiedBy "Schor"
ModifiedDateFormat "%<Auto>"
LastModifiedDate "Fri Oct 15 14:07:10 2004"
ModelVersionFormat "1.%<AutoIncrement:41>"
ConfigurationManager "None"
SimParamPage "Solver"
LinearizationMsg "none"
Profile off
ParamWorkspaceSource "MATLABWorkspace"
AccelSystemTargetFile "accel.tlc"
AccelTemplateMakefile "accel_default_tmf"
AccelMakeCommand "make_rtw"
TryForcingSFcnDF off
ExtModeMexFile "ext_comm"
ExtModeBatchMode off
ExtModeTrigType "manual"
ExtModeTrigMode "normal"
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 on
ExtModeSkipDownloadWhenConnect off
ExtModeLogAll on
ExtModeAutoUpdateStatusClock on
BufferReuse on
RTWExpressionDepthLimit 5
SimulationMode "normal"
Solver "ode45"
SolverMode "Auto"
StartTime "0.0"
StopTime "4"
MaxOrder 5
MaxStep "auto"
MinStep "auto"
MaxNumMinSteps "-1"
InitialStep "auto"
FixedStep "auto"
RelTol "1e-3"
AbsTol "auto"
OutputOption "RefineOutputTimes"
OutputTimes "[]"
Refine "1"
LoadExternalInput off
ExternalInput "[t, u]"
LoadInitialState off
InitialState "xInitial"
SaveTime on
TimeSaveName "tout"
SaveState off
StateSaveName "xout"
SaveOutput on
OutputSaveName "yout"
SaveFinalState off
FinalStateName "xFinal"
SaveFormat "Array"
Decimation "1"
LimitDataPoints on
MaxDataPoints "1000"
SignalLoggingName "sigsOut"
ConsistencyChecking "none"
ArrayBoundsChecking "none"
AlgebraicLoopMsg "warning"
BlockPriorityViolationMsg "warning"
MinStepSizeMsg "warning"
InheritedTsInSrcMsg "warning"
DiscreteInheritContinuousMsg "warning"
MultiTaskRateTransMsg "error"
SingleTaskRateTransMsg "none"
CheckForMatrixSingularity "none"
IntegerOverflowMsg "warning"
Int32ToFloatConvMsg "warning"
ParameterDowncastMsg "error"
ParameterOverflowMsg "error"
ParameterPrecisionLossMsg "warning"
UnderSpecifiedDataTypeMsg "none"
UnnecessaryDatatypeConvMsg "none"
VectorMatrixConversionMsg "none"
InvalidFcnCallConnMsg "error"
SignalLabelMismatchMsg "none"
UnconnectedInputMsg "warning"
UnconnectedOutputMsg "warning"
UnconnectedLineMsg "warning"
SfunCompatibilityCheckMsg "none"
RTWInlineParameters off
BlockReductionOpt on
BooleanDataType on
ConditionallyExecuteInputs on
ParameterPooling on
OptimizeBlockIOStorage on
ZeroCross on
AssertionControl "UseLocalSettings"
ProdHWDeviceType "Microprocessor"
ProdHWWordLengths "8,16,32,32"
RTWSystemTargetFile "grt.tlc"
RTWTemplateMakefile "grt_default_tmf"
RTWMakeCommand "make_rtw"
RTWGenerateCodeOnly off
RTWRetainRTWFile off
TLCProfiler off
TLCDebug off
TLCCoverage off
TLCAssertion off
BlockDefaults {
Orientation "right"
ForegroundColor "black"
BackgroundColor "white"
DropShadow off
NamePlacement "normal"
FontName "Helvetica"
FontSize 10
FontWeight "normal"
FontAngle "normal"
ShowName on
}
BlockParameterDefaults {
Block {
BlockType Constant
Value "1"
VectorParams1D on
ShowAdditionalParam off
OutDataTypeMode "Inherit from 'Constant value'"
OutDataType "sfix(16)"
ConRadixGroup "Use specified scaling"
OutScaling "2^0"
}
Block {
BlockType DeadZone
LowerValue "-0.5"
UpperValue "0.5"
SaturateOnIntegerOverflow on
LinearizeAsGain on
ZeroCross on
}
Block {
BlockType Derivative
}
Block {
BlockType DigitalClock
SampleTime "1"
}
Block {
BlockType Display
Format "short"
Decimation "10"
Floating off
SampleTime "-1"
}
Block {
BlockType Gain
Gain "1"
Multiplication "Element-wise(K.*u)"
ShowAdditionalParam off
ParameterDataTypeMode "Same as input"
ParameterDataType "sfix(16)"
ParameterScalingMode "Best Precision: Matrix-wise"
ParameterScaling "2^0"
OutDataTypeMode "Same as input"
OutDataType "sfix(16)"
OutScaling "2^0"
LockScale off
RndMeth "Floor"
SaturateOnIntegerOverflow on
}
Block {
BlockType Relay
OnSwitchValue "eps"
OffSwitchValue "eps"
OnOutputValue "1"
OffOutputValue "0"
ShowAdditionalParam off
OutputDataTypeScalingMode "All ports same datatype"
OutDataType "sfix(16)"
OutScaling "2^0"
ConRadixGroup "Use specified scaling"
ZeroCross on
}
Block {
BlockType Saturate
UpperLimit "0.5"
LowerLimit "-0.5"
LinearizeAsGain on
ZeroCross on
}
Block {
BlockType Scope
Floating off
ModelBased off
TickLabels "OneTimeTick"
ZoomMode "on"
Grid "on"
TimeRange "auto"
YMin "-5"
YMax "5"
SaveToWorkspace off
SaveName "ScopeData"
LimitDataPoints on
MaxDataPoints "5000"
Decimation "1"
SampleInput off
SampleTime "0"
}
Block {
BlockType Step
Time "1"
Before "0"
After "1"
SampleTime "-1"
VectorParams1D on
ZeroCross on
}
Block {
BlockType Sum
IconShape "rectangular"
Inputs "++"
ShowAdditionalParam off
InputSameDT on
OutDataTypeMode "Same as first input"
OutDataType "sfix(16)"
OutScaling "2^0"
LockScale off
RndMeth "Floor"
SaturateOnIntegerOverflow on
}
Block {
BlockType Switch
Criteria "u2 >= Threshold"
Threshold "0"
ShowAdditionalParam off
InputSameDT on
OutDataTypeMode "Inherit via internal rule"
RndMeth "Floor"
SaturateOnIntegerOverflow on
ZeroCross on
}
Block {
BlockType TransferFcn
Numerator "[1]"
Denominator "[1 2 1]"
AbsoluteTolerance "auto"
Realization "auto"
}
Block {
BlockType TransportDelay
DelayTime "1"
InitialInput "0"
BufferSize "1024"
PadeOrder "0"
TransDelayFeedthrough off
}
}
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 "simplified25"
Location [2, 70, 1022, 720]
Open on
ModelBrowserVisibility off
ModelBrowserWidth 200
ScreenColor "white"
PaperOrientation "landscape"
PaperPositionMode "auto"
PaperType "usletter"
PaperUnits "inches"
ZoomFactor "100"
ReportName "simulink-default.rpt"
Block {
BlockType Scope
Name "Accommodation\nresponse scope"
Ports [1]
Position [880, 328, 910, 362]
FontName "Arial"
Location [5, 72, 1029, 762]
Open off
NumInputPorts "1"
TickLabels "on"
ZoomMode "xonly"
List {
ListType AxesTitles
axes1 "%<SignalLabel>"
}
List {
ListType SelectedSignals
axes1 ""
}
TimeRange "4"
YMin "0"
YMax "3"
DataFormat "StructureWithTime"
}
Block {
BlockType TransportDelay
Name "F-B latency"
Position [205, 415, 235, 445]
Orientation "up"
NamePlacement "alternate"
FontName "Arial"
DelayTime ".300"
}
Block {
BlockType Sum
Name "Internal F-B"
Ports [2, 1]
Position [250, 375, 270, 395]
ShowName off
FontName "Arial"
IconShape "round"
Inputs "|+-"
}
Block {
BlockType Display
Name "Response meter"
Ports [1]
Position [875, 183, 945, 217]
FontName "Arial"
Decimation "1"
}
Block {
BlockType Sum
Name "Sum"
Ports [2, 1]
Position [70, 375, 90, 395]
ShowName off
FontSize 12
IconShape "round"
Inputs "|+-"
}
Block {
BlockType Sum
Name "Sum1"
Ports [2, 1]
Position [705, 265, 725, 285]
ShowName off
FontSize 12
IconShape "round"
Inputs "|++"
}
Block {
BlockType Sum
Name "Sum2"
Ports [2, 1]
Position [525, 300, 545, 320]
ShowName off
FontSize 12
IconShape "round"
Inputs "|++"
}
Block {
BlockType Sum
Name "Sum3"
Ports [2, 1]
Position [210, 375, 230, 395]
ShowName off
FontSize 12
IconShape "round"
Inputs "|++"
}
Block {
BlockType Sum
Name "Sum5"
Ports [2, 1]
Position [600, 265, 620, 285]
NamePlacement "alternate"
ShowName off
FontSize 12
IconShape "round"
Inputs "|++"
}
Block {
BlockType Derivative
Name "accel"
Position [960, 260, 990, 290]
}
Block {
BlockType Scope
Name "accel\nscope"
Ports [1]
Position [955, 184, 985, 216]
Orientation "left"
FontName "Arial"
Location [151, 166, 923, 762]
Open off
NumInputPorts "1"
ZoomMode "xonly"
List {
ListType AxesTitles
axes1 "%<SignalLabel>"
}
List {
ListType SelectedSignals
axes1 ""
}
TimeRange "0.04096000000000002"
YMin "-100"
YMax "90"
SaveName "ScopeData1"
DataFormat "StructureWithTime"
}
Block {
BlockType TransferFcn
Name "acceleration pulse \nintegrator"
Position [400, 217, 460, 253]
NamePlacement "alternate"
FontName "Arial"
Numerator "[3.75]"
Denominator "[1 5]"
}
Block {
BlockType Gain
Name "choroid compliance"
Position [650, 260, 680, 290]
NamePlacement "alternate"
Gain ".25"
}
Block {
BlockType Switch
Name "closed-loop\nstep switch"
Position [315, 289, 360, 331]
Threshold ".435"
}
Block {
BlockType DeadZone
Name "dead zone"
Position [110, 370, 140, 400]
LowerValue "-0.14"
UpperValue "0.14"
SaturateOnIntegerOverflow off
}
Block {
BlockType Derivative
Name "derivative"
Position [65, 150, 95, 180]
}
Block {
BlockType DigitalClock
Name "digital Clock"
Position [210, 258, 275, 282]
NamePlacement "alternate"
SampleTime ".005"
}
Block {
BlockType TransferFcn
Name "fast step\n integrator"
Position [380, 292, 440, 328]
FontName "Arial"
Numerator "[7.5]"
Denominator "[1 .5]"
}
Block {
BlockType TransportDelay
Name "latency"
Position [160, 370, 190, 400]
FontName "Arial"
DelayTime ".300"
}
Block {
BlockType TransferFcn
Name "lens\ntime constant"
Position [745, 257, 805, 293]
NamePlacement "alternate"
FontName "Arial"
Numerator "[1.75]"
Denominator "[1 5.2]"
}
Block {
BlockType Gain
Name "negative\nfeedback gain"
Position [140, 448, 180, 482]
Orientation "left"
FontName "Arial"
}
Block {
BlockType Gain
Name "phasic-velocity\nsignal gain"
Position [500, 219, 535, 251]
NamePlacement "alternate"
FontName "Arial"
Gain "4"
}
Block {
BlockType Switch
Name "plateau width\nswitch"
Position [315, 214, 360, 256]
NamePlacement "alternate"
Threshold ".435"
}
Block {
BlockType DeadZone
Name "position \nthreshold"
Position [20, 150, 50, 180]
NamePlacement "alternate"
LowerValue "-.75"
UpperValue ".75"
SaturateOnIntegerOverflow off
}
Block {
BlockType Gain
Name "posterior zonule\ncompliance"
Position [650, 295, 680, 325]
Gain ".4"
}
Block {
BlockType Saturate
Name "pulse plateau\n height"
Position [240, 150, 270, 180]
NamePlacement "alternate"
UpperLimit "7.5"
LowerLimit "0"
}
Block {
BlockType TransferFcn
Name "pulse ramp \nslope"
Position [165, 147, 225, 183]
FontName "Arial"
Numerator "[15]"
Denominator "[1 .05]"
}
Block {
BlockType TransferFcn
Name "slow step\nintegrator"
Position [465, 337, 525, 373]
FontName "Arial"
Numerator "[.5]"
Denominator "[1 .01]"
}
Block {
BlockType Step
Name "step Input"
Position [20, 370, 50, 400]
Time "0"
After "2"
SampleTime "0.005"
}
Block {
BlockType Constant
Name "switch\noff-level"
Position [285, 365, 315, 395]
Orientation "up"
Value "0"
}
Block {
BlockType Constant
Name "switch \noff level"
Position [285, 150, 315, 180]
Orientation "down"
Value "0"
}
Block {
BlockType Gain
Name "tonic-position\nsignal gain"
Position [565, 295, 605, 325]
Gain "4"
}
Block {
BlockType Derivative
Name "velocity"
Position [885, 260, 915, 290]
}
Block {
BlockType Scope
Name "velocity \nscope"
Ports [1]
Position [960, 328, 990, 362]
FontName "Arial"
Location [5, 59, 1029, 749]
Open off
NumInputPorts "1"
ZoomMode "yonly"
List {
ListType AxesTitles
axes1 "%<SignalLabel>"
}
List {
ListType SelectedSignals
axes1 ""
}
TimeRange "0.1024000000000001"
YMin "5.17269"
YMax "5.17342"
SaveName "1Dvelocity"
DataFormat "StructureWithTime"
LimitDataPoints off
}
Block {
BlockType Relay
Name "velocity thresh"
Position [115, 150, 145, 180]
NamePlacement "alternate"
OnSwitchValue "600"
OffSwitchValue "-.1"
OnOutputValue "4.86"
}
Line {
SrcBlock "negative\nfeedback gain"
SrcPort 1
Points [-55, 0]
DstBlock "Sum"
DstPort 2
}
Line {
SrcBlock "accel"
SrcPort 1
Points [5, 0]
DstBlock "accel\nscope"
DstPort 1
}
Line {
SrcBlock "velocity"
SrcPort 1
Points [10, 0; 10, 0]
Branch {
DstBlock "accel"
DstPort 1
}
Branch {
Points [0, 70]
DstBlock "velocity \nscope"
DstPort 1
}
}
Line {
SrcBlock "Sum"
SrcPort 1
DstBlock "dead zone"
DstPort 1
}
Line {
SrcBlock "dead zone"
SrcPort 1
DstBlock "latency"
DstPort 1
}
Line {
SrcBlock "latency"
SrcPort 1
Points [0, 0; 0, -35]
Branch {
Points [-190, 0]
DstBlock "position \nthreshold"
DstPort 1
}
Branch {
DstBlock "Sum3"
DstPort 1
}
}
Line {
SrcBlock "step Input"
SrcPort 1
DstBlock "Sum"
DstPort 1
}
Line {
SrcBlock "position \nthreshold"
SrcPort 1
DstBlock "derivative"
DstPort 1
}
Line {
SrcBlock "derivative"
SrcPort 1
DstBlock "velocity thresh"
DstPort 1
}
Line {
SrcBlock "velocity thresh"
SrcPort 1
DstBlock "pulse ramp \nslope"
DstPort 1
}
Line {
SrcBlock "Sum1"
SrcPort 1
DstBlock "lens\ntime constant"
DstPort 1
}
Line {
SrcBlock "pulse ramp \nslope"
SrcPort 1
DstBlock "pulse plateau\n height"
DstPort 1
}
Line {
Labels [2, 0]
SrcBlock "fast step\n integrator"
SrcPort 1
Points [0, 0; 5, 0]
Branch {
DstBlock "Sum2"
DstPort 1
}
Branch {
DstBlock "slow step\nintegrator"
DstPort 1
}
}
Line {
SrcBlock "slow step\nintegrator"
SrcPort 1
Points [5, 0]
DstBlock "Sum2"
DstPort 2
}
Line {
SrcBlock "digital Clock"
SrcPort 1
Points [15, 0]
Branch {
Points [0, -35]
DstBlock "plateau width\nswitch"
DstPort 2
}
Branch {
Points [0, 40]
DstBlock "closed-loop\nstep switch"
DstPort 2
}
}
Line {
SrcBlock "switch \noff level"
SrcPort 1
DstBlock "plateau width\nswitch"
DstPort 1
}
Line {
SrcBlock "switch\noff-level"
SrcPort 1
DstBlock "closed-loop\nstep switch"
DstPort 3
}
Line {
SrcBlock "plateau width\nswitch"
SrcPort 1
DstBlock "acceleration pulse \nintegrator"
DstPort 1
}
Line {
SrcBlock "Sum2"
SrcPort 1
DstBlock "tonic-position\nsignal gain"
DstPort 1
}
Line {
Labels [1, 0]
SrcBlock "tonic-position\nsignal gain"
SrcPort 1
DstBlock "Sum5"
DstPort 2
}
Line {
SrcBlock "choroid compliance"
SrcPort 1
DstBlock "Sum1"
DstPort 1
}
Line {
SrcBlock "posterior zonule\ncompliance"
SrcPort 1
Points [30, 0]
DstBlock "Sum1"
DstPort 2
}
Line {
SrcBlock "Sum5"
SrcPort 1
Points [0, 0]
Branch {
DstBlock "choroid compliance"
DstPort 1
}
Branch {
Points [5, 0; 0, 35]
DstBlock "posterior zonule\ncompliance"
DstPort 1
}
}
Line {
SrcBlock "closed-loop\nstep switch"
SrcPort 1
DstBlock "fast step\n integrator"
DstPort 1
}
Line {
SrcBlock "F-B latency"
SrcPort 1
DstBlock "Sum3"
DstPort 2
}
Line {
SrcBlock "Sum3"
SrcPort 1
DstBlock "Internal F-B"
DstPort 1
}
Line {
SrcBlock "Internal F-B"
SrcPort 1
Points [0, -90]
DstBlock "closed-loop\nstep switch"
DstPort 1
}
Line {
SrcBlock "pulse plateau\n height"
SrcPort 1
Points [0, 85]
DstBlock "plateau width\nswitch"
DstPort 3
}
Line {
SrcBlock "acceleration pulse \nintegrator"
SrcPort 1
DstBlock "phasic-velocity\nsignal gain"
DstPort 1
}
Line {
SrcBlock "phasic-velocity\nsignal gain"
SrcPort 1
Points [15, 0; 0, 40]
DstBlock "Sum5"
DstPort 1
}
Line {
SrcBlock "lens\ntime constant"
SrcPort 1
Points [10, 0]
Branch {
Points [25, 0]
Branch {
DstBlock "velocity"
DstPort 1
}
Branch {
Points [0, 70]
DstBlock "Accommodation\nresponse scope"
DstPort 1
}
Branch {
Points [0, -75]
DstBlock "Response meter"
DstPort 1
}
}
Branch {
Points [0, 0; 0, 190; -560, 0]
Branch {
Points [0, -65]
DstBlock "Internal F-B"
DstPort 2
}
Branch {
Points [-40, 0]
Branch {
DstBlock "negative\nfeedback gain"
DstPort 1
}
Branch {
DstBlock "F-B latency"
DstPort 1
}
}
}
}
}
}
--
EYE-MOVEMENT mailing list ([log in to unmask])
N.B. Replies are sent to the list, not the sender
To unsubscribe, etc. see http://www.jiscmail.ac.uk/files/eye-movement/introduction.html
Other queries to list owner at [log in to unmask]
|