第二章测量误差分布优秀PPT.ppt

上传人:石*** 文档编号:78959555 上传时间:2023-03-19 格式:PPT 页数:61 大小:4.88MB
返回 下载 相关 举报
第二章测量误差分布优秀PPT.ppt_第1页
第1页 / 共61页
第二章测量误差分布优秀PPT.ppt_第2页
第2页 / 共61页
点击查看更多>>
资源描述

《第二章测量误差分布优秀PPT.ppt》由会员分享,可在线阅读,更多相关《第二章测量误差分布优秀PPT.ppt(61页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、第二章测量误差分布第一页,本课件共有61页主要内容 熟悉误差分布的基本概念、常见误差分布特征与处理熟悉误差分布的基本概念、常见误差分布特征与处理方法方法直方图的绘制直方图的绘制概率密度分布图概率密度分布图误差分布的特征值误差分布的特征值常见的误差分布常见的误差分布常用的统计量分布常用的统计量分布误差分布的统计检验误差分布的统计检验 第二页,本课件共有61页主主主主 要要要要 内内内内 容容容容2.1 测量误差的统计特性 2.2 常见测量误差分布2.3常见的统计量分布2.4误差分布的分析与检验1234第三页,本课件共有61页2.1测量误差的统计特性测量误差的统计特性一、测量点列图一、测量点列图某

2、钢球工件直径重复测量某钢球工件直径重复测量150次,次,得到一个测量样本得到一个测量样本7.0857.3357.585单 峰 性:数 据 集 中 在7.335附近有界性:数据分布在7.085至7.585之间对称性:正负误差的数目大致相同;抵偿性:误差的总和大致趋于零第四页,本课件共有61页(1 1)分分 组组 数数=11=11,组组 距距=0.05mm=0.05mm;(2 2)依依次次定定各各组组的的频频数数、频率和频率密度;频率和频率密度;(3 3)以以数数据据为为横横坐坐标标,频频率率密密度度为为纵纵坐坐标标,在在横横坐坐标标上上划划出出等等分分的的子子区区间间,划划出出各各子子区区间间的

3、的直直方方柱柱,即即为为所所求求统统计计直直方图。方图。77.17.27.37.47.57.60510152025二、统计直方图二、统计直方图第五页,本课件共有61页绘制统计直方图注意事项绘制统计直方图注意事项(1 1)样本大小:)样本大小:确定误差的分布范围时,取确定误差的分布范围时,取 n=50 n=50200200 确定误差分布规律时,最好取确定误差分布规律时,最好取n=200n=20010001000()子区间个数、间距:子区间个数、间距:当当n=50n=50100100时,时,个数个数=6=61010当当n=100n=100200200时,个数时,个数=9=91212当当n=200n

4、=200500500时,个数时,个数=12=121717当当n=500n=500以上时,个数以上时,个数=20=20可用下列两个公式之一来计算分组数 或间距 或第六页,本课件共有61页 把把各各直直方方柱柱顶顶部部中中点点用用直直线线连连接接起起来来,便便得得到到一一条条由由许许多多折折线线连连接接起起来来的的曲曲线线。当当测测量量样样本本数数n n无无限限增增加加,分分组组间间隔隔趋趋于于零零,图图中中直直方方图图折折线线变变成成一一条条光光滑滑的的曲曲线线,即即测测量量总总体体的的概概率率(分分布布)密密 度度 曲曲 线线,记记 为为 。这这就就是是用用实实验验方方法法由由样样本本得得到到

5、的的概概率率密密度度分分布曲线。布曲线。77.17.27.37.47.57.60510152025三、概率密度(分布)图三、概率密度(分布)图(测量总体)第七页,本课件共有61页概率密度曲线完好的描述了随机误差的统计规律。概率密度曲线完好的描述了随机误差的统计规律。概率密度函数的几何意义 置信区间置信区间 显著性水平(又称显著度或危险率)显著性水平(又称显著度或危险率)置信概率(或置信水置信概率(或置信水平),简记为符号平),简记为符号概率密度的性质概率密度的性质有两个性质有两个性质第八页,本课件共有61页四、统计分布特征值四、统计分布特征值尽管误差分布反映了该误差的全貌,但在实际使用中更尽管

