分数阶梯形法的收敛性与步长选择_第1页
分数阶梯形法的收敛性与步长选择_第2页
分数阶梯形法的收敛性与步长选择_第3页
分数阶梯形法的收敛性与步长选择_第4页
分数阶梯形法的收敛性与步长选择_第5页
已阅读5页,还剩19页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

19/24分数阶梯形法的收敛性与步长选择第一部分分数阶梯形法的收敛性条件 2第二部分稳定性因子对收敛性的影响 4第三部分步长选择对收敛速度的影响 7第四部分局部收敛与全局收敛的区分 10第五部分不同步长选择算法的比较 12第六部分自适应步长选择策略 14第七部分多步长收敛性分析 16第八部分初始条件对收敛性的影响 19

第一部分分数阶梯形法的收敛性条件关键词关键要点主题名称:分数阶梯形法的收敛性条件

1.初始一致性:分数阶梯形法的收敛性取决于初始状态的选取。如果初始状态满足一定条件,称为一致性条件,则该方法收敛到方程式的解。

2.数值稳定性:分数阶梯形法的数值稳定性受步长选择的影响。合适的步长可以提高计算稳定性,避免累积误差过大。

3.步长条件:分数阶梯形法的收敛性与步长大小密切相关。存在一个临界步长值,当步长超过该值时,该方法发散,而当步长较小时,该方法缓慢收敛。

主题名称:收敛阶数

分数阶梯形法的收敛性条件

分数阶梯形法是一种数值计算方法,用于求解分数阶微分方程组。其收敛性由以下条件保证:

一致性条件:

对于分数阶微分方程组:

```

```

其中,$y(t)=[y_1(t),y_2(t),...,y_m(t)]^T$是待求解的未知函数向量,$0<α_i<1$是分数阶导数阶数,$f_i(t,y(t))$是连续函数。

分数阶梯形法的收敛性要求一致性条件满足:

```

f_i(t,y(t))\inC[a,b],i=1,2,...,m

```

这意味着函数$f_i(t,y(t))$在区间$[a,b]$上连续。

稳定性条件:

分数阶梯形法的稳定性条件由以下不等式给出:

```

h^α≤cT(α),i=1,2,...,m

```

其中,$h$是时间步长,$T(α)$是Mittag-Leffler函数,定义为:

```

```

对于$\alpha>0$,Mittag-Leffler函数满足:

```

```

常数$c$由下式确定:

```

```

收敛性定理:

如果一致性条件和稳定性条件都满足,那么分数阶梯形法的误差满足以下估计:

```

|y(t)-y_h(t)|≤Ch^α,t∈[a,b]

```

其中,$y(t)$是分数阶微分方程组的精确解,$y_h(t)$是分数阶梯形法的数值解,$C$是一个常数。

步长选择:

为了确保分数阶梯形法的收敛性,时间步长$h$的选择至关重要。根据稳定性条件,时间步长需要满足:

```

```

在实际应用中,通常选择接近这个界限的值作为时间步长,例如:

```

```

这种选择有助于平衡计算精度和效率。第二部分稳定性因子对收敛性的影响关键词关键要点【稳定性因子对收敛性的影响】

1.稳定性因子的大小决定了梯形方法的稳定性。稳定性因子越小,方法越稳定,解的误差越小。

2.当稳定性因子为1时,即显式梯形方法,该方法对于常系数线性常微分方程是A级稳定的。但对于非线性常微分方程,稳定性因子应该取小于1的值,以避免解发散。

3.对于具有刚性特征的常微分方程,需要选择较小的稳定性因子才能保持收敛性。

1.稳定性因子与步长密切相关。一般来说,较小的步长需要较大的稳定性因子,以保证方法的稳定性。

2.当步长足够小时,稳定性因子可以取更大的值,从而加快收敛速度。

3.因此,在选择稳定性因子时,需要综合考虑步长和常微分方程的特性。

1.对于高阶分数阶常微分方程,稳定性因子也起着重要的作用。

2.在选择稳定性因子时,需要考虑分数阶导数的阶数。

