股票指数动态的弹簧块类比

对于这两个收益值，概率分布函数的特征是一个单一的最大值，尽管最优投资期限发生了变化，这在文献中被称为损益不对称[37]。当传送带以u速度启动时，整个系统与传送带一起移动，直到第一块开始滑动。当作用在块体上的弹性力之和超过静摩擦力（Fst）的最大值时，就会产生滑动力矩。滑动运动由动力摩擦力调节。当相对速度vr=vi时，挡块再次粘附在皮带上- u变成零。之后，该过程可能会再次开始（有关弹簧块系统动力学的amovie可参考[48]）。在这里进行的类比中，单个股票的价格可以建模为区块的xi（t）位置，从外部静态点测量。例如，如果我们想对代表30只大型美国股票经济状态的道琼斯工业平均指数（DJIA index）进行类比，那么模型中将考虑由N=30个区块组成的链。由于人们可以根据股票在指数中的价格在股票之间建立一个排名，因此链中各区块之间预先存在的顺序并非完全无关。当然，股票的价格排名会随着时间的推移而变化。相应的事件是，块体相互穿过，但在这个简单的模型中没有考虑这种非物理现象。在我们的模型框架中，X（t）指数的值被定义为系统的质心，在静态坐标系中测量：X（t）=PNi=1xi（t）N.（5）通过弹簧作用于块之间的弹性力对应于股票相关性。传送带通过摩擦力作用于系统。

皮带施加的拉力模拟了增长的愿望，这是所有公司和社会都期待的普遍趋势。对于较小的皮带速度，粘滑动力学将导致不同尺寸的“雪崩”。因此，链的长度（由最后一个区块的位置xN确定）随时间变化，振幅变化很大，如图8所示。a、 参考文献[25]。当所有股票价格以相关方式下跌时，这种自组织行为会导致指数价值急剧下降[42]。为了解释指数持续上升时的随机跳跃，以及观察到的正对数回报的肥尾分布，皮带的u速度是随机变化的（更多细节见附录A）。使所提出的类比起作用的一个重要因素是在X（t）的动力学中选择采样时间。金融时间序列（在我们的例子中是道琼斯工业平均指数）代表了大约80年的每日收盘价格。与文献中的许多其他研究类似，本研究未考虑该指数的日内变异性。因此，在弹簧块类比中，除了动力学的显著连续时间外，我们还需要定义一个离散时间步。为此，模拟的时间序列需要以明确的频率进行采样，以这种方式产生一个离散的时间序列，该时间序列应对应于一个交易日间隔。该采样频率将影响Spring block模型生成的人工指数中测量的波动率。

由于道琼斯工业平均指数的对数收益波动率已知，我们可以通过施加该实验值来确定模型中的抽样区间（见附录A）。该模型中引入的类比的力量在于，X（t）动力学与观察到的股指动力学有许多相似之处。附录A.4中给出了弹簧块链的动力学、计算机模拟和使用参数的详细信息。模型结果与历史股市数据对比考虑到附录A中规定的模拟参数和附录B中描述的积分方法，我们现在比较SB训练模型和DJIA指数测量的统计结果。选择参数值是为了以最佳方式再现DJIA测量的所有程式化事实。当然，人们可以选择其他参数，这些参数将为某些统计指标带来更好的结果，但在这种情况下，其他指标将获得更大的差异。然而，我们记得，我们在这里声明的目的并不是对股票指数的动态给出区域性的描述，因此，用这样一个简单的力学方法无法重现股票指数动态的所有观察到的统计特征也就不足为奇了。9.29.39.4DJIA17700 17800 17900 18000 18100 18200t[天]8.28.38.4弹簧块图4：DJIA对数指数动态和弹簧块模型得出的对数指数（ln[X（t）]的相同大小（500天）的样本。在图4中，我们展示了DJIA的对数指数（s（t））与SB模型生成的对数指数（ln[X（t）]）之间的首次视觉比较。尽管人们可以观察到视觉上可察觉的统计差异，但结果还是很有希望的。

