CN1682482A - 确定数字调制信号的相位噪声谱和/或振幅噪声谱的方法 - Google Patents

确定数字调制信号的相位噪声谱和/或振幅噪声谱的方法 Download PDF

Info

Publication number
CN1682482A
CN1682482A CN03822303.1A CN03822303A CN1682482A CN 1682482 A CN1682482 A CN 1682482A CN 03822303 A CN03822303 A CN 03822303A CN 1682482 A CN1682482 A CN 1682482A
Authority
CN
China
Prior art keywords
complex
real
sampling
ideal
phase
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN03822303.1A
Other languages
English (en)
Other versions
CN100338901C (zh
Inventor
克里斯托弗·巴尔兹
马丁·霍夫迈斯特
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Rohde and Schwarz GmbH and Co KG
Original Assignee
Rohde and Schwarz GmbH and Co KG
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Rohde and Schwarz GmbH and Co KG filed Critical Rohde and Schwarz GmbH and Co KG
Publication of CN1682482A publication Critical patent/CN1682482A/zh
Application granted granted Critical
Publication of CN100338901C publication Critical patent/CN100338901C/zh
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L1/00Arrangements for detecting or preventing errors in the information received
    • H04L1/20Arrangements for detecting or preventing errors in the information received using signal quality detector
    • H04L1/206Arrangements for detecting or preventing errors in the information received using signal quality detector for modulated signals

Abstract

本发明涉及一种用于确定数字调制输入信号(S)的相位噪声谱和/或振幅噪声谱的方法。所述方法包括以下步骤:通过对基带中的所述输入信号(S)的相位分量(I)和相位正交分量(Q)进行数字采样,产生实际复数采样(Areal [n]) (S100);根据产生的实际采样(Areal [n])确定理想复数采样(Aideal [n])(S101);根据实际复数采样(Areal [n])和理想复数采样(Aideal [n])形成复数商(ΔA1 [n]=Areal [n]/Aideal [n]);通过对所述复数商指定值1,产生经修正复数商(B[n])(S106);以及对经修正复数商(B[n])进行傅立叶变换(S107)。本发明还涉及一种类似的用于确定所述数字调制输入信号的振幅噪声谱的方法。

Description

确定数字调制信号的相位噪声谱和/或振幅噪声谱的方法
技术领域
本发明涉及一种用于确定数字调制信号的相位噪声谱和/或振幅噪声谱的方法。
背景技术
对于数字调制信号的分析和技术测量评估,信号处理中涉及的振荡器的相位噪声谱和振幅噪声谱的图形表示是重要的测量量。在对例如为QAM(正交振幅调制)调制或mVSB(残留边带)调制的电视信号进行数字传输的情况下,这种测量尤其重要。
通常以没有残留载波或只有小残留载波的方式来处理数字调制信号。调制信号的有效频谱分布在相对大的带宽上。然而,要测量的相位噪声或者要测量的振幅噪声的频谱也位于调制信号的该有效频谱内。为了测量相位噪声或振幅噪声的频谱,迄今通常切断调制并传送连续的CW(连续波)信号。然后可以用频谱分析仪来测试该CW信号,并可对相位噪声谱或振幅噪声谱进行检测,这种将相位噪声与振幅噪声分离的方法仍然存在困难。在这种用于测量相位或振幅噪声的操作状态中,不可能同时传送数据。这非常不利,因为为了进行测量必须中断正常的操作,这在操作传输模式中的业务测量过程中是不可能的。
US 6,366,621 B1提出了用于确定8VSB或16VSB信号的基准相位的方法。该文献中提出了用计算机重建导频信号(pilot signal)。使用这种方法,对短期相位波动(相位抖动)的测量、特别是对分布于有效数据频谱上的频谱的测量是不可能的。
发明内容
本发明的主要目的是使得可以在正常调制操作过程中不需切断调制就可以确定数字调制信号的相位噪声谱和/或振幅噪声谱。
关于对相位噪声谱的确定,通过权利要求1的特征来实现所述目的;关于对振幅噪声谱的确定,通过权利要求2的特征来实现所述目的。
本发明的主要原理是,由于测量到的实际复数采样(分别具有同相分量(I)和正交相位分量(Q))与理想复数采样相关联,所以可以对叠加在要测量的相位噪声或振幅噪声的频谱上的有效数据的频谱进行计算。在测量到的实际复数采样与由于调制而产生的理想复数采样之间因此出现的相位差或由此出现的振幅比,是仍然存在的经修正调制的相位波动或振幅波动,它们形成经修正调制的测量量。
因为给定理想基带信号的情况下相位波动或者振幅波动与受控调制相关联,所以由此检测到的相位波动或振幅波动完全独立于刚才传送的调制信号。因此无须中断此操作。例如,无需中断电视发射机传送的节目就可以测量该电视发射机上的相位噪声谱或振幅噪声谱。
可以通过根据测量到的实际复数采样和从其提取的理想复数采样来形成商而以简单的方式产生理想基带信号的基准。通过形成所述商,一方面产生了实际复数采样与理想复数采样之间的相位差,另一方面,产生了实际复数采样的值与理想复数采样的值的振幅比。在对相位噪声谱进行确定的情况下,应该将商的值设置为1。在对振幅噪声谱进行确定的情况下,应该将系数的相位设置为0。在进行傅立叶变换后,可以得到对应的频谱。
权利要求3至7包括本发明的有利发展。
如果待测试信号是mVSB信号,则根据实际采样(实际上根据其同相分量)只确定理想采样的同相分量是有效的。理想采样的正交相位分量是根据理想采样的同相分量、然后通过对这种单边带调制类型常用的希尔伯特变换而产生的。
特别地,当评估mVSB信号时,如果脱离了容许值域,特别是当实际采样的值低于第一阈值或者实际采样的虚部大于第二阈值或小于第三阈值时,用插值来代替实际采样与理想采样的复数商是有利的。
权利要求8至11涉及用于实现根据本发明的方法的数字存储介质、计算机程序或者计算机程序产品。
附图说明
接下来参照附图进一步详细说明本发明的实施例。在附图中示出:
图1是用于说明根据本发明的用于确定相位噪声谱的方法的流程图;
图2是用于说明根据本发明的用于确定振幅噪声谱的方法的流程图;
图3是受相位抖动干扰的8VSB信号的实际星座图(constellationdiagram);
图4是与图3相关联的理想星座图;
图5A是作为采样指标(sample index)n的函数的相位误差Δ;
图5B是实际采样和理想采样的复数系数值(已设置了值=1)的I/Q图;
图5C是由根据本发明的方法确定的相位噪声谱;
图6A是图5A的放大部分;
图6B是与图6A相关联的I/Q图;
图6C是与图6A相关联的相位噪声谱,表示图5C的放大部分;以及
图7是用于说明实现根据本发明的方法的装置的框图。
具体实施方式
接下来参照图1和图7,对用于确定数字调制信号的相位噪声谱的根据本发明的方法、或用于实现根据本发明的方法的装置进行说明。
在图7所示的根据本发明的装置1的情况下,待分析的数字调制高频信号S首先被提供给高频单元2。以普通的方式,该信号通过第一混频器3(其与本机振荡器或可变振荡器4相连)被下混频至中频,并在中频单元5中受到进一步处理。所述中频信号由第二混频器6和第三混频器7变换到基带。为此,第一混频器6直接与第二本机振荡器8相连,第三混频器7通过90°移相器9与本机振荡器8相连。因此,提供给混频器6或7的振荡器信号相对于彼此相移90°。在第二混频器6的输出中,出现基带信号的同相分量I,所述分量通过第一低通滤波器10提供给第一模/数转换器11。在第三混频器7的输出中,可以得到基带信号的正交相位分量Q,所述分量通过第二低通滤波器12提供给第二模/数转换器13。因此,在模/数转换器11和13的输出中,可以得到复数采样Areal[n],其表示输入信号S的复数基带信号。n是采样指标。第一模/数转换器11的输出中的同相分量I表示这些复数采样Areal[n]的实部,而第二模/数转换器13的输出中的正交相位分量Q表示这些复数采样Areal[n]的虚部。
还应该注意,一个前提是基带信号相对于频率和时间同步地出现。
上述的产生实际复数采样Areal[n]对应于图1的流程图中的步骤S100。在方法步骤S101中,根据实际复数采样Areal[n]产生理想复数采样Aideal[n]。为此,在星座图中建立入口区(entry region),其将一个特定的实际I/Q值精确地指定给一个理想I/Q值。参照图3和4,通过8VSB调制的示例对此进行说明,所述8VSB调制例如用于传送数字电视的视频信号。
图3表示实际复数采样Areal[n]的星座图。根本上,该处理可以在将实际复数采样Areal[n]指定给理想复数采样Aideal[n]的过程中按如下方式发生:即,在该指定过程中对I值和Q值都加以考虑,从而使每个理想采样在实际采样上都具有平坦的入口区。这种操作模式例如适合于QAM调制的情况。在图3中出现mVSB调制的情况下,不同的操作模式是有效的:只对实际采样Areal[n]的实部(即,同相分量I)进行评估,并对实际采样Areal[n]的每个同相分量I分别指定理想采样Aideal[n]的一个同相分量I。在图3中,间隔15至22表示实际复数采样Areal[n]的同相分量I的入口区,所述入口区分别精确地指向理想采样Aideal[n]的一个同相分量。
因为在mVSB调制的情况下的每个正交相位分量Q都可以根据同相分量I的时间连续序列通过希尔伯特变换来计算出,所以对应于根据本发明的进展,提出了不根据实际采样Areal[n]的正交相位分量Q来获得理想采样Aideal[n]的正交相位分量Q,而通过对理想采样Aideal[n]的同相分量序列进行希尔伯特变换的计算来获得理想采样Aideal[n]的正交相位分量Q。
由此获得的理想采样Aideal[n]示于图4中。因此很明显,图4中不再存在图3中存在的以标号14表示的值域极限。这里,如果有必要,则仍然必需实行对应的插值措施。这将进一步说明。
上述的根据实际复数采样Areal[n]来产生理想复数采样Aideal[n]是在图7所示的指定装置23中实现的。在商形成装置24中,商
ΔA1[n]=Areal[n]/Aideal[n]    (1)
即,复数商△A1[n]是根据实际复数采样Areal[n]和理想复数采样Aideal[n]来计算的。这在图1的流程图中由步骤S102示出。
如果复数商超过特定值范围而不可靠,则可以在插值器25中实现的可选方法步骤S103中可选地进行对复数商的插值。例如,如果实际采样Areal[n]的虚部Im{Areal[n])大于预定最大值(即大于阈值Amax)、或者小于预定最小值(即小于预定阈值Amin),那么商△A1[n]就不能再由数字格式进行数字表示,并且绝不能考虑这些受限值。应该根据前面和/或后继的值以插值代替这些值。
用于确定商△A1[n]的I/Q值的分辨率是由Areal[n]的量化步骤的数量确定的。因此,△A1[n]的相对误差越大,Areal[n]的值就越小。为了最小化这类随机误差的效应,优选地,如果△A1[n]的值相对小,则同样应该舍弃△A1[n]的值并以插值代替,整个结果不会由此受到篡改。因此,优选地,当实际复数采样Areal[n]的值小于阈值Minvalue时也应该进行插值。
在方法步骤S104中实现上述的插值标准的确定,受到插值作用的采样以标记(标志)U[n]标记。可以在方法步骤S103中针对所有商值△A1[n]计算插值△A2[n],只有当设置了插值标记U[n]时才在方法步骤S105中接收所述商值。插值后(如果有必要进行插值)产生的复数商△A3[n]可以写为如下的极坐标形式:
根据本发明,通过在图1中的步骤S106中或在图7中的修正装置26中将复数商△A3[n]的值|△A3[n]|设为1,产生了经修正的复数系数B[n]以表示相位噪声谱:
当确定相位噪声谱时,不考虑振幅波动而只考虑相位波动的频谱。相位波动由相位差△3[n]确定,因为通过在步骤S102中形成商从而出现了相位差△1[n]=real-ideal(即,实际采样Areal[n]的相位real与理想采样Aideal[n]的相位ideal之间的差)。△3[n]与△1[n]的不同仅在于必要时实施了插值。根据本发明的本质发明点在于:如果对应于根据本发明的方法,通过对理想采样的重构来对受控调制的瞬时相位进行重构,并从测量的实际相位real[n]提出由此重构的基准相位 ideal[n],那么就可以独立于调制所指定的瞬时相位来评估相位波动。
在方法步骤S107或在傅立叶变换单元27中进行傅立叶变换之后,得到相位噪声谱,并可用显示装置28(例如显示器)显示该相位噪声谱。
为了说明本发明,图5A示出了作为采样指标n的函数的相位波动△[n]的示例。在图5B中,示出了经修正的复数系数B[n]的关联I/Q图。检测到值B[n]在具有单位半径的圆上移动。在图5c中,示出了由根据本发明的方法确定的关联相位噪声谱。图6A示出了图5A的放大部分,图6B示出与这部分相关联的对应I/Q图。图6c示出对应地更精确分辨的相位噪声谱。
也可以以类似的方式来估计振幅噪声谱。为此所需的方法步骤示于图2所示的流程图。方法步骤S100到S105与已经参照图1进行了说明的步骤S100到S105相同。在图2中的方法步骤S108中,与图1中的方法步骤S106不同,通过将复数商△A3[n]的相位△3[n]设为0从而产生经修正的复数商B[n]。
B[n]=|△A3[n]|·ej·0    (4)
这样,相位波动不影响通过步骤S107中的傅立叶变换而生成的频谱。相反地,频谱的特征在于(必要时经插值的)商△A3[n]的值|△A3[n]|的波动。在图7中的修正装置29中产生用于振幅噪声谱的经修正的复数商B[n]。可以通过切换装置30对用于傅立叶变换装置27的输入信号在修正装置26和29之间进行切换。
有利地,可以将每滤波器带宽功率密度计算为另一单位,例如dBc/Hz(也就是功率密度=每1Hz带宽的功率),所述滤波器带宽是所使用的FFT(快速傅立叶变换)的设置基数,并且取决于原始I/Q值的时间间隔。这在对噪声干扰的评估中特别有用。在对窄带干扰(CW干扰)进行估计的情况下,明智的做法是保持水平轴(level axis)的单位不变。如果有必要,可以用切换器来选择希望的单位或标度。
本发明不限于所述实施例。相反,在本发明范围内可以进行许多修改和改进。例如,当根据实际采样Areal[n]产生理想采样Aideal[n]时,也可以对通常存在的纠错编码进行估计,结果,因为对错误理想采样Aideal[n]的缺陷分配而造成的不稳定相位和/或振幅波动实际上不存在,所以进一步提高了精确度。

Claims (11)

1、用于确定数字调制输入信号(S)的相位噪声谱的方法,具有以下方法步骤:
通过对基带中的输入信号(S)的同相分量(I)和正交相位分量(Q)进行数字采样,生成实际复数采样(Areal[n])(S100),
根据实际采样(Areal[n])确定理想复数采样(Aideal[n])(S101),
根据实际复数采样(Areal[n])和理想复数采样(Aideal[n]),形成复数商(ΔA1[n]=Areal[n]/Aideal[n])(S102),
通过将复数商的值设为1,产生经修正复数商(B[n])(S106),以及
对经修正复数商(B[n])进行傅立叶变换(S107)。
2、用于确定数字调制输入信号(S)的振幅噪声谱的方法,具有以下步骤:
通过对基带中的输入信号(S)的同相分量(I)和正交相位分量(Q)进行数字采样,生成实际复数采样(Areal[n])(S100),
根据实际采样(Areal[n])确定理想复数采样(Aideal[n])(S101),
根据实际复数采样(Areal[n])和理想复数采样(Aideal[n]),形成复数商(ΔA1[n]=Areal[n]/Aideal[n])(S102),
通过将复数商的相位设为0,产生经修正复数商(B[n])(S106),以及
对经修正复数商(B[n])进行傅立叶变换(S107)。
3、根据权利要求1或2所述的方法,
其特征在于,
输入信号(S)是根据mVSB方法、特别是8VSB方法进行了数字调制的。
4、根据权利要求3所述的方法,
其特征在于,
仅理想采样(Aideal[n])的同相分量(I)是根据实际采样(Areal[n])的同相分量(I)确定的,并且
理想采样(Areal[n])的正交相位分量(Q)是通过对理想采样(Aideal[n])的同相分量(I)进行希尔伯特变换而产生的。
5、根据权利要求1至4中的任何一项所述的方法,
其特征在于,
当相关联的实际采样的值(|Areal[n]|)小于第一阈值(Minvalue)时,以插值(ΔA2[n])代替复数商(ΔA1[n])。
6、根据权利要求1至5中的任何一项所述的方法,
其特征在于,
当相关联的实际采样的虚部(Im{Areal[n]})大于第二阈值(Amax)时,以插值(ΔA2[n])代替复数商(ΔA1[n])。
7、根据权利要求1至6中的任何一项所述的方法,
其特征在于,
当相关联的实际采样的虚部(Im{Areal[n]})小于第三阈值(Amin)时,以插值(ΔA2[n])代替复数商(ΔA1[n])。
8、具有控制信号的数字存储介质,所述控制信号可以被电读出,并且可以与可编程计算机或数字信号处理器协作以实现根据权利要求1至7之一的方法。
9、具有程序代码装置的计算机程序产品,所述程序代码装置存储在机器可读的载体上,以便当在计算机上或在数字信号处理器上运行该程序时能够实现根据权利要求1至7之一的所有步骤。
10、计算机程序,其具有程序代码装置,以便当在计算机上或在数字信号处理器上运行该程序时能够实现根据权利要求1至7之一的所有步骤。
11、计算机程序,其具有程序代码装置,以便当该程序存储在机器可读的载体上时能够实现根据权利要求1至7之一的所有步骤。
CNB038223031A 2002-10-04 2003-09-17 确定数字调制信号的相位噪声谱和/或振幅噪声谱的方法 Expired - Fee Related CN100338901C (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
DE10246316.6 2002-10-04
DE10246316A DE10246316A1 (de) 2002-10-04 2002-10-04 Verfahren zum Bestimmen des Phasen- und/oder des Amplituden-Rauschspektrums eines digital modulierten Signals

Publications (2)

Publication Number Publication Date
CN1682482A true CN1682482A (zh) 2005-10-12
CN100338901C CN100338901C (zh) 2007-09-19

Family

ID=32010210

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB038223031A Expired - Fee Related CN100338901C (zh) 2002-10-04 2003-09-17 确定数字调制信号的相位噪声谱和/或振幅噪声谱的方法

Country Status (8)

Country Link
EP (1) EP1547293B1 (zh)
JP (1) JP4291777B2 (zh)
CN (1) CN100338901C (zh)
AT (1) ATE458323T1 (zh)
DE (2) DE10246316A1 (zh)
DK (1) DK1547293T3 (zh)
ES (1) ES2338225T3 (zh)
WO (1) WO2004034630A1 (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102013213657A1 (de) * 2013-07-12 2015-01-15 Rohde & Schwarz Gmbh & Co. Kg Verfahren und Vorrichtung zur Messung des Amplituden- und/oder Phasenrauschens eines Signals
RU2756024C1 (ru) * 2020-09-14 2021-09-24 Акционерное общество "ОКБ-Планета" АО "ОКБ-Планета" Цифровое устройство определения спектра принимаемых сигналов с высоким разрешением по частоте
CN116298513B (zh) * 2023-03-01 2024-03-12 扬州宇安电子科技有限公司 一种多扫描模式的宽带频谱监测系统

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4953186A (en) * 1988-02-19 1990-08-28 Silicon Systems, Inc. Phase jitter tracker
US5001724A (en) * 1989-01-13 1991-03-19 Hewlett-Packard Company Method and apparatus for measuring phase accuracy and amplitude profile of a continuous-phase-modulated signal
US5187719A (en) * 1989-01-13 1993-02-16 Hewlett-Packard Company Method and apparatus for measuring modulation accuracy
US6246717B1 (en) * 1998-11-03 2001-06-12 Tektronix, Inc. Measurement test set and method for in-service measurements of phase noise
US6366621B1 (en) * 1998-11-03 2002-04-02 Tektronix, Inc. Method of estimating pilot signal phase in a digitally modulated RF signal
US6275523B1 (en) * 1998-11-03 2001-08-14 Textronic, Inc. In-service measurement of transmitter nonlinearities
JP3388346B2 (ja) * 1999-05-10 2003-03-17 アンリツ株式会社 ディジタル変調信号測定装置
US6249180B1 (en) * 1999-09-08 2001-06-19 Atmel Corporation Phase noise and additive noise estimation in a QAM demodulator
US20010037189A1 (en) * 2000-01-20 2001-11-01 Dan Onu Method of estimating phase noise spectral density and jitter in a periodic signal

Also Published As

Publication number Publication date
EP1547293B1 (de) 2010-02-17
EP1547293A1 (de) 2005-06-29
ES2338225T3 (es) 2010-05-05
DK1547293T3 (da) 2010-06-14
JP2006502633A (ja) 2006-01-19
WO2004034630A1 (de) 2004-04-22
ATE458323T1 (de) 2010-03-15
CN100338901C (zh) 2007-09-19
JP4291777B2 (ja) 2009-07-08
DE50312426D1 (de) 2010-04-01
DE10246316A1 (de) 2004-04-15

Similar Documents

Publication Publication Date Title
US6246717B1 (en) Measurement test set and method for in-service measurements of phase noise
JP3611862B2 (ja) デジタル直角変調及び復調方法、並びにデジタル直角変調器及び復調器
US8116363B2 (en) Method of determining true error vector magnitude in a wireless LAN
US7010278B2 (en) Sideband suppression method and apparatus for quadrature modulator using magnitude measurements
US6272172B1 (en) Measurement acquisition and display apparatus
JP6389355B2 (ja) 入力チャネルのs/n比評価方法及び試験測定装置
EP1518342B1 (en) Satellite travelling-wave tube amplifier on-line non-linearity measurement
US20090185037A1 (en) Measuring hum in a digital tv signal
EP1509015A1 (en) Quadrature modulator carrier quadrature error detection method and quadrature modulation device
CN1682482A (zh) 确定数字调制信号的相位噪声谱和/或振幅噪声谱的方法
US9350446B2 (en) Optical signal power selection and control
US20050175127A1 (en) Method for determining the phase-and/or amplitude-noise spectrum of a digitally modulated signal
JP3701851B2 (ja) ディジタル変調信号受信装置
US20050094746A1 (en) Characteristic measurement system for a digital modulation signal transmission circuit
JP2004120789A (ja) 妨害波検出装置及び妨害波検出方法
JP6358679B2 (ja) 剰余数系アナログ−デジタル変換を使用した位相角測定
Nolan et al. Modulation scheme recognition techniques for software radio on a general purpose processor platform
EP0762700B1 (en) Method and device for correcting a phase and an amplitude error of a carrier signal
EP1209839B1 (en) System and method for automatic detection of impairments in a digital quadrature amplitude modulated signal
US20060155493A1 (en) Method for determining the envelope curve of a modulated signal
JPH07245634A (ja) 位相被変調信号の特性を記述する方法とシステム
US7054383B2 (en) Method for the determination of the average long-term power of digital modulated signal
CN116319234A (zh) Iq不平衡校正方法、装置及电子设备
Law et al. Identifying RF-related impairments in full-service digital networks
CN117949941A (zh) 雷达信号的接收侧去失真的方法、存储介质、电子控制设备和雷达传感器

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20070919

Termination date: 20210917