3.一般来说,分数阶导数的阶数越高,所需的稳定性因子越小。

1.除了常系数常微分方程,稳定性因子对非线性常微分方程的收敛性也有影响。

2.对于非线性常微分方程,稳定性因子需要满足一定的条件才能保证收敛性。

3.这些条件通常与常微分方程的局部Lipschitz连续性和非线性项的增长率有关。

1.稳定性因子对收敛性影响的最新研究表明,在某些情况下,较大的稳定性因子可以加速收敛速度。

2.这种现象被称为"反稳定现象"。

3.反稳定现象的出现与常微分方程的刚性特征和步长有关。

1.在实际应用中,选择合适的稳定性因子至关重要。

2.对于一般常微分方程,推荐使用稳定性因子介于0.5和1之间的值。

3.对于具有刚性特征的常微分方程,需要使用较小的稳定性因子,如0.1或更小。稳定性因子对收敛性的影响

分数阶梯形法是一种数值积分方法,其收敛性受稳定性因子的影响。稳定性因子是一个常数,反映积分区域函数的特性对数值积分结果的影响。

稳定性因子的定义

对于分数阶梯形法,稳定性因子定义为:

```

γ=h^α/Γ(1+α)

```

其中:

*h是步长

*α是分数阶(0<α<1)

*Γ是Γ函数

稳定性因子对收敛性的影响

稳定性因子对分数阶梯形法的收敛性具有至关重要的影响:

*当γ<1时,方法是收敛的。在稳定区域内,数值积分结果随着步长h的减小而收敛到真实值。

*当γ>1时,方法是发散的。在不稳定区域内,数值积分结果随着步长h的减小而远离真实值。

*当γ=1时,方法的收敛性取决于被积函数的性质。对于某些函数,方法可能收敛,而对于其他函数,方法可能发散。

稳定区域

稳定区域由满足以下条件的步长范围定义:

```

0<h<h_max=Γ(1+α)/sup(|f''(x)|)^(1/α)

```

其中:

*f(x)是被积函数

*sup(|f''(x)|)是f''(x)的上确界

步长选择的策略

为了确保分数阶梯形法的收敛性,应选择一个位于稳定区域内的步长h。以下是选择步长的几种策略:

*基于稳定性因子的策略:根据稳定性因子γ的选择步长。例如,可以将γ设置为接近但不超过1的值。

*基于误差估计的策略:通过计算不同步长下的数值积分结果的误差估计来选择步长。选择使误差估计最小的步长。

*自适应步长控制策略:使用算法自动调整步长以保持稳定性。

例子

考虑积分函数f(x)=x^2,其f''(x)=2。

*当α=0.5时,稳定区域为0<h<1.414。

*当α=0.8时,稳定区域为0<h<0.707。

如果选择步长h=0.5,则对于α=0.5,方法是收敛的,而对于α=0.8,方法是发散的。

结论

稳定性因子对于分数阶梯形法的收敛性至关重要。通过选择一个位于稳定区域内的步长,可以确保方法的收敛性。有几种策略可以用来选择步长,包括基于稳定性因子、误差估计和自适应步长控制。第三部分步长选择对收敛速度的影响分数阶梯形法的收敛速度与步长选择

步长选择对收敛速度的影响

步长选择是影响分数阶梯形法收敛速度的关键因素。步长太小,会导致计算过程缓慢;步长太大,则会导致数值不稳定。因此,选择合适的步长对于提高计算效率和精度至关重要。

收敛速度的度量

收敛速度一般用局部截断误差(LTE)来衡量。LTE是分数阶导数的数值逼近值与真值之间的误差。步长越小,LTE也越小。

步长与LTE的关系

对于分数阶梯形法,LTE与步长的关系可以表示为:

```

LTE=O(h^p)

```

其中:

*h为步长

*p为分数阶导数的阶数

从该式可以看出,LTE与步长的p次方成正比。因此,为了获得更高的精度,需要缩小步长。

最优步长

理论上,最优步长可以最大限度地减少LTE,从而提高收敛速度。最优步长取决于分数阶导数的阶数和被积分函数的平滑度。

自适应步长策略