如果我们分析这些股票的对数收益率，也可以得出类似的结果t=1天（图5）。如果我们绘制对数指数的分布函数（图6），则通过目视检查对数回报所检测到的差异和相似性更好地说明。从对数图可以得出结论，在传送带的随机驱动下，SB模型成功地再现了高返回值下观察到的幂律尾。然而，在低回报限额方面存在显著差异。与我们的SB方法相比，DJIA的分布函数似乎更窄，不对称性更小。为了进一步进行统计比较，可以构造DJIA和SBIA模型的雪崩大小分布。雪崩被定义为时间矩与t之间指数的连续单调下降。雪崩指数变化的分布( = S（t）-S（t）代表DJIA和 = X（t）-X（t）对于B模型）生成雪崩大小分布。测量结果如图7所示。正如预期的那样-0.050.05DJIA11400 11600 11800 12000t[天]-0.050.05spring-Block图5：对数回报率之间800天时间段的视觉比较（对于在DJIA数据和spring区块模型中，t=1天）。±0.001 0.01 0.1rt=1>0弹簧块-0.01-0.1rt=1<00.00010.0010.010.1图6：正对数收益和负对数收益的分布(t=1天），用于DJIA指数和SBA模型生成的艺术指数。在这两个系统中都有许多小的雪崩，其中一个会发现一些非常大的雪崩。这种情况是自组织临界性的特征，这似乎是两个系统都具有的特征。

这两个分布函数非常相似，尽管SB模型的趋势比DJIA的情况更复杂。最后，我们考虑反向统计，即固定对数回报值的投资期限分布。我们考虑了相同的ρ=+5σ和ρ=-DJIA指数（图3）的5σ对数返回值和SB模型的结果绘制在图8中。SB模型成功地再现了真实股指数据中观察到的损益不对称性。虽然这些曲线的最大值位置移动了几天，但曲线的形状和投资范围内的不对称性是相似的。因此，SB模型也能够揭示股票指数的这种微妙的统计特征。5.结论考虑了一种简单的力学类比来模拟股票指数的动力学特征。一组由弹簧连接的块体被用作模型10 100 1000 10000，这些块体受到运行传送带的连续驱动0.00010.0010.010.1spring-Blockdjia图7:DJIA指数动态和弹簧块模型中的雪崩大小分布。1100 1000τ0.0050.010.0150.020.0250.03ρ=+5σρ=-5σ图8：考虑|ρ|=5σ，弹簧块模型生成的艺术指数的投资期限分布≈ 0.05（即5%的回报率）。将结果与DJIA指数的结果进行比较（图3）。股票指数系统。与最近报道的大量模型相反，我们的目标是截然不同的。我们没有考虑与经济相关的复杂描述，而是测试了一个普遍而简单的物理模型：弹簧块系统的适用性。

我们所做的类比是有充分动机的，强调了个别股票价格与其集体行为之间的联系。我们已经确定了模型的驱动（传送带的速度）和自由参数，以便DJIAstock指数已知的大多数程式化事实都能合理地重现。在判断我们的方法的成功或失败时，我们还必须考虑这样一个事实，即我们并没有像最近文献中提出的大多数模型那样，只关注少数几个和选定的程式化事实。相反，我们考虑了道琼斯工业平均指数（DJIA stockindex）的所有已使用和可测量的统计特征，试图用固定的参数集将它们作为一个整体重现。引入的模型的一个重要结果是，它解释了在逆静力学中观察到的令人困惑的增益-损耗不对称现象，并且该模型也很好地描述了雪崩尺寸分布曲线。尽管在定性地再现指数动态的已知统计特征方面取得了成功，但内部统计数据显示出明显的差异。我们认为这是正常的，因为人们不能指望从这样一个简单的机械类比中对一个复杂的社会经济现象有更深的理解。总之，我们认为，本文介绍的模型主要对经济物理学的跨学科领域具有教学价值。它提供了一个吸引人的力学类比，让物理学家在股票市场中研究集体现象。确认该研究得到了欧盟和匈牙利州的支持，由欧洲社会基金会在T\'AMOP 4.2.4的框架内共同资助。A/2-11-1-2012-0001国家卓越计划。附录A.弹簧块系统的动力学细节我们考虑由N个相同质量块m形成的链。

块体由相同的弹簧连接，弹簧常数为k，长度为l。计算机建模面临的挑战是摩擦和弹簧力的量化，以及运动方程的数值积分。m=1、k=1和l=50时使用无量纲单位。选择l的值是为了便于图形显示（弹簧长度对应于50像素）。为了使这些无量纲值符合实际实验情况，单位选择为：[m]=0.1158 kg（质量），[k]=19.8 N/m（弹簧常数），[l]=1.4·10-3米长。其他数量的单位遵循尺寸考虑。因此，时间、速度和力的单位是[t]=√[m] /[k]=0.0765秒，[u]=[l]/[t]=0.0183米/秒，[F]=[k].[l]=0.0277牛顿。链的第i块的运动方程是¨xi=Fe(L-) - Fe(l++Ff弗瑞(L-) - Fe(l+）, （6） 在哪里L-= xi- xi-1.- 50和l+=xi+1- xi- 分别为50和Vr表示块体相对于传送带的相对速度。弹性力Fe和摩擦力Ff定义如下。用于积分运动方程的数值方法见附录B，时间步长选择为dt=0.01时间单位。弹簧中的弹力是线性的，达到一定的变形值，lmax。对于更高的变形，这种依赖关系被假定为指数，负变形的指数（模量）更大（见图9）。因此，Fe(l）=lmax+beb(|l|-lmax）-b、 如果l<-lmax，-l、 如果- lmax6l 6lmax，-[lmax+beb(|l|-lmax）-b] ，如果l>lmax，（7）我们选择的地方lmax=20，b=0.2l<0，b=0.01l>0。

由于考虑了弹簧力的非线性，避免了块体之间的碰撞，因此有了这些选择，该模型对于实验实现来说是非常现实的（参见示例[25]）。同时，我们还没有对这些参数进行系统的更改。对于与速度有关的摩擦，使用简单的库仑摩擦定律。静摩擦力和动摩擦力都与速度模量无关。在合力fexe超过静摩擦力Fst的值之前，块体保持在斗杆状态。对于较高的外力值，块体在存在动摩擦力Fk时开始滑动。我们假设，静摩擦力和动摩擦力之比Fk/Fst=Fst是常数。作用在滑车上的摩擦力F既取决于S（vr），其中S（·）是符号函数，vr是滑车相对于传送带的速度，也取决于作用在滑车上的合力值。在我们的1D设置中，摩擦力方向仅由符号Ff（vr，Fex）给出=-如果vr=0，|Fex |<Fst，-S（vr）fsFst，如果vr，0，（8），其中vr=v-u和v是试块相对于实验室框架的速度。为了使用与我们之前实验[25]中相同的摩擦力值，在无量纲单位中，静态摩擦力取Fst=71.4，我们认为fs=0.45-150-100-50-40 -20  0  20  40  60  80  100lFe(l） 理想弹簧现实弹簧图9：本研究中使用的非线性弹簧力模型，与理想弹簧力模型进行比较。为了再现对数收益的经验测量厚尾（幂律）分布，传送带的速度随机变化，呈幂律分布，在u∈ [5100]速度间隔。

幂函数的指数被选择为α=-3.0. 通过这种方式，该模型再现了实证测量的正收益和负收益的幂律分布[47]。假设DJIA对数收益率的日波动率与我们的模型给出的日波动率相同，我们确定了X（t）动态的δt采样时间。在我们的模拟单元中，我们得到了δt=5（500个积分步）。在这个采样频率下，艺术指数的对数收益率波动率变为σa=0.011，对应于高达三位小数的经验值。附录B——运动方程的数值积分此处简要描述了用于积分运动方程（6）和处理块体不连续粘滑动力学的方法。当一个块体粘在传送带上时，它会在恒定的ELOCITY u处与之一起移动。因此，第i个块体相对于地面的位置通过simplexi（t+dt）=xi（t）+u·dt（9）方程计算。当块体相对于皮带滑动时，基本的Verlet方法xi（t+dt）=2xi（t）- xi（t）- dt）+ai（t）dt+O（dt）（10）用于更新其位置。可以看出，这是一种三阶方法，它也可以扩展到速度空间：vi（t+dt）=xi（t）- xi（t）- dt）dt++dt[11ai（t）- 2ai（t- dt）]+O（dt）。（11） 当相对速度Vrich改变其符号时，发现第i个块粘在皮带上，而当作用在其上的外力总和超过静摩擦力的最大值时，即当-Fst>0。还开发了一种更复杂的随机数值方法来处理粘滑动力学，但发现该方法不会显著改变所给出的结果[25]。参考文献[1]R.Burridge和L.Knopo Off，Bull。地震。Soc。Am，第57卷，341371（1967年）。[2] P。

Bak，C.Tang和K.Wiesenfeld，Phys。牧师。莱特。59, 381 (1987).[3] J.M.卡尔森和J.S.兰格，物理系。牧师。A 406470（1989）。[4] 安徒生，物理系。牧师。B 499981（1994年）。[5] 梁国泰和内达，物理系。牧师。莱特。85, 662 (2000).[6] F.Jarai Szabo，S.Astilean和Z.Neda，Chem。菲斯。莱特。408, 241 (2005).[7] F.Jarai Szabo，E.-A.Horvat，R.Vajtai和Z.Neda，Chem。菲斯。莱特。511, 378 (2011).[8] E-A.霍瓦特、F.贾莱·萨博、Y.布雷谢和Z.内达，分。欧元。J.Phys。10, 926 (2012).[9] M.A.Lebyodkin、Y.Brechet、Y.Estrin和L.P.Kubin，Phys。牧师。莱特。74, 4758 (1995).[10] K.科瓦奇、Y.布雷谢和Z.内达，模特。Simul。马特。Sc.131341（2005年）。[11] F.Jarai Szabo和Z.Neda，Physica A 3915727（2012）。[12] G.Mate，Z.Neda和J.Benedek，《公共科学图书馆综合》第6卷，e16518（2011）[13]M.de Sousa Vieira，Phys。牧师。A 466288（1992年）。[14] M.de Sousa Vieira和H.J.Herrmann，Phys。牧师。E 494534（1994年）。[15] T.鲍姆伯格、C.卡罗利、B.佩林和O.朗辛，Phys。牧师。E 514005（1995年）。[16] G.L.瓦康塞洛斯，物理系。牧师。莱特。76, 4865 (1996).[17] F-J.埃尔默，J.菲斯。A-数学。Gen.30，6057（1997年）。[18] H.Sakaguchi，J.Phys。Soc。Jpn。72, 69 (2003).[19] A.约翰森、P.戴蒙、C.埃勒加德、J.S.拉森和H.H.鲁格，Phys。牧师。E 484779（1993年）。[20] M.de Sousa Vieira和A.J.Lichtenberg，Phys。牧师。E 531441（1996）。[21]C.V.Chianca，J.S.Sa Martins和P.M.C.de Oliveira，欧洲。菲斯。J.B 68549（2009）。[22]苏萨·维埃拉先生，物理系。莱特。A 198407（1995年）。[23]J.Szkutnik和K.Kulakowski，国际国防部J。菲斯。C 13、41（2002年）。[24]B.Erickson，B.Birnir和D.Lavalley，非线性过程。Geoph。15, 1 (2008).【25】B.Sandor，F.Jarai Szabo，T.Tel和Z.Neda，Phys。牧师。E 87042920（2013年）。[26]R.Mantegna和H.E。