6、误差分布反映了该误差的全貌,但在实际使用中更关心代表该误差分布的若干数字特征量。关心代表该误差分布的若干数字特征量。p数学期望数学期望p标准偏差标准偏差p偏态系数偏态系数p峰态系数峰态系数p协方差协方差p相关系数相关系数第九页,本课件共有61页p 数学期望(加权平均)数学期望(加权平均)定义定义一阶原点矩,它表示随机变一阶原点矩,它表示随机变量分布的位置特征。它与真值量分布的位置特征。它与真值之差即为系统误差,如果系统之差即为系统误差,如果系统误差可以忽略,则误差可以忽略,则 就是被测就是被测量的真值量的真值 三条测量值分布曲线的精密度相同,三条测量值分布曲线的精密度相同,但正确度不同。但正确

7、度不同。数学期望代表了测量的最佳数学期望代表了测量的最佳估计值,或相对真值的估计值,或相对真值的系统误系统误差大小差大小第十页,本课件共有61页p 标准偏差标准偏差二阶中心矩,称为二阶中心矩,称为X X的标准(偏)的标准(偏)差,差,的大小表征了随机,的大小表征了随机误差的误差的分散程度分散程度,即大部分分布在,即大部分分布在 范围内,可作为随机误差的评定范围内,可作为随机误差的评定尺度尺度 定义定义三条误差分布曲线的正确度相同,三条误差分布曲线的正确度相同,但精密度不同但精密度不同 标准差代表了该测量条件下的标准差代表了该测量条件下的测量结果分散性的大小,或是测量结果分散性的大小,或是该测量

8、分布的该测量分布的随机误差随机误差大小大小 第十一页,本课件共有61页p 偏态系数偏态系数定义定义三阶中心矩,三阶中心矩,将将 无量纲无量纲化,化,称为偏态系数,称为偏态系数,描述描述了测量总体及其误差分布的了测量总体及其误差分布的非对称程度非对称程度 曲线曲线具有正具有正(右右)偏态,曲偏态,曲线线具有负具有负(左左)偏态偏态第十二页,本课件共有61页p 峰态系数峰态系数定义定义 表征了测量总体表征了测量总体及其误差分布的峰凸程度。及其误差分布的峰凸程度。是将是将 无量纲化,也称峰度,无量纲化,也称峰度,而而 是按标准正态分布归零,是按标准正态分布归零,即对于正态分布超越系数即对于正态分布超

9、越系数 视为零视为零 较尖峭的分布有较尖峭的分布有 ,较平,较平坦的分布有坦的分布有 第十三页,本课件共有61页p 协方差协方差定义定义式中式中协方差协方差 表示了两变量间的相关程度表示了两变量间的相关程度 第十四页,本课件共有61页p 相关系数相关系数定义定义表示了两个变量间线性相关表示了两个变量间线性相关的程度的程度 越小,越小,X X,Y Y之间线性相之间线性相关程度越小,关程度越小,取值越大,取值越大,X X,Y Y之间线性相关程度越大之间线性相关程度越大 当当 ,X X与与Y Y正相关,当正相关,当 ,X X与与Y Y负相关负相关 线性相关正相关负相关线性不相关第十五页,本课件共有6

10、1页数学期望名称定义方差几何意义误差意义偏态系数峰态系数协方差位置特征实际值正确度弥散分散性,精密度不对称误差分布不对称性尖峭误差分布尖峭程度两误差关联程度统计分布常用的特征值统计分布常用的特征值第十六页,本课件共有61页2.2常见测量误差分布常见测量误差分布p正态分布正态分布p标准偏差标准偏差p均匀分布均匀分布p三角分布三角分布p瑞利分布瑞利分布p反正弦分布反正弦分布p分布分布几种常见的误差分布几种常见的误差分布第十七页,本课件共有61页一、正态分布一、正态分布服从正态分布的条件服从正态分布的条件 误误差差因因素素多多而而小小,无无一一个个占占优优,彼彼此此相相互互独独立立(中心极限定理)。

