首页 / 专利库 / 电信 / 波峰因数 / 生成用于通信系统之中或之上具有预定质量标准的信号的方法和装置

生成用于通信系统之中或之上具有预定质量标准的信号的方法和装置

阅读:218发布:2020-09-14

专利汇可以提供生成用于通信系统之中或之上具有预定质量标准的信号的方法和装置专利检索,专利查询,专利分析的服务。并且这里公开了一种具有预定 质量 标准并用于通信系统的 信号 ,一种生成这种信号的方法和系统,一种利用这种信号测试通信系统操作的方法以及配置为运用这种方法的(电信)通信系统。用于生成具有预定质量标准的方法包括以下步骤:表示第一信号(10),所述第一信号具有多个 频率 分量,每个频率分量具有 频谱 幅度性质和 相位 性质,以及-通过配置(11)信号的频谱幅度性质来处理所表示信号,以及-通过根据预定质量标准配置(11)信号的频谱幅度性质来处理所表示信号。,下面是生成用于通信系统之中或之上具有预定质量标准的信号的方法和装置专利的具体信息内容。

1.一种配置适合用于通信系统之中或之上具有至少一种预定质 量标准的信号的方法,所述方法包括以下步骤:
-表示第一信号,所述第一信号包括多个频率分量,所述每个 频率分量具有频谱幅度特性和相位特性,以及
-通过根据所述或每种预定质量标准配置所述频谱幅度特性并 配置随机相位特性来处理所述表示的第一信号。
2.如权利要求1所述的方法,其特征在于所述第一信号由指定 所述各频率分量的频谱幅度和相位的第一组数字表示。
3.如权利要求1所述的方法,其特征在于所述第一信号由具有 实部虚部的第二组复数表示,所述实部和虚部的组合指定所述各 频率分量的频谱幅度和相位。
4.如权利要求1所述的方法,其特征在于所述第一信号由第三 组数字表示,其中的每个数字在时域中指定所述第一信号的幅度。
5.如权利要求4所述的方法,其特征在于还包括以下步骤:将 所述第三组数字从时域变换到频域,从而由第四组数字表示所述第 一信号,所述第四组数字指定所述各频率分量的频谱幅度和相位。
6.如权利要求4所述的方法,其特征在于还包括以下步骤:将 所述第三组数字从时域变换到频域,从而由具有实部和虚部的第五 组复数表示所述第一信号,所述实部和虚部的组合指定所述各频率 分量的频谱幅度和相位。
7.如以上权利要求中的任何一项所述的方法,其特征在于还包 括以下步骤:对所述已处理的所表示第一信号进行后处理,从而实 现与所述或各种预定质量标准的进一步匹配。
8.如以上权利要求中的任何一项所述的方法,其特征在于还包 括以下步骤:将所述已处理的所表示信号从频域变换到时域。
9.如权利要求8所述的方法,其特征在于具有所述预定质量标 准的所述信号在时域中由第六组数字表示。
10.如以上权利要求中的任何一项所述的方法,其特征在于所 述或各种预定质量标准包括至少一个已调载波,所述或各个已调载 波包括由载频、载波幅度、调制深度和调制宽度组成的组中的任一 些。
11.如权利要求1-9中的任何一项所述的方法,其特征在于所 述或各种预定质量标准包括由预定时域幅度分布和预定频谱幅度包 络组成的组中的任一些。
12.如权利要求11所述的方法,其特征在于还包括以下步骤: 根据预定时域幅度分布配置所述已处理的所表示第一信号。
13.如权利要求11或12所述的方法,其特征在于还包括根据 预定频谱幅度包络配置所述已处理的所表示第一信号。
14.如权利要求12或13所述的方法,其特征在于由迭代过程 进行所述时域幅度分布和频谱幅度包络中的至少一个。
15.如权利要求14所述的方法,其特征在于所述迭代过程包括 将所述已处理的所表示第一信号的所述时域幅度分布和频谱幅度包 络中的任一些与预定时域幅度分布和预定频谱幅度包络进行比较。
16.如以上权利要求中的任何一项所述的方法,其特征在于通 过对根据权利要求10所处理的以及根据权利要求11-15中的任何 一项所处理的多个所表示信号进行组合而提供具有所述或各种预定 质量标准的所述信号。
17.如以上权利要求中的任何一项所述的方法,其特征在于具 有所述或各种预定质量标准的所述信号是噪声信号。
18.如以上权利要求中的任何一项所述的方法,其特征在于由 一组指令提供具有所述或各种预定质量标准的所述信号,其中所述 指令组为代码格式并可按预定顺序在处理装置上执行。
19.一组指令,所述指令组为代码格式并可按预定顺序在处理 装置上执行,所述指令组被配置为用于从依照根据权利要求18的方 法的第一信号表示生成具有所述或各种预定质量标准和随机相位特 性的信号。
20.一种装置,它包括处理单元、存储单元以及任意波形发生 器单元,所述装置被配置为根据以上权利要求中的任何一项生成具 有至少一种预定信号质量标准和随机相位特性的信号。
21.一种信号,它具有至少一种预定信号质量标准和随机相位 特性,并根据以上权利要求中的任何一项生成。
22.一种数据载体装置,它包括用于以上权利要求中的任何一 项的信号表示库。
23.一种测试通信系统的操作的方法,所述方法包括以下步骤:
-根据以上权利要求中的任何一项生成具有至少一种预定质量 标准的信号,以及
-通过所述通信系统传送所述信号。
24.一种系统,它包括用于根据以上权利要求中的任何一项生 成具有至少一种预定信号质量标准的信号的单元、调制解调器单元、 电缆单元以及处理器单元,其中为了进行自动测量和/或监控,所述 处理器单元被配置为控制所述生成单元、调制解调器单元以及电缆 单元。
25.一种电信系统,它被配置为运用权利要求24所述的方法。
26.一种配置适合用于通信系统之中或之上的信号的方法,所 述方法包括以下步骤:
-在时域中表示具有时域幅度分布的第一信号,所述信号在频 域中具有频谱密度,从而得到所表示信号;
-按照非线性变换处理所述所表示信号,所述非线性变换得到 至少一种预定质量标准,
-至少用预定时域幅度分布的反函数处理所述所表示信号的所 述时域幅度分布。
27.如权利要求26所述的方法,其特征在于还包括以下步骤: 将所述所表示信号的所述时域幅度分布与所述预定时域幅度分布进 行比较,并在其后配置所述非线性变换以便得到已处理的所表示信 号,所述已处理的所表示信号具有接近所述预定时域幅度分布的时 域幅度分布。
28.如权利要求26或27所述的方法,其特征在于所述已处理 的所表示信号g(t)是所述所表示第一信号f(t)的函数Q{f(t)},并且其 中所述函数Q被定义为:
Q(x)=sign(x)·G-1(F/x/))
其中sign(x)=x/|x|其中x<>0;sign(x)=0其中x=0;
F是所述所表示信号的所述时域幅度分布;而G是所述预定时 域幅度分布函数。
29.如权利要求28所述的方法,其特征在于还包括以下步骤: 根据预定频谱密度质量标准得到具有频谱密度的所述所表示信号。
30.如权利要求26所述的方法,其特征在于还包括以下步骤:
-将所述第一信号变换到频域;以及
-在所述频域中将所述第一信号乘以谱包络以得到相乘信号; 以及
-将所述相乘信号变换到时域。
31.如权利要求30所述的方法,其特征在于还包括以下步骤: 将所述所表示信号的所述时域幅度分布与所述预定时域幅度分布进 行比较,并在其后配置所述非线性变换以便得到已处理的所表示信 号,所述已处理的所表示信号具有接近所述预定时域幅度分布的时 域幅度分布。
32.如权利要求31所述的方法,其特征在于还包括以下步骤: 根据预定频谱密度质量标准得到具有频谱密度的所述所表示相乘信 号。
33.如权利要求32所述的方法,其特征在于所述步骤中的至少 两个步骤是迭代地执行的。
34.如权利要求32所述的方法,其特征在于所述步骤中的至少 两个步骤是迭代地执行的,直到得到预定波峰因数
35.如以上权利要求中的任何一项所述的方法,其特征在于所 述信号是噪声信号。
36.如权利要求35所述的方法,其特征在于将在频域中的表示 中的所述第一信号生成为一组随机数,所述随机数最好是复数,所 述复数的模表征幅度,所述复数的幅表征相位。
37.如权利要求36所述的方法,其特征在于基本上每个所述复 数的实部和/或虚部都是根据高斯分布选择的。
38.如权利要求36所述的方法,其特征在于基本上每个所述复 数的所述模基本上都等于所述预定频谱密度的所述幅度。
39.如权利要求36所述的方法,其特征在于基本上每个所述复 数的所述幅角都是随机的。
40.一种配置适合用于通信系统之上或之中的信号的方法,所 述方法包括以下步骤:
-在时域中表示具有时域幅度分布的第一信号,所述信号在频 域中具有频谱密度,从而得到所表示信号;以及
-处理所述所表示信号,直到得到具有根据预定频谱密度质量 标准的频谱密度的信号。
41.一种配置适合用于通信系统之上或之中的信号的方法,所 述方法包括以下步骤:
-在时域中表示具有幅度分布的第一信号,并且所述信号在频 域中具有频谱密度,从而得到所表示信号;以及
-对频域中的所述所表示信号进行滤波,包括在频域中估计至 少部分所述所表示信号,并在其后在频域中处理所表示信号的步骤。
42.如权利要求40或41所述的方法,其特征在于所述处理步 骤包括迭代处理。
43.如权利要求40或41所述的方法,其特征在于所述处理步 骤包括迭代处理,直到得到预定波峰因数。
44.一种信号,它包括至少一个随机噪声信号,所述随机噪声 信号在时域中具有根据预定质量标准的幅度分布并在频域中具有根 据预定质量标准的频谱密度,所述随机信号由随机数阵列构成。
45.如权利要求44所述的信号,其特征在于还包括离散频谱。
46.如权利要求44所述的信号,其特征在于利用一组指令生成 所述噪声信号,所述指令组为代码格式并按预定顺序执行。
47.一种生成随机信号的方法,所述随机信号包括至少一个随 机噪声信号,所述随机信号在时域中具有根据预定质量标准的幅度 分布并在频域中具有根据预定质量标准的频谱密度,所述随机信号 由随机数阵列构成,所述方法包括以下步骤:利用一组指令生成随 机数组,所述指令组为代码格式并按预定顺序执行。
48.如权利要求47所述的方法,其特征在于还包括以下步骤: 生成离散频谱,所述离散频谱利用角度测定函数并用噪声特征调制 基本上每个所述离散频率。
49.如权利要求48所述的信号,其特征在于还包括以下步骤: 利用一组指令组合所述随机噪声信号和所述离散频谱,所述指令组 为代码格式并按预定顺序执行。
50.一组指令,所述指令组为代码格式并可按预定顺序执行, 所述指令组被配置为生成随机噪声信号和离散频谱,所述随机信号 在时域中具有根据预定质量标准的幅度分布并在频域中具有根据预 定质量标准的频谱密度。
51.一种用于测试通信系统的操作的测试系统,所述测试系统 包括一组指令,所述指令组为代码格式并可按预定顺序在装置上执 行和编译,所述指令组被配置为生成噪声信号,所述噪声信号包括 随机噪声信号和离散频谱中的至少一种,所述随机信号在时域中具 有根据预定质量标准的幅度分布并在频域中具有根据预定质量标准 的频谱密度。
52.一种测试具有调制解调器的通信系统的操作的方法,所述 方法包括以下步骤:在由所述调制解调器收发的信号上叠加一个信 号,所述一个信号包括随机噪声信号和离散频谱中的至少一种,所 述随机信号在时域中具有根据预定质量标准的幅度分布并在频域中 具有根据预定质量标准的频谱密度,另外,所述噪声信号由随机数 阵列构成。
53.一种测试具有调制解调器的通信系统的操作质量的方法, 所述方法包括以下步骤:
-在由所述调制解调器收发的信号上叠加一个信号,所述一个 信号包括随机噪声信号和离散频谱中的至少一种,所述随机信号在 时域中具有根据预定质量标准的幅度分布并在频域中具有根据预定 质量标准的频谱密度,另外,所述噪声信号由随机数阵列构成;以 及
-根据预定质量标准估计所述收发信号。
54.一种改进通信系统的设计和/生产的方法,所述方法包括以 下步骤:
-在由调制解调器收发的信号上叠加一个信号,所述一个信号 包括随机噪声信号和离散频谱中的至少一种,所述随机信号在时域 中具有根据预定质量标准的幅度分布并在频域中具有根据预定质量 标准的频谱密度,另外,所述噪声信号由随机数阵列构成;
-根据预定质量标准估计所述收发信号;以及
-迭代地配置所述调制解调器的所述设计,从而更加接近用于 估计所述收发信号的所述质量标准。
55.一种包括这样一种信号的电信网络,所述一种信号包括随 机噪声信号和离散频谱中的至少一种,所述随机信号在时域中具有 根据预定质量标准的幅度分布并在频域中具有根据预定质量标准的 频谱密度,另外,所述噪声信号由随机数阵列构成。

