期刊文献+
共找到349篇文章
< 1 2 18 >
每页显示 20 50 100
A New Technique for Constructing Higher-order Iterative Methods to Solve Nonlinear Systems
1
作者 XIAO Xiaoyong 《应用数学》 北大核心 2025年第3期762-774,共13页
In this paper,a new technique is introduced to construct higher-order iterative methods for solving nonlinear systems.The order of convergence of some iterative methods can be improved by three at the cost of introduc... In this paper,a new technique is introduced to construct higher-order iterative methods for solving nonlinear systems.The order of convergence of some iterative methods can be improved by three at the cost of introducing only one additional evaluation of the function in each step.Furthermore,some new efficient methods with a higher-order of convergence are obtained by using only a single matrix inversion in each iteration.Analyses of convergence properties and computational efficiency of these new methods are made and testified by several numerical problems.By comparison,the new schemes are more efficient than the corresponding existing ones,particularly for large problem sizes. 展开更多
关键词 systems of nonlinear equation Order of convergence Higher-order method Extended newton iteration Computational efficiency
在线阅读 下载PDF
基于Newton/Gauss-Seidel迭代的DGM隐式方法 被引量:8
2
作者 刘伟 张来平 +2 位作者 赫新 贺立新 张涵信 《力学学报》 EI CSCD 北大核心 2012年第4期792-796,共5页
在Newton迭代方法的基础上,对高阶精度间断Galerkin有限元方法 (DGM)的时间隐式格式进行了研究.Newton迭代法的优势在于收敛效率高效,并且定常和非定常问题能够统一处理,对于非定常问题无需引入双时间步策略.为了避免大型矩阵的求逆,采... 在Newton迭代方法的基础上,对高阶精度间断Galerkin有限元方法 (DGM)的时间隐式格式进行了研究.Newton迭代法的优势在于收敛效率高效,并且定常和非定常问题能够统一处理,对于非定常问题无需引入双时间步策略.为了避免大型矩阵的求逆,采用一步Gauss-Seidel迭代和Matrix-free技术消去残值Jacobi矩阵的上、下三角矩阵,从而只需计算和存储对角(块)矩阵.对角(块)矩阵采用数值方法计算.空间离散采用Taylor基,其优势在于对于任意形状的网格,基函数的形式是一致的,有利于在混合网格上推广.利用该方法,数值模拟了Bump绕流和NACA0012翼型绕流.计算结果表明,与显式的Runge-Kutta时间格式相比,隐式格式所需的迭代步数和CPU时间均在很大程度上得到减少,计算效率能够提高1~2个量级. 展开更多
关键词 间断Galerkin有限元 Taylor基函数 newton迭代 Gauss-seidel迭代 时间隐式方法
在线阅读 下载PDF
基于Newton-Raphson迭代与PSO数字图像相关方法 被引量:15
3
作者 杜亚志 王学滨 《计算机工程与应用》 CSCD 2012年第34期184-189,共6页
鉴于基于Newton-Raphson(N-R)迭代的数字图像相关方法对迭代初值的敏感性问题,提出了一种基于N-R迭代与粒子群优化(PSO)算法的数字图像相关方法。该方法利用了PSO算法中的全局搜索能力与N-R迭代中的局部搜索能力,通过它们的交替迭代以... 鉴于基于Newton-Raphson(N-R)迭代的数字图像相关方法对迭代初值的敏感性问题,提出了一种基于N-R迭代与粒子群优化(PSO)算法的数字图像相关方法。该方法利用了PSO算法中的全局搜索能力与N-R迭代中的局部搜索能力,通过它们的交替迭代以改善初值;以改善后的初值进行最后的N-R迭代。经检验,该方法比N-R迭代对初值的要求宽松。这是由于在N-R迭代与PSO算法多次交替迭代之后,位移和应变的初值都得到了一定的改善。对相似材料模型受载破坏后表面上的点的位移和应变进行了计算,获得了较好的结果。 展开更多
关键词 数字图像相关方法 newton-Raphson迭代 粒子群优化算法 初值
在线阅读 下载PDF
基于Newton-Raphson迭代的动压气体止推箔片轴承特性研究 被引量:5
4
作者 闫佳佳 刘占生 +1 位作者 王铮 高佳颖 《汽轮机技术》 北大核心 2017年第2期116-120,124,共6页
以波箔型气体止推箔片轴承为研究对象,给出一种快速收敛的静特性数值求解方法,结合Newton-Raphson迭代法和有限差分法,耦合求解Reynolds方程和润滑膜厚度方程,研究气体止推箔片轴承的静态性能,分析结构参数对气体止推箔片轴承轴向承载... 以波箔型气体止推箔片轴承为研究对象,给出一种快速收敛的静特性数值求解方法,结合Newton-Raphson迭代法和有限差分法,耦合求解Reynolds方程和润滑膜厚度方程,研究气体止推箔片轴承的静态性能,分析结构参数对气体止推箔片轴承轴向承载力、黏性摩擦力矩的影响规律。建立转子止推盘倾斜的润滑膜厚度方程,求解获得轴承压力分布规律,分析转子止推盘倾斜对气体止推箔片轴承静态性能的影响。结果表明:采用基于NewtonRaphson迭代的数值求解方法收敛速度较快;气体止推箔片轴承的轴向承载力和黏性摩擦力矩均随箔片变形刚度的增大而增大;止推盘倾斜时气体止推箔片轴承的静态性能增大,且倾斜角度越大静态特性结果也越大。 展开更多
关键词 气体止推箔片轴承 newton—Raphson迭代 轴向承载力 黏性摩擦力矩 润滑膜厚度 止推盘倾斜
在线阅读 下载PDF
修正Newton-Raphson方法和弧长法在结构推覆分析中的联合应用
5
作者 宋少沪 曾玉山 张佩军 《上海大学学报(自然科学版)》 CAS CSCD 北大核心 2006年第1期101-106,共6页
基于性能抗震设计理论将成为未来建筑抗震设计的主流,推覆分析和非线性时程分析是基于性能抗震设计的主要手段.该文在前人的理论基础上,为了提高分析的精度,搜索结构的软化段,采用了Newton-Raphson方法和弧长法联合应用的非线性方程组... 基于性能抗震设计理论将成为未来建筑抗震设计的主流,推覆分析和非线性时程分析是基于性能抗震设计的主要手段.该文在前人的理论基础上,为了提高分析的精度,搜索结构的软化段,采用了Newton-Raphson方法和弧长法联合应用的非线性方程组的求解技术,对高层结构推覆分析进行了较为深入的探讨. 展开更多
关键词 推覆分析 非线性方程 newton-Raphson方法 弧长法
在线阅读 下载PDF
ITER馈线S弯结构的应力分析和优化设计
6
作者 徐薇薇 宋云涛 +3 位作者 路建武 王忠伟 程勇 刘旭峰 《核聚变与等离子体物理》 CAS CSCD 北大核心 2012年第3期250-254,共5页
根据能量方法原理,通过单位力法和卡氏定理分析方法,理论计算得出了给定位移变形条件下沿S弯的弯矩、应力分布。分析表明,沿S弯结构的应力主要由热收缩变形产生,但重力影响也不可忽略,考虑重力影响后沿S弯最大应力较不考虑重力时增大了... 根据能量方法原理,通过单位力法和卡氏定理分析方法,理论计算得出了给定位移变形条件下沿S弯的弯矩、应力分布。分析表明,沿S弯结构的应力主要由热收缩变形产生,但重力影响也不可忽略,考虑重力影响后沿S弯最大应力较不考虑重力时增大了10%。根据该理论原理编程计算发现,在S弯高度受限的条件下,圆弧段半径越小,结构最大应力越小,实际设计中可根据强度要求结合加工工艺选择合适的圆弧半径。该分析方法为有效优化结构设计提供了理论指导。 展开更多
关键词 iter馈线s弯结构 能量方法 应力分析 结构优化
在线阅读 下载PDF
Gauss-Newton法的半局部收敛性
7
作者 张文红 李冲 《东南大学学报(自然科学版)》 EI CAS CSCD 北大核心 2001年第5期135-139,共5页
设f:Rn→Rm 是Frechet可微的 ,m≥n .则非线性最小二乘问题可描述为下面的极小化问题 :minF(x) :=12 f(x) Tf(x) .Gauss Newton法是求解非线性最小二乘问题的最基本的方法之一 ,其n + 1步迭代定义为 :xn + 1=xn - f′(xn) Tf′(x) -1f′... 设f:Rn→Rm 是Frechet可微的 ,m≥n .则非线性最小二乘问题可描述为下面的极小化问题 :minF(x) :=12 f(x) Tf(x) .Gauss Newton法是求解非线性最小二乘问题的最基本的方法之一 ,其n + 1步迭代定义为 :xn + 1=xn - f′(xn) Tf′(x) -1f′(xn) Tf(xn) .本文主要研究解非线性最小二乘问题的Gauss Newton法的半局部收敛性 .假设f(x)在B(x0 ,r)内连续可导且f′(x0 )满秩 ,若f的导数满足Lipschitz连续F′(x) -f′(x′)≤γx -x′ , x ,x′∈B(x0 ,r) .在一个关于初始点x0 的判断准则c =f(x0 ) ,β =f′T(x0 )f′(x0 ) -1f′(x0 ) T ,β2 cγ <1 1 0下 ,Gauss Newton法产生的序列 {xn}收敛到一个驻点x ,从而给出了Gauss Newton法的半局部收敛性 . 展开更多
关键词 非线性最小二乘问题 Garuss-newton 半局部收敛性 最优化方法
在线阅读 下载PDF
解非光滑方程组的Newton-GMRES算法
8
作者 孟泽红 张建军 《华东地质学院学报》 2003年第2期121-122,共2页
给出了求解非光滑方程组的Newton GMRES迭代法。该方法在求解半光滑方程组时,不需要计算广义Jacobi矩阵,同时使求解相应广义Newton方程组也变得容易。尤其对于大型问题,该方法特别适用。数值例子显示了这种方法的有效性。
关键词 非光滑方程组 newton—GMREs算法 迭代法
在线阅读 下载PDF
重力异常界面反演的Gauss-Newton方法及其隐式迭代实现 被引量:2
9
作者 李会元 许梦杰 +1 位作者 彭瑞仁 陈基明 《上海大学学报(自然科学版)》 CAS CSCD 1999年第1期29-32,共4页
首先根据重力异常积分公式及重力异常实测数据将确定异常源的密度界面问题归结为求非线性泛函极小化的变分问题,并在Hilbert空间中应用Gauss-Newton迭代方法进行求解.在实际计算中采用有限维试验函数空间对模型进... 首先根据重力异常积分公式及重力异常实测数据将确定异常源的密度界面问题归结为求非线性泛函极小化的变分问题,并在Hilbert空间中应用Gauss-Newton迭代方法进行求解.在实际计算中采用有限维试验函数空间对模型进行离散,每次迭代的增量由隐式迭代方法给出。 展开更多
关键词 界面反演 隐式迭代 高斯-牛顿法 重力异常界面
在线阅读 下载PDF
一类非线性代数方程组的Newton-Triangle Splitting迭代法 被引量:3
10
作者 胡纪洋 王川龙 温瑞萍 《工程数学学报》 CSCD 北大核心 2015年第1期29-38,共10页
Triangle Splitting迭代方法是求解大型稀疏非Hermitian正定线性代数方程组的一种有效迭代算法.为了有效求解大型稀疏且Jacobi矩阵为非Hermitian正定的非线性代数方程组,本文将Triangle Splitting迭代方法作为不精确Newton方法的内迭代... Triangle Splitting迭代方法是求解大型稀疏非Hermitian正定线性代数方程组的一种有效迭代算法.为了有效求解大型稀疏且Jacobi矩阵为非Hermitian正定的非线性代数方程组,本文将Triangle Splitting迭代方法作为不精确Newton方法的内迭代求解器,构造了不精确Newton-Triangle Splitting迭代方法.在适当的约束条件下,给出了该方法的两类局部收敛性定理.通过数值实验结果验证了该方法的可行性和有效性,并说明了该方法在计算时间和迭代次数方面比Newton-BTSS迭代方法更有优势. 展开更多
关键词 TRIANGLE splitting迭代方法 非线性代数方程组 不精确newton方法 局部收敛性
在线阅读 下载PDF
Novel Newton’s learning algorithm of neural networks 被引量:2
11
作者 Long Ning Zhang Fengli 《Journal of Systems Engineering and Electronics》 SCIE EI CSCD 2006年第2期450-454,共5页
Newton's learning algorithm of NN is presented and realized. In theory, the convergence rate of learning algorithm of NN based on Newton's method must be faster than BP's and other learning algorithms, because the ... Newton's learning algorithm of NN is presented and realized. In theory, the convergence rate of learning algorithm of NN based on Newton's method must be faster than BP's and other learning algorithms, because the gradient method is linearly convergent while Newton's method has second order convergence rate. The fast computing algorithm of Hesse matrix of the cost function of NN is proposed and it is the theory basis of the improvement of Newton's learning algorithm. Simulation results show that the convergence rate of Newton's learning algorithm is high and apparently faster than the traditional BP method's, and the robustness of Newton's learning algorithm is also better than BP method' s. 展开更多
关键词 newton's method Hesse matrix fast learning BP method neural network.
在线阅读 下载PDF
考虑圆弧修形和交错角的人字齿轮啮合行为分析 被引量:1
12
作者 侯晓燕 王曦 周越 《西安交通大学学报》 北大核心 2025年第4期180-192,共13页
针对人字齿轮制造加工导致左、右齿面存在不同程度的对称性偏差,以及传统轴向位移迭代法受误差尺度限制的问题,以某型高速列车人字齿轮为研究对象,提出了一种新的轴向位移迭代方法,并建立人字齿轮承载接触分析(LTCA)模型。该模型以有限... 针对人字齿轮制造加工导致左、右齿面存在不同程度的对称性偏差,以及传统轴向位移迭代法受误差尺度限制的问题,以某型高速列车人字齿轮为研究对象,提出了一种新的轴向位移迭代方法,并建立人字齿轮承载接触分析(LTCA)模型。该模型以有限元法和有限长线接触理论为基础,在载荷-变形协调方程中考虑不同啮合点处螺旋角的差异,采用牛顿迭代法计算主动齿轮轴向位移,通过LTCA模型分析了轮齿圆弧修形与齿面交错角对啮合激励和载荷分布的影响。研究结果表明:牛顿迭代法可适用于含不同误差类型的人字齿轮LTCA模型,且小误差下牛顿迭代法比传统迭代法的计算时间减少70%,圆弧修形量从2μm增至15μm,时变啮合刚度均值降低38%,啮合误差均值提高39%;当输入扭矩为500 N·m时,两侧齿面中部最大法向啮合力由100 N增大到370 N;随着齿面交错角的增大,时变啮合刚度均值提高,啮合误差均值减小,交错角相位差从0增加到π,时变啮合刚度、啮合误差的幅值分别降低了54%和79%。该研究结果可为人字齿轮啮合稳定性的提高提供新思路。 展开更多
关键词 人字齿轮 承载接触分析 圆弧修形 齿面交错角 轴向位移 牛顿迭代法
在线阅读 下载PDF
改进FastICA算法在谐波检测中的应用 被引量:21
13
作者 汪斌 王年 +3 位作者 蒋云志 程志友 王继 鲍文霞 《电力自动化设备》 EI CSCD 北大核心 2011年第3期135-138,共4页
介绍了一种基于改进的快速独立分量分析(FastICA)算法,并将其引入到谐波检测中。该算法在FastICA算法的基础上对牛顿迭代法进行了改进,使其满足三阶收敛,同时依据负熵极大的独立性准则实现了谐波信号的盲分离。为了更好地逼近真实信号,... 介绍了一种基于改进的快速独立分量分析(FastICA)算法,并将其引入到谐波检测中。该算法在FastICA算法的基础上对牛顿迭代法进行了改进,使其满足三阶收敛,同时依据负熵极大的独立性准则实现了谐波信号的盲分离。为了更好地逼近真实信号,对分离后的信号进行幅值修正,从而完成谐波的检测。仿真实验结果表明了该算法相对于FastICA算法的优点在于减少了迭代次数和加快了收敛速度,同时在谐波分离准确性方面也明显优于FastICA算法。 展开更多
关键词 电力系统 独立分量分析 谐波检测 FAsTICA 牛顿迭代法
在线阅读 下载PDF
改进的高阶收敛FastICA算法 被引量:13
14
作者 季策 胡祥楠 +1 位作者 朱丽春 张志伟 《东北大学学报(自然科学版)》 EI CAS CSCD 北大核心 2011年第10期1390-1393,共4页
高阶收敛的FastICA具有形式简单、收敛速度快的特点,但其对初始值的选择比较敏感,若初始值选择不当很容易影响收敛的效果,甚至造成不收敛的结果.针对这一问题,采用最速下降法对三阶和五阶收敛的FastICA算法进行改进.首先,应用最速下降... 高阶收敛的FastICA具有形式简单、收敛速度快的特点,但其对初始值的选择比较敏感,若初始值选择不当很容易影响收敛的效果,甚至造成不收敛的结果.针对这一问题,采用最速下降法对三阶和五阶收敛的FastICA算法进行改进.首先,应用最速下降法求出初值,再用高阶收敛的FastICA算法求出最优解.语音信号的分离实验表明:改进后的算法对混合信号进行了较好的分离,并且有效地克服了初值敏感性的问题. 展开更多
关键词 独立分量分析 牛顿迭代法 FAsTICA 最速下降法 初值敏感性
在线阅读 下载PDF
Rolls-Royce发动机基线挖掘方法 被引量:30
15
作者 钟诗胜 崔智全 付旭云 《计算机集成制造系统》 EI CSCD 北大核心 2010年第10期2265-2270,共6页
采用基于核函数的多元非线性回归分析方法,对Rolls-Royce公司的发动机基线方程进行了求解,建立了标准化过程中的核函数和性能参数的偏差值模型。把厂家系统计算出的偏差值作为目标值,其他相关参数的测量值作为输入,采用高斯牛顿迭代法... 采用基于核函数的多元非线性回归分析方法,对Rolls-Royce公司的发动机基线方程进行了求解,建立了标准化过程中的核函数和性能参数的偏差值模型。把厂家系统计算出的偏差值作为目标值,其他相关参数的测量值作为输入,采用高斯牛顿迭代法对该多元非线性方程组进行回归分析,能够迅速收敛,准确地求解出模型未知系数,经整理可得出准确的基线方程。通过比较自主基线求解的偏差值与厂家系统求解的偏差值可知,采用上述基线挖掘方法得到的基线方程模型精度能够满足实际应用的要求。 展开更多
关键词 基线 回归分析 核函数 偏差值 高斯迭代法 故障诊断 航空发动机
在线阅读 下载PDF
基于改进FastICA算法的入侵检测样本数据优化方法 被引量:14
16
作者 杜晔 张亚丹 +1 位作者 黎妹红 张大伟 《通信学报》 EI CSCD 北大核心 2016年第1期42-48,共7页
为更好实现对入侵检测样本数据的优化处理,提出了一种改进的快速独立成分分析(Fast ICA)算法,采用基于加权相关系数进行白化处理以减少信息损失,并优化牛顿迭代法使其满足三阶收敛。对算法进行了细致描述,分析了算法的时间复杂度。实验... 为更好实现对入侵检测样本数据的优化处理,提出了一种改进的快速独立成分分析(Fast ICA)算法,采用基于加权相关系数进行白化处理以减少信息损失,并优化牛顿迭代法使其满足三阶收敛。对算法进行了细致描述,分析了算法的时间复杂度。实验结果表明,该方法可有效减少数据信息损失,具有迭代次数少、收敛速度快等优点,可有效提高入侵检测样本数据的优化效率。 展开更多
关键词 入侵检测 快速独立成分分析 数据优化 牛顿迭代法
在线阅读 下载PDF
NURBS插补中的速度规划与参数计算 被引量:16
17
作者 王允森 杨东升 +1 位作者 刘荫忠 孙一兰 《计算机集成制造系统》 EI CSCD 北大核心 2014年第8期1896-1902,共7页
为满足数控机床高速度、高质量加工的需求,提出一种新的非均匀有理B样条曲线插补算法。该算法包括速度规划和参数计算两部分。速度规划部分采用五段S曲线加减速控制方法,能够保证高速运行过程中加速度的连续,使机床运行平稳,避免产生激... 为满足数控机床高速度、高质量加工的需求,提出一种新的非均匀有理B样条曲线插补算法。该算法包括速度规划和参数计算两部分。速度规划部分采用五段S曲线加减速控制方法,能够保证高速运行过程中加速度的连续,使机床运行平稳,避免产生激烈的震颤;参数计算部分应用抛物线插值结合牛顿迭代的方法计算插补参数,将实时插补时产生的进给速度波动控制到理想水平,从而进一步减小机床震颤。仿真实验表明,该算法能够减小机床振动,实现高质量加工。 展开更多
关键词 数控 非均匀有理B样条插补 五段s曲线加减速控制 抛物线插值 牛顿迭代法
在线阅读 下载PDF
基于牛顿迭代法的NURBS曲线插补算法 被引量:17
18
作者 王允森 盖荣丽 +2 位作者 孙一兰 杨东升 徐明子 《组合机床与自动化加工技术》 北大核心 2013年第4期13-17,共5页
通过分析数控机床加工中样条曲线插补原理,提出了一种基于牛顿迭代法的非均匀有理B样条(NURBS)曲线插补算法。该算法首先进行速度控制,保证了加工精度,然后利用牛顿迭代法来计算插补参数,避免了传统方法的大量求导运算。通过将二阶泰勒... 通过分析数控机床加工中样条曲线插补原理,提出了一种基于牛顿迭代法的非均匀有理B样条(NURBS)曲线插补算法。该算法首先进行速度控制,保证了加工精度,然后利用牛顿迭代法来计算插补参数,避免了传统方法的大量求导运算。通过将二阶泰勒展开式化成简单的线性递推方程来预测插补参数的迭代初值,进一步减小了运算开销。仿真对比实验表明,该算法稳定性好,运算精度高,能够满足实时插补的要求,且能够将实时插补产生的速度波动限制到理想的水平。 展开更多
关键词 数控插补 NURBs曲线 速度控制 牛顿迭代法 线性递推方程
在线阅读 下载PDF
基于超松弛因子的高阶复值FastICA算法 被引量:6
19
作者 季策 陈雷 +2 位作者 王艳茹 王丽娟 沙毅 《系统仿真学报》 CAS CSCD 北大核心 2014年第12期2834-2838,共5页
高阶收敛的复值FastICA(CFICA)算法是一种有效的独立分量分析算法,具有收敛速度快、形式简单的特点。其对初始值的选择要求较高,如果初始值选择不当,不仅会影响收敛的效果,甚至有可能造成不收敛的结果。针对这一问题,在高阶收敛的CFICA... 高阶收敛的复值FastICA(CFICA)算法是一种有效的独立分量分析算法,具有收敛速度快、形式简单的特点。其对初始值的选择要求较高,如果初始值选择不当,不仅会影响收敛的效果,甚至有可能造成不收敛的结果。针对这一问题,在高阶收敛的CFICA的基础之上,通过引入超松弛因子对随机产生的初始值进行处理,放宽算法对初始值的要求。在保证收敛速度的前提下,获得了能有效克服初值敏感性的高阶CFICA改进算法。仿真实验结果表明:改进后的算法不依赖于初始值的选择,避免了收敛速度不均衡的现象,提高了算法的收敛性能。 展开更多
关键词 快速独立分量分析 牛顿迭代法 超松弛因子 复值 初值敏感性
在线阅读 下载PDF
大容量链式STATCOM的三角载波移相-开关频率最优PWM法控制策略的研究 被引量:16
20
作者 王立杰 孙玉坤 +1 位作者 孙运全 刘永敏 《继电器》 CSCD 北大核心 2007年第3期41-44,共4页
提出了适用于大容量链式静止无功发生器(STATCOM)的一种新的脉冲控制策略,即三角载波移相-开关频率最优PWM法(PS-SFO-PWM)。这种新的PWM方法兼具三角载波移相PWM法(PSPWM)法和开关频率优化PWM法(SFOPWM)的优点,在不提高开关频率的情况下... 提出了适用于大容量链式静止无功发生器(STATCOM)的一种新的脉冲控制策略,即三角载波移相-开关频率最优PWM法(PS-SFO-PWM)。这种新的PWM方法兼具三角载波移相PWM法(PSPWM)法和开关频率优化PWM法(SFOPWM)的优点,在不提高开关频率的情况下,实现等效高频载波的效果,并提高了直流电压的利用率。研究了链式多电平逆变器的PS-SFO-PWM技术,给出了PS-SFO-PWM控制方法的基本算法,并以2个H桥串联的5电平链式多电平变换器为例,通过Matlab仿真,验证了方法的正确性。 展开更多
关键词 链式多电平逆变器 sTATCOM Ps-sFO-PWM法 牛顿迭代法
在线阅读 下载PDF
上一页 1 2 18 下一页 到第
使用帮助 返回顶部