11、(中心极限定理)。一般认为,当影响测量的因素在一般认为,当影响测量的因素在1515个以上,且相互个以上,且相互独立,其影响程度相当,可以认为测量值服从正态分独立,其影响程度相当,可以认为测量值服从正态分布;若要求不高,影响因素则应在布;若要求不高,影响因素则应在5 5个(至少个(至少3 3个)以个)以上,也可视为正态分布。上,也可视为正态分布。第十八页,本课件共有61页概率密度函数概率密度函数正态分布的密度函数:正态分布的密度函数:为测量总体的数学期望,如不计系统误为测量总体的数学期望,如不计系统误差,则差,则 即为随机误差即为随机误差 为测量总体的标准差,也是为测量总体的标准差,也是 随机误

12、随机误差的标准差差的标准差 第十九页,本课件共有61页(1 1)单单峰峰性性:小小误误差差出出现现的的概概率率比大误差出现的概率大。比大误差出现的概率大。(2 2)对称性:正误差出现的概)对称性:正误差出现的概率与负误差出现的概率相等。率与负误差出现的概率相等。(3 3)抵偿性:随测量次数增)抵偿性:随测量次数增加,算术平均值趋于零。加,算术平均值趋于零。分布的误差特性分布的误差特性正态分布的这三个特点与误差大样本下的统计特性相符。但在理论上,正态分布无界,这也是正态分布与实际误差有界性不相符之处。第二十页,本课件共有61页正态分布的置信概率 误差在分布区间 的置信概率 式中68.26%95.

13、45%99.73%置信概率置信概率 正态积分函数,已制成正态积分表 置信因子第二十一页,本课件共有61页正态分布的某些正态分布的某些k值的置信概率值的置信概率3.33.02.582.01.961.6451.00.67450.9990.99730.990.9540.950.900.6830.50.0010.00270.010.0460.050.100.3170.5第二十二页,本课件共有61页 (1)(1)经典误差理论都是建立在正态分布的经典误差理论都是建立在正态分布的基础基础上。凡上。凡是有是有3 3、5 5个以上的、差不多微小的、独立影响的个以上的、差不多微小的、独立影响的合成分布都趋近正态分

14、布。这是被前人早已证明合成分布都趋近正态分布。这是被前人早已证明了的中心极限定理告诉我们的一个事实。了的中心极限定理告诉我们的一个事实。正态分布在误差理论和实践中的地位正态分布在误差理论和实践中的地位(2)(2)许多非正态分布可以用正态分布来表示。许多非正态分布可以用正态分布来表示。(3)(3)正态分布的概率密度函数具有简单的数学形式正态分布的概率密度函数具有简单的数学形式和优良的性质。和优良的性质。当然,也有不少的误差分布并不能简单地用正态分布当然,也有不少的误差分布并不能简单地用正态分布来描述。因而,现代误差理论及其实践需要进一步研来描述。因而,现代误差理论及其实践需要进一步研究非正态分布

15、的问题。究非正态分布的问题。第二十三页,本课件共有61页二、均匀分布二、均匀分布若误差在某一范围中出现的概率相等,称其服从若误差在某一范围中出现的概率相等,称其服从均匀分布,也称为等概率分布。均匀分布,也称为等概率分布。概率密度函数概率密度函数 数学期望数学期望方差方差标准方差标准方差置信因子置信因子 o-a a 第二十四页,本课件共有61页服从均匀分布的可能情形服从均匀分布的可能情形(1)(1)数据截尾引起的舍入误差数据截尾引起的舍入误差;(2)(2)数字显示末位的截断误差数字显示末位的截断误差(3)(3)瞄准误差瞄准误差;(4)(4)数字仪器的量化误差数字仪器的量化误差;(5)(5)齿齿轮

16、轮回回程程所所产产生生的的误误差差以以及及基基线线尺尺滑滑轮轮摩摩擦擦引起的误差;引起的误差;(6)(6)多中心值不同的正态误差总和服从均匀分布。多中心值不同的正态误差总和服从均匀分布。第二十五页,本课件共有61页三、三角分布三、三角分布概率密度函数概率密度函数 数学期望数学期望标准方差标准方差当两个分布范围相等的均匀分布,其合成误差就当两个分布范围相等的均匀分布,其合成误差就是三角分布。是三角分布。第二十六页,本课件共有61页四、反正弦分布四、反正弦分布概率密度函数概率密度函数 数学期望数学期望标准方差标准方差a a-a-a o o 服从反正弦分布的可能情形服从反正弦分布的可能情形 度盘偏心

