信号分析基础

上传人:cl****1 文档编号:568223742 上传时间:2024-07-23 格式:PPT 页数:84 大小:2.18MB
返回 下载 相关 举报
信号分析基础_第1页
第1页 / 共84页
信号分析基础_第2页
第2页 / 共84页
信号分析基础_第3页
第3页 / 共84页
信号分析基础_第4页
第4页 / 共84页
信号分析基础_第5页
第5页 / 共84页
点击查看更多>>
资源描述

《信号分析基础》由会员分享,可在线阅读,更多相关《信号分析基础(84页珍藏版)》请在金锄头文库上搜索。

1、第二章第二章 信号分析基础信号分析基础 2.1 信号信号一、概述一、概述 信号可以分为确定性信号和非确定性信号两大类。信号可以分为确定性信号和非确定性信号两大类。确定性信号确定性信号确定性信号确定性信号 是指可以用数学关系式或图表来明确描述是指可以用数学关系式或图表来明确描述其关系的信号。其关系的信号。非确定性信号非确定性信号非确定性信号非确定性信号具有随机特点,每次观察的结果都不相具有随机特点,每次观察的结果都不相同,无法用数学式或图表描述其关系,更不能确切预同,无法用数学式或图表描述其关系,更不能确切预测,只能用概率统计方法由过去估计未来,因此也叫测,只能用概率统计方法由过去估计未来,因此

2、也叫随机信号。随机信号。 确定性信号确定性信号又可分为又可分为周期信号周期信号和和非周期信号非周期信号随机信号随机信号又可分又可分平稳平稳和和非平稳非平稳的信号两种的信号两种周期信号周期信号是经过一定时间可以重复出现的信号,是经过一定时间可以重复出现的信号,满足条件:满足条件: x ( t ) = x ( t + N t )式中:式中:T周期,周期,T20; 0基频基频 N0,十十1确定信号与随机信号确定信号与随机信号当信号是一确定的时间函数时,给定某一时间值,就可以确定一相应的函数值。这样的信号称为确定信号。随机信号不是确定的时间函数,只知道该信号取某一数值的概率。带有信息的信号往往具有不可

3、预知的不确定性,是一种随机信号。除实验室发生的有规律的信号外,通常的信号都是随机的,因为确定信号对受信者不可能载有信息。 连续信号与离散信号连续信号与离散信号如果在某一时间间隔内,对于一切时间如果在某一时间间隔内,对于一切时间值,除若干不连续点外,该函数都能给值,除若干不连续点外,该函数都能给出确定的函数值,此信号称为连续信号。出确定的函数值,此信号称为连续信号。和连续信号相对应的是离散信号。代表和连续信号相对应的是离散信号。代表离散信号的时间函数只在某些不连续的离散信号的时间函数只在某些不连续的时间值上给定函数值。时间值上给定函数值。一般而言,模拟信号是连续的(时间和一般而言,模拟信号是连续

4、的(时间和幅值都是连续的),数字信号是离散的幅值都是连续的),数字信号是离散的。连续信号模拟信号连续信号f(t)0t0tf(t)f0f1f2离散信号01234-1tf(tk)(3)(2)(4.5)(1.5)(6)(-1) 周期信号与非周期信号周期信号与非周期信号用确定的时间函数表示的信号,可以分为周期信号和非周期信号。当且仅当 则信号f(t)是周期信号,式中常数T 是信号的周期。换言之,周期信号是每隔固定的时间又重现本身的信号,该固定的时间间隔称为周期。非周期信号无此固定时间长度的循环周期。严格数学意义上的周期信号,是无始严格数学意义上的周期信号,是无始无终地重复着某一变化规律的信号。无终地重

5、复着某一变化规律的信号。实际应用中,周期信号只是指在较长实际应用中,周期信号只是指在较长时间内按照某一规律重复变化的信号。时间内按照某一规律重复变化的信号。实际上周期信号与非周期信号之间没实际上周期信号与非周期信号之间没有绝对的差别,当周期信号有绝对的差别,当周期信号f fT T( (t t) )的周的周期期T T 无限增大时,则此信号就转化为无限增大时,则此信号就转化为非周期信号非周期信号f f( (t t) )。即即 周期方波的描述周期方波的描述二、能量信号与功率信号二、能量信号与功率信号1、能量信号、能量信号 在所分析的区间在所分析的区间(-,)内,能量为内,能量为有限值的信号称为能量信

