期刊文献+
共找到61篇文章
< 1 2 4 >
每页显示 20 50 100
狄拉克矩阵在超复数系下的关联性
1
作者 张玉贤 冯晓丽 +1 位作者 黄志祥 冯乃星 《安徽大学学报(自然科学版)》 北大核心 2025年第1期8-18,共11页
在超复数系下,重点研究狄拉克矩阵的导出以及其固有的关联性,通过计算机来进行数值验证.先从超复数中的四元数定义出发,结合实虚矩阵的表达方式,给出一种获取狄拉克矩阵的广义表示方法.基于对四元数情况下的狄拉克矩阵进行重点研究分析... 在超复数系下,重点研究狄拉克矩阵的导出以及其固有的关联性,通过计算机来进行数值验证.先从超复数中的四元数定义出发,结合实虚矩阵的表达方式,给出一种获取狄拉克矩阵的广义表示方法.基于对四元数情况下的狄拉克矩阵进行重点研究分析,给出一组实矩阵与超复数系的对应关系,探索性地对超复数系中的其他情形给予定性研究.通过二元数、八元数以及十六元数等情形的计算实例,进一步验证狄拉克矩阵在超复数系下关联性的可行性和准确性. 展开更多
关键词 超复数 狄拉克矩阵 广义表示方法 计算实例
在线阅读 下载PDF
内充气吹式玉米排种器工作性能EDEM-CFD模拟与试验 被引量:29
2
作者 韩丹丹 张东兴 +4 位作者 杨丽 李克鸿 张天亮 王云霞 崔涛 《农业工程学报》 EI CAS CSCD 北大核心 2017年第13期23-31,共9页
针对内充气吹式排种器组合气嘴低位安装对供气系统气流稳定性要求高、田间作业适应性差的问题,为优化内充气吹式排种器结构、改善工作性能,运用EDEM-CFD耦合分析方法,采用玉米籽粒粘结颗粒模型,以气嘴安装位置、工作压强和前进速度为试... 针对内充气吹式排种器组合气嘴低位安装对供气系统气流稳定性要求高、田间作业适应性差的问题,为优化内充气吹式排种器结构、改善工作性能,运用EDEM-CFD耦合分析方法,采用玉米籽粒粘结颗粒模型,以气嘴安装位置、工作压强和前进速度为试验因素进行正交试验优化仿真,提取关键参数充种极限速度和充种时长为评价指标进行分析,结果表明气嘴安装位置、工作压强和前进速度对充种极限速度影响显著而对充种时长影响不显著,充种极限速度和充种时长在0.05水平上呈显著负相关关系,较大充种极限速度和较短充种时长有利于提高排种器充种性能;对不同作业状态下的合格率、重播率、漏播率等作业指标进行分析表明气嘴安装位置在中位时合格率和充种极限速度呈显著正相关,低位和高位时因清种效果不佳导致合格率和充种极限速度呈负相关,中位时排种效果较佳。为验证仿真结果,进行气嘴安装位置和工作压强的全因素试验,结果表明,气嘴安装位置为中位时作业效果最佳,合格率大于90%,明显优于低位、高位状态,与仿真结果一致;在确定组合气嘴中位安装结构参数下,进行工作压强、前进速度的全因素试验,结果表明,前进速度为5~10 km/h、工作压强为4.5~5.5 k Pa时,排种器合格率较高,均为90%以上,该区间各工作参数下对合格率影响不显著。 展开更多
关键词 农业机械 计算机仿真 优化 精量排种器 充种性能 离散单元法
在线阅读 下载PDF
考虑颗粒碰撞的多重Monte Carlo算法 被引量:7
3
作者 赵海波 郑楚光 陈胤密 《力学学报》 EI CSCD 北大核心 2005年第5期564-572,共9页
从减少计算代价和改进碰撞算法出发,提出了考虑颗粒碰撞的多重MonteCarlo算法,它采用直接模拟MonteCarlo算法来考虑颗粒碰撞,并与求解颗粒拉氏Langevin方程的MonteCarlo算法耦合起来,跟踪比实际颗粒数目小得多的虚拟颗粒.提出了时间步... 从减少计算代价和改进碰撞算法出发,提出了考虑颗粒碰撞的多重MonteCarlo算法,它采用直接模拟MonteCarlo算法来考虑颗粒碰撞,并与求解颗粒拉氏Langevin方程的MonteCarlo算法耦合起来,跟踪比实际颗粒数目小得多的虚拟颗粒.提出了时间步长选定标准、虚拟碰撞伙伴所在控制容积的判断准则、颗粒碰撞发生的判断准则、虚拟碰撞伙伴的选择、基于随机碰撞角度的碰撞动力学,构成了考虑颗粒碰撞的完整多重MonteCarlo算法.对理想工况的细微颗粒流和粗重颗粒流进行了数值模拟,颗粒碰撞率的模拟结果与理论分析解和DNS结果均符合很好,颗粒场演变的细节信息,如时间平均和特定时刻的颗粒数密度,速度和颗粒湍动能等,均与DNS结果符合很好.数值模拟结果证明该算法不仅具有较低的计算代价,而且能够达到足够的计算精度. 展开更多
关键词 颗粒 碰撞 多重Monte Carlo算法 计算代价 计算精度
在线阅读 下载PDF
复杂非均匀介质伪谱法波场数值模拟 被引量:10
4
作者 程冰洁 李小凡 徐天吉 《石油物探》 EI CSCD 2007年第1期16-19,共4页
在地震波场数值正演模拟方法的研究中,计算精度和计算效率是评价方法的有效性及优越性的2个关键问题。从—阶速度—应力弹性波动方程出发,利用伪谱法求解波动方程,对复杂非均匀介质模型中的波场进行了正演模拟,并利用经典的Mamousi速度... 在地震波场数值正演模拟方法的研究中,计算精度和计算效率是评价方法的有效性及优越性的2个关键问题。从—阶速度—应力弹性波动方程出发,利用伪谱法求解波动方程,对复杂非均匀介质模型中的波场进行了正演模拟,并利用经典的Mamousi速度模型验证了该方法所具有的优势及存在的问题。将伪谱法模拟结果与交错网格高阶有限差分法的模拟结果比较可知:对于较为简单的非均匀模型,伪谱法和交错网格高阶有限差分法生成了几乎相同的波场;而当模型非常复杂且存在变化较剧烈的速度间断面时,伪谱法的模拟结果比较差。尽管如此,伪谱法计算速度快,计算效率高,能够直观、高效地反映介质中波场的传播规律。因而仍不失为一种很好的地震波模拟方法。 展开更多
关键词 非均匀介质 伪谱法 数值模拟 波动方程 计算精度 计算效率
在线阅读 下载PDF
船体零部件补偿量的计算方法研究 被引量:4
5
作者 刘玉君 胡日强 +1 位作者 田丰增 李瑞 《哈尔滨工程大学学报》 EI CAS CSCD 北大核心 2007年第7期723-726,共4页
船体建造精度控制是当代造船生产的关键技术之一,也是现代化造船科学管理的重要组成部分,其目的是最大限度地减少建造过程中的现场重复和修正作业量,缩短建造周期,降低造船成本.从船体建造精度控制的技术层面上看,如何能合理准确地给出... 船体建造精度控制是当代造船生产的关键技术之一,也是现代化造船科学管理的重要组成部分,其目的是最大限度地减少建造过程中的现场重复和修正作业量,缩短建造周期,降低造船成本.从船体建造精度控制的技术层面上看,如何能合理准确地给出船体零部件的补偿量,是船体建造精度控制方法需要研究的重要课题之一.船体零部件加工尺寸补偿就是在各零部件的基本尺寸上增加一个确定的量,以代替其工艺余量.文中开发了船体零部件加工测量数据库,应用统计质量控制理论建立了船体零部件加工补偿量数据统计的计算方法,应用尺寸链理论研究了船体零部件公差补偿量的计算方法. 展开更多
关键词 造船工艺 精度控制 补偿量 计算方法
在线阅读 下载PDF
4种计算自然常数e的方法及精度比较 被引量:4
6
作者 李鑫 王璐 +2 位作者 林金花 韩冬 谷德山 《东北师大学报(自然科学版)》 CAS CSCD 北大核心 2010年第4期57-61,共5页
用Visual C++6.0语言设计了3种计算方法:利用极限、级数、积分计算自然常数e的值.经过不断改进,让计算机按照人的思维方式对e的值进行了计算.计算结果表明:利用极限计算的e值能精确到小数点后8位;利用级数和定积分计算的e值精确到小数点... 用Visual C++6.0语言设计了3种计算方法:利用极限、级数、积分计算自然常数e的值.经过不断改进,让计算机按照人的思维方式对e的值进行了计算.计算结果表明:利用极限计算的e值能精确到小数点后8位;利用级数和定积分计算的e值精确到小数点后16位;利用人的思维方式计算e的值能精确到小数点后800位. 展开更多
关键词 自然常数e VISUAL C++6.0 精度 人的计算方式
在线阅读 下载PDF
基于单元的子区间摄动有限元方法研究 被引量:5
7
作者 苏静波 邵国建 刘宁 《计算力学学报》 EI CAS CSCD 北大核心 2007年第4期524-528,共5页
针对区间相关性导致区间扩张的问题,探讨了区间数之间的相关性并给出了降低区间扩张的子区间摄动方法。文中给出了基于单元的子区间摄动有限元计算公式和子区间划分数目的近似计算公式,同时文中讨论了区间有限元计算精度问题,给出了可... 针对区间相关性导致区间扩张的问题,探讨了区间数之间的相关性并给出了降低区间扩张的子区间摄动方法。文中给出了基于单元的子区间摄动有限元计算公式和子区间划分数目的近似计算公式,同时文中讨论了区间有限元计算精度问题,给出了可提高计算效率的一些措施。对桁架结构和平面应力问题梁结构算例分析结果表明文中方法可以达到一定的计算精度,并且是合理可行的。 展开更多
关键词 区间相关性 区间扩张 计算精度 区间有限元 摄动方法
在线阅读 下载PDF
阶梯矩阵及其一般化在迭代法中的应用 被引量:2
8
作者 邵新慧 沈海龙 李长军 《应用数学和力学》 CSCD 北大核心 2006年第8期971-977,共7页
Lu Hao首先给出了阶梯矩阵及其一般性的定义和性质.这类矩阵为迭代法提供了新矩阵分裂的基础.基于此新矩阵类的迭代方法的显著特征是它对于并行计算很容易被实现.应用这一新的分解方法,给出了一般的加速松弛方法(GAOR),而关于AOR方法的... Lu Hao首先给出了阶梯矩阵及其一般性的定义和性质.这类矩阵为迭代法提供了新矩阵分裂的基础.基于此新矩阵类的迭代方法的显著特征是它对于并行计算很容易被实现.应用这一新的分解方法,给出了一般的加速松弛方法(GAOR),而关于AOR方法的一些性质可以被延伸到该新方法中,并针对Hermite正定矩阵进行了新方法收敛性的分析.最后,给出了一些例子来表明新方法的优越性. 展开更多
关键词 阶梯矩阵 迭代法 平行计算 一般加速松弛方法(GAOR)
在线阅读 下载PDF
变电站空载母线波过程的精细积分计算方法 被引量:4
9
作者 赵进全 马西奎 邱关源 《电力系统自动化》 EI CSCD 北大核心 2002年第3期52-55,共4页
精细积分计算方法能够对具有任意负载的有损耦合输电线的波过程进行有效的分析。文中将该方法扩展到具有分支结构的输电线路的波过程分析之中 ,并通过对变电站空载母线加电的波过程进行数值分析 。
关键词 变电站 空载母线 波过程 精细积分计算方法 输电线路
在线阅读 下载PDF
舍入误差对大气环流模式模拟结果的影响 被引量:12
10
作者 王鹏飞 王在志 黄刚 《大气科学》 CSCD 北大核心 2007年第5期815-825,共11页
 此文旨在研究气候数值模式的长期计算时受舍入误差的影响。通过对大气环流谱模式SAMIL采用不同CPU数计算时获得的长时间积分结果进行分析,发现使用不同CPU数进行单精度计算时,其十年平均月平均500hPa高度场随机误差在正负6~8gpm范...  此文旨在研究气候数值模式的长期计算时受舍入误差的影响。通过对大气环流谱模式SAMIL采用不同CPU数计算时获得的长时间积分结果进行分析,发现使用不同CPU数进行单精度计算时,其十年平均月平均500hPa高度场随机误差在正负6~8gpm范围内,而使用双精度计算时相应的误差为正负3~4gpm。对于气候平均场而言,作者的试验表明SAMIL在并行计算时由于计算顺序改变而引起的误差在可接受范围之内。然而,虽然舍入误差的全球平均值不大,但其误差分布的差别范围并不小。数值试验得到的不同模拟结果之间误差大小与模拟结果的自身年际变化大小在同样的量级,因此对于“年际变化”这样的问题来说,其影响是不可忽略的,必须要使用集合预报的办法来减小误差的影响。文中列出了3种研究复杂数值模式舍入误差的实验方法,指出其一定条件下的等效性和不同适用范围,对于其他模式的舍入误差影响研究有一定的参考价值。在舍入误差分析的基础上,介绍了一种新型的专门针对舍入误差的集合预报方法(舍入误差平均集合,RME),指出了其在气候模拟研究中的应用价值。 展开更多