在实际应用中,由于被积分函数的平滑度通常未知,因此使用自适应步长策略可以动态调整步长,以适应不同的函数特性。自适应步长策略通过监测LTE来调整步长:当LTE较大时,减小步长;当LTE较小时,增大步长。

步长选择的经验原则

虽然没有通用的步长选择规则,但以下经验原则可以作为指导:

*对于光滑函数,可以使用较大的步长。

*对于非光滑函数,需要使用较小的步长。

*对于高阶分数阶导数,需要使用较小的步长。

步长选择的例子

考虑如下分数阶积分方程:

```

_0D_x^αy(x)=f(x),0≤x≤1

```

其中:

*α为分数阶导数的阶数

*f(x)为已知函数

对于不同的α值,最优步长如下:

|α|最优步长|

|||

|0.5|0.01-0.05|

|1.0|0.001-0.01|

|1.5|0.0001-0.001|

从表中可以看出,随着分数阶导数阶数的增加,最优步长相应地减小。

结论

步长选择是影响分数阶梯形法收敛速度的重要因素。通过选择合适的步长,可以提高计算效率和精度。自适应步长策略可以动态调整步长,以适应不同的函数特性。通过遵循经验原则和考虑具体的积分方程,可以为分数阶梯形法选择最优步长。第四部分局部收敛与全局收敛的区分关键词关键要点局部收敛与全局收敛的区分

主题名称:局部收敛

1.局部收敛是算法在初始点附近某个区域内收敛,但不能保证在整个定义域内收敛。

2.局部收敛受到初始点的选择影响,不同的初始点可能导致不同的局部极值。

3.局部收敛算法可能陷入局部极值或鞍点,无法找到全局最优解。

主题名称:全局收敛

局部收敛与全局收敛的区分

分数阶梯形法是一种用于求解分数阶微分方程组的数值方法。收敛性是指分数阶梯形法在特定条件下逼近真实解的程度。

局部收敛和全局收敛是收敛性的两个不同概念。

局部收敛

如果分数阶梯形法在某些初始条件和步长范围内逼近真实解,则称其局部收敛。局部收敛意味着该方法在这些条件下产生了可接受的近似值,但并不保证它在所有可能的条件下都会收敛。

全局收敛

如果分数阶梯形法在所有可能的初始条件和步长范围内都逼近真实解,则称其全局收敛。全局收敛是一个更强的收敛性标准,因为它确保该方法在任何情况下都能产生可接受的近似值。

收敛性分析

分数阶梯形法的收敛性分析涉及以下因素:

*阶数:分数阶阶数决定了方程的复杂性,并且会影响收敛性。

*步长:步长是分数阶梯形法中使用的增量,选择合适的步长至关重要。

*初始条件:初始条件是方程解的起点。

*稳定性:稳定性是指分数阶梯形法是否产生有界的解,这与步长选择有关。

步长选择对收敛性的影响

步长选择在分数阶梯形法的收敛性中起着至关重要的作用。如果步长过小,则计算成本太高,而如果步长过大,则可能会导致不稳定或发散。

局部收敛的步长选择

对于局部收敛性,可以采用以下准则选择步长:

*经验法则:根据经验,步长通常选为阶数和步数的函数。

*自适应步长选择:根据先前的计算结果动态调整步长,以优化收敛性。

全局收敛的步长选择

对于全局收敛性,需要使用更严格的步长选择准则:

*稳定性分析:通过分析分数阶梯形法的稳定性条件确定步长范围。

*一致收敛准则:确保在所有可能的情况下都满足收敛条件。

总结

局部收敛性表示分数阶梯形法在某些条件下逼近真实解,而全局收敛性则表示该方法在所有条件下都收敛。步长选择对收敛性至关重要,必须根据方程的性质和收敛性目标进行选择。通过适当的步长选择,分数阶梯形法可以有效地求解各种分数阶微分方程组。第五部分不同步长选择算法的比较不同步长选择算法的比较

引言

分数阶梯形法的步长选择对于其收敛性至关重要。在不同的应用场景中,不同的步长选择算法可以产生显著不同的收敛速度和精度。本文将比较常用的步长选择算法,包括固定步长、自适应步长和混合步长方法。