17、引起的测角误差;度盘偏心引起的测角误差;正弦(或余弦)振动引起的位移误差;正弦(或余弦)振动引起的位移误差;无线电中失配引起的误差。无线电中失配引起的误差。第二十七页,本课件共有61页五、瑞利分布五、瑞利分布概率密度函数概率密度函数 数学期望数学期望标准方差标准方差服从瑞利分布的可能情形服从瑞利分布的可能情形 偏心值偏心值在非负值的单向误差中,由于偏在非负值的单向误差中,由于偏心因素所引起的轴的径向跳动心因素所引起的轴的径向跳动 刻度盘、圆光栅盘的最大分度误差刻度盘、圆光栅盘的最大分度误差 齿轮和分度盘的最大齿距累积误差齿轮和分度盘的最大齿距累积误差 第二十八页,本课件共有61页六、贝塔分布六

18、、贝塔分布概率密度函数概率密度函数 数学期望数学期望标准方差标准方差第二十九页,本课件共有61页在给定分布界限下通在给定分布界限下通过参数取不同值,贝塔过参数取不同值,贝塔分布可呈对称分布、非对称分布可呈对称分布、非对称分布、单峰分布、递增或递分布、单峰分布、递增或递减分布等,可逼近常见的正减分布等,可逼近常见的正态、三角、均匀、反正弦、态、三角、均匀、反正弦、瑞利等各种典型分布。贝塔瑞利等各种典型分布。贝塔分布具有可逼近各种实际误分布具有可逼近各种实际误差分布的多态性。差分布的多态性。贝塔分布在理论上就是有界的。贝塔分布在理论上就是有界的。不像正态、瑞利等呈拖尾型分布,不像正态、瑞利等呈拖尾

19、型分布,完全符合误差的基本特性即有界性。完全符合误差的基本特性即有界性。贝塔分布的性质与密度函数图贝塔分布的性质与密度函数图第三十页,本课件共有61页常见分布的数字特征量常见分布的数字特征量名称正态分布区间半宽度标准差期望等价均匀分布三角分布反正弦分布瑞利分布第三十一页,本课件共有61页2.3常见的统计量分布常见的统计量分布本节介绍常用的统计量分布,包括本节介绍常用的统计量分布,包括t t分布分布 F F分分布,分布。布,分布。前边介绍主要是单个统计量分布,实际中要前边介绍主要是单个统计量分布,实际中要常用到变量间组合也就是函数的统计分布。常用到变量间组合也就是函数的统计分布。第三十二页,本课

20、件共有61页一、分布定义定义若为独立服从同分布 的随机误差,则称服从为自由度为的分布。概率密度函数概率密度函数 数学期望数学期望标准方差标准方差第三十三页,本课件共有61页二、t分布定义定义若随机误差,随机误差,且和相互独立,则服从的分布称为自由度为的t分布。概率密度函数概率密度函数 数学期望数学期望标准方差标准方差o 第三十四页,本课件共有61页 当自由度足够大时,t分布趋近于正态分布。t t分布在误差理论和实践中的应用分布在误差理论和实践中的应用t分布在研究正态小子样(测量次数较少时),是一个严密而有效的理论分布。正态样本的算术平均值构成的如下统计量服从自由度为的t分布。其测量算术平均值满

21、足 t分布的临界值,满足第三十五页,本课件共有61页三、F分布定义定义若,则称服从为自由度为的F分布。概率密度函数概率密度函数 数学期望数学期望标准方差标准方差第三十六页,本课件共有61页2.4误差分布的分析与检验误差分布的分析与检验本节介绍确定误差分布规律的几种方法,包括物本节介绍确定误差分布规律的几种方法,包括物理来源法,函数关系法以及图形判断法。最后介绍理来源法,函数关系法以及图形判断法。最后介绍有关分布检验的知识,有关分布检验的知识,包括正态分布统计检验(夏皮罗包括正态分布统计检验(夏皮罗-威尔克检验、威尔克检验、偏态系数和峰态系数检验)和一般分布检验(皮偏态系数和峰态系数检验)和一般

22、分布检验(皮尔逊检验)尔逊检验)。第三十七页,本课件共有61页物理来源判断法物理来源判断法 根据测量误差产生的来源,可以判断其根据测量误差产生的来源,可以判断其属于何种类型属于何种类型 如其测量受到至少有如其测量受到至少有三个以上三个以上独立的、微小独立的、微小而大小相近的因素的影响,则可认为它服从或而大小相近的因素的影响,则可认为它服从或接近正态分布。接近正态分布。测量值在某范围内各处出现的机会相等,测量值在某范围内各处出现的机会相等,则可认为它服从均匀分布。则可认为它服从均匀分布。一、误差分布的分析与判断一、误差分布的分析与判断第三十八页,本课件共有61页函数关系法函数关系法 利用随机变量