关键词 大气环流模式 舍入误差 计算精度 误差研究方法 集合预报
在线阅读 下载PDF
刚塑性无网格方法中精度影响因素的研究 被引量:3
11
作者 赵国群 吴欣 管延锦 《塑性工程学报》 EI CAS CSCD 北大核心 2006年第3期13-18,共6页
将无网格方法引入金属塑性成形过程分析,基于刚塑性材料假设,采用移动最小二乘法对速度场进行逼近,利用变换法施加本质边界条件,实现了金属塑性成形过程的无网格方法模拟。在此基础上,结合典型锻造过程分析实例,重点研究了应用无网格方... 将无网格方法引入金属塑性成形过程分析,基于刚塑性材料假设,采用移动最小二乘法对速度场进行逼近,利用变换法施加本质边界条件,实现了金属塑性成形过程的无网格方法模拟。在此基础上,结合典型锻造过程分析实例,重点研究了应用无网格方法分析金属塑性成形问题时,诸如:基函数、正交化与否、权函数、权函数影响域大小、积分方案、高斯积分阶次、节点分布密度以及针对体积闭锁问题而提出的体积应变率映射方法等对于模拟精度和效率的影响规律,给出了提高无网格模拟精度和效率的影响因素的选择方法。 展开更多
关键词 刚塑性 无网格方法 金属成形 精度
在线阅读 下载PDF
非均匀传输线时域响应的精细计算 被引量:5
12
作者 赵进全 马西奎 邱关源 《西安交通大学学报》 EI CAS CSCD 北大核心 1998年第8期5-7,共3页
应用精细计算法对非均匀传输线时域响应进行了分析.该方法考虑了非均匀传输线的非均匀性,避免了使用突变的级联模型及频域、时域转换,直接时域分析,极大地提高了分析的精度和效率.计算结果证实了该方法的可靠性.
关键词 非均匀传输线 时域响应 精细计算法
在线阅读 下载PDF
一种求解地震波方程的高效并行谱元格式 被引量:5
13
作者 林灯 崔涛 +1 位作者 冷伟 张林波 《计算机研究与发展》 EI CSCD 北大核心 2016年第5期1147-1155,共9页
地震波数值模拟在地震学和地震勘探中扮演着非常重要角色.在已有工作的基础上,提出1种高效并行的地震波PML方程谱元格式.PML被引入地震波方程以吸收外向波进而模拟无界区域.进一步,为了适应复杂地形同时允许时间显式推进,谱元方法被用... 地震波数值模拟在地震学和地震勘探中扮演着非常重要角色.在已有工作的基础上,提出1种高效并行的地震波PML方程谱元格式.PML被引入地震波方程以吸收外向波进而模拟无界区域.进一步,为了适应复杂地形同时允许时间显式推进,谱元方法被用来离散地震波PML方程.由此得到地震波PML方程谱元格式.在此基础上,阐述了单元刚度矩阵分解性质,并说明了利用单元刚度矩阵分解可以大幅减少刚度矩阵存储量同时显著加速刚度矩阵与向量乘积,进而显著减少格式的计算量和存储量.此外,算法复杂性分析表明格式无论在计算量上还是在存储量上都优于几种已知的1阶地震波PML方程谱元格式.结合并行技术,给出了高效并行的地震波PML方程谱元格式.数值实验验证了格式的正确性、良好的强弱并行可扩展性以及对复杂地形的适应性. 展开更多
关键词 地震波方程 数值模拟 完美匹配层 谱元方法 单元刚度矩阵分解 并行计算
在线阅读 下载PDF
精细积分法应用于地震碰撞力反应谱计算研究 被引量:3
14
作者 张瑞杰 李青宁 尹俊红 《计算力学学报》 CAS CSCD 北大核心 2015年第6期733-738,共6页
相邻结构的碰撞时程分析是计算地震碰撞力反应谱的基础。结构碰撞时程分析要求采用稳定性好、精度高及计算效率高的数值分析方法。精细积分法将二阶动力微分方程通过增元降阶的方式转换成Hamilton对偶变量体系,得到了动力微分方程的精... 相邻结构的碰撞时程分析是计算地震碰撞力反应谱的基础。结构碰撞时程分析要求采用稳定性好、精度高及计算效率高的数值分析方法。精细积分法将二阶动力微分方程通过增元降阶的方式转换成Hamilton对偶变量体系,得到了动力微分方程的精确解。基于此,本文将精细积分法引入结构碰撞时程分析及地震碰撞力反应谱计算中。在推导精细积分法公式的基础上,在MATLAB环境下编制了结构碰撞时程分析程序和碰撞力反应谱计算程序,并实现了碰撞力反应谱程序的并行化。经算例验证,精细积分法应用于结构碰撞时程分析及地震碰撞力反应谱计算是可行的,程序计算结果准确。 展开更多
关键词 地震 结构碰撞 精细积分法 反应谱 并行计算
在线阅读 下载PDF
结构动力响应的精细时程积分并行算法 被引量:6
15
作者 李红云 沈为平 《上海交通大学学报》 EI CAS CSCD 北大核心 1998年第8期81-85,共5页
大型结构动力响应问题采用传统的时间步长直接积分方法求解是非常耗时的,因此发展相应的并行计算方法成为必然.传统的直接积分方法是极度串行而不适于并行计算,而精细时程积分方法可以使用大步长单步计算出求解区间任意时间点的响应... 大型结构动力响应问题采用传统的时间步长直接积分方法求解是非常耗时的,因此发展相应的并行计算方法成为必然.传统的直接积分方法是极度串行而不适于并行计算,而精细时程积分方法可以使用大步长单步计算出求解区间任意时间点的响应值,为并行计算提供了极大的可能性.结构动力响应方程通过变量变换可以转化为一阶线性常微分方程,该方程组的解由表示初值影响的齐次方程解和反映荷载作用的积分之和组成.其中第一项用矩阵指数函数计算;第二项在文中采用精细时程积分傅里叶展开方法计算(设计了3种相应的并行算法),并在TRANSPUTER并行机上实现.结果表明,3种HPD-F并行算法有很高的加速比和并行效率. 展开更多
关键词 动力响应力 精细时程积分 并行算法 结构
在线阅读 下载PDF
计算机控制光学表面成形中大规模驻留时间求解 被引量:5
16
作者 罗丽丽 何建国 +3 位作者 王亚军 张云飞 黄文 吉方 《强激光与粒子束》 EI CAS CSCD 北大核心 2011年第12期3207-3212,共6页
采用基于稀疏矩阵的大规模非负最小二乘法,对大口径、微浮雕结构光学元件加工中的驻留时间进行了分析与求解,并对该算法开展了正则化研究。仿真结果表明:与传统非负最小二乘法相比,基于稀疏矩阵的大规模非负最小二乘法精度高、效率快。... 采用基于稀疏矩阵的大规模非负最小二乘法,对大口径、微浮雕结构光学元件加工中的驻留时间进行了分析与求解,并对该算法开展了正则化研究。仿真结果表明:与传统非负最小二乘法相比,基于稀疏矩阵的大规模非负最小二乘法精度高、效率快。采用该算法仿真加工平均振幅为1.177 6倍波长的大口径、微浮雕结构光学元件,误差面形均方根收敛至0.067倍波长。 展开更多
关键词 计算机控制光学表面成形 驻留时间 大规模非负最小二乘法 正则化 面形精度
在线阅读 下载PDF
农业机器人轨迹优化自动控制研究——基于BP神经网络与计算力矩 被引量:13
17
作者 袁铸 申一歌 《农机化研究》 北大核心 2017年第6期33-37,共5页
以农业机器人精密轨迹优化自动控制为目标,在优化算法中引入BP神经网络与计算力矩法结合的自动控制器,旨在减少作业过程中的运动误差,提高其工作效率。首先,建立农业机器人数学模型,分析其运动学和动力学原理;然后,设计了农业机器人运... 以农业机器人精密轨迹优化自动控制为目标,在优化算法中引入BP神经网络与计算力矩法结合的自动控制器,旨在减少作业过程中的运动误差,提高其工作效率。首先,建立农业机器人数学模型,分析其运动学和动力学原理;然后,设计了农业机器人运动控制系统,引入BP神经网络对不确定动力学因素进行判断,并提出解决该因素的自适应学习法;最后,对该系统运用Mat Lab进行了仿真。试验表明:以BP神经网络与计算力矩法结合的自动控制器可以有效优化机器人运动路径,提高机器人整体作业效率,系统运行稳定、可靠性强,且对外部环境的干扰因素具有较强的自适应学习能力。 展开更多
关键词 农业机器人 精密轨迹优化 BP神经网络 计算力矩法
在线阅读 下载PDF
利用积分算子方程的小波-矩量法分析频率选择反射面 被引量:1
18
作者 宋玉明 方大纲 《南京理工大学学报》 CAS CSCD 1996年第6期537-540,共4页
小波-矩量法以小波基作为矩量法中的基函数,通过Galerkin法可以获得很稀疏的阻抗矩阵,用迭代法求解时可减少计算机内存和计算时间。该文以有限尺寸的频率选择反射面为例,用小波-矩量法分析散射特性和电流分布,着重讨论了... 小波-矩量法以小波基作为矩量法中的基函数,通过Galerkin法可以获得很稀疏的阻抗矩阵,用迭代法求解时可减少计算机内存和计算时间。该文以有限尺寸的频率选择反射面为例,用小波-矩量法分析散射特性和电流分布,着重讨论了阻抗矩阵的稀疏化程度对计算精度的影响。结果表明,当非零元素仅为10%时,仍能达到满意的精度。 展开更多
关键词 稀疏矩阵 电磁场 积分算子方程 小波-矩量法
在线阅读 下载PDF
线弹性结构的高阶摄动随机有限元法 被引量:3
19
作者 廖光明 陈兵 王遂泸 《四川大学学报(工程科学版)》 EI CSCD 北大核心 2014年第3期44-48,共5页
计算精度和误差控制是摄动随机有限元法在数值计算中的重要问题。根据随机理论和摄动有限元方法,推导了随机变量均值和n阶中心矩的计算公式;提出了摄动随机有限元实施过程中精度的控制方法,能够在计算结构响应均值和方差之前就分别确定... 计算精度和误差控制是摄动随机有限元法在数值计算中的重要问题。根据随机理论和摄动有限元方法,推导了随机变量均值和n阶中心矩的计算公式;提出了摄动随机有限元实施过程中精度的控制方法,能够在计算结构响应均值和方差之前就分别确定计算结果的精度;研究了摄动随机有限元0~N阶递推方程组计算列式;运用高阶摄动随机有限元方法对受轴向拉力的1维线弹性杆进行分析,结果表明,与2阶摄动法相比,N阶摄动技术能有效改进计算精度。 展开更多
关键词 随机有限元 摄动技术 随机参数 计算精度
在线阅读 下载PDF
航空摄影测量技术的数学方法探讨 被引量:3
20
作者 李海滨 申家双 王耿峰 《海洋测绘》 CSCD 2018年第2期21-24,共4页
航空摄影测量技术经历了模拟、解析、数字摄影测量三个阶段,随着计算机技术及数学算法更加完善,全自动数字摄影技术即将到来,以数字航空摄影测量处理流程为基础,梳理了涉及的数学理论及算法,并从学科融合的角度将控制学、云计算等理论... 航空摄影测量技术经历了模拟、解析、数字摄影测量三个阶段,随着计算机技术及数学算法更加完善,全自动数字摄影技术即将到来,以数字航空摄影测量处理流程为基础,梳理了涉及的数学理论及算法,并从学科融合的角度将控制学、云计算等理论引入算法流程中,以提高航空摄影测量处理自动化程度,为拓展摄影测量领域应用,提高自动化处理能力提供理论依据。 展开更多
关键词 数字摄影测量 傅里叶变换 数值计算方法 云计算 矩阵论 概率论 多视图几何 控制系统理论
在线阅读 下载PDF
上一页 1 2 4 下一页 到第
使用帮助 返回顶部