固定步长方法

固定步长方法使用一个在整个积分过程中保持不变的步长大小。这种方法简单且易于实现,但可能难以平衡收敛速度和精度。对于平滑函数,固定步长方法可能收敛较快,但对于非平滑函数或快速变化的函数,它可能会导致不稳定的计算或过度振荡。

自适应步长方法

自适应步长方法通过根据函数的局部特征动态调整步长大小来提高收敛性。这些方法通常基于局部误差估计,其中当前计算的误差被用于确定下一个步长的最佳大小。自适应步长方法可以自动适应函数的复杂性和变化率,从而提高收敛效率。

混合步长方法

混合步长方法结合了固定步长和自适应步长的优点。它们从一个固定步长开始,然后根据局部误差估计动态调整步长大小。这种方法可以平衡收敛速度和稳定性,同时避免自适应步长方法的过度调整。

不同步长选择算法的比较

表1:不同步长选择算法的比较

|算法类型|优点|缺点|

||||

|固定步长|简单、易于实现|对于非平滑函数收敛缓慢|

|自适应步长|收敛速度快、精度高|可能过度调整步长大小|

|混合步长|平衡收敛速度和稳定性|比固定步长方法复杂|

具体比较

在实际应用中,不同步长选择算法的性能会根据具体问题而有所不同。对于平滑函数,固定步长方法通常足以获得合理的精度。然而,对于非平滑函数或快速变化的函数,自适应步长方法可以显著提高收敛速度。混合步长方法通常在各种问题上提供良好的折衷方案。

下表提供了不同步长选择算法在特定问题上的定量比较:

表2:不同步长选择算法在特定问题上的比较

|问题|固定步长|自适应步长|混合步长|

|||||

|平滑函数|快速收敛|收敛速度稍慢|收敛速度中等|

|非平滑函数|收敛缓慢|快速收敛|收敛速度良好|

|快速变化的函数|不稳定|稳定|稳定|

结论

分数阶梯形法的步长选择对于收敛性至关重要。不同的步长选择算法在不同情况下表现出不同的性能。固定步长方法简单但对于非平滑函数收敛缓慢,自适应步长方法收敛速度快但可能过度调整步长大小,混合步长方法则平衡了收敛速度和稳定性。具体选择哪种算法取决于所解决问题的特征。第六部分自适应步长选择策略关键词关键要点【自适应步长选择策略】,

1.基于误差估计:通过估计局部截断误差或全局截断误差,动态调整步长以控制误差在可接受范围内。

2.基于收敛率:监测收敛率,当收敛缓慢时,增大步长以加速收敛;当收敛过快时,减小步长以提高精度。

3.基于稳定性:步长选择还应考虑方法的稳定性,过大的步长可能导致数值解发散或不稳定。

【趋势和前沿】,自适应步长选择策略

分数阶梯形法是一种用于求解分数阶微分方程组的数值方法。自适应步长选择策略旨在动态调整步长大小,以平衡收敛性、稳定性和计算效率之间的关系。以下是分数阶梯形法中常用的几种自适应步长选择策略:

局部截断误差估计

这种策略通过估计局部的截断误差来调整步长大小。在每一步中,使用较高阶的梯形法计算解的近似值,并将其与较低阶的梯形法计算的值进行比较。如果局部截断误差超过某个容差,则减小步长大小;如果误差小于容差,则增大步长大小。

Romberg准则

Romberg准则使用一系列梯形法近似值来估计解。在每一步中,使用不同的步长大小计算解的近似值,并使用插值技术获得更准确的近似值。如果相邻两次迭代的近似值之间的差值超过某个容差,则减小步长大小。

PID控制器

PID控制器(比例-积分-微分)是一种反馈控制系统,用于调整步长大小。控制器监控解的误差并根据误差的比例、积分和微分值生成一个控制信号。控制信号用于调整步长大小,以将误差保持在某个预定的范围内。

步长倍增策略

