CN1674009A - 用于s参数计算的设备和方法及其程序和记录介质 - Google Patents
用于s参数计算的设备和方法及其程序和记录介质 Download PDFInfo
- Publication number
- CN1674009A CN1674009A CNA2005100067267A CN200510006726A CN1674009A CN 1674009 A CN1674009 A CN 1674009A CN A2005100067267 A CNA2005100067267 A CN A2005100067267A CN 200510006726 A CN200510006726 A CN 200510006726A CN 1674009 A CN1674009 A CN 1674009A
- Authority
- CN
- China
- Prior art keywords
- parameter
- length
- individual devices
- parameters transformation
- calculation
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Algebra (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Measurement Of Resistance Or Impedance (AREA)
- Monitoring And Testing Of Transmission In General (AREA)
Abstract
用于S参数计算的设备和方法及其程序和记录介质。根据本发明,S参数文件输入部输入针对构成经受传输分析的系统的每个器件的S参数,端口位置对应关系转换部转换所输入的S参数,以使其与待连接的预设端口位置相对应,S参数/T参数转换部把S参数转换为T参数,连接计算部执行T参数的连接计算以获得整个系统的T参数,T参数/S参数转换部把整个系统的T参数转换为整个系统的S参数,而S参数文件输出部输出整个系统的S参数。
Description
技术领域
本发明涉及用于设计高频电路等的S参数计算技术,更具体来说,涉及用于S参数计算的设备和方法及其程序和记录介质,其能够准确地计算由多个器件构成的整个系统的S参数,并计算诸如缆线或布线的具有大纵横比的单个器件的任意长度的S参数。
背景技术
首先,解释S参数。S参数是表示电路的输入与输出之间的关系的参数。如图13中所示,当从带有4个端口端子的电路的端口1进行输入时,在从端口1到端口4中的各端口处将出现与特定特性相对应的输出。这些特性是针对端口1的反射、针对端口3的透射以及针对端口2和端口4的串扰。
图13中所示电路的输入与透射、反射及串扰的输出之间的关系可表示为图14中所示的矩阵计算的形式。图14中的矩阵(复矩阵)的矩阵元素(S11和S44)即为S参数,它们表示了电路的输入与输出之间的关系。准确地说,S参数表示多个端口中的电压和相位的绝对值之间的关系。
下面把图15中所示的带有4个端口端子的电路的输入/输出作为示例进行说明。在该图中,假定附于相应端口的标号1到4表示端口编号,ai(1≤i≤4)表示从每个端口的输入,而bi(1≤i≤4)表示每个端口的输出。
图15中所示的带有4个端口端子的电路的输入与输出之间的关系,可通过矩阵元素由S参数S11到S44组成的复矩阵表示为:
可以通过三维电磁场分析和实际测量来提取S参数。使用通用电路仿真工具SPICE来执行利用S参数的传输分析(transmission analysis)。SPICE通过合并输入S参数来执行传输分析。
当计算连接有多个装置的整个系统的S参数时,利用基于三维电磁场分析的技术通过把整个系统视为仿真对象来提取S参数的尝试将使得仿真规模变得太大,以至于不能计算S参数。因此,基于三维电磁场分析的技术存在无法同时提取整个系统的S参数的问题。
例如,如图16中所示,当包括LSI 110、BGA(球栅阵列)111、联接器112、BP(底板)113、为BP 113设置的通孔(未示出)、PCB(印刷电路板)114以及布线115的系统经受三维电磁场分析时,网格数和计算时间增加,因此无法同时提取整个系统的S参数。
而且,不能通过三维电磁场分析以足够的准确度来提取具有大纵横比的器件(例如,如图17中所示的布线或缆线)的S参数,并且每当布线或缆线等的长度改变时,都必须对所述S参数进行实际测量。
另一方面,对于基于实际测量的技术,例如,如图18中所示,当将诸如探针(probe)的测量系统210接合在连接器112的上面和下面以测量连接器112的S参数时,被测量的S参数不仅包括连接器112的S参数,而且包括测量系统210的S参数,这导致在提取连接器112本身的S参数的过程中的准确度降低的问题。
此外,还存在另一问题,即,难以测量诸如缆线或布线的具有大纵横比的器件的S参数。
如上所述,传统上一直难以计算由多个器件构成的整个系统的S参数,因此,当对整个系统进行传输分析时,常规技术通过逐个地提取单个器件的S参数并在SPICE模型上连接所提取的S参数,来执行传输分析。
例如,常规技术对图16中所示的LSI 110、BGA 111、布线115以及连接器112独立地执行三维电磁场分析,逐个装置地提取S参数,在如图19中所示的SPICE模型上连接各个装置的S参数,并执行对如图16中所示的整个系统的传输分析。
下面通过利用图20A、图20B以及图21,对常规技术进行更具体的阐述。假定提取了如图20A中所示的附连试验板(test coupon)211的S参数,并且执行了使用SPICE的传输分析。在图20A中,标号212表示探针,115表示布线。图20B是图20A中所示的附连试验板211的截面图。图20B中所示的附连试验板211的标记“SUS”表示一不锈钢部分。
如图20A中所示,S参数提取目标范围与长度为4cm的部分相对应,该长度对应于布线115的导线长度。由于S参数提取目标范围较长,所以无法作为单个模型来实现三维电磁场分析。因此,把附连试验板211划分成两个区域:左部和右部,并提取各个区域的S参数。
图21示出了常规技术的原理。通过三维电磁场分析来单独地提取所划分的附连试验板211的左部的S参数和右部的S参数。然后,根据提取的左部的S参数和右部的S参数来分别创建SPICE模型,并且这两个SPICE模型被连接起来,并经受传输分析。
没有文献描述用于执行整个系统的传输分析的技术,这种技术基于多个S参数使用连接代码来自动计算整个系统的S参数,获取整个系统的S参数,然后创建SPICE模型。
当使用常规技术来执行传输分析时,对以下公式1中所示的逆FFT(快速傅立叶变换)执行与在SPICE模型上所连接的S参数的个数相同的次数。
公式1可按离散区表示为公式2:
由于常规技术重复地在SPICE模型上连接S参数并按公式2执行积分计算,误差可能累积,从而产生影响传输分析的误差。因此,难以执行准确的传输分析。
图22示出了对图20A中所示的附连试验板211所施行的传输分析的结果,通过利用常规技术来提取如图21中所示的左部和右部上的附连试验板的S参数,并通过根据各个S参数来创建SPICE模型,来施行该传输分析。该图中的细线表示测量值,而粗线表示使用SPICE的传输分析的结果。如图22中由虚线所包围的区域所示,当利用常规技术来执行SPICE传输分析时,可以观察到产生了误差(波形的扰动)。
而且,常规技术无法根据诸如缆线或布线的具有大纵横比的单个器件的某些S参数来计算任意长度的S参数。
发明内容
本发明的一个目的是提供用于计算S参数的装置,其能够解决上述常规技术的问题,准确地计算由多个器件所组成的整个系统的S参数,并计算诸如缆线或布线的具有大纵横比的单个器件的任意长度的S参数。
本发明采用了以下所述的技术。图1示出了本发明的原理。例如,当执行对图20A中所示的附连试验板211的传输分析时,本发明执行附连试验板211的左部的S参数和右部的S参数(这些S参数已通过三维电磁场分析或实际测量而被单独地提取出来了)的连接计算,并计算整个附连试验板211的S参数。
利用由本发明所计算出的整个附连试验板211的S参数,由SPICE进行的传输分析只需要执行一次公式2所示的逆FFT,从而能够以最小的误差执行传输分析,并以高准确度执行传输分析。
即,本发明是一种用于计算由多个相同或不同器件组成的整个系统的S参数的S参数计算设备。该设备包括:输入装置,用于输入针对每个所述器件的S参数;端口位置对应关系转换装置,用于转换所输入的S参数,以使其与每个所述器件的预设端口位置相对应S参数/T参数转换装置,用于把所述转换后的S参数转换为针对每个所述器件的T参数;计算装置,用于利用针对每个所述器件的T参数来执行矩阵计算,并且计算所述整个系统的T参数;T参数/S参数转换装置,用于把所述整个系统的所述计算出的T参数转换为所述整个系统的S参数;以及输出装置,用于输出由所述T参数/S参数转换装置所获得的所述整个系统的S参数。
而且,本发明是一种用于计算长度为一定长度L的N(N表示值为2或更大的任意整数)倍的单个器件的S参数的S参数计算设备。该设备包括:输入装置,用于输入所述长度为L的单个器件的S参数;端口位置对应关系转换装置,用于转换所输入的S参数,以使其与连接所述单个器件时的预设端口位置相对应;S参数/T参数转换装置,用于把所述转换后的S参数转换为T参数;计算装置,用于利用所述T参数执行矩阵计算,将对被顺序连接的长度与所述长度为L的单个器件的长度相对应的单个器件的T参数的计算过程重复(N-1)次,并且计算长度为L×N的单个器件的T参数;T参数/S参数转换装置,用于把所述长度为L×N的单个器件的所述计算出的T参数转换为长度为L×N的单个器件的S参数;以及输出装置,用于输出由所述T参数/S参数转换装置所获得的所述长度为L×N的单个器件的S参数。
而且,本发明是一种用于计算由多个相同或不同器件组成的整个系统的S参数的方法。该方法包括以下步骤:输入步骤,输入针对每个所述器件的S参数;端口位置对应关系转换步骤,转换所输入的S参数,以使其与每个所述器件的预设端口位置相对应S参数/T参数转换步骤,把所述转换后的S参数转换为针对每个所述器件的T参数;计算步骤,利用针对每个所述器件的T参数来执行矩阵计算,并且计算所述整个系统的T参数;T参数/S参数转换步骤,把所述整个系统的所述计算出的T参数转换为所述整个系统的S参数;以及输出步骤,输出由所述T参数/S参数转换步骤所获得的所述整个系统的S参数。
而且,本发明是一种用于计算由多个相同或不同器件组成的整个系统的S参数的S参数计算程序。该程序使得计算机执行以下步骤:输入步骤,输入针对每个所述器件的S参数;端口位置对应关系转换步骤,转换所输入的S参数,以使其与每个所述器件的预设端口位置相对应;S参数/T参数转换步骤,把所述转换后的S参数转换为针对每个所述器件的T参数;计算步骤,利用针对每个所述器件的T参数来执行矩阵计算,并且计算所述整个系统的T参数;T参数/S参数转换步骤,把所述整个系统的所述计算出的T参数转换为所述整个系统的S参数;以及输出步骤,输出由所述T参数/S参数转换步骤所获得的所述整个系统的S参数。
而且,本发明是一种记录S参数计算程序的计算机可读记录介质,该S参数计算程序用于计算由多个相同或不同器件组成的整个系统的S参数。该程序使得计算机执行以下步骤:输入步骤,输入针对每个所述器件的S参数;端口位置对应关系转换步骤,转换所输入的S参数,以使其与每个所述器件的预设端口位置相对应;S参数/T参数转换步骤,把所述转换后的S参数转换为针对每个所述器件的T参数;计算步骤,利用针对每个所述器件的T参数来执行矩阵计算,并且计算所述整个系统的T参数;T参数/S参数转换步骤,把所述整个系统的所述计算出的T参数转换为所述整个系统的S参数;以及输出步骤,输出由所述T参数/S参数转换步骤所获得的所述整个系统的S参数。
本发明以高准确度计算由多个器件组成的整个系统的S参数。因此,根据本发明,与通过把每个器件的S参数转换为SPICE模型来逐个地执行传输分析的常规技术相比,可以更快且更准确地执行分析。
而且,根据本发明,甚至可以计算针对诸如缆线或布线的具有大纵横比的单个器件的任意长度的S参数。
附图说明
图1示出本发明的原理。
图2示出S参数计算设备的结构的示例。
图3示出S参数文件的数据结构的示例。
图4示出从S参数到T参数的转换。
图5示出两个4端口电路的连接。
图6A示出输入文件端口位置和基于输入文件端口位置的S参数的示例。
图6B示出代码设置端口位置和基于代码设置端口位置的S参数的示例。
图7示出端口编号输入屏面的示例。
图8示出不同S参数的连接过程流程。
图9示出任意长度的S参数的计算。
图10示出任意长度的缆线的S参数计算过程流程的示例。
图11示出本发明的应用示例。
图12A示出计算微带线(micro strip line)的S参数时的反射特性。
图12B示出计算微带线的S参数时的传输特性。
图13示出电路的输入/输出关系。
图14示出S参数。
图15示出带有4个端口端子的电路的输入/输出关系。
图16示出一个系统的示例。
图17示出布线或缆线。
图18示出连接器的S参数的测量示例。
图19示出在SPICE模型上连接S参数。
图20A示出附连试验板。
图20B是附连试验板的截面图。
图21示出常规技术的原理。
图22示出使用常规技术的传输分析结果。
具体实施方式
图2示出本发明的S参数计算设备的结构的示例。在该图中,标号1表示S参数计算设备,其包括诸如CPU和存储器的硬件以及软件程序,10表示为每个所预先计算的S参数文件。在S参数计算设备1中,标号11表示S参数文件输入部,12表示端口位置对应关系转换部,13表示S参数/T参数转换部,14表示连接计算部,15表示T参数/S参数转换部,以及16表示S参数文件输出部。
S参数文件输入部11从预先设置的S参数文件10为构成传输分析目标系统的每个器件输入S参数。而且,S参数文件输入部11在计算缆线或布线等的任意长度的S参数时,还可以通过指定针对缆线或布线等的一部分所提取的S参数文件10来输入S参数。
图3示出S参数文件的数据结构的示例。例如,如图3中所示,S参数文件10针对每个频率存储S参数(S11到S44)。典型的数据格式是检验标准(Touchstone)格式。这里,某个频率处的S参数表示为MA=(幅值,幅角)。除此之外,还有DB=(dB,幅角)和RI=(实部,虚部)的记法,并且本发明还适用于其他记法的S参数。
端口位置对应关系转换部12执行端口位置对应关系转换处理,其使得从S参数文件10输入的S参数的端口位置与按为连接计算所预设的顺序的端口位置相对应,并对S参数进行重排。后面将详细描述该端口位置对应关系转换处理。
S参数/T参数转换部13将S参数转换为T参数,以使得可以进行矩阵计算。例如,如图4中所示,将S参数S11到S44转换为T参数T11到T44。例如,图4中所示的T参数是表示图15中所示的电路的左侧的输入/输出与右侧的输入/输出之间的关系的矩阵元素。
连接计算部14执行通过为每个器件转换S参数而获得的T参数的连接计算,并获得待经受传输分析的整个系统的T参数。T参数/S参数转换部15将整个系统的T参数转换为整个系统的S参数。S参数文件输出部16输出所计算出的整个系统的S参数。
下面将对根据本发明的S参数的连接计算进行详细阐述。S参数/T参数转换部13如图4中所示将每个器件的S参数转换为T参数,从而使得可以连接电路。这是因为S参数表示电路的输入与输出之间的关系,而T参数表示电路的左侧的输入/输出与右侧的输入/输出之间的关系,因此可以使用矩阵计算来执行连接/分离操作。
图5示出两个4端口电路的连接。在图5中,aiA(1≤i≤4)和biA(1≤i≤4)表示电路A的从端口1到端口4中的每个端口的输入和输出,而aiB(1≤i≤4)和biB(1≤i≤4)表示电路B的从端口1到端口4中的每个端口的输入和输出。电路A的T参数为TA,而电路B的T参数为TB。
利用TA将电路A的左侧的输入/输出与右侧的输入/输出之间的关系表示为:
而利用TB将电路B的左侧的输入/输出与右侧的输入/输出之间的关系表示为:
当电路A的端口4连接到电路B的端口1并且电路A的端口2连接到电路B的端口3时有:
由此,根据图5中所示的连接,整个电路的左侧的输入/输出与右侧的输入/输出之间的关系可通过以下公式3中所示的矩阵计算来表示:
公式3中的两个电路的T参数TA·TB的相乘结果变成连接后的整个电路的T参数。当把连接后的该整个电路的T参数转换为S参数时,可以获得连接后的整个电路的S参数。
接下来,利用图6A、图6B以及图7对由端口位置对应关系转换部12进行的处理进行详细阐述。为测量S参数,对每个端子(端口)进行编号。编号顺序在实际测量的情况下为探针设置顺序,而在三维电磁场分析的情况下为端口设置顺序,并且编号顺序可以随测量者而变化。
端口编号被与S参数的矩阵元素编号相关联,并且连接要求相同的端口编号排列。然而,实际测量与三维电磁场分析之间的端口位置不同。为解决此问题,在从S参数文件输入S参数之后,本发明根据由代码(程序)所预设的端口位置来对输入S参数的矩阵元素进行重排。
图6A示出输出文件端口位置和基于该输入文件端口位置的S参数的示例,图6B示出代码设置端口位置和基于该代码设置端口位置的S参数的示例。图6A中的左图示出输入文件端口位置(即,与输入文件的S参数相对应的端口位置),而图6A中的右图示出基于输入文件的左图中的端口位置的S参数。此外,图6B中的左图示出代码设置端口位置(即,由代码所预设的端口位置),而图6B中的右图示出基于代码设置端口位置的S参数。
在执行连接计算的程序中,根据依照图6B中的端口位置(编号)而排列的S参数来执行计算处理,因此有必要根据代码设置端口位置,将从S参数文件输入的图6A中的右图中的S参数的元素转换为如图6B中的左图中所示,并使用转换后的S参数(S11到S44)。端口位置对应关系转换部12显示图7中所示的端口编号输入屏面,输入来自用户的端口位置对应关系信息,并执行S参数元素的位置的转换处理。
图7示出用于端口位置对应关系转换处理的端口编号输入屏面的示例。代码设置端口位置(其为与图6B中的左图上相同的端口位置)示出在图7中的输入屏面中的右侧。
在由图7中的虚线所包围的输入文件端口位置设置屏面中,用户在四个端口位置处输入端口编号1到4。例如,在图7中,通过输入端口编号,设置了与图6A中的左图中相同的端口位置。图6A中的右图示出了对应于这些输入文件端口位置的S参数。如果用户在设置了输入文件端口位置之后点击“确定”按钮,则将如图6A中的右图所示的S参数的排列转换为图6B中的右图中的S参数(S11到S44)。
在本发明中,根据代码设置位置来转换输入文件的S参数。因此,即使输入的S参数文件对应于不同的端口位置,也可以转换S参数。在本发明中,通过执行上述端口位置对应关系转换处理,可以连接S参数的实际测量值和分析值。
图8示出根据本发明的S参数的连接过程流程。首先,输入通过实际测量或三维电磁场分析所获得的构成系统的多个器件的S参数文件(步骤S1)。接着,为了使得可以进行对实际测量值和分析值组合的连接,根据代码设置端口位置来转换S参数(步骤S2)。
将S参数转换为可对其应用矩阵计算的T参数(步骤S3)。然后,使用公式3执行连接计算(步骤S4)。把作为连接计算的计算结果的T参数转换为S参数(步骤S5),并把所计算出的S参数输出为一文件(步骤S6),然后处理结束。
在上述处理中,计算具有两个器件的整个系统的S参数的计算时间约为10秒到20秒。因此,显然,即使器件数目增加,也可以在相对短的时间中计算S参数。
如上所述,本发明通过复矩阵(由T参数组成的矩阵)的相乘来连接针对每个器件的S参数。本发明的范围不限于带有四个端子或八个端子的电路,而是也可应用于带有更多端口的电路。
接下来,对利用本发明的对任意长度的布线/缆线的S参数计算处理进行阐述。如图9中所示,本发明对作为对例如约5mm长的缆线执行三维电磁场分析的结果而获得的S参数文件进行连接,以计算具有任意长度的缆线的S参数。根据本发明,可以在较短的时间内通过利用缆线或布线等的一部分的长度的S参数,来计算任意长度的S参数。
图10示出任意长度的缆线的S参数计算过程流程的示例。首先,指定一S参数文件,其存储具有缆线的一部分的长度的S参数(步骤S11),输入连接计数N(步骤S12),并根据代码设置端口位置来转换S参数(步骤S13)。
例如,在图7中所示的端口编号输入屏面中,在“重复计数(整数)”中输入任意连接计数N(图7中为“5”),并在由虚线包围的区域所指示的输入文件端口位置设置屏面上输入端口编号1到4,以执行步骤S12、S13。
接下来,把S参数转换为T参数(步骤S14)。接下来,对由公式3所表示的矩阵计算执行N-1次(步骤S15),以获得具有希望长度的T参数。将所获得的T参数转换为S参数(步骤S16),将所获得的任意长度的S参数输出为文件(步骤S17),然后结束处理。
如上所述,本发明可以通过对由T参数所组成的复矩阵的重复计算来计算任意长度的缆线等的S参数,并且即使端口数目增加时(正如上述器件的连接的情况),本发明也是适用的。
图11示出了通过利用根据本发明所计算出的S参数的SPICE的传输分析结果与实际测量值之间的对比。在此实验中,把图20A中所示的附连试验板211划分成如图1中所示的左部和右部,利用根据本发明的S参数计算方法连接由对两个部分进行实际测量所测量到的S参数,并利用为整个附连试验板211计算出的S参数来执行SPICE传输分析。
在该图中,虚线表示整个附连试验板211的实际测量结果,而实线表示SPICE的传输分析结果。SPICE传输分析的结果与实际测量符合得很好,并且没有观察波形扰动。根据与图22中所示的使用常规技术的传输分析结果(即,把附连试验板211的左部和右部的S参数转换为SPICE模型并执行传输分析的结果)的比较,可以显见,通过利用由本发明所计算出的整个系统的S参数的传输分析,可以执行更准确的传输分析。
图12A和图12B示出在对利用三维电磁场分析而获得的5cm的微带线的S参数根据本发明施行连接计算并计算了具有双倍长度(10cm)和三倍长度(15cm)的微带线的S参数之后的瞬态分析结果。
图12A中的电压降100表示在利用三维电磁场分析获得5cm微带线的S参数时的反射特性,电压降101表示在连接S参数并获得具有双倍长度的微带线的S参数时的反射特性,而电压降102表示在连接S参数并获得具有三倍长度的微带线的S参数时的反射特性。
另一方面,图12B中的电压阶跃200表示在利用三维电磁场分析获得5cm微带线的S参数时的传输特性,电压阶跃201表示在连接S参数并获得具有双倍长度的微带线的S参数时的传输特性,而电压阶跃202表示在连接S参数并获得具有三倍长度的微带线的S参数时的传输特性。
对于反射,如从图12A显见,电压降的定时随着微带线的长度增加而延迟,而对于透射,如从图12B显见,电压阶跃的定时随着微带线的长度增加而延迟。
此外,电压下降定时100的与电压下降定时101的之间的时间间隔等于电压下降定时101与电压下降定时102之间的时间间隔,而电压阶跃定时200的与电压阶跃定时201之间的时间间隔等于电压阶跃定时201与电压阶跃定时202之间的时间间隔。
因此,显然,根据本发明正确地连接了S参数。由此,可以说根据本发明的S参数连接使得可以准确地计算缆线的任意长度的S参数。
可以利用计算机和软件程序来实现至此所述的根据本发明的S参数计算处理,并且还可以提供记录在计算机可读记录介质中的程序,或者在网络上提供该程序。
Claims (12)
1、一种S参数计算设备,用于计算由多个相同或不同器件组成的整个系统的S参数,该S参数计算设备包括:
输入装置,用于输入针对每个所述器件的S参数;
端口位置对应关系转换装置,用于转换所输入的S参数,以使其与每个所述器件的预设端口位置相对应;
S参数/T参数转换装置,用于把所述转换后的S参数转换为针对每个所述器件的T参数;
计算装置,用于利用针对每个所述器件的T参数来执行矩阵计算,并且计算所述整个系统的T参数;
T参数/S参数转换装置,用于把所述整个系统的所述计算出的T参数转换为所述整个系统的S参数;以及
输出装置,用于输出由所述T参数/S参数转换装置所获得的所述整个系统的S参数。
2、如权利要求1所述的S参数计算设备,其中,所述端口位置对应关系转换装置显示一屏面,以输入关于与所输入的S参数对应的器件的端口位置与用于连接计算的该器件的预设端口位置之间的对应关系的信息,并基于从该输入屏面输入的对应关系信息对所输入的S参数进行重排。
3、一种S参数计算设备,用于计算长度为一定长度L的N倍的单个器件的S参数,其中N表示值为2或更大的任意整数,该S参数计算设备包括:
输入装置,用于输入所述长度为L的单个器件的S参数;
端口位置对应关系转换装置,用于转换所输入的S参数,以使其与连接所述单个器件时的预设端口位置相对应;
S参数/T参数转换装置,用于把所述转换后的S参数转换为T参数;
计算装置,用于利用所述T参数执行矩阵计算,将对被顺序连接的长度与所述长度为L的单个器件的长度相对应的单个器件的T参数进行计算的过程重复(N-1)次,并且计算长度为L×N的单个器件的T参数;
T参数/S参数转换装置,用于把所述长度为L×N的单个器件的所述计算出的T参数转换为长度为L×N的单个器件的S参数;以及
输出装置,用于输出由所述T参数/S参数转换装置所获得的所述长度为L×N的单个器件的S参数。
4、如权利要求3所述的S参数计算设备,其中,所述端口位置对应关系转换装置显示一屏面,以输入关于与所输入的S参数对应的器件的端口位置与用于连接计算的该器件的预设端口位置之间的对应关系的信息,并基于从该输入屏面输入的对应关系信息对所输入的S参数进行重排。
5、一种用于利用计算机来计算由多个相同或不同器件构成的整个系统的S参数的方法,该方法包括以下步骤:
输入步骤,输入针对每个所述器件的S参数;
端口位置对应关系转换步骤,转换所输入的S参数,以使其与每个所述器件的预设端口位置相对应;
S参数/T参数转换步骤,把所述转换后的S参数转换为针对每个所述器件的T参数;
计算步骤,利用针对每个所述器件的T参数来执行矩阵计算,并且计算所述整个系统的T参数;
T参数/S参数转换步骤,把所述整个系统的所述计算出的T参数转换为所述整个系统的S参数;以及
输出步骤,输出由所述T参数/S参数转换步骤所获得的所述整个系统的S参数。
6、一种用于利用计算机来计算长度为一定长度L的N倍的单个器件的S参数的方法,其中N表示值为2或更大的任意整数,该方法包括以下步骤:
输入步骤,输入所述长度为L的单个器件的S参数;
端口位置对应关系转换步骤,转换所输入的S参数,以使其与连接所述单个器件时的预设端口位置相对应;
S参数/T参数转换步骤,把所述转换后的S参数转换为T参数;
计算步骤,利用所述T参数执行矩阵计算,将对被顺序连接的长度与所述长度为L的单个器件的长度相对应的单个器件的T参数进行计算的过程重复(N-1)次,并且计算长度为L×N的单个器件的T参数;
T参数/S参数转换步骤,把所述长度为L×N的单个器件的所述计算出的T参数转换为长度为L×N的单个器件的S参数;以及
输出步骤,输出由所述T参数/S参数转换步骤所获得的所述长度为L×N的单个器件的S参数。
7、一种用于计算由多个相同或不同器件构成的整个系统的S参数的S参数计算程序,该S参数计算程序使得计算机执行以下步骤:
输入步骤,输入针对每个所述器件的S参数;
端口位置对应关系转换步骤,转换所输入的S参数,以使其与每个所述器件的预设端口位置相对应;
S参数/T参数转换步骤,把所述转换后的S参数转换为针对每个所述器件的T参数;
计算步骤,利用针对每个所述器件的T参数来执行矩阵计算,并且计算所述整个系统的T参数;
T参数/S参数转换步骤,把所述整个系统的所述计算出的T参数转换为所述整个系统的S参数;以及
输出步骤,输出由所述T参数/S参数转换步骤所获得的所述整个系统的S参数。
8、一种用于利用计算机来计算长度为一定长度L的N倍的单个器件的S参数的S参数计算程序,其中N表示值为2或更大的任意整数,该S参数计算程序使得计算机执行以下步骤:
输入步骤,输入所述长度为L的单个器件的S参数;
端口位置对应关系转换步骤,转换所输入的S参数,以使其与连接所述单个器件时的预设端口位置相对应;
S参数/T参数转换步骤,把所述转换后的S参数转换为T参数;
计算步骤,利用所述T参数执行矩阵计算,将对被顺序连接的长度与所述长度为L的单个器件的长度相对应的单个器件的T参数进行计算的过程重复(N-1)次,并且计算长度为L×N的单个器件的T参数;
T参数/S参数转换步骤,把所述长度为L×N的单个器件的所述计算出的T参数转换为长度为L×N的单个器件的S参数;以及
输出步骤,输出由所述T参数/S参数转换步骤所获得的所述长度为L×N的单个器件的S参数。
9、一种记录S参数计算程序的计算机可读记录介质,该S参数计算程序用于计算由多个相同或不同器件构成的整个系统的S参数,该S参数计算程序使计算机执行以下步骤:
输入步骤,输入针对每个所述器件的S参数;
端口位置对应关系转换步骤,转换所输入的S参数,以使其与每个所述器件的预设端口位置相对应;
S参数/T参数转换步骤,把所述转换后的S参数转换为针对每个所述器件的T参数;
计算步骤,利用针对每个所述器件的T参数来执行矩阵计算,并且计算所述整个系统的T参数;
T参数/S参数转换步骤,把所述整个系统的所述计算出的T参数转换为所述整个系统的S参数;以及
输出步骤,输出由所述T参数/S参数转换步骤所获得的所述整个系统的S参数。
10、一种记录S参数计算程序的计算机可读介质,该S参数计算程序用于利用计算机来计算长度为一定长度L的N倍的单个器件的S参数,其中N表示值为2或更大的任意整数,该S参数计算程序使计算机执行以下步骤:
输入步骤,输入所述长度为L的单个器件的S参数;
端口位置对应关系转换步骤,转换所输入的S参数,以使其与连接所述单个器件时的预设端口位置相对应;
S参数/T参数转换步骤,把所述转换后的S参数转换为T参数;
计算步骤,利用所述T参数执行矩阵计算,将对被顺序连接的长度与所述长度为L的单个器件的长度相对应的单个器件的T参数进行计算的过程重复(N-1)次,并且计算长度为L×N的单个器件的T参数;
T参数/S参数转换步骤,把所述长度为L×N的单个器件的所述计算出的T参数转换为长度为L×N的单个器件的S参数;以及
输出步骤,输出由所述T参数/S参数转换步骤所获得的所述长度为L×N的单个器件的S参数。
11、一种输入/输出特性分析设备,用于利用S参数来分析由多个相同或不同器件构成的整个系统的输入/输出特性,该输入/输出特性分析设备包括:
用于利用所述S参数计算设备根据单个器件的S参数来计算整个系统的S参数的装置;和
用于利用所述整个系统的计算出的S参数通过SPICE来分析输入/输出特性的装置,
其中,所述S参数计算设备还包括:
输入装置,用于输入针对每个所述器件的S参数;
端口位置对应关系转换装置,用于转换所输入的S参数,以使其与每个所述器件的预设端口位置相对应;
S参数/T参数转换装置,用于把所述转换后的S参数转换为针对每个所述器件的T参数;
计算装置,用于利用针对每个所述器件的T参数来执行矩阵计算,并且计算所述整个系统的T参数;
T参数/S参数转换装置,用于把所述整个系统的所述计算出的T参数转换为所述整个系统的S参数;以及
输出装置,用于输出由所述T参数/S参数转换装置所获得的所述整个系统的S参数。
12、一种输入/输出特性分析设备,用于利用S参数来分析由多个相同或不同器件构成的整个系统的输入/输出特性,该输入/输出特性分析设备包括:
用于利用所述S参数计算设备根据单个器件的S参数来计算所述整个系统的S参数的装置;和
用于利用所述整个系统的计算出的S参数通过SPICE来分析输入/输出特性的装置,
其中,所述S参数计算设备用于计算长度为一定长度L的N倍的单个器件的S参数,其中N表示值为2或更大的任意整数,该S参数计算设备还包括:
输入装置,用于输入所述长度为L的单个器件的S参数;
端口位置对应关系转换装置,用于转换所输入的S参数,以使其与连接所述单个器件时的预设端口位置相对应;
S参数/T参数转换装置,用于把所述转换后的S参数转换为T参数;
计算装置,用于利用所述T参数执行矩阵计算,将对被顺序连接的长度与所述长度为L的单个器件的长度相对应的单个器件的T参数的计算过程重复(N-1)次,并且计算长度为L×N的单个器件的T参数;
T参数/S参数转换装置,用于把所述长度为L×N的单个器件的所述计算出的T参数转换为长度为L×N的单个器件的S参数;以及
输出装置,用于输出由所述T参数/S参数转换装置所获得的所述长度为L×N的单个器件的S参数。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2004088478A JP2005274373A (ja) | 2004-03-25 | 2004-03-25 | Sパラメータ算出装置,sパラメータ算出方法,sパラメータ算出プログラムおよびそのプログラムを記録したコンピュータ読み取り可能な記録媒体 |
JP2004088478 | 2004-03-25 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1674009A true CN1674009A (zh) | 2005-09-28 |
CN100410951C CN100410951C (zh) | 2008-08-13 |
Family
ID=34991197
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2005100067267A Expired - Fee Related CN100410951C (zh) | 2004-03-25 | 2005-01-31 | 用于s参数计算的设备和方法 |
Country Status (5)
Country | Link |
---|---|
US (1) | US7389191B2 (zh) |
JP (1) | JP2005274373A (zh) |
KR (1) | KR100682005B1 (zh) |
CN (1) | CN100410951C (zh) |
TW (1) | TWI292544B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107490732A (zh) * | 2016-06-13 | 2017-12-19 | 株式会社村田制作所 | 传输特性测定方法以及传输特性测定装置 |
CN107943641A (zh) * | 2017-11-16 | 2018-04-20 | 郑州云海信息技术有限公司 | 一种生成多端口s参数文档的方法及装置 |
CN108804775A (zh) * | 2018-05-22 | 2018-11-13 | 郑州云海信息技术有限公司 | 一种基于Matlab的高速信号散射参数计算方法及装置 |
CN109376342A (zh) * | 2018-11-01 | 2019-02-22 | 郑州云海信息技术有限公司 | 一种s参数端口转换方法与装置 |
CN112182997A (zh) * | 2020-09-04 | 2021-01-05 | 曙光信息产业(北京)有限公司 | 散射参数确定方法、装置、计算机设备和存储介质 |
WO2021031424A1 (zh) * | 2019-08-19 | 2021-02-25 | 苏州华太电子技术有限公司 | 散射参数测量方法及器件校准方法 |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8990280B2 (en) * | 2005-09-30 | 2015-03-24 | Nvidia Corporation | Configurable system for performing repetitive actions |
US7433981B1 (en) * | 2005-09-30 | 2008-10-07 | Nvidia Corporation | System and method for using co-processor hardware to accelerate highly repetitive actions |
CN101087148B (zh) * | 2006-06-07 | 2011-01-05 | 鸿富锦精密工业(深圳)有限公司 | 高速差分信号模拟分析系统及方法 |
JP2008116369A (ja) * | 2006-11-06 | 2008-05-22 | Fujitsu Ltd | パラメータ算出装置、方法、及びプログラム |
US7865319B1 (en) * | 2006-11-30 | 2011-01-04 | Lecroy Corporation | Fixture de-embedding method and system for removing test fixture characteristics when calibrating measurement systems |
TWI407115B (zh) * | 2007-12-31 | 2013-09-01 | Hon Hai Prec Ind Co Ltd | 一種測量無源電路s參數之方法和實施該方法之測量裝置 |
TWI659322B (zh) * | 2010-02-12 | 2019-05-11 | 新思科技股份有限公司 | 擴展具有不同傳輸線長度的模擬方法與裝置 |
US8892414B1 (en) | 2010-02-26 | 2014-11-18 | Sas Ip, Inc. | Transmission-line simulators and methods |
US8938372B1 (en) * | 2010-09-09 | 2015-01-20 | Sas Ip, Inc. | Simulating signal integrity structures |
US9063882B1 (en) | 2010-09-09 | 2015-06-23 | Sas Ip, Inc. | Matrix preconditioners for simulations of physical fields |
US8924186B1 (en) | 2010-09-09 | 2014-12-30 | Sas Ip, Inc. | Simulations of physical systems for multiple excitations |
US9477792B1 (en) | 2012-02-09 | 2016-10-25 | Sas Ip, Inc. | Enhancements to parameter fitting and passivity enforcement |
KR101363259B1 (ko) * | 2012-11-21 | 2014-02-12 | 경북대학교 산학협력단 | 저항 추출 장치, 저항 추출 방법 및 컴퓨터 판독가능 기록매체 |
KR102541065B1 (ko) * | 2021-06-28 | 2023-06-08 | 한국자동차연구원 | S-파라미터 생성 장치 및 방법 |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3016252B2 (ja) * | 1990-07-26 | 2000-03-06 | アジレント・テクノロジー株式会社 | ネットワーク・アナライザ |
JPH04104506A (ja) * | 1990-08-23 | 1992-04-07 | Fujitsu Ltd | 高周波電力増幅装置 |
DE4122730A1 (de) * | 1991-07-05 | 1993-01-07 | Siemens Ag | Schaltungsanordnung fuer eine leitungsnachbildungseinrichtung |
JP3147987B2 (ja) * | 1992-04-19 | 2001-03-19 | アジレント・テクノロジー株式会社 | 平衡系特性パラメータ測定装置および測定方法 |
JP3017367B2 (ja) * | 1992-09-04 | 2000-03-06 | 京セラ株式会社 | マルチポート回路のsパラメータ測定方法 |
JPH08160084A (ja) * | 1994-12-01 | 1996-06-21 | Nippon Steel Corp | 高周波用半導体デバイスの特性測定方法およびそのための測定用パターン |
US6188968B1 (en) * | 1998-05-18 | 2001-02-13 | Agilent Technologies Inc. | Removing effects of adapters present during vector network analyzer calibration |
JP2001153904A (ja) * | 1999-11-25 | 2001-06-08 | Advantest Corp | ネットワークアナライザ、ネットワーク解析方法および記録媒体 |
JP2001242207A (ja) * | 2000-02-25 | 2001-09-07 | Anritsu Corp | 高周波デバイス測定プログラムを記録した記憶媒体並びに高周波デバイス測定装置及び方法 |
JP2001244452A (ja) * | 2000-02-28 | 2001-09-07 | Japan Radio Co Ltd | 高周波半導体デバイスモデルパラメータ抽出方法、高周波半導体デバイスモデルパラメータ抽出装置及び記録媒体 |
US6744262B2 (en) * | 2002-03-14 | 2004-06-01 | Agilent Technologies, Inc. | Method, apparatus, and article of manufacture for characterizing a device and predicting electrical behavior of the device in a circuit |
US6757625B2 (en) * | 2002-04-22 | 2004-06-29 | Agilent Technologies, Inc. | Method, apparatus, and article of manufacture for predicting electrical behavior of a multiport device having balanced device ports |
US7113891B2 (en) * | 2004-05-12 | 2006-09-26 | Agilent Technologies, Inc. | Multi-port scattering parameter calibration system and method |
-
2004
- 2004-03-25 JP JP2004088478A patent/JP2005274373A/ja active Pending
-
2005
- 2005-01-24 TW TW094102006A patent/TWI292544B/zh not_active IP Right Cessation
- 2005-01-25 US US11/041,246 patent/US7389191B2/en not_active Expired - Fee Related
- 2005-01-28 KR KR1020050008027A patent/KR100682005B1/ko not_active IP Right Cessation
- 2005-01-31 CN CNB2005100067267A patent/CN100410951C/zh not_active Expired - Fee Related
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107490732A (zh) * | 2016-06-13 | 2017-12-19 | 株式会社村田制作所 | 传输特性测定方法以及传输特性测定装置 |
CN107490732B (zh) * | 2016-06-13 | 2019-11-29 | 株式会社村田制作所 | 传输特性测定方法以及传输特性测定装置 |
CN107943641A (zh) * | 2017-11-16 | 2018-04-20 | 郑州云海信息技术有限公司 | 一种生成多端口s参数文档的方法及装置 |
CN108804775A (zh) * | 2018-05-22 | 2018-11-13 | 郑州云海信息技术有限公司 | 一种基于Matlab的高速信号散射参数计算方法及装置 |
CN108804775B (zh) * | 2018-05-22 | 2022-02-18 | 郑州云海信息技术有限公司 | 一种基于Matlab的高速信号散射参数计算方法及装置 |
CN109376342A (zh) * | 2018-11-01 | 2019-02-22 | 郑州云海信息技术有限公司 | 一种s参数端口转换方法与装置 |
WO2021031424A1 (zh) * | 2019-08-19 | 2021-02-25 | 苏州华太电子技术有限公司 | 散射参数测量方法及器件校准方法 |
CN112182997A (zh) * | 2020-09-04 | 2021-01-05 | 曙光信息产业(北京)有限公司 | 散射参数确定方法、装置、计算机设备和存储介质 |
CN112182997B (zh) * | 2020-09-04 | 2024-04-12 | 曙光信息产业(北京)有限公司 | 散射参数确定方法、装置、计算机设备和存储介质 |
Also Published As
Publication number | Publication date |
---|---|
US20050216239A1 (en) | 2005-09-29 |
TWI292544B (en) | 2008-01-11 |
TW200532528A (en) | 2005-10-01 |
KR20050095553A (ko) | 2005-09-29 |
CN100410951C (zh) | 2008-08-13 |
KR100682005B1 (ko) | 2007-02-15 |
US7389191B2 (en) | 2008-06-17 |
JP2005274373A (ja) | 2005-10-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1674009A (zh) | 用于s参数计算的设备和方法及其程序和记录介质 | |
CN101413973B (zh) | 一种线路板特性阻抗测试系统和方法 | |
CN101201389B (zh) | 利用模拟信号进行装置中扫描链测试的系统与方法 | |
JP2009223885A5 (zh) | ||
CN102723950A (zh) | 一种模数转换器非线性参数测试适配器及测试方法 | |
CN101655882B (zh) | 基于统计模型最差情况的建模方法 | |
CN107345986A (zh) | 一种去嵌入方式的阻抗测试方法 | |
JP4500338B2 (ja) | テストパターン評価方法及びテストパターン評価装置 | |
Cao et al. | A new training approach for robust recurrent neural-network modeling of nonlinear circuits | |
WO2004055529A1 (fr) | Procede de mesure des interferences electromagnetiques et systeme associe | |
CN103809022A (zh) | 使用增强递归离散傅里叶变换的片上专用分析 | |
CN203811126U (zh) | 一种应力环变形检测装置 | |
CN203164407U (zh) | 一种超高频局放状态检测仪器校验装置 | |
CN102062834A (zh) | 半导体射频器件的射频测试方法 | |
CN205749621U (zh) | 一种频谱分析及波形显示的数字示波器 | |
Rumyantsev et al. | Development of a Measurement System-on-Chip and Simulation on FPGA | |
CN111079250B (zh) | 一种电子产品疲劳寿命评估、评估模型建立方法及装置 | |
US9213051B2 (en) | Impedance measurement system, impedance measurement method and program | |
CN102012955B (zh) | 获取背板连接器针头钻阵列的散射矩阵的方法和装置 | |
TW200905208A (en) | System and method for analyzing non-monotonicity of a waveform curve of signals | |
EP2726887B1 (en) | Device and method for emi source location | |
CN109977522B (zh) | 一种可伸缩去嵌模型的形成及去嵌方法 | |
Rodzi et al. | Design and development of a human height measurement system for sports health screening | |
KR101478625B1 (ko) | 헬리칼 구조 전류 프로브, 이를 포함하는 반도체 칩 및 전류 산출 방법 | |
Chu et al. | Linear system DFT deconvolution algorithm based on frequency spectrum tailoring. |
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 |
Granted publication date: 20080813 Termination date: 20180131 |
|
CF01 | Termination of patent right due to non-payment of annual fee |