6、号,满足条件有限值的信号称为能量信号,满足条件信号能量的解释信号能量的解释:对于电信号,通常是电压或电流,电压:对于电信号,通常是电压或电流,电压和电流在己知区间和电流在己知区间 (t1, t2) 内消耗在电阻内消耗在电阻R上的能量为上的能量为R R1 1时。上述两式具有相同形式。时。上述两式具有相同形式。 定义定义:当区间:当区间( (t1, t2)t1, t2)为为(-(-,)时,时,能量为有限值的信号称为能量信号,或能量为有限值的信号称为能量信号,或称为能量有限信号,称为能量有限信号, 2、功率信号、功率信号 有许多信号,如周期信号、随机信有许多信号,如周期信号、随机信号等,它们在区间号

7、等,它们在区间(-(-,)内能量不是内能量不是有限值在这种情况下,研究信号的平均有限值在这种情况下,研究信号的平均功率更为合适。功率更为合适。显而易见,一个能量信号具有显而易见,一个能量信号具有0平均功率,平均功率,而一个功率信号具有无限大能量。而一个功率信号具有无限大能量。2.2 信号的时域统计分折信号的时域统计分折一、均值一、均值、方差方差2 和均方值和均方值2各态历经信号的均值各态历经信号的均值为为方差方差方差的正平方根叫方差的正平方根叫标准差标准差 ,是随机数据是随机数据分析的重要参数。分析的重要参数。均方值均方值2描述随机信号的强度描述随机信号的强度 均方值的正平方根称为均方值的正平

8、方根称为均方根值均方根值,可表,可表示为示为xrms。均值、方差和均方值的相互关系是均值、方差和均方值的相互关系是二、概率密度函数二、概率密度函数 随机信号的概率密度函数表示信号随机信号的概率密度函数表示信号幅值落在指定区间内的概率。幅值落在指定区间内的概率。定义幅值概率密度函数定义幅值概率密度函数概率密度函数分析仪原理方框图概率密度函数分析仪原理方框图典型信号的概率密度函数典型信号的概率密度函数 含正弦波随机信号的概率密度函效含正弦波随机信号的概率密度函效23 信号的相关分析信号的相关分析变量变量x和和y之间的相关程度常用相关系数表示之间的相关程度常用相关系数表示式中:式中:E数学期望数学期

9、望 xEx随机变量随机变量x的均值的均值 yEy随机变量随机变量y的均值的均值 x,y随机变量随机变量x,y的标准差:的标准差:可以证明,可以证明,xy1。当当xy 1时,则所有的点都落在时,则所有的点都落在y-y = m (x-x) 的直线上,说明的直线上,说明x,y两变量两变量是理想的线性相关。是理想的线性相关。xy -1也是理想的线性相关,只是直线的也是理想的线性相关,只是直线的斜率为负。斜率为负。xy 0 表示表示x,y两变量之间完全无关。两变量之间完全无关。用用Rx()表示自相关函数:表示自相关函数: 则则例例1 求正弦函数求正弦函数 的自相关的自相关函数函数作业:作业:1、证明均值

10、、方差和均方值的相互关、证明均值、方差和均方值的相互关系系2*、根据图、根据图2-4方框图设计一个概率密方框图设计一个概率密度函数分析仪电路,画出电路原理图。度函数分析仪电路,画出电路原理图。 二、互相关函数二、互相关函数对于各态历经过程,两个随机信号对于各态历经过程,两个随机信号x(t)和和y(t)的互相关函的互相关函数及数及Rxy()定义为定义为互相关函数不是偶函数:互相关函数不是偶函数:四、相关函数的性质四、相关函数的性质根据定义,相关函数有如下性质:根据定义,相关函数有如下性质:1、自相关函数是偶函数、自相关函数是偶函数 互相关函数不是偶函数,也不是奇函数,而满足下互相关函数不是偶函数

11、,也不是奇函数,而满足下式式 2、自相关函数在、自相关函数在=0处取得最大值处取得最大值 这个性质极为重要,它是相关技术这个性质极为重要,它是相关技术确定同名点的依据确定同名点的依据 两边取时间两边取时间T的平均值并取极限的平均值并取极限 3、 周期信号的自相关函数仍然是同周期信号的自相关函数仍然是同 频率的周期信号,但不具有原信号频率的周期信号,但不具有原信号的相位信息。的相位信息。4、随机信号的自相关函数将随、随机信号的自相关函数将随值值增大而很快趋于零。增大而很快趋于零。互相关函数具有以下性质: 两周期信号具有相同的频率,才有互相关函数,即两个非同频的周期信号是不相关的。 两个相同周期的

12、信号的互相关函数仍是周期函数,其周期与原信号的周期相同,并不丢失相位信息。 两信号错开一个时间间隔0处相关程度有可能最高,它反映两信号x(t)、y(t)之间主传输通道的滞后时间。 五、相关分析应用五、相关分析应用1、影像相关影像相关原理原理 影像相关是利用互相关函数,评价两块影像的相似性以确定同名点 。示意图示意图 互相互相关函关函数数相似程相似程度度同名点同名点目目标标区区搜搜索索区区影像匹配影像匹配-同名点寻找同名点寻找2、电子相关、电子相关n 电子相关就是采用电子线路构成的相关器电子相关就是采用电子线路构成的相关器来实现相关的功能来实现相关的功能 图5-1-1电子相关3、光学相关、光学相

13、关n光的干涉和衍射光的干涉和衍射-傅立叶变换特性傅立叶变换特性 相干光学计算机相干光学计算机 n相干光学相关系统相干光学相关系统 三个傅立叶透镜三个傅立叶透镜L1,L2,L3及激光及激光源与光电倍增管等器件组成源与光电倍增管等器件组成 3、数字相关、数字相关n二维相关二维相关 数字相关是利用计算机对数字影像进数字相关是利用计算机对数字影像进行数值计算的方式完成影像的相关行数值计算的方式完成影像的相关 目标区目标区 搜搜索索区区相相似似性性测测度度4工程应用工程应用2.4 信号的频域分析信号的频域分析确定信号的时间特性确定信号的时间特性表示信号的时间函数,包含了信号的全部表示信号的时间函数,包含

14、了信号的全部信息量,信号的特性首先表现为它的时间信息量,信号的特性首先表现为它的时间特性。特性。时间特性主要指信号随时间变化快慢、幅时间特性主要指信号随时间变化快慢、幅度变化的特性。度变化的特性。同一形状的波形重复出现的周期长短同一形状的波形重复出现的周期长短信号波形本身变化的速率(如脉冲信号的脉信号波形本身变化的速率(如脉冲信号的脉冲持续时间及脉冲上升和下降边沿陡直的程冲持续时间及脉冲上升和下降边沿陡直的程度)度)以时间函数描述信号的图象称为时域图,以时间函数描述信号的图象称为时域图,在时域上分析信号称为时域分析。在时域上分析信号称为时域分析。确定信号的频率特性确定信号的频率特性信号还具有频

15、率特性,可用信号的频谱函数来表示。在频谱信号还具有频率特性,可用信号的频谱函数来表示。在频谱函数中,也包含了信号的全部信息量。函数中,也包含了信号的全部信息量。频谱函数表征信号的各频率成分,以及各频率成分的振幅和频谱函数表征信号的各频率成分,以及各频率成分的振幅和相位。相位。频谱:对于一个复杂信号,可用傅立叶分析将它分解为许频谱:对于一个复杂信号,可用傅立叶分析将它分解为许多不同频率的正弦分量,而每一正弦分量则以它的振幅和多不同频率的正弦分量,而每一正弦分量则以它的振幅和相位来表征。将各正弦分量的振幅与相位分别按频率高低相位来表征。将各正弦分量的振幅与相位分别按频率高低次序排列成频谱。次序排列

16、成频谱。频带:复杂信号频谱中各分量的频率理论上可扩展至无限,频带:复杂信号频谱中各分量的频率理论上可扩展至无限,但因原始信号的能量一般集中在频率较低范围内,在工程但因原始信号的能量一般集中在频率较低范围内,在工程应用上一般忽略高于某一频率的分量。频谱中该有效频率应用上一般忽略高于某一频率的分量。频谱中该有效频率范围称为该信号的频带。范围称为该信号的频带。以频谱描述信号的图象称为频域图,在频域上分析信号称为以频谱描述信号的图象称为频域图,在频域上分析信号称为频域分析。频域分析。时域和频域 时域特性与频域特性的联系时域特性与频域特性的联系信号的频谱函数和信号的时间函数既然都包含了信号的全部信息量,

17、都能表示出信号的特点,那么,信号的时间特性与频率特性必然具有密切联系。例:周期性脉冲信号的重复周期的倒数就是该信号的基波频率,周期的大或小分别对应着低的或高的基波和谐波频率;信号分析中将进一步揭示两者的关系。不同频率信号的时域图和频域图不同频率信号的时域图和频域图信号分析时域分析信号时域分析(线性系统叠加原理)卷积积分的应用及其数学描述频域分析周期信号的频域分析(三角与指数傅立叶级数)非周期信号的频域分析(傅立叶积分)信号在频域与时域之间的变换(正反傅立叶变换式)频谱与时间函数的关系0000tttS(t)r(kt)r(t)ktktkts(kt)时时域域分分析析法法示示意意图图r(kt)激励函数

18、激励函数(输入输入信号信号)的分解的分解第第k个脉冲的个脉冲的冲激响应冲激响应(输输出信号出信号)波形波形冲激响应叠加冲激响应叠加后的总响应后的总响应(输输出信号出信号)波形波形第第k个脉冲函数之面积个脉冲函数之面积(当(当t 0,脉冲函脉冲函数可近似表示为冲激函数可近似表示为冲激函数)数)系统对第系统对第k个冲激函数个冲激函数的冲激响应函数的冲激响应函数 时域分析的方法时域分析的方法式中h(t)是单位冲激函数(t)对应的响应,称为单位冲激响应函数。单位冲激函数(t) 也称狄拉克函数或函数,其定义是:在t0时,函数值均为0;在t=0处,函数值为无穷大,而脉冲面积为1,即当t无限趋小而成为d时,

19、上式中不连续变量kt成了连续变量,对各项求和就成了求积分。于是有 这种叠加积分称为卷积积分卷积积分。频域分析频域分析作为时间函数的激励和响应,可通过傅立叶变换将时间变量变换为频率变量去进行分析,这种利用信号频率特性的方法称为频域分析法。频域是最常用的一种变换域。如同时域分析把信号始终看成是时间的函数一样,在频域分析中,任何信号又可看成是频率函数。频域分析的基本工具是傅立叶分析,包括傅立叶级数和傅立叶变换。 周期信号的频域分析方法周期信号的频域分析方法考察信号考察信号 式中式中1 1=2=2f f1 1。1 1称为称为基波频率基波频率,简称基频,简称基频,1 1的倍数称为的倍数称为谐波谐波。该信

20、号的波形图和其频谱图见下图。该信号的波形图和其频谱图见下图。对于周期信号而言,其频谱由离散的频率成分,对于周期信号而言,其频谱由离散的频率成分,即基波与谐波构成。图中,每一条谱线代表一个即基波与谐波构成。图中,每一条谱线代表一个正弦分量,谱线的位置代表这一正弦分量的角频正弦分量,谱线的位置代表这一正弦分量的角频率,谱线的高度代表该正弦分量的振幅。信号率,谱线的高度代表该正弦分量的振幅。信号f f( (t t) )的成分正好是角频率为的成分正好是角频率为1 1、331 1 、551 1和和7 71 1的正弦波。的正弦波。复杂周期信号波形复杂周期信号波形数字信号的谐波数字信号的谐波 分解周期信号的

21、条件狄利希莱条件 要将一周期信号分解为谐波分量,代表这一周期信号的函数f(t)应当满足下列条件:在一周期内,函数是绝对可积的,即在一周期内,函数是绝对可积的,即 应为有限值;应为有限值;在一周期内,函数的极值数目为有限;在一周期内,函数的极值数目为有限;在一周期内,函数在一周期内,函数f(t)或者为连续的,或者具有有限或者为连续的,或者具有有限个这样的间断点,即当个这样的间断点,即当t从较大的时间值和较小的时从较大的时间值和较小的时间值分别趋向间断点时,函数具有两个不同的有限的间值分别趋向间断点时,函数具有两个不同的有限的函数值。函数值。测试技术中的周期信号,大都满足该条件。测试技术中的周期信

22、号,大都满足该条件。 周期信号的频域分析方法周期信号的频域分析方法根据傅立叶变换原理,通常任何信号都可表示成各种频率成分的正弦波之和。对于任何一个周期为T、且定义在区间(- T/2, T/2)内的周期信号f(t),都可以用上述区间内的三角傅立叶级数表示:a0是频率为零的直流分量(如图),式中系数值为傅立叶级数的这种形式称为三三角角函函数数展展开开式式或称正弦-余弦表示,是用正交函数集来表示周期信号的一种常用方法。傅立叶级数还可以改写成:傅立叶级数还可以改写成:An-,n-分别称为幅值谱和相位谱,统称为频谱。带有直流分量的信号带有直流分量的信号 指数傅立叶级数指数傅立叶级数用正交函数集来表示周期

23、信号另一种更常用的方法是傅立叶级数的指数表示法,称为指数傅立叶级数。三角傅立叶级数与指数傅立叶级数并不是两种不同类型的级数,而只是同一级数的两种不同的表示方法。指数级数形式比三角级数形式更简化更便于计算。根据欧拉公式当当当当n n n n取取取取-和和和和+之间包括之间包括之间包括之间包括0 0 0 0在内在内在内在内的所有整数,则函数集的所有整数,则函数集的所有整数,则函数集的所有整数,则函数集e e e ejnjnjnjntttt(其其其其中中中中n n n n=0,1,2,=0,1,2,=0,1,2,=0,1,2,)为一完备为一完备为一完备为一完备的正交函数集。任意周期信号的正交函数集。

24、任意周期信号的正交函数集。任意周期信号的正交函数集。任意周期信号f f f f(t)(t)(t)(t)可在时间区间(可在时间区间(可在时间区间(可在时间区间(- - - -T/2T/2T/2T/2, , , ,T/2T/2T/2T/2)内用此函数集表示为内用此函数集表示为内用此函数集表示为内用此函数集表示为 求出求出求出求出C C C Cn n n n,信号分解的任务就完成了。信号分解的任务就完成了。信号分解的任务就完成了。信号分解的任务就完成了。 非周期信号的频域分析方法非周期信号的频域分析方法对于定义于区间(-,+)上的非周期函数,也能分解成许多正弦波的叠加。(也要满足狄利希莱条件)如果在

25、表示周期信号f(t)的傅立叶级数中令周期T,则在整个时间内表示f(t)的傅立叶级数也能在整个时间内表示非周期信号。f (t)的指数傅立叶级数可写为 式中Fn是复数振幅,将其代入f(t),得到 非周期信号的频域分析方法非周期信号的频域分析方法当T 增加时,基频1变小,频谱线变密,且各分量的振幅也减小,但频谱的形状不变。在T的极限情况下,每个频率分量的幅度变为无穷小,而频率分量有无穷多个,离散频谱变成了连续频谱。这时,f(t)已不是n1的离散函数,而是的连续函数。以上过程可以用计算式说明。由于相邻频率分量间隔为=(=(n n+1)+1)1 1- -nn1 1= =1 1 周期T 可写为 于是,有

26、非周期信号的频域分析方法非周期信号的频域分析方法当T 时,求和变成了取积分,变成d d ,n1用表示。因此有式中方括号是原函数f(t)的频谱密度函数,简称频谱函数,它具有单位频带振幅的量纲,记作F() 。即将原函数写成这就是非周期信号f(t)的傅立叶积分表示式,它与周期信号的傅立叶级数相当。 和傅立叶级数中的复数振幅相当,是无穷小量,频谱密度函数反映了各分量振幅间的相对比例关系。傅立叶变换傅立叶变换通过非周期信号的频谱分析得知,时域上的原函数中含有包含全部信息量的频谱函数,而频谱函数中也含有原函数。因此我们可以在时域与频域之间对信号进行相互变换。这种变换通过称之为傅立叶变换式的公式来实现。即我

27、们前面已经推导出的一对傅立叶积分表示式:前者称为傅立叶正变换式,它将时域内t t 的函数变换为频域内的函数;后者称为傅立叶逆变换式或反变换式,可把的函数变换为t t 的函数。傅立叶变换式简记为 傅立叶变换的应用傅立叶变换的应用傅立叶变换可将时域上较复杂的运算简化为相对简单的频域运算。作为时域上卷积积分例子的函数r(t)对应的频域函数为上式即卷积定理卷积定理,激励s(t)通过频率特性为H()的系统时,响应r(t)的频谱函数R()等于s(t)的频谱函数S() 和H()的乘积运算。 频谱与时间函数的关系频谱与时间函数的关系通过时域与频谱分析的讨论,可总结为两个关系式 R() = S() H()r(t

28、) = s(t)*h(t) 其中两个关系式的意义是:两个频谱相乘,其乘积的时间函两个频谱相乘,其乘积的时间函数就是相应的两个时间函数相卷积。数就是相应的两个时间函数相卷积。反之,两个时间函两个时间函数相卷积,其频谱就是相应的两个频谱相乘数相卷积,其频谱就是相应的两个频谱相乘。从滤波角度看,该两关系式的意义是:滤波可以两种方式实现。一是在频域上实现,将频谱H()与 S() 相乘得到R(),再由R()作傅立叶反变换得到r(t)。 二是在时域上直接实现,将时间函数h(t) 与s(t) 相卷积得到r(t)。几种典型信号的傅立叶变换几种典型信号的傅立叶变换数字信号中典型的波形是矩形窗函数(矩形脉冲函数)

29、。矩形脉冲g(t)及其对应的频域函数为G()分别如图和下面两式: 当当=0=0=0=0时,时,时,时, G G( ()=)=A A ; = = = =2 2kk/ / 时,时,时,时, G G( ()=0)=0。 (t)函数的性质函数的性质:1.抽样性抽样性2. 单位脉冲函数的积分等于阶跃函数单位脉冲函数的积分等于阶跃函数3. 函数与其他函数的卷积函数与其他函数的卷积 4. 函数的频谱函数的频谱 功率谱密度和带宽功率谱密度和带宽对于一个矩形脉冲信号,其能量主要集中在频谱中零频率到第一个过零点之间( ) ,所含能量达到信号全部能量的90%以上,故可将其定义为矩形脉冲信号的有效带宽。一般而言,任何

30、一个有限时间的信号之频谱宽度是无限的。然而,信号的大部分功率实际上只集中在某个有限的频谱宽度内。所谓信号的有效带宽就是指包含信号大部分功率的这部分频谱的宽度。见图。为了精确地说明以上概念,需要定义信号的功率谱密度。实际频谱与有效频谱(有效带宽)实际频谱与有效频谱(有效带宽) 信号的能量谱与功率谱信号的能量谱与功率谱除时域和频域的关系外,时间信号的另一个重要特征是能量和功率随时间分布的关系,即能量谱密度和功率谱密度。信号f(t)在1电阻上所消耗的能量定义为信号的归一化能量,简称能量,表示为只有在上式给出的积分值为有限时信号能量的概念才有意义。当信号能量趋于于无穷大时,存在其平均功率,简称功率,即

31、上式可理解为信号f(t)在1电阻上所消耗的平均功率。该平均功率也就是f(t)的均方值,记作 。 信号的能量谱与功率谱信号的能量谱与功率谱帕什瓦尔定理若f(t)为能量信号,且其傅立叶变换为F(),则有如下关系:若f(t)为周期性功率信号,则有: 式中,T 为信号f(t)的周期,Fn为f(t)的傅立叶级数系数。前式说明时域内能量信号的总能量等于频时域内能量信号的总能量等于频域内各个频率分量能量的连续和域内各个频率分量能量的连续和。后式说明周期信号的功率等于该信号在完周期信号的功率等于该信号在完备正交函数集中各分量功率之和。备正交函数集中各分量功率之和。 信号的能量谱与功率谱信号的能量谱与功率谱设能量以E表示,功率以P表示,如果在频域内有 则称E()为能量谱密度函数, P()为功率谱密度函数。能量谱密度和功率谱密度简称能量谱和功率谱。能量谱的单位为J/Hz,功率谱的单位为W/Hz。对于能量信号f(t),其能量谱E()当然一定存在,将前式与帕什瓦尔定理前式对照,可得由于 ,故能量谱是的一个实偶函数,此时信号能量E可简化为fT(t)f(t)ttT2T2-信号信号f(t)及其截断函数及其截断函数

展开阅读全文
相关资源
正为您匹配相似的精品文档
相关搜索

最新文档


当前位置:首页 > 商业/管理/HR > 营销创新

电脑版 |金锄头文库版权所有
经营许可证:蜀ICP备13022795号 | 川公网安备 51140202000112号