步长倍增策略是一种简单的自适应步长选择策略,其中步长大小在每一步中增加或减少一个倍数。倍数因子根据解的收敛行为进行调整。如果解收敛良好,则增加倍数因子;如果解不收敛,则减小倍数因子。

基于梯度的策略

基于梯度的策略使用梯度信息来选择步长大小。梯度信息表示解对自变量的导数,它可以用来估计解的变化率。如果梯度较大,则减小步长大小;如果梯度较小,则增大步长大小。

这些自适应步长选择策略的目的是在保证收敛性的同时,尽可能高效地求解分数阶微分方程组。通过动态调整步长大小,可以减少计算时间,同时保持解的准确性。

选择策略的因素

选择最佳的自适应步长选择策略取决于以下因素:

*方程组的非线性度

*解的收敛特性

*计算资源的可用性

对于非线性方程组或收敛速度较慢的解,局部截断误差估计或Romberg准则可能是更好的选择。对于线性方程组或收敛速度较快的解,步长倍增策略或基于梯度的策略可能是更有效率的。

试验和误差

由于没有一种自适应步长选择策略适用于所有情况,因此通常需要进行试验和误差才能确定针对特定方程组最有效的策略。通过比较不同策略的性能,可以找到最佳的策略,以在收敛性、稳定性和计算效率之间取得适当的平衡。第七部分多步长收敛性分析多步长收敛性分析

分数阶梯形法的多步长收敛性分析是指研究在使用多个步长的情况下,数值解的误差随步长的变化规律。其主要目的是确定最佳步长选择,以在保证数值稳定性和收敛性前提下获得较高的精度。

基本原理

多步长收敛性分析基于以下原理:

*数值解的误差由局部截断误差和全局截断误差组成。

*局部截断误差与步长成正比,而全局截断误差与步长成高阶幂次方(通常为步长的平方)。

因此,在步长较小时,局部截断误差起主要作用,导致误差随步长减小而减小。然而,当步长增大时,全局截断误差变得显着,导致误差随步长增大而增大。

收敛性阶数

分数阶梯形法的多步长收敛性阶数是指在步长较小时,误差随步长减小的速率。收敛性阶数越高,数值解的精度就越高。

分数阶梯形法的收敛性阶数取决于所使用的公式。例如,对于二阶梯形法,收敛性阶数为2;对于四阶梯形法,收敛性阶数为4。

最佳步长选择

最佳步长选择取决于具体问题和使用的分数阶梯形公式。一般来说,以下准则可以指导步长选择:

*局部截断误差控制:选择一个步长,使得局部截断误差在可接受范围内。

*稳定性考虑:选择一个步长,使得数值解稳定,避免出现发散或振荡。

*计算效率:在满足精度和稳定性要求的情况下,选择一个能够最大化计算效率的步长。

自适应步长选择

自适应步长选择是一种动态调整步长的技术,以在不同的求解阶段满足收敛性和计算效率要求。其基本思想是根据局部截断误差的估计值来调整步长。

自适应步长选择算法有多种,例如:

*控制局部截断误差:根据局部截断误差的估计值调整步长,以将其保持在可接受范围内。

*预测局部截断误差:使用高阶公式预测局部截断误差,并根据预测结果调整步长。

*基于稳定性的步长选择:使用稳定性指标(如增长因子)来评估数值解的稳定性,并根据稳定性指标调整步长。

数值实验

通过数值实验可以验证分数阶梯形法的多步长收敛性。通常,以下步骤可以帮助确定最佳步长选择:

1.使用不同的步长求解方程。

2.计算误差(与精确解的误差)。

3.绘制误差随步长的关系图。

4.根据收敛性阶数和计算效率选择最佳步长。

结论

分数阶梯形法的多步长收敛性分析对于选择最佳步长以获得准确和有效的数值解至关重要。通过理解收敛性阶数、局部截断误差控制和自适应步长选择,可以优化分数阶梯形法的性能,以满足特定问题的要求。第八部分初始条件对收敛性的影响关键词关键要点【初始条件对收敛性的影响】

1.初始条件对分数阶梯形法的收敛性有显著影响,选择合适的初始条件至关重要。