23、的函数关系,来判断误差属于何种分布。若与都在-a,a内服从均匀分布,则服从三角分布 若与都服从正态分布,则 服从偏心分布(瑞利分布)若服从均匀分布,则 服从反正弦分布第三十九页,本课件共有61页图形判断法图形判断法对重复测量获得的样本数据绘出频率密度直方图,并与各种常见的概率密度分布曲线相比较,判断它与何种分布相接近。第四十页,本课件共有61页统计检验的步骤统计检验的步骤 1 1、概念、概念事先对分布形式作出某种假设事先对分布形式作出某种假设然后利用样本信息来判断原假设是否成立然后利用样本信息来判断原假设是否成立2 2、类型、类型正态分布统计检验正态分布统计检验一般分布检验一般分布检验v夏皮罗

24、夏皮罗-威尔克检验威尔克检验v偏态系数检验偏态系数检验v峰态系数检验峰态系数检验v皮尔逊检验皮尔逊检验 二、误差分布的统计检验二、误差分布的统计检验第四十一页,本课件共有61页皮尔逊 检验(且已知 )1、提出原假设总体 的分布函数 未知 某个已知的分布函数 2、计算统计量总体中抽取出一个容量为 的样本 把整个数轴分成 个区间 频数,样本的观察值落在第 个区间的个数 由 计算出总体 在各区间内取值的概率 第四十二页,本课件共有61页检验(续)3、在给定显著性水平 下,由分布表查得临界值 。4、作出决策。若,拒绝,则认为 。反之,思路是当样本个数充分大时,频率和概率应当相思路是当样本个数充分大时,

25、频率和概率应当相差不会太大,如果超出某种限度,则假设就会推差不会太大,如果超出某种限度,则假设就会推翻。翻。第四十三页,本课件共有61页 皮尔逊检验(分布中含有未知参数)1、提出原假设总体 的分布函数 未知 某个已知形式的分布函数,未知参数 2、计算统计量总体中抽取出一个容量为 的样本 第四十四页,本课件共有61页在 下利用样本给出 的极大似然估计 把整个数轴分成 个区间 频数,样本的观察值落在第 个区间的个数 由 计算出总体 在各区间内取值的概率 3、在给定显著性水平 下,由分布表查得临界值 。4、作出决策。若 ,拒绝 皮尔逊检验(续)第四十五页,本课件共有61页【例【例2-12-1】用阿贝

26、比较仪测量某轴承直径 100次,依次测得 ,的数据见下所列,的单位0.1 。检验 是否服从正态分布。0 -5 11 -10 17 -3 -13 6 4 7 1 -5 -6 -3 13 -1 -1 5 9 7 -3 9 -8 3 -2 -24-30 -2 1 -2 4 2 -5 -13 1 -7 -1 0 -4 -7 0 7 17 5 10 0 -2 6 3 8 6 -3 -3-10 0 5 2 -8 0 4 2 2 6-11 5 2 7 -1 12 0 -19 10 -1 7 9 2 -5 14 -6 -5 8 3 8 -9 4 -5 -8 8 -8 4 -13 -9-10 -10 2 13

27、2 -4 6 -7 第四十六页,本课件共有61页计算步骤【解】【解】检验 由于 中含有未知参数,故需先进行参数估计。在正态分布下,和 的极大似然估计为 将 取值分成8组,然后计算概率 第四十七页,本课件共有61页计算结果频数70.10710.75-3.751.31150.16016.01-1.010.06130.13313.37-0.370.0890.0989.87-0.870.08100.0989.870.130160.13313.372.630.52210.16016.014.991.5690.10710.75-1.750.281003.82第四十八页,本课件共有61页结论 给定显著性水平

