专利名称：通信系统中基于差别相关性的框架同步装置与方法
技术领域：
本发明涉及一种通信系统，且特别涉及一种通信系统中基于差别相关性 (differential correlation)的框架同步(frame synchronization)装置与方法。
背景技术：
通信系统(例如DVB-S或DVB-S2)可以使用多种不同的框架型态(frametype)。通信系统的框架同步(frame synchronization)装置用于侦测出实体层框架标头(Physical Layer Frame Header ；以下简称PLHeader)的位置，以使后续的频偏修正(frequency offset correlation)与前向纠错解码器(FEC decoder)能够运作。不同于DVB-S的是， DVB-S2的框架长度更长，同时讯杂比(Signal to NoiseRatio ；以下简称SNR)更低。例如，以带有前导信号(pilot)四相相移键控(Quaternary Phase Shift Keying;以下简称 QPSK)的框架型态为例，其框架长度为33282个符元(symbol)，而讯杂比为-2. 35dB (for 1/4 code rate)。因此，有效利用DVB-S2的框架结构以快速锁定PLHeader的位置是非常重要的。例如，SuruJiang 与 Lee 等发表的论文"Frame Synchronization and PilotStructure for Second Generation DVB via Satellite，， (International Journal ofSatellite Communications and Networking，2004 ;22 :319_339)，此先前技术在搜寻窗 (search window)中计算出多个相关值(correlation)，并从这些相关值中找出最大值。假设第K种框架型态的第T个搜寻窗WK,τ长度为Lk。此先前技术利用同步字元(sync word)对每个符元进行差别检测(differential detection)，以在搜寻窗WK, τ当中多个相关值中找出最大值的索引值(index) Ικ,τ。重复此法在连续多个搜寻窗WK,T、WK,T+1、WK,T+2、...中各自找到最大相关值之索引值为IK,T、IK,T+1、IK,T+2、...。若Ικ,τ = Ικ,τ+ι = Ικ,τ+2 =...，即可确认前述索引值〗!^、、,^、、^、. · ·的对应位置为标头位置，也就是PLHeader位置。然而，此先前技术未利用前导信号的资讯，而且在低杂讯(low noise)情况下侦测漏失(mis-detection) 机率过高，导致同步时间过久。美国专利公告号第US 7，308，064号专利案利用“sync word correlation”步骤以及 “multi-step correlation test” 两大步骤搜寻 PLHeader 的位置。在 “syncword correlation”步骤中，只要PLHeader的相关值大于临界值(threshold) 即进入 "multi-step correlation test，，步骤。在"multi-step correlation test，，步骤中，计算前导信号的相关值大于临界值Γ 2，则测试成功次数I累进加1(即I = 1+1)。当I = Pn 即判定同步成功。然而，“sync word correlation”步骤使用临界值进行判断，在低讯杂比条件下易造成假警报(false-alarm)。“multi-st印correlation test」”步骤针对每个前导信号进行判断，然而Pn大易造成侦测漏失，反之易造成假警报。再者，此现有技术未考虑到卫星通道中的频谱反转(spectrum inversion)现象，导致此现象产生时“multi-st印 correlation test”步骤会失效。另外，此现有技术未考虑无前导信号的框架(non-pilot frame)0
发明内容
本发明提供一种通信系统中基于差别相关性的框架同步装置与方法，以快速搜寻标头位置。本发明实施例提出一种通信系统中基于差别相关性的框架同步方法，包括下述步骤。使用同步字元(sync word)于前搜寻窗与后搜寻窗中各自进行相关性分析 (correlation analysis)，以分别获得多个前窗相关值(correlation)与多个后窗相关值。 于这些前窗相关值中选取N个最大相关值，其中这些最大相关值的索引值分别为I1 IN， 而N为大于0的整数。将这些后窗相关值与第一临界值(threshold)进行比较。若这些后窗相关值其中一个相关值大于第一临界值，则将此相关值于该后搜寻窗的索引值定义为J。 若所述索引值J等于所述索引值I1 In其中一者，则使用该索 引值J进行一检查步骤，以确认该索引值J的对应位置为一标头位置。本发明实施例提出一种通信系统中基于差别相关性的框架同步装置，包括差别检测器以及框架搜寻单元。差别检测器接收基频信号，而于前搜寻窗与后搜寻窗中使用同步字元进行相关性分析，以分别获得多个前窗相关值与多个后窗相关值。框架搜寻单元于这些前窗相关值中选取N个最大相关值，其中这些最大相关值的索引值分别为I1 In，而N为大于0的整数。框架搜寻单元将这些后窗相关值与第一临界值进行比较。若这些后窗相关值其中一个相关值大于第一临界值，则将此相关值于后搜寻窗的索引值定义为J。框架搜寻单元判断所述索引值J是否等于所述索引值I1 In其中一者。其中，所述框架搜寻单元包含一检查功能单元，当所述索引值J等于所述索引值I1 In其中一者时，所述检查功能单元使用索引值J以确认该索引值J的对应位置是否为标头位置。基于上述，本发明实施例于前搜寻窗中搜寻一个或多个最大相关值，以及在后搜寻窗中找出大于第一临界值的相关值。然后，比较前搜寻窗的最大相关值之索引位置I1 In与后搜寻窗中大于第一临界值的相关值的索引位置J。依据前述索引位置J与I1 In的关系，可以快速搜寻标头位置。然后，使用索引值J进行检查步骤，以确认该索引值J的对应位置为标头位置。为让本发明的上述特征和优点能更明显易懂，下文特举实施例，并配合附图作详细说明如下。
图1为本发明实施例说明一种通信系统中基于差别相关性的框架同步装置。图2为本发明实施例说明一种通信系统中基于差别相关性的框架同步方法。图3为本发明实施例说明相关值与时间索引值的关系示意图。图4为本发明实施例说明图2中检查步骤S20的细部流程。图5为本发明另一实施例说明图2中检查步骤S20的细部流程。图6为检查功能单元依照本发明实施例说明图5中步骤S510的示意图。图7为本发明另一实施例说明一种通信系统中基于差别相关性的框架同步方法。主要元件符号说明10:信号处理单元100:框架同步装置
110差别检测器120 框架搜寻单元
605第-一延迟器610 共轭器
640A-Ap · 弟-二延迟器615 第一乘法器
625A-Ap · 弟-二乘法器645 第三乘法器
620第-一相关分析单元660 第二相关分析单元
630加法器635 第一多路复用器
650A-Ap · 弟-二多路复用器670 第三多路复用器
Cthi第--临界值SlO 获取步骤
S20检查步骤WT_i、WT 搜寻窗
Sl 第—-控制信号S2 第二控制信号
S3 第三控制信号665 乘法器
Ιτ-ι,ι、！τ-1,2Λ Jt,ι、JT,2、JT^^JTj4 索引值
S210 S240、S710 S760 获取步骤的子步骤
S410 S470、S510 S520 检查步骤的子步骤
具体实施例方式图1为本发明实施例说明一种通信系统(例如DVB-S或DVB-S2)中基于差别相关性的框架同步装置100。框架同步装置100于基频信号中搜寻框架标头(frame header), 然后将标头位置告知信号处理单元10以完成框架同步。此信号处理单元10可以依据框架同步装置100所提供的标头位置，对基频信号进行后续的频偏修正与/或前向纠错解码等操作。框架同步装置100包括差别检测器110以及框架搜寻单元120。图2为本发明实施例说明一种通信系统中基于差别相关性的框架同步方法。请参照图1与图2，本框架同步方法可分为两部分获取步骤(AcquisitionSt印)SlO与检查步骤(Check St印)S20。获取步骤SlO定义第T个搜寻窗为Wt，每个搜寻窗Wt的长度为 L，并定义搜寻窗中的时间索引值(time index)为I，而第T个搜寻窗PLHeader的相关值 (correlation)为Γτ>1。在搜寻窗W^1找出各相关值中N个最大值，并储存其索引值为It+
......It-!,^设搜寻窗Wt中相关值超过第一临界值Cthi的索引值为J。若J符合以下条
件(J == IT_ia)or(J == IT_lj2)or. ..or(J== It^1j N)则检查步骤 S20 开始。例如，图3为本发明实施例说明相关值与时间索引值的关系示意图。差别检测器110接收基频信号，于前搜寻窗WT_i与后搜寻窗Wt中使用同步字元(sync word)对基频信号进行相关性分析，以分别获得多个前窗相关值与多个后窗相关值(步骤S210)。在某些实施例中，前搜寻窗Wp1与后搜寻窗Wt可能是在时间上不相邻接的两个搜寻窗。应用本实施例者可以视其设计需求而以任何方式实现差别检测器110，例如可以使用ETSI TR 102 376 VL 1. 1 标准规格书关于 second generation system of Digital Video Broadcasting(DVB-S2)中 User guidelines C3. 1 章节(page 77 page 78)所述技术来实现差别检测器110。框架搜寻单元120耦接至差别检测器110，以接收多个相关值。框架搜寻单元120 于前搜寻窗Wh期间的多个相关值(即前窗相关值)中选取并储存N个最大相关值，其中这些最大相关值的索引值I1 In，而N为大于0的整数(步骤S220)。此N值须视设计需求而定，例如N为2。图3的前搜寻窗Wp1期间仅绘出最大的2个相关值，且此2个最大相关值的索引值分别为注记为It-U与It-U,表示此2个索引值I1与I2是属于第T-I个搜寻窗W”的索引值。接下来，框架搜寻单元120将后搜寻窗Wt期间的多个相关值(即后窗相关值)与第一临界值Cthi进行比较(步骤S230)。若该些后窗相关值其中一个相关值大于第一临界值Cthi，则将所述相关值于后搜寻窗Wt的索引值定义为J。若所述索引值J等于所述索引值 I1 In其中一者(步骤S240)，则框架搜寻单元120使用该索引值J进行检查步骤S20，以确认该索引值J的对应位置为一标头位置。反之，若后搜寻窗Wt的所有索引值J皆不符合步骤S240的判断条件，则框架搜寻单元120将以上述相同方法搜寻下一个搜寻窗WT+1期间的多个相关值。以图3为例，差别检测器110依序输出相关值串，因而框架搜寻单元120可以在前搜寻窗Wp1期间储存2个最大相关值的索引值Ipu与Ip1,随着差别检测器110依序输出相关值串，框架搜寻单元120会在搜寻窗Wt期间找到第1个大于第一临界值Cthi的相关值(其索引值注记为Ju)。由于索引值Ju既不等于索引值It-U，也不等于索引值It-U,因此判定索引值Ju的对应位置不是标头位置。所以，框架搜寻单元120舍弃索引值Ju而找到第2个大于第一临界值Cthi的相关值(其索引值注记为Jt,2)。由于索引值JT,2既等于索引值It-U，表示索引值JT,2的对应位置是标头位置的可能性很高。因此，框架搜寻单元120 使用该索引值JT,2进行检查步骤S20，以确认该索引值JT,2的对应位置是否为标头位置。检查步骤S20在低讯杂比的条件下可提升同步速度，降底侦测漏失 (mis-detection)机率。检查步骤S20可以针对带有前导信号的框架以及没有前导信号的框架两种不同情况进行检查。以下将以基频信号框架不含前导信号的情况下，说明检查步骤S20的实现范例。图4为本发明实施例说明图2中检查步骤S20的细部流程。检查步骤 S20包括步骤S410 S470。框架搜寻单元120进行步骤S410，以便依据索引值J对下一个标头(例如PLHeader)进行相关性分析，而获得标头相关值Γτ+1。详细而言，步骤S410 自索引值J位置始，经L个符元(symbol)后使用同步字元进行相关性分析，以获得标头相关值Γτ+1。前述L为搜寻窗的长度。若沿用前述图3的范例，则步骤S410会从索引值Jt, 2位置开始，经L个符元(symbol)后进行相关性分析而获得标头相关值Γτ+1。若标头相关值Γ T+1大于第四临界值Cu (步骤S420)，则第二计数值Nscs累进加1 (即Nscs = Nscs+1，步骤 S440)。第二计数值Ns。s表示检查成功次数。若步骤S420判断结果为否，则判断标头相关值Γ T+1是否小于第六临界值CJ步骤 S430)。前述第四临界值Cu大于第六临界值Q。若标头相关值Γ T+1是介于第四临界值Cu 与第六临界值Q之间，则回到步骤S410以对在下一个标头进行相关性分析而获得下一个标头相关值Γ τ+2，并重复上述操作。若步骤S430判断结果为是，则第三计数值累进加 1 (即Nfa1L = NFAIL+1，步骤S450)。第三计数值Nfa1L表示检查失败次数。在步骤S460判断第二计数值Nscs是否等于第五临界值NSCS_TH。若第二计数值Nscs 达第五临界值NS。S_TH，则确认索引值J的对应位置为该标头位置。若步骤S460判断结果为否，则判断第三计数值是否达第七临界值NFAtt_TH (步骤S470)。若第三计数值Npm达第七临界值NFm_TH，则确认该索引值J的对应位置不是标头位置。若第二计数值Ns。s未达第五临界值Ns。s_th，且第三计数值未达第七临界值NFAtt_TH，则回到步骤S410以对在下一个标头进行相关性分析而获得下一个标头相关值Γ τ+2，并重复上述操作。经过L个符元(symbol)后，若相关值Γτ>1大于上限门槛(第四临界值Cu),则成功次数累进加1 (即Nscs = Nscs+1)。反之，若相关值Γ u小于下限门槛(第六临界值⑦， 则失败次数累进加1 (即= Nfail+1)。设定上下两个门槛Cu与Q的功用是保留未决区域(non-decision area)，使得在低讯杂比的情况下可同时降低侦测漏失(mis-detection) 与假警报(false alarm)的机率。以下将以基频信号框架含有前导信号的情况下，说明检查步骤S20的实现范例。 图5为本发明另一实施例说明图2中检查步骤S20的细部流程。检查步骤S20包括步骤 S510 S520。框架搜寻单元120进行步骤S510，以便依据索引值J于后搜寻窗Wt中对化 个前导信号进行相关分析，以获得前导信号相关累积值NPU_S。S。若前导信号相关累积值化『％ 达第三临界值Nm_ses_TH (步骤S520)，则确认该索引值J的对应位置为该标头位置。例如，沿用前述图3的范例，则该索引值J为Jt,2。由于前导信号在框架中的位置，以及前导信号的内容均属已知，因此框架搜寻单元120在步骤S510中可以依据索引值JT,2的位置往后找到一个或多个前导信号，并且对这些前导信号进行相关分析，然后累计分析成功次数(即前导信号相关累积值NPU_S。S)。若前导信号相关累积值Nm_s。s大于第三临界值Npu_s。s_th则判定为同步成功。Npu是依据设计需求而决定的可调整参数。Npu增加则判定的准确度越高，但同步所需时间也越久。图6为检查功能单元，位于框架搜寻单元120之内，为依照本发明实施例说明图5 中步骤S510的示意图。第一延迟器605将基频信号延迟后输出至第一乘法器615。共轭器 610将基频信号的共轭复数输出给第一乘法器615。第一乘法器615接收第一延迟器605与共轭器610的输出并进行乘法运算。第一相关分析单元620与第二相关分析单元660的实施方式相仿。考虑到卫星通道中的频谱反转(spectrum inversion)现象，第二相关分析单元660被用来分析频谱反转的前导信号。以下将说明第一相关分析单元620的实施细节， 而第二相关分析单元660可以参照第一相关分析单元620的说明。在第一相关分析单元620中，第二乘法器625将第一乘法器615的输出乘上系数 Cp (即前导信号的样本)，以将相关值输出给加法器630。加法器630将第二延迟器640所输出的相关累积值以及第二乘法器625所输出的相关值进行加总，然后将新的相关累积值传送给多路复用器635。第一多路复用器635依据第一控制信号Sl的控制而决定是否将加法器630所输出的新相关累积值传送给第二延迟器640。因此，当前导信号来临时，藉由将控制信号Sl设为1，而使前导信号段(pilot segment)的相关值可以进行加总。当不在前导信号期间内时，藉由将控制信号Sl设为0，而避免将非前导信号的相关值加入相关累积值中。第三乘法器645计算前述相关累积值的平方(或是求取前述相关累积值的绝对值)。 第二多路复用器650依据第二控制信号S2的控制而决定是否将第三乘法器645所输出的前导信号相关累积值传送给第三多路复用器670。此控制信号S2与前述前导信号个数Npu 有关。当完成Nm个前导信号的相关分析后，控制信号S2会被设为1而将前导信号相关累积值Npm传送给第三多路复用器670。第二相关分析单元660相似于第一相关分析单元620，不同之处在于第二相关分析单元660内乘法器665是接收系数Cp ’(即频谱反转后的前导信号样本)，而不是系数Cp。因此相关分析单元660可以分析频谱反转的前导信号，然后将前导信号相关累积值Npu2传送给第三多路复用器670。当前导信号相关累积值Npm大于前导信号相关累积值Nm2时，控制信号S3会被设为1，否则控制信号S3会被设为0。因此第三多路复用器670可以选择前导信号相关累积值Npm与化 2 二者较大者，并将较大者输出为前导信号相关累积值Nm_scs， 以供图5步骤S520进行判断。所以，不论通道是否发生频谱反转现象，步骤S510均可以对这些前导信号进行相关分析，然后累计分析成功次数(即前导信号相关累积值NPU_SCS)。图7为本发明另一实施例说明一种通信系统中基于差别相关性的框架同步方法。 图7所示方法大致上相同于图2所示方法，不同之处在于图7所示方法更于获取步骤SlO 内包含步骤S710 S760。为提高搜寻效率，在进入检查步骤S20的同时，获取步骤SlO继续进行步骤S710。请参照图7，若步骤S240判断所述索引值J等于索引值I1 In其中一者，则进行步骤S710以判断第一计数值Nakj是否为初始值0。若第一计数值Nakj为0，则P = Jt (步骤 S720)，以及第一计数值累进加1 (即Naqj = Nacq+1，步骤S730)。接下来进行步骤S760以判断第一计数值Naqj是否达第二临界值CA。Q FEM。若步骤S760的判断结果为否，则框架搜寻单元120将回到步骤S210以搜寻下一个搜寻窗WT+1期间的多个相关值。若步骤S710判断第一计数值Nakj不等于初始值0，则进行步骤S740以判断索引值Jt是否等于P。若步骤S740的判断结果为是，则进行步骤S730，即Na。q = Naqj+1。若步骤 S740的判断结果为否，则进行步骤S750，即P = JT。完成步骤S750后接着进行步骤S760， 以判断第一计数值Nakj是否达第二临界值CA。Q—FEM。若第一计数值Naqj达第二临界值Caqj fkm， 则确认该索引值J的对应位置为标头位置，因此判定为同步成功。综上所述，本实施例于前搜寻窗中搜寻一个或多个最大相关值，以及在后搜寻窗中找出大于第一临界值的相关值。然后，比较前搜寻窗的最大相关值之索引位置I1 In与后搜寻窗中大于第一临界值的相关值之索引位置J。依据前述索引位置J与I1 In的关系，可以快速搜寻标头位置。然后，使用索引值J进行检查步骤S20，以确认该索引值J的对应位置为标头位置。虽然本发明已以实施例揭示如上，但是并非用以限定本发明，任何本技术领域人员，在不脱离本发明的精神和范围内，当可作些许的更动与润饰，故本发明额保护范围当视所附的权利要求所界定者为准。
权利要求
1.一种通信系统中基于差别相关性的框架同步方法，包括于一前搜寻窗与一后搜寻窗中使用一同步字元进行一相关性分析，以分别获得多个前窗相关值与多个后窗相关值；于所述多个前窗相关值中选取N个最大相关值的索引值I1 In，其中N为大于O的整数；将所述多个后窗相关值与一第一临界值进行比较；若所述多个后窗相关值其中一个相关值大于所述第一临界值，则将所述相关值于所述后搜寻窗的索引值定义为J ；以及若所述索引值J等于所述索引值I1 In其中一者，则使用所述索引值J进行一检查步骤，以确认所述索引值J的对应位置为一标头位置。
2.根据权利要求1所述通信系统中基于差别相关性的框架同步方法，还包括 若所述索引值J等于所述索引值I1 In其中一者，则一第一计数值累进加1 ；以及若所述第一计数值达一第二临界值，则确认所述索引值J的对应位置为所述标头位置。
3.根据权利要求1所述通信系统中基于差别相关性的框架同步方法，其中所述检查步骤包括依据所述索引值J于所述后搜寻窗中进行至少一个前导信号分析，以获得一前导信号相关累积值；以及若所述前导信号相关累积值达一第三临界值，则确认所述索引值J的对应位置为所述标头位置。
4.根据权利要求1所述通信系统中基于差别相关性的框架同步方法，其中所述检查步骤包括自所述索引值J位置始，经L个符元后使用所述同步字元进行所述相关性分析，以获得一标头相关值，其中L为所述前搜寻窗或所述后搜寻窗的长度；若所述标头相关值大于一第四临界值，则一第二计数值累进加1 ；以及若所述第二计数值达一第五临界值，则确认所述索引值J的对应位置为所述标头位置。
5.根据权利要求4所述通信系统中基于差别相关性的框架同步方法，其中所述检查步骤还包括若所述标头相关值小于一第六临界值，则一第三计数值累进加1，其中所述第六临界值小于所述第四临界值；以及若所述第三计数值达一第七临界值，则确认所述索引值J的对应位置不是所述标头位置。
6.一种通信系统中基于差别相关性的框架同步装置，包括一差别检测器，接收一基频信号，于一前搜寻窗与一后搜寻窗中使用一同步字元进行一相关性分析，以分别获得多个前窗相关值与多个后窗相关值；一框架搜寻单元，于所述多个前窗相关值中选取N个最大相关值的索引值I1 In，其中 N为大于O的整数；将所述多个后窗相关值与一第一临界值进行比较；若所述多个后窗相关值其中一个相关值大于所述第一临界值，则将所述相关值于所述后搜寻窗的索引值定义为J，并判断所述索引值J是否等于所述索引值I1 In其中一者；其中，所述框架搜寻单元包含一检查功能单元，当所述索引值J等于所述索引值I1 In其中一者时，使用所述索引值 J以确认所述索引值J的对应位置为一标头位置。
7.根据权利要求6所述通信系统中基于差别相关性的框架同步装置，其中若所述索引值J等于所述索引值I1 In其中一者，则所述框架搜寻单元将一第一计数值累进加1 ；以及若所述第一计数值达一第二临界值，则所述框架搜寻单元确认所述索引值J的对应位置为所述标头位置。
8.根据权利要求6所述通信系统中基于差别相关性的框架同步装置，其中该检查功能单元包括一第一延迟器，将一基频信号延迟并输出；一共轭器，计算所述基频信号的共轭复数并输出；一第一乘法器，接收所述第一延迟器与所述共轭器的输出并进行乘法运算；以及一第一相关分析单元，接收所述第一乘法器的输出并进行多个前导信号的相关分析， 以输出一第一前导信号相关累积值。
9.根据权利要求8所述通信系统中基于差别相关性的框架同步装置，其中所述第一相关分析单元包括一第二乘法器，将所述第一乘法器的输出乘上一前导信号样本；一加法器，将一相关累积值以及所述第二乘法器的输出进行加总；一第一多路复用器，依据一第一控制信号而决定是否将所述加法器的输出通过所述第一多路复用器；一第二延迟器，所述第二延迟器的输入端连接至所述第一多路复用器的输出端，所述第二延迟器的的输出端提供所述相关累积值给所述加法器；一第三乘法器，计算前述第二延迟器的输出的平方；以及一第二多路复用器，依据一第二控制信号而决定是否将所述第三乘法器的输出通过所述第二多路复用器，其中所述第二多路复用器输出所述第一前导信号相关累积值。
10.根据权利要求8所述通信系统中基于差别相关性的框架同步装置，其中该检查功能单元还包括一第二相关分析单元，接收所述第一乘法器的输出并进行多个频谱反转后的前导信号的相关分析，以输出一第二前导信号相关累积值；以及一第三多路复用器，选择所述第一前导信号相关累积值与所述第二前导信号相关累积值二者较大者，并将较大者输出。
全文摘要
本发明提供了一种通信系统中基于差别相关性的框架同步装置与方法。使用同步字元于前搜寻窗与后搜寻窗中各自进行相关性分析，以分别获得多个前窗相关值与多个后窗相关值。于这些前窗相关值中选取N个最大相关值，其中这些最大相关值的索引值分别为I1～IN，而N为大于0的整数。将这些后窗相关值与第一临界值进行比较。若这些后窗相关值其中一个相关值大于第一临界值，则将此相关值于该后搜寻窗的索引值定义为J。若所述索引值J等于所述索引值I1～IN其中一者，则使用该索引值J进行一检查步骤，以确认该索引值J的对应位置为一标头位置。
文档编号H04L1/00GK102223224SQ20101015191
公开日2011年10月19日 申请日期2010年4月19日 优先权日2010年4月19日
发明者卢耿志 申请人:联咏科技股份有限公司