2.对于求解非线性方程组,初始条件通常影响迭代过程的稳定性和收敛速度。

3.对于求解分数微分方程,初始条件决定了分数解的全局行为和渐近性质。

【初始条件选择原则】

初始条件对分数阶梯形法的收敛性的影响

在分数阶梯形法中,初始条件的选择会显著影响收敛性。理想情况下,对于给定的分数阶微分方程,初始条件应满足以下条件:

-一致性:初始条件应与微分方程的阶数一致。

-连续性:积分阶数小于或等于微分阶数的初始条件应连续。

#初始条件与收敛性之间的关系

研究表明,当初始条件不满足一致性和连续性条件时,分数阶梯形法可能会出现以下收敛性问题:

-发散:初始条件不满足一致性条件时,数值解可能会发散到无穷大。

-振荡:初始条件不满足连续性条件时,数值解可能会出现振荡现象,无法收敛到准确解。

#选择合适初始条件的策略

为了确保分数阶梯形法的收敛性,建议采用以下策略选择合适的初始条件:

1.分析微分方程:确定微分方程的阶数和积分阶数,并确保初始条件与之相一致。

2.使用辅助方程:对于积分阶数小于或等于微分阶数的初始条件,可以求解以下辅助方程:

```

D^αu(t)=f(t),t>0,α∈(0,1]

```

其中,u(t)是辅助变量,f(t)是微分方程的右端项。辅助方程的解可以作为初始条件。

3.使用插值方法:如果辅助方程难以求解,可以使用插值方法来估计初始条件。例如,可以利用已知解的离散数据进行拉格朗日插值或样条插值。

4.利用物理意义:如果微分方程来自一个物理模型,可以利用模型的物理意义来推断合理的初始条件。例如,在弹簧-质量系统中,初始速度和位移是合理的初始条件。

#步长选择的注意事项

除了初始条件外,步长选择也是影响分数阶梯形法收敛性的重要因素。以下注意事项应予以考虑:

-稳定性条件:对于显式分数阶梯形法,步长h必须满足稳定性条件:

```

h<c*(T^(-α)+h^(-α))^(-1/α)

```

其中,c是一个常数,通常取值为0.25或0.5,T是最终时间。

-收敛性:步长越小,数值解的精度越高。但过小的步长可能会导致计算成本过高。因此,需要在精度和效率之间进行权衡。

-自适应步长方法:为了提高效率和适应性,可以使用自适应步长方法,根据局部误差调整步长。这可以避免不必要的计算,同时确保收敛性。关键词关键要点主题名称:步长选择对收敛速度的影响

关键要点:

1.较小步长有利于提高精度,但会减慢收敛速度。较小步长可以更精确地逼近函数的局部行为,但也会导致迭代次数增加,从而减慢整体收敛速度。

2.较大步长可以加速收敛速度,但可能导致不稳定或发散。较大步长可以减少迭代次数,但如果步长过大,可能会导致数值积分不稳定或发散,从而无法得到有意义的结果。

3.自适应步长选择方法可以优化收敛速度。自适应步长选择方法通过根据每次迭代的误差或稳定性来调整步长,可以在提高精度和收敛速度之间取得平衡。

主题名称:步长选择对准确性的影响

关键要点:

1.步长大小与近似误差成正比。较大的步长会导致积分近似值与实际值之间更大的误差,而较小的步长可以减少误差。

2.误差估计可以指导步长选择。通过估计积分过程中的误差,可以确定适当的步长以平衡精度和计算成本。

3.自适应误差控制方法可以保证精度。自适应误差控制方法通过动态调整步长以满足预定的精度要求,从而确保近似值满足可接受的误差范围。

主题名称:步长选择对计算成本的影响

关键要点:

1.较大的步长减少了迭代次数,但增加了每次迭代的计算成本。较大步长虽然可以加速收敛速度,但可能需要更复杂的积分算法,从而增加每次迭代的计算成本。

2.较小的步长增加了迭代次数,但降低了每次迭代的计算成本。较小的步长虽然可以提高精度,但会增加迭代次数,从而增加整体计算成本。

3.

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论