28、 ,自由度8-2-1=5,由 分布表查得临界值 因为所以,接受 ,故可认为这些测量服从正态分布 组数未知数个数第四十九页,本课件共有61页夏皮罗威尔克检验夏皮罗-威尔克检验又称W检验 时检验效果最佳,并且计算简便。只能用于正态性检验第五十页,本课件共有61页W检验的实施步骤从总体中抽取出一个容量为 的样本(1)将样本的观测值按由小到大排列成为其次序统计量(2)计算检验统计量(3)查表。由夏皮罗-威尔克值表查出 ,为给定的显著性水平;(4)判断。若 ,则拒绝正态性假设夏皮罗夏皮罗威尔克威尔克当当n为偶数时取为偶数时取n/2,当,当n奇数时取(奇数时取(n1)/2第五十一页,本课件共有61页【例【

29、例2-22-2】用夏皮罗-威尔克法检验该组数据是否来自正态分布。将某量独立测得结果按从小到大排列成(n=10)108,109,110,110,110,112,112,116,119,124 【解】【解】查夏皮罗-威尔克系数 表得出 第五十二页,本课件共有61页计算结果计算 给定显著性水平 ,查表得因为,故拒绝正态性假设 第五十三页,本课件共有61页偏态系数检验(1)给出备择假设 (正偏)或 (负偏)(2)计算检验统计量(3)查表。根据显著性水平 和样本容量 ,由偏态统计量的分位数表查出 (4)判断。当备择假设为 时,若 ,则拒绝正态性假设;当备择假设为 时,若 ,则拒绝正态性假设第五十四页,本

30、课件共有61页【例【例2-32-3】有下列一组测量数据,确定这批数据是否来自正态分布-0.40-1.80-2.14 0.40-1.40 0.67-1.40-1.51 1.40 -1.40-1.38-1.40 1.20-2.14-0.60-2.33 1.24-0.40 -0.32-0.22-1.60-1.40-0.51-0.20-1.40-1.72-1.60 -1.20-1.80 1.20-1.40-0.80-1.72-0.71-1.40-1.20 -1.91-0.69-1.60-1.39-2.20-1.40-0.40 0.40-1.80 -1.80 -1.60 0 -1.95 1.20 第五十

31、五页,本课件共有61页计算结果计算统计量 由得因此,选择备择假设 给定显著性水平 ,当n=50时,查表得因为,故拒绝正态性假设 第五十六页,本课件共有61页峰态系数检验(1)给出备择假设 (正偏)或 (负偏)(2)计算检验统计量(3)查表。根据显著性水平 和样本容量 ,由峰态统计量的分位数表查出 或 (4)判断。当备择假设为 时,若 ,则拒绝正态性假设;当备择假设为 时,若 ,则拒绝正态性假设第五十七页,本课件共有61页【例【例2-42-4】利用某测量仪器进行40次测量,测得与理论值的如下一系列偏差数据,确定这批数据是否来自正态分布 0.038 0.240 0.124 0.054 -0.061

32、 -0.004 -0.004 -0.006 0.007 0.001 0.061 0.043 0.035 0.163 -0.008 -0.010 0.006 -0.008 -0.024 0.007 0.028 0.108 0.155 -0.159 -0.032 0.003 -0.007 -0.018 -0.008 -0.011 0.060 0.067 -0.025 -0.096 -0.223 0.004 -0.007 -0.007 -0.010 0.014 第五十八页,本课件共有61页计算结果计算统计量 由得因此,选择备择假设 给定显著性水平 ,当n=40时,查表得因为,故拒绝正态性假设 第五十九页,本课件共有61页思考与练习题 第六十页,本课件共有61页第六十一页,本课件共有61页

展开阅读全文
相关资源
相关搜索

当前位置:首页 > 生活休闲 > 资格考试

本站为文档C TO C交易模式,本站只提供存储空间、用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。本站仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知淘文阁网,我们立即给予删除!客服QQ:136780468 微信:18945177775 电话:18904686070

工信部备案号:黑ICP备15003705号© 2020-2023 www.taowenge.com 淘文阁