说明书全文

技术领域

发明一般涉及通信系统,更具体地说,涉及用于通信系统的 信号,用于生成这种信号的方法和系统,利用这种信号测试通信系 统的操作的方法,测试系统以及被配置为运用这种方法的(电信) 通信系统。

发明背景

为了测试诸如xDSL收发信机以及电缆或网络的通信系统和通信 设备,特别需要这样的测试信号,它们以代表每条电缆上具有大量 的系统或系统装置的实际配置情况的方式加载到通信系统和通信装 置。
通过在实际(有噪声)测试条件下测量系统或系统装置的传输 性能,人们可以改进系统或者装置的设计,和/或证明它们的性能符 合诸如ETSI、ITU或ANSI或其它(电信)通信团体发布的标准。
一种执行这种性能测试的方法是生成被称为减损(impairment) 的信号。更具体地说,可以将减损细分为:
(i)串音噪声,它具有由例如来自相邻xDSL系统的谱包络和频 谱幅度分布表征的噪声分布;
(ii)干扰噪声,它由离散频率分量组成,也被称为射频干扰音 (rfi-tone),具有由起源于例如无线电或业余广播的多个离散频 率分量以及频谱幅度、调制深度和调制宽度参数表征的噪声分布, 以及
(iii)脉冲噪声,由例如开关操作和开关元件引起的信号脉冲表 征。
对于干扰噪声来说,频率可以随时间变化(摆动)。
一种用于生成减损的装置被称为减损发生器,它被配置为生成 所述串音噪声和干扰噪声中的至少一种,特别用于通信系统之中或 之上。
在实践中,为了测试通信系统和通信装置是否符合标准,已定 义了各种噪声分布,其中这些噪声分布根据如诸如通信电缆中线对 的长度和数量以及传输数据速率系统参数而变化。
另外,诸如电缆、电缆或光纤或其它类型的电缆等各种不同 类型或长度的传输媒体要求不同的噪声信号。
用于生成噪声分布的方法和装置在本领域中是已知的。具体地 说,众所周知,滤波技术和滤波器被用于从输入信号生成噪声,提 供具有特定谱包络和频谱幅度分布的输出信号
但是,通过利用滤波技术和滤波器,在输入信号和输出信号间 建立了因果关系。本领域的技术人员将意识到,这种类型的信号不 太适合现实模拟实际操作的通信系统和通信装置。
WO 00/16181公开了一种用于生成近似于预定幅度分布图的随机 时域信号的方法和装置。在第一步骤中,通过对诸如白噪声信号的 噪声信号进行滤波来创建信号,进而产生具有预定谱包络的信号。 在下一步骤中,将非线性函数应用于这种滤波后的噪声信号,以产 生近似于预定幅度分布图的所需时域信号。再下一步,将脉冲响应 滤波应用于时域信号,以校正其谱包络并得到具有所需谱包络的输 出信号。非线性函数和脉冲响应滤波函数是根据提供的谱包络而选 择的特殊函数。
WO 00/16181限于这样一种情况:它仅提供一种只具有预定谱包 络的时域信号。WO 00/16181未对要提供的时域信号所具有的其它质 量标准作出说明,特别未对相位特性作出说明。
发明概述
本发明的目的是提供一种改进信号,它被用于通信系统和通信 装置,特别是用于根据预定(标准化)噪声分布测试这些系统和装 置。
在本发明的第一个方面中,公开了一种配置适合用在通信系统 之中或之上具有预定质量标准的信号的方法,所述方法包括以下步 骤:
-表示第一信号,所述第一信号包括多个频率分量,每个分量 具有频谱幅度特性以及相位特性,以及
-通过根据所述或各种预定质量标准配置频谱幅度特性以及配 置随机相位特性来处理所表示的第一信号。
修改谱包络的传统方法是利用数字滤波器组。对本发明来说, 所表示的第一信号和要提供的信号之间不必建立因果关系,所以这 种传统方法决非理想。这种根据本发明对问题的理解大大简化了频 率整形的完成。
例如,从具有随机相位特性的第一信号开始,所表示的第一信 号的频率整形操作足以提供满足预定质量标准的信号。可以以几种 方式执行根据本发明的频率整形。
在本发明的实施例中,第一信号由指定各频率分量的频谱幅度 和相位的第一组数字表示。在保持信号的随机相位特性的同时,各 频率分量的频谱幅度的调整足以实现频域中所表示信号的频率整 形。
在本发明的另一个实施例中,第一信号由具有实部虚部的第 二组复数表示,其中实部和虚部的组合指定各频率分量的频谱幅度 和相位。通过充分地调整这些复数来实现频率整形,但是,在对所 表示的第一信号进行调整之后要保持随机相位特性。
在本发明的又一个实施例中,第一信号由第三组数字表示,各 数字在时域中指定第一信号的幅度。通过利用例如快速傅里叶变换 (FFT)算法将第三组数字从时域变换到频域,第一信号由指定各频 率分量的频谱幅度和相位的第四组数字表示。可以由频率整形操作 对第四组数字作进一步处理,如以上针对第一组数字所公开的那样。
但是根据本发明,也可以将第三组数字从时域变换到频域,从 而由第五组具有实部和虚部的复数表示第一信号。如上所述,为了 频率整形,必须充分地调整第五组复数。
在所表示的第一信号具有非随机相位特性的情况下,通过适当 地配置第二、第四和第五组数字,可以得到随机相位特性。
利用实数或复数的比例因子,可以由乘法运算调用频域中的缩 放。将频率分量的频谱幅度的期望值除以其实际频谱幅度得到频率 分量的频谱幅度的乘法比例因子。
依照根据本发明的方法的另一个实施例,为了达到与所述或各 种预定质量标准的进一步匹配,提供了已处理的所表示第一信号的 后处理。
但是,为了用于根据本发明的通信系统之中或之上,必须利用 例如反快速傅里叶变换(IFFT)算法,将在频域中配置的所表示第 一信号变换到时域。
另外,以上公开的处理步骤还可以包括诸如信号的卷积或去卷 积或者乘法或加法运算。在时域中,尤其可以由时域中第六组数字 表示满足所述或各种预定质量标准的已处理的所表示第一信号。
但是,通过以上方法,所提供的满足频率中诸如预定频谱幅度 包络以及随机相位特性的质量标准的信号可能还不满足时域中诸如 预定时域幅度分布的质量标准。
在根据本发明的方法的再一个实施例中,所述或各种预定质量 标准包括包含预定时域幅度分布以及预定频谱幅度包络的组中的任 何一些。
因此,在根据本发明的方法的另一个实施例中,根据预定时域 幅度分布配置已处理的所表示第一信号。
在根据本发明的方法的再一个实施例中,根据预定频谱幅度包 络配置已处理的所表示第一信号。
根据本发明,为了提供准确满足频域和时域中的预定质量标准 的信号,时域幅度分布和频谱幅度包络中的至少一种是由迭代过程 完成的。可以按照所需次数重复幅度和频率整形,直到两个形状都 在合理精确度之内满足要求。
在本发明的实施例中,迭代过程包括在任一迭代步骤后,将已 处理的所表示第一信号的时域幅度分布和频谱幅度包络中的任一些 与预定时域幅度分布与预定频谱幅度包络进行比较。
已经发现,不需要在频率整形后执行全时域特征检查,以弄清 时域特征是否足够接近要求。已证明,在实践中,波峰因数要求的 简单检查足以能够决定是停止还是继续迭代。信号的波峰因数被定 义为信号音平均或rms(均方根)值与信号音最大或峰值幅度的关系。
以上公开的根据本发明的方法特别适合生成串音噪声等。
在根据本发明的方法的第二方面中,如果要生成具有干扰噪声 特征的信号,所述或各种预定质量标准包括至少一个已调载波,所 述已调载波包括由载频、载波幅度、调制深度以及调制宽度组成的 组中的任一些。
通过根据一种质量标准或上述质量标准对所表示第一信号进行 整形,可以容易并且非常有效地提供表示特定类型的干扰噪声、具 有特定时域幅度分布以及预定频谱幅度包络的信号。
根据本发明的方法,可以通过组合多个如上文所公开的内容处 理的信号来提供满足所述或各种预定质量标准的信号。
为了在例如通信网络或通信装置的测试中利用这种信号,必须 利用例如FFT算法等将已处理的所表示信号从频域变换到时域。
本发明还用于组合根据以上公开的第一和第二方面生成的信 号。但也可以包括其它信号分量。
具体地说,根据本发明的方法,具有所述或各种预定质量标准 的信号是噪声信号。
在本发明的第三方面,公开了一种测试通信系统的操作的方法, 其中所述方法包括以下步骤:
-根据以上公开的本发明的方法,生成具有预定质量标准的信 号,以及
-通过测试中的通信系统传送这种信号。
可以利用一组为代码格式并按预定顺序在装置上执行的指令来 生成和存储所述信号。该组指令可以是软件代码,它在计算机上编 译并存储在计算机或计算机网络软盘或CD-ROM上,或通过因特 网存储。软件和/或产生的信号还可以被存储在任意波形发生器 (AWG)卡上,并且AWG可被用于生成这种信号或者从存储器再 现存储信号。因此有可能存在一种可用信号库,其中这些信号存储 在数据载体上,可被用于执行或使用根据本发明的方法。
通信系统可以是诸如xDSL调制解调器、或这种调制解调器内的 或用于这种调制解调器的芯片、或网络中的电缆、或用于(电信) 通信的网络的装置。
在本发明的第四方面中,公开了配置用于通信系统之上或之中 的信号的另一种方法。信号最好是噪声信号。信号可包括串音噪声, 它是在频域和时域中具有预定特性的随机信号。信号还可以包括射 频干扰音,它具有离散频谱。这种信号中也可以包括其它信号分量。
所述方法包括以下步骤:
-在时域中表示具有时域幅度分布的第一信号,所述信号在频 域中具有频谱密度,从而得到所表示信号;
-根据非线性变换处理所表示信号,所述非线性变换得到至少 一种预定质量标准;
-至少通过预定时域幅度分布的反函数处理所表示信号的时域 幅度分布。
这种方法还可以包括以下步骤:将所表示信号的时域幅度分布 与预定时域幅度分布进行比较,并在其后配置非线性变换从而得到 已处理的所表示信号,该信号具有近似于预定时域幅度分布的时域 幅度分布。
在本发明的第五方面中,公开了一种方法,这种方法还包括以 下步骤:将所表示信号的时域幅度分布与预定时域幅度分布进行比 较,并在其后配置非线性变换从而得到已处理的所表示信号,该信 号具有近似于预定时域幅度分布的时域幅度分布。
根据本发明的第五方面,所述方法还可以包括以下步骤:在时 域中表示具有幅度分布的第一信号,并且该信号在频域中具有频谱 密度,从而得到所表示信号,以及在频域中对所表示信号进行滤波 包括以下步骤:在频域中估计至少部分信号表示,并在其后在频域 中处理所表示信号。
可以组合本发明方法的第四和第五方面。本发明方法的第四和 第五方面允许在不同的迭代步骤产生信号,其中这种信号具有预定 幅度分布和/或具有预定频谱密度或具有根据预定质量标准的幅度分 布和/或具有根据预定质量标准的频谱密度。预定质量标准可以是信 号的波峰因数,它是信号音的最大值或峰值与信号音的平均值或rms 值的关系。上文所列举的处理步骤可以包括以下步骤:快速傅里叶 变换(FFT)或反快速傅里叶变换(IFFT)。处理步骤还可以包括诸 如信号的卷积或反卷积或者乘法或加法的运算。
在本方法的第四方面中,对所表示信号的幅度分布进行处理包 括预定幅度分布的函数,该函数可包括预定幅度分布的反函数。
本发明的第四和第五方面所列举的方法还可以包括以下步骤: 在频域中变换第一信号;将频域中的第一信号乘以谱包络从而得到 相乘信号;以及之后在时域中表示相乘信号。
在这些方法中,可将第一信号在频域中的表示生成为一组随机 数、最好是复数,其中复数的模表示幅度,复数的幅表示相位, 以及可以根据高斯分布选择基本上每个复数的实部和/或虚部。各个 复数可以基本上等于预定频谱密度的幅度。
在本发明的第六个方面中,公开了一种信号,它包括至少一个 随机噪声信号,这种随机信号具有根据预定质量标准的时域幅度分 布以及根据预定质量标准的频域频谱密度,随机信号由随机数阵列 组成。信号还可以包括离散频谱。可以利用一组为代码格式并按预 定顺序执行的指令来生成噪声信号。该组指令可以是在计算机上编 译并存储在计算机或计算机网络或软盘或CD-ROM上的、或通过 因特网存储的软件代码。还可以将这种软件存储在任意波形发生器 (AWG)卡上,并且AWG可被用于生成这种信号或从存储器再现 存储信号。因此,有可能存在一种可用信号库,其中这些信号可被 用于执行或使用本发明方法的第四和第五方面。
在本发明的第七方面中,公开了一种生成包括至少一个随机噪 声信号的信号,这种随机信号在时域中具有根据预定质量标准的幅 度分布并在频域中具有根据预定质量标准的频谱密度。随机信号由 随机数阵列构成,所述方法包括以下步骤:利用一组为代码格式并 按预定顺序执行的指令来生成随机数组。所述方法还可以包括以下 步骤:生成离散频谱,所述离散频谱利用角度测定函数并用噪声特 性对基本上每个离散频率进行调制。可以利用一组为代码格式并按 预定顺序执行的指令来组合随机噪声信号和离散频谱。
在本发明的第八个方面中,公开了一组为代码格式并可按预定 顺序执行的指令,这组指令被配置为生成随机噪声信号和离散频谱, 这种随机信号在时域中具有根据预定质量标准的幅度分布并在频域 中具有根据预定质量标准的频谱密度。这组指令可以是在计算机上 编译并存储在计算机或计算机网络或软盘或CD-ROM上或通过因 特网存储的软件代码。软件还可以被存储在任意波形发生器(AWG) 卡上,并且这种AWG可被用于生成信号或从存储器再现存储信号。 因此有可能存在一种可被利用的可用信号库。软件可以是C语言代 码或可在MATLAB环境下编译。
在本发明的第九方面中,公开了一种用于测试通信系统的操作 的系统,它包括一组为代码格式的指令,该组指令可按预定顺序在 装置上执行和编译,这组指令被配置为生成包括随机噪声信号和离 散频谱中的至少一种的噪声信号,随机信号在时域中具有根据预定 质量标准的幅度分布并在频域中具有根据预定质量标准的频谱密 度。根据本发明的这个方面的测试系统可以包括用于生成噪声信号 的减损发生器。
将减损发生器连接到被测试的通信系统的连接部件(变压器、 有源器件、衰减器等)可以具有不需要的频率从属响应。例如可以 通过在减损发生器中生成特定测试信号来测量这种不需要的频率从 属响应。可以通过乘以信号的期望频谱密度除以连接部件的不需要 频率从属响应来补偿这种不需要的频率从属响应。
在本发明的第十个方面中,公开了一种测试诸如xDSL调制解调 器的通信系统的操作的方法,这种方法包括以下步骤:在经调制解 调器收发的信号上叠加一个信号,其中所述一个信号包括随机噪声 信号和离散频谱中的至少一种,所述随机信号在时域中具有根据预 定质量标准的幅度分布并在频域中具有根据预定质量标准的频谱密 度,此外,所述噪声信号由随机数阵列构成。
在本发明的第十一方面中,公开了一种测试通信系统的操作质 量的方法。这种方法包括以下步骤:在经调制解调器收发的信号上 叠加一个信号,其中所述一个信号包括随机噪声信号和离散频谱中 的至少一种,所述随机信号在时域中具有根据预定质量标准的幅度 分布并在频域中具有根据预定质量标准的频谱密度,此外,所述噪 声信号由随机数阵列构成,并根据预定质量标准估计所收发的信号。
另外,在本发明的第十二方面中,公开了一种改进通信系统的 设计和/或生产的方法,这种方法包括以下步骤:在经调制解调器收 发的信号上叠加一个信号,其中所述一个信号包括随机噪声信号和 离散频谱中的至少一种,所述随机信号在时域中具有根据预定质量 标准的幅度分布并在频域中具有根据预定质量标准的频谱密度,此 外,所述噪声信号由随机数阵列构成;根据预定质量标准估计所收 发的信号;以及迭代地配置调制解调器的设计,以便得到更接近用 于估计所收发信号的质量标准。
在本发明的第十三方面中,公开了包括这样一种信号的电信网 络,这种信号包括随机噪声信号和离散频谱中的至少一种,这种随 机信号在时域中具有根据预定质量标准的幅度分布并在频域中具有 根据预定质量标准的频谱密度,另外,所述噪声信号由随机数阵列 构成。
可以组合上述各方面的特性和本发明的实施例。
上述信号、方法和指令组将提供在诸如电话电缆或无线媒体等 媒体上的更高质量信号传输。更好的信号传输允许为通信系统的用 户在更大范围提供更多服务。
附图说明
图1示意了利用根据本发明的方法操作的减损发生器在通信系 统中进行性能测试的设置方框图
图2示意了根据本发明方法的流程图类型的实施例。
图3在图示中示意了根据本发明用于生成干扰噪声的方法的实 施例。
图4示意了根据本发明方法的另一些流程图类型的实施例。
图5示意了放大信号中高幅度峰值或信号音的幅度失真(非线 性变换)函数Q(x)。
图6示意了本发明的例示实施例的流程图。
图7-11示意了根据本发明的实施例得到的结果。
实施例详述
为了说明本发明,以下对本发明的信号和方法以及系统的各方 面和实施例进行描述。本领域的技术人员将理解,可以设想并实现 本发明的其它备选和等效实施例,并不背离本发明的真实精神。本 发明的范围仅由所附权利说明书进行限制。
在本发明的实施例中,公开了一种用于测试诸如xDSL收发信机 的通信系统的操作的系统。对运行存取网络中的高穿透性系统情况 的测试设备的设置进行描述。
公开了一种配置用于通信系统之中或之上的信号的方法。
传输性能测试的目的是以这样一种方式加载xDSL收发信机,这 种方式表示运行存取网络中的高穿透性系统情况。这种高穿透性方 法允许:
(i)组件和系统设计者量化性能,并用此改进他们的设计,从而证 明符合标准;以及
(ii)操作员定义应用于大多数运行情况的配置规则。
图1示意了可能的测试设置1的功能描述。它包括:
●测试回路2,它是真电缆或电缆模拟器
●添加部件3,它向测试回路2注入减损噪声;
●高抗阻,以及非常均衡的差动电压探针4,它连接到诸如频
谱分析仪或rms电压表(未示出)等电平检测器5,以及
●测试中的xDSL收发信机(调制解调器)6、7。 当测试中的xDSL系统需要外部分流器时(用于POTS或ISDN信号), 可以将这些分流器包括在测试中的调制解调器6、7中。
流经测试设备设置1的信号是从端口Tx到端口Rx,这意味着 测量上游和下游性能需要交换收发信机位置和测试“电缆”端。当 端口Tx以及端口Rx的末端是测试中的xDSL收发信机(调制解调 器)6、7时,端口Rx的接收信号电平是节点A2和B2之间测量的 电平。在这种测量期间断开减损发生器8。在相同情况下,端口Tx 的发送信号电平是节点A1和B1之间测量的电平。
减损发生器8应注入到测试设置1中的噪声与频率有关。减损 发生器8注入到测试设置1中的噪声应该是实际(频谱受影响)存 取网络的实际表示,并且它是:
(a)取决于测试回路2的长度,以及
(b)在下游性能测试和上游性能测试中不同。 在节点A2和B2之间测量的这种减损噪声通常是随机噪声、冲击噪 声和谐波噪声(射频干扰音)的混合。一组特性被标识为“噪声分 布”。
用相当均衡的差动电压探针4测量信号和噪声电平。
在全自动测试设置1中,测试回路2,3、电压探针4以及电平 检测器5、测试中的调制解调器6,7以及减损发生器8可以连接到 中央处理单元(CPU)9,如图中的虚线所示。本领域的技术人员将 理解,与CPU9的连接可涉及CPU9进行的远程测试的数据链路。
以下是与测试设备的使用有关的说明: ●在满信号带上测量这种设置中的rms电压Urms[V],意味着P[dBm]   的功率电平等于P=10×log10(Urms2/Rv×1000)[dBm]; ●在小频带Δf(以赫兹(Hz)为单位)内测量这种设置中的rms电   压Urms[V],意味着该滤波带内的功率频谱密度电平P[dBm/Hz]等   于P=10×log10(Urms2/Rv×1000/Δf))[dBm/Hz]; ●带宽Δf标识滤波器的噪声带宽,并不是-3dB带宽。
图2示意性地显示了用于根据本发明配置信号Ui(t)i=1,2,3,… 以用于通信系统之上或之中的方法的实施例。信号可包括串音噪声, 它是在频域中和在时域中具有预定特性的随机信号。
如图2的流程[1]所示,这种方法包括以下步骤:表示第一信号, 所述第一信号包括多个频率分量,每个频率分量具有频谱幅度和相 位特性,框10“第一信号”,并通过根据至少一种预定质量标准配 置频谱幅度特性来处理所表示信号,以及配置随机相位特性,框11 “配置”,因此得到已处理的所表示信号。
第一信号可由第一组数字表示,这组数字指定各频率分量的频 谱幅度和相位。另外,第一信号可以由具有实部和虚部的第二组复 数表示,实部和虚部的组合指定各频率分量的频谱幅度和相位。也 就是说,复数的模表征频率分量的频谱幅度而复数的幅角表征频率 分量的相位。
根据本发明,对所表示第一信号10进行处理,以配置随机相位 特性。但是,从具有随机相位特性的所表示第一信号10开始,对于 根据预定质量标准对信号的频率进行整形,它足以对频率分量的频 谱幅度进行整形。
这种方法还可以包括以下步骤:在时域中表示第一信号,其中 第一信号由第三组数字表示,其中每个数字指定第一信号在时域中 的幅度。通过例如利用FFT算法将第三组数字从时域变换到频域, 得到指定各频率分量的频谱幅度和相位的第四组数字。同样地,通 过根据至少一种预定质量标准配置频谱幅度特性以及配置其随机相 位特性来处理第四组数字。
通过例如利用反FFT算法(IFFT)将已处理的所表示信号从频 域变换到时域,框12“变换”,最终生成具有至少一种预定质量标 准的信号U1(t)。
尽管信号U1(t)满足频域中的至少一种预定质量标准,比如谱包 络和/或预加重特性,还可能需要提供在时域中具有预定质量标准的 信号。
如图2中流程[2]和[3]所公开,即框13“幅度形状”和框14“频 率形状”,时域中的质量标准可包括预定幅度分布和/或预定频谱幅 度包络。
这种方法还可以包括以下步骤:在不同的迭代步骤中产生信号, 见图2[4],框15“测试形状”和反馈回路16。因此,信号可以具有 预定时域幅度分布和/或预定频谱幅度包络和/或根据预定质量标准的 频谱密度。
至少一种预定质量标准可以是信号的波峰因数,也就是信号音 的最大值或峰值与信号音的平均值或rms值的关系。
可以利用一组为代码格式并可按预定顺序在装置上执行并编译 的指令来生成并存储信号。这组指令可以是在计算机上编译并存储 在计算机或计算机网络或软盘或CD-ROM上、或通过因特网存储 的软件代码。这种软件还可以存储在任意波形发生器(AWG)卡上, 并且AWG可被用于生成信号或从存储器再现存储信号。因此,有可 能存在可用信号库,其中这些信号可被用于执行或使用本发明的方 法的第一和第二方面。通信系统可以是诸如xDSL调制解调器6、7 的装置,或这种调制解调器6、7内或用于这种调制解调器6、7的 芯片或电信网络。
进行图2的框13“幅度形状”中的处理,从而得到对时域特性 的影响或控制。选择幅度失真(变换)函数Q(x)来放大信号的高幅 度峰值或信号音。可以从信号的实际幅度分布函数和预定幅度分布 函数重构非线性变换函数Q(x)。
对于0到T之间时段t中的噪声信号f(t),信号的幅度分布F(a) 被定义为噪声f的绝对值大于a的时间部分。如果G(a)是预定幅 度分布(如增强高斯函数,见下文),并且G-1(a)是其反函数,则从 噪声信号f(t)产生中间或最终信号g(t)的变换函数Q(x)可被定义为: Q(x)=sign(x)·G-1(F(|x|))   (1) g(t)=Q{f(t)}                 (2) sign(x)=x/|x|其中x<>0;sign(x)=0其中x=0; 因此,g(t)将具有预定幅度分布G(a)。在许多情况下,Q(x)可以是解 析函数但也可以用数字构造。增强高斯函数的例示如下:
高斯类型噪声的幅度分布是: G ( x ) = 1 - erf ( x 2 σ ) . . . . . . . ( 3 ) 其中 erf ( x ) = 2 π 0 x dtexp ( - t 2 ) . . . . . . . ( 4 ) 并且其中σ是信号的RMS值。 “增强”高斯分布被定义为:
如果VRMS是噪声采样的期望RMS值,并且Cf是期望波峰因数, 选择:
A=Cf×VRMS,并且       (6) σ = ( ( 1 + a ) V 2 RMS - A 2 · a / 3 ) . . . . . . . . . . . ( 7 )
已证明有用的a的典型值是0.001和0.01之间的大小,并且它表 示增强高斯分布与实际高斯分布的偏差。
在图2的框14频率形状中,作为后处理步骤,改进了信号的频 域特性,从而得到与所述或各种质量标准的进一步匹配。例如,通 过用(中间)信号U2(t)的测量频谱密度与预定频谱密度进行比较 (除),可得到校正频率曲线。在下面对FFT函数的卷积的描述中 的最佳模式实施例中给出其例示。
图2的部分[4]中,示意了上文详细描述的步骤的迭代过程如何 可以导致对用于通信系统之中或之上的完成或最终信号的进一步改 进。执行迭代过程、即框15测试形状的频率形状测试以及反馈回路 16,直到得到预定质量标准。
图3在示意性图解中说明利用根据本发明的方法提供干扰噪声 信号。
如在前序中所公开的,干扰噪声可以由离散载频fci,i=1,2,3… 上的多个频率分量表征。载频fci上的频率分量各自具有载波幅度 Aci,i=1,2,3…,以及(如果可适用)具有调制宽度、即相关载频fci 的左手边和右手边的多个离散频率,以及具有调制深度,这种调制 深度是与相应载频fci有关的边频幅度。
图3在具有平或频率轴f和垂直或幅度轴A的图形表示中, 仅作为例示显示了由两个载频fc1和fc2组成的信号,这两个载频分 别具有载波幅度Ac1和Ac2。
载频fc1的两侧旁配置了三个边频分量,各边频具有幅度A1。 对于载频fc2两侧的频率分量,配置了两个边频分量,各边频分量具 有幅度A2。
根据本发明,为了提供具有至少一种预定质量标准的信号,必 须对频率分量的幅度进行整形,如在图3中用虚线1和2所示。
从具有随机相位特性的所表示第一信号开始,根据本发明,通 过对幅度进行整形,在要提供的具有预定质量标准的信号中保持随 机相位特性。
图4示意性地显示用于配置在通信系统之中或之上使用的信号 的方法的另一实施例,具体地说,这种方法在已经提供随机相位特 性时使用。该信号包括串音信号,它是在频域中和在时域中具有预 定特性的随机信号。该信号还可以包括具有离散频谱的射频干扰音。 另外,这种信号中可以包括其它信号分量。
这种方法还可以包括以下步骤:在时域中表示第一信号,所述 第一信号具有幅度分布,并且该信号在频域中具有频谱密度,进而 得到所表示信号;以及根据非线性变换处理所表示信号,所述非线 性变换实现预定质量标准。这在图4中被示为幅度整形,流程[2-4]。
这种方法还可以包括以下步骤:在时域中表示第一信号,所述 第一信号具有幅度分布,并且该信号在频域中具有频谱密度,进而 得到所表示信号;根据预定频谱密度质量标准处理所表示信号,直 到得到具有频谱密度的信号。这在图4中被示为频率整形,流程[2-4]。 频率整形步骤还可以包括以下步骤:在频域中对所表示信号进行滤 波,其中包括在频域中估计至少部分的信号表示,以及其后处理频 域中的信号表示的步骤。
这种方法还可以包括以下步骤:在不同的迭代步骤中产生信号, 见图4的流程[4]。因此,信号可以具有预定时域幅度分布和/或预定 频谱密度或根据预定质量标准的时域幅度分布和/或频谱密度。预定 质量标准可以是信号的波峰因数,也就是信号音的最大值或峰值与 信号音的平均值比较的关系。可以利用一组为代码格式并可按预定 顺序在装置上执行和编译的指令来生成并存储这种信号。
同样地,这组信号可以是在计算机上编译并存储在计算机或计 算机网络或软盘或CD-ROM上、或通过因特网存储的软件代码。 这种软件还可以被存储在任意波形发生器(AWG)卡上,并且AWG 可被用于生成信号或从存储器再现存储信号。因此,有可能存在可 用信号库,它可以用于执行或使用本发明的方法。通信系统可以是 诸如xDSL调制解调器、或这种调制解调器内或用于这种调制解调器 的芯片、或电信网络的装置。在图4中详细说明了以下实施例。
利用软件生成随机数,框16“创建噪声”。在硬件中可以生成 白噪声。对随机数进行过滤直到得到预定频谱密度。所生成的随机 数各自表示一个频率分量。通过调整复数的幅度来执行得到预定频 谱密度的必要处理,之后进行IFFT处理,从而产生期望噪声信号。 执行这种方法的另一种方式是生成表示各频率分量相位的随机数, 并在其后复数的幅度被配置为近似于或等于预定频谱密度。
为得到对时域特性的影响或控制,进行框13“幅度形状”的处 理。图5示意了选择幅度失真(变换)函数Q(x)来放大信号中的高 幅度峰值或信号音。如以上针对公式(1-7)所述的,可以从信号 的实际幅度分布函数和预定幅度分布函数重构非线性变换函数 Q(x)。
在图4的框14“频率形状”中,改进了信号的频域特性。通过 用(中间)信号的预定频谱密度与测量频谱密度进行比较(除)得 到校正频率曲线。在后续对FFT函数的卷积中描述的最佳模式实施 例中给出其例示。
在图4的流程[4]中,就象图2中流程[2]一样,再次示意了上文 详细描述的步骤的迭代过程如何可以导致用于通信系统之上或之中 的最终进一步改进。执行迭代过程,直到得到预定质量标准。
如上所述,通过根据本发明的方法,利用诸如减损发生器8的 装置,可以生成表示串音噪声和干扰噪声的信号,见图1,其中的减 损发生器可被配置为提供由串音噪声信号和干扰噪声信号构成的信 号,同时如果需要,可以在要提供的输出信号中加入其它信号分量。
在本发明的实施例中,可以有利地将要提供的信号作为时域中 的第六组数字提供,例如数字阵列。
图6在流程类型的图中示意了运行在个人计算机20上的本发明 的例示实施例。被称为SPOCS的框21生成减损噪声,它包括框22 “白噪声”、FFT提供的框23“频谱整形”、IFFT从框23的输出 产生的框24“期望噪声信号”,框24的结果信号被存储在AWG卡 25中。在串音情况中,即框26,建立噪声PSD(功率频谱密度), 框27,由框23对噪声PSD进一步处理。
下文公开了本发明的一组指令的最佳模式实施例。下文给出的 代码在MATLAB环境下编译。%符号后面给出了关于代码功能的注 释。对于本领域的技术人员来说,所提供的代码是不解自明的。
图7-9示意了通过最佳模式实施例得到的结果。图8示意了生 成噪声样本以及噪声分布的PSD的频谱曲线。图9示意了生成噪声 样本在时域中的曲线。图10示意了生成噪声样本的分布函数曲线。 图11示意了生成噪声样本的累积分布函数曲线。图7示意了图形用 户界面以及AWG控制的设置。
%------------------------------------------------------------- 函数DemoImpair2; %------------------------------------------------------------- %DemoImpair2 % Matlab编程语音编写的代码,演示了减损发生器的基本算法。 % % 所演示的算法全面控制预定质量标准,例如: % -在用连续频谱产生噪声时的频域和时域特性(幅度;概率分布) % -在用离散频谱生成噪声时的载波幅度、载频、调制深度和调制 % 宽度 % 独立计算两种类型的噪声,并且在时域中作为数字阵列表示。 % 可以通过增加这些阵列元素的方式同时得到两种类型的噪声。 % % % % % (c)2000-2001 KPN Research; % %演示函数 % DemoXtalkNoise -表示产生连续噪声的过程 % DemoIngressNoise -表示产生离散噪声的过程 % %主函数:噪声被表示为随机数阵列 % DefineShape -初始化所有用户可定义参数 % CreateNoiseCont -产生连续噪声 % CreateNoiseDiscr_Fast -产生离散噪声,快速算法 % CreateNoiseDiscr_Slow -产生离散噪声,慢速算法 % PrequencyShape -修改连续噪声的频谱密度 % AmplitudeShape -修改连续噪声的幅度分布 % %支持函数: % CalcSpec -计算噪声的频谱密度 % CalcNBSV -计算噪声的窄带信号电压 % % CalcCrest -计算噪声的波峰因数 % CalcDistrib -计算噪声的概率分布 % CalcCumDistrib -计算噪声的累积分布 % CalcSmooth -平滑频谱,如在实际频谱分析仪中那样 % % CalcEnhancedGaussDistribution -近高斯分布的示例 % CalcDemodulation -计算载波上调制的噪声 %---------------------------------------------------------- Shape=DefineShape; DemoXtalkNoise(Shape); DemoIngressNoise(Shape); %---------------------------------------------------------- function [U,t]=DemoXtalkNoise(Shape); %---------------------------------------------------------- %演示用连续频谱生成噪声 %如对于串音测试 R =Shape.R; CF_min= Shape.Xtalk.CF_min; [U,t]=CreateNoiseCont(Shape);plot(t,U); title(′Xtalk method 1′);shg;pause [X,f]=CalcSpec(U,t); plot(f,x); title(′Xtalk method 1′);shg;pause [X,f]=CalcSpec(U,t); plot(f,dBm(X,R));title(′Xtalk method 1′);shg;pause [P,u]=CalcCumDistrib(U); plot(u,P); title(′Xtalk method 1′);shg;pause % U=AmplitudeShape(U,Shape); plot(t,U); title(′Xtalk method 2′);shg;pause [X,f]=CalcSpec(U,t); plot(f,X) title(′Xtalk method 2′);shg;pause [X,f]=CalcSpec(U,t); plot(f,dBm(X,R));title(′Xtalk method 2′);shg;pause [P,u]=CalcCumDistrib(U); plot(u P); title(′Xtalk method 2′);shg;pause % U=FrequencyShape(U,Shape); plot(t,U); title(′Xtalk method 3′);shg;pause [X,f]=CalcSpec(U,t); plot(f,X); title(′Xtalk method 3′);shg;pause [X,f]=Calcspec(U,t); plot(f,dBm(X,R));title(′Xtalk method 3′);shg;pause [P,u]=CalcCumDistrib(U); plot(u,P); title(′Xtalk method 3′);shg;pause % for i=2∶10 i U=AmplitudeShape(U,Shape); [X,f]=CalcSpec(U,t); % U=FrequencyShape(U,Shape); [X,f]=CalcSpec(U,t); % if CalcCrest(U)>CF_min,break;end; end; [P,u]=CalcCumDistrib(U); plot(t,U); title(′Xtalk method 4′);shg;%暂停 plot(f,dBm(X,R)); title(′Xtalk method 4′);shg;%暂停 plot(u,P); title(′Xtalk method 4′);shg;%暂停 %---------------------------------------------------------- function[U,t]=DenoIngressNoise(Shape); %---------------------------------------------------------- %演示用离散频谱生成噪声,如对于干扰测试 R=Shape.R; [U,t]=CreateNoiseDiscr_Fast(Shape); %[U,t]=CreateNoiseDiscr_Slow(Shape);%给出同一结果 [X,f]=CalcNBSV(U,t);plot(f,dBm(X,R));title(′Ingress method′);shg;pause for ToneNr=[1∶2] [Uac,Uac_rms]=CalcDemodulation(U,t,Shape,ToneNr); plot(t,Uac);title(′demodulated ingress noise ofone carrier′);shg;pause [P,u]=CalcDistrib(Uac/Uac_rms); plot(u,P);title(′distribution of demod noise′);shg,pause end; %-------------------------------------------------------------------------------- function[Shape]=DefineShape; %-------------------------------------------------------------------------------- % 创建应该生成的噪声的噪声分布,按照 % -频谱密度(本例示中实际上是矩形) % -概率分布(本例示中是近高斯分布) % -信号音和调制 % 以伏特每平方根赫兹为单位的频谱 %--------------------------------------------------------------------------------- Fmax=4E6;Fl=300E3;Fh=700E3;N=2^18;R=135; cf=5.5; %期望波峰因数 cf min=5.1;%期望波峰因数 m=N/2; % Shape.N=N; %时间取样数 Shape.m=m; %频率取样数 Shape.dF=Fmax./(m-1); %频率间隔 Shape.dT=l/(N*Shape.dF); %时间间隔 Shape.R=R; %期望噪声源的阻抗 %定义串音噪声目标(频谱密度&幅度分布) Shape.Xtalk.freq=[0:M-1]′*Shape.dF; Shape.Xtalk.spec=(Shape.Xtalk.freq>=Fl).*(Shape.Xtalk.freq<=Fh)*(1/300); Shape.Xtalk.DistU=0:cf/1000:cf; Shape.Xtalk.DistP=CalcEnbancedGaussDistribution(Shape.Xtalk.DistU,cf);%P Shape.Xtalk.CF_min=of_min; %定义干扰噪声目标(RFI音) P_dBm=[-70;-50;-60;-60;-40;-60;-60;-40;-70;-40];%dBm@135 ohm P=(10).^(P_dBm/10)*1E-3; Shape.Ingress.ToneU=sqrt(P*135); %U=sqrt(P*R);有效值 Shape.Ingress.ToneF=[99;207;333;387;531;603;711;801;909;981]*1E3; Shape.Ingress.ModDepth=0.32*ones(10,1);%=mod index>0.8,at CF>2.5) Shape.Ingress.ModWidth=2*4.5E3*ones(10,1);%=-10KHz..+10KHz) %--------------------------------------------------------------------------- function[U,t]=CreateNoiseCont(Shape); %--------------------------------------------------------------------------- %创建噪声电压U(t),具有预定频域特性(频谱),但具有无控制的时域特性 %(分布) N=Shape.N; %要产生的取样数 %U=rand(N,1); %均匀分布的白噪声 U=randn(N,1); %高斯分布的白噪声 U=FrequencyShape(U,Shape); %整形噪声 t=[0:N-1]′*Shape-dT; %相关时间轴 %---------------------------------------------------------------------------- function[U,t]=CreateNoiseDiscr_Fast(Shape); %---------------------------------------------------------------------------- %创建电压U(t),具有AM已调载波(RFI音);每个具有单独的预定频率、 %幅度、调制宽度以及调制深度。 %每个载波音下边带的随机相位被镜像,以将任意QAM调制转换到 %(不是镜像为)更严格的AM调制(全镜像) %标记:本算法中X指的是(近谐波)干扰噪声信号的傅里叶级数的分量, %但是它在(伪随机)串音噪声信号的情况下指的是频谱密度 % % %计算时间随着取样数量大约线性增长 %全部计算时间中的大约80%由反傅里叶变换造成 %----------------------------------------------------------------- N=Shape.N; %取样数量 m=Shape.m; %该数量的一半 Nc=round(Shape.Ingress.ToneF/Shape.dF)+1; %载频的指数(只有正) Nm=round(Shape.Ingress.ModWidth/Shape.dF/2); %调制次数 %分量 Xc=0.5*Shape.Ingress.ToneU; %载波幅度的幅度 Xm=Shape.Ingress.ModDepth.*Xc./sqrt(2*Nm); %调制频带的幅度 X=zeros(N,1); %初始化 Xc=Xc.*exp(j*1000*rand(size(Xc))); %随机载波相位 Xcc=(Xc.*Xc)./abs(Xc.*Xc); for k=1:length(Nc) %对于所有已调载波,进行: Nmp=Nc(k)+[1:Nm(k)]′; %定位上边带 %频率 Nmn=Nc(k)-[1:Nm(k)]′; %定位上边带 %频率 Xmp=Xm(k).*exp(j*1000*rand(size(Nmp))); %创建上边带 Xmn=conj(Xmp)*Xcc(k); %镜像下边带 X(Nmp)=Xmp; %插入上边带 X(Nmn)=Xmn; %插入下边带 end; X(Nc)=Xc; %插入所有载波 X(N:-1:m+2)=conj(X(2:ceil(m))); %附加频谱 %(负频率) %U=real(ifft(x))*N; %变换到时域 U=real(fft(X)); %变换到时域 %(快了10%) t=Shape.dT*[0:N-1]′; %相关时间轴 %------------------------------------------------------------------------------- function[U,t]=CreateNoiseDiscr_Slow(Shape); %------------------------------------------------------------------------------- %创建电压U(t),具有预定频率上的RFI音、幅度和调制宽度以及调制 %深度。 %本算法是直接的,非常低效的,并且只用于演示 %可以证明CreateNoiseDiscr_Fast返回相同的结果 %------------------------------------------------------------------------------- N=Shape.N; %要生成的取样数量 m=Shape.m; f=[0:N-1]′*Shape.dF; t=[0:N-1]′*Shape.dT; Fc=Shape.Ingress-ToneF; %载波频率列表 Fc=Shape.dF*round(Fc/Shape.dF); %强制整数周期 % U=0; for k=1:length(Shape.Ingress.ToneF); %-产生噪声调整,具有U_avg=0并且U_rms=1 Nm=round(Shape.Ingress.ModWidth(k)/Shape.dF/2); Xn0=([1:N]<=(Nm+1))′; %形状调制噪声幅度 Xn=Xn0.*exp(j*1000*rand(N,1)); %形状调制噪声相位 Xn(1)=0; %消去DC分量 Xn(N:-1:m+2)=conj(Xn(2:ceil(m)));%附加频谱(负频率) Noise=real(ifft(Xn)); %变换到时域 Noise=Noise/sqrt(sum(Noise.*Noise)/N); %强制rms=1 %--执行调制 Carrier=Shape.Ingress.ToneU(k)*cos(2*pi*Fc(k)*t+1000*rand); Modulate=Shape.Ingress.ModDepth(k)*Noise; U=U+Carrier.*(1+Modulate); end; %---------------------------------------------------------------- function [U]=FrequcyShape(U,Shape) %---------------------------------------------------------------- %重新整形取样U的频谱,如目标形状指定的 %输入: % U: 取样的连续值 % Fs: 取样频率 % 频谱: 期望PSD(以伏特每平方根赫兹为单位) %---------------------------------------------------------------- N=length(U); m=length(Shape.Xtalk.spec); %m=N/2 t=[0:N-1]′*Shape.dT; %执行频率调整 Scaling=Shape.Xtalk.spec./CalcSpec(U,t); X=fft(U); %变换到频域 X(1)=0; %消去DC分量 X(2:m+1)= X(2:m+1).*Scaling; %调整频谱(正频率) X(N:-1:m+2)=conj(X(2:ceil(m))); %附加频谱(负频率) U=real(ifft(X)); %变换到时域 %------------------------------------------------------------------------------------- function[U]=AmplitudeShape(U,Shape) %------------------------------------------------------------------------------------- %本函数利用幅度相关(非线性)失真函数Q(x)对函数U的幅度分布进行 %整形。 %结果是U(t)=Q{U(t)} % %令FF为取样的实际累积分布函数,以及 %令GG为期望累积分布函数, %则给出的失真函数为: % Q(X)=GG^{-1}FF(X) %------------------------------------------------------------------------------------- U0=sqrt(sum(U.*U)/length(U)); %缩放因子(用于标准化) %计算失真函数Q [DistP1,DistU1]=CalcCumDistrib(U/U0); %实际分布 Q=interp1(Shape.Xtalk.DistP,Shape.Xtalk.DistU,DistP1);%失真函数 U=U0*interp1(DistU1,Q,abs(U/U0)).*sign(U); %执行失真 %plot(DistU1,Q);shg;%暂停 %------------------------------------------------------------------------------------- function[X,f]=CalcSpec(U,t); %------------------------------------------------------------------------------------- %当要在指定分辨率带宽上‘测量’时,计算信号的频谱密度 % RBW=1E3; %RBW: U的频谱的期望分辨率 N=length(U);m=N/2; dT=t(2)-t(1); %时间间隔 dF=1/dT/N; %频率间隔 f=[0:m-1]′*dF; %所有正频率 X=fft(U)*dT; %到频域 X=abs(X(2:m+1)); %无DC并且无负频率 X=sqrt(CalcSmooth(X.*X,f,RBW)); %在整个带宽RBW上对它进行平均 %------------------------------------------------------------------ function [X,f]=CalcNBSV(U,t); %------------------------------------------------------------------ %当要在指定分辨率带宽上‘测量’时,计算信号的频谱密度 % RBW=1E3; %RBW: U的频谱的期望分辨率 N=length(U);m=N/2; dT=t(2)-t(1); %时间间隔 dF=1/dT/N; %频率间隔 f=[0:m-1]′*dF; %所有正频率 %X=fft(U)*dT*sqrt(dF); %到频域 %X=fft(U)/N; %到频域 X=fft(U)/N*2; %到频域 X=abs(X(2:m+1)); %无DC并且无负频率 %X=sqrt(CalcSmooth(X.*X,f,RBW));%在整个带宽RBW上对它进行平均 %-------------------------------------------------------------------- function[CF]=CalcCrest(U) %-------------------------------------------------------------------- %计算信号U(t)的波峰因数,它是峰值除以rms值 % Urms=sqrt(sum(U.^2)/length(U)); Upeak=max(abs(U)); CF=Upeak/Urms; %-------------------------------------------------- function[F]=CalcEnhancedGaussDistribution(x,Cf); %-------------------------------------------------- %生成累积分布函数F(x),它被标识为‘增强高斯分布’ % %Cf=波峰因数 Alpha=le-3; Sigma=sqrt((1+Alpha)-cf^2*Alpha/3); x=x.*(x>0).*(x<Cf)+Cf*(x>=Cf); denominator=Alpha+erf(Cf(sqrt(2)*Sigma)); F=1-(Alpha*x/Cf+erf(x/(sqrt(2)*Sigma)))/denominator; %--------------------------------------------------- function[DistP,DistU,P]=CalcDistrib(U) %--------------------------------------------------- %计算信号U的幅度分布 N=length(U); Nbins=100; [cumbin,xx]=hist(U,Nbins); dX=xx(3)-xx(2); DistP=cumbin(:)N/dX; %强制sum(DistP)*dX DistU=xx(:); Urms=sqrt(sum(U.*U)/N); P=exp(-0.5*(DistU/Urms).^2);P=P/sum(P)/dX; DistU=[DistU,DistU]; DistP=[DistP,P]; %------------------------------------------------------- function[DistP,DisU]=CalcCumDistib(U) %------------------------------------------------------- %计算信号U的(向后)累积幅度分布 len=length(U); Ueff=sqrt(sum(U.*U)/length(U)); U=abs(U/Ueff); %---估计分布函数 Nbins=min([50,floor(len/10)]); [cumbin,xx]=hist(U,Nbins); BinWidth=xx(2)-xx(1); DistU=xx-BinWidth/2;%偏移 for n=[Nbins-1:-1:1];cumbin(n)=cumbin(n)+cumbin(n+1);end DistP=cumbin/len; %---当它们利用这个结果时,改进其它程序的数值稳定性 DistU=[0,DistU(2:end)]; %从x=0开始 DistP=[DistP,1/len]; DistU=[DistU,xx(Nbins)+0.999*BinWidth/2];%增加最终(信号)点 DistP=[DistP,1e-100]; DistU=[DistU,xx(Nbins)+(1.001)*BinWid/2];%因数1.001得到稳定性 %---------------------------------------------------------------- function,[PSD,freq]=CalcSmooth(PSD,freq,RBW) %---------------------------------------------------------------- %模拟实际频谱分析仪,具有有限分辨率的带宽以及高斯整形带通滤波器 %PSD=“功率频谱密度”,它是“频谱密度”的去平方(de square), %单位为伏特每平方赫兹 % %---------------------------------------------------------------- N=length(PSD); df=freq(2)-freq(1); br=3*floor(RBW/df); factor=2*br+1; if(factor>1) ..ff=df*(-br:br); %平滑间隔 ..mask=exp(-ff*ff/(2*RBW^2)); ..mask=mask/sum(mask); %分辨率带通滤波器的高斯掩蔽 ..xhelp=[PSD;zeros(2*br,1)]; lp=filter(mask,1,xhulp); %智能卷积 PSD=yhelp(br+1:end-br); end; %----------------------------------------------------------------- function[Uac,Uac_rms]=CalcDemodulation(U,t,Shape,ToneNr); %----------------------------------------------------------------- %解调已在离散噪声的载波上调制的噪声,只是为了演示,以证明离散噪声 %满足预定参数 % %解调器利用同步检测,不是相的 %结果是全解调频带上的未知衰减 %通过测量DC电平进行后向校正,并且放大解调信号直到此DC电平被 %归一化为1伏特 % %证明:(psi未知) % %令 Uff=cos(w*t+psi)*(1+Uac); %=用“1+Uac”调制载波 % Uc=cos(w*t); %=载波 % Ud=Urf*Uc; %=同步检测信号 % %则 Ud=1/2*(cos(psi+2*w*t)+cos(psi))*(1+Uac); % Ulf=cos(psi)/2*(1+Uac); %低通滤波之后 % Udc=cos(psi)/2; %=通过平均U1f % Uac=(Ulf/Udc)-1; %----------------------------------------------------------- N=Shape.N; Fc=Shape.Ingress.ToneF(ToneNr); %选择载频 Fc=Shape.dF*round(Fc/Shape.dF); %强制整数周期 ModWidth=Shape.Ingress.ModWidth(ToneNr); ModDepth=Shape.Ingress.ModDepth(ToneNr); % Ud=U.*cos(2*pi*Fc*t); %已调载波的同步检测 Nm=round(1.1*ModWidth/Shape.dF/2); %计算滤波频率 mask=zeros(N,1);mask([1:Nm,N-Nm:N]=1;%创建滤波器 Ulf=real(ifft(fft(Ud).*mask; %执行低通滤波 Udc=sum(Ulf)/N; %找到未归一化DC电平 Uac=Ulf/Udc-1; %归一化整个电平,并且去掉DC % Uac_rms=sqrt(sum(Uac.*Uac/N)); %必须等于ModDepth,因为Udc=1 Scale=Uac_rms/ModDepth; %必须为“1”
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