CN102346233B - System and method for analyzing scattering parameter passivity - Google Patents

Abstract

The invention provides a system and a method for analyzing scattering parameter passivity, which can be applied to a calculating device. The calculating device is connected with a measuring instrument which is used for measuring a signal transmitted on a circuit to obtain a scattering parameter file. The system comprises a parameter reading module, a vector fitting module, a matrix conversion module and a passivity analysis module, wherein the system reads the scattering parameter file by the modules; each scattering parameter is respectively subjected to vector fitting to generate a non-common extreme-value form rational function matrix of the scattering parameter; according to the rational function matrix, an equivalent circuit is generated, and the rational function matrix is converted into a state space matrix; finally, the state space matrix obtained by conversion is substituted into a Hamilton matrix by the system; and whether the characteristic value of the Hamilton matrix comprises a pure imaginary numerical value is analyzed to judge whether the rational function matrix satisfies the passivity or not so as to judge whether the equivalent circuit satisfies the passivity requirement or not.

Description

Scattering parameter passivity analysis system and method

Technical field

The present invention relates to a kind of breadboardin system and method, especially about a kind of scattering parameter passivity analysis system and method.

Background technology

In low-frequency channel, the size of electronic component (for example transmission line) can be ignored for the wavelength of signal.But in high-frequency microwave circuit, because wavelength is shorter, the size of electronic devices and components just cannot be ignored.Microwave Net method is equivalent to reactance or resistance device by microwave component, actual guided wave transmission system is equivalent to transmission line, thereby actual microwave system is reduced to Microwave Net.The electronic component that can not produce power supply such as reactance, resistance etc. is called passivity electronic component, and the circuit that only includes passive electronic components is called passivity circuit.

It is one of basic means of Microwave circuit analysis that scattering parameter (scattering parameters is called for short S parameter) is measured.S Parametric Representation be the signal relation between each port in circuit, as reflection, loss, crosstalk etc., and be used to the behavior of simulation electronic components and parts under different frequency.At present, in high-frequency microwave circuit design, designer is to the S parameter measuring from circuit, utilize vectorial matching (vector fitting) algorithm to produce the rational function (rational function approximation model) of S parameter, according to this rational function, produce equivalent-circuit model, then carry out time domain transient analysis after this circuit model is connected with other Circuits System.If equivalent-circuit model does not meet passivity, can cause simulation software cannot reach the problem of convergence in time domain solution procedure.In order to guarantee that equivalent-circuit model also meets passivity, the rational function of S parameter must meet passivity.

Summary of the invention

In view of above content, be necessary to propose a kind of scattering parameter passivity analysis system, can carry out passivity analysis to the equivalent electrical circuit producing according to scattering parameter, to judge whether equivalent electrical circuit meets the requirements.

In addition, be also necessary to propose a kind of scattering parameter passivity analysis method, can carry out passivity analysis to the equivalent electrical circuit producing according to scattering parameter, to judge whether equivalent electrical circuit meets the requirements.

A scattering parameter passivity analysis system, is applied to calculation element, and this calculation element is connected with measuring instrument, and the signal transmitting in this measuring instrument measurement circuit, obtains scattering parameter file.This system comprises parameter read module, Vector Quasi compound module, matrix conversion module and passivity analysis module.Parameter read module, for reading scattering parameter file, this scattering parameter file is included in the scattering parameter measuring from each port of circuit under different frequency.Vector Quasi compound module, for each scattering parameter being carried out respectively to the Rational Matrix of the non-common extreme value form of vectorial matching generation scattering parameter, produces equivalent electrical circuit according to this Rational Matrix.Matrix conversion module, for being converted to state space matrices by the Rational Matrix of the non-common extreme value form of scattering parameter.Passivity analysis module, for the state space matrices substitution Hamiltonian Matrices being converted to is analyzed, according to the eigenwert of Hamiltonian Matrices, whether comprise pure imaginary number value judges whether the Rational Matrix of the non-common extreme value form of scattering parameter meets passivity, thereby judgement equivalent electrical circuit meets passivity requirement.

A scattering parameter passivity analysis method, is applied to calculation element, and this calculation element is connected with measuring instrument, and the signal transmitting in this measuring instrument measurement circuit, obtains scattering parameter file.The method comprises: (A) read scattering parameter file, this scattering parameter file is included in the scattering parameter measuring from each port of circuit under different frequency; (B) each scattering parameter is carried out respectively to the Rational Matrix that vectorial matching produces the non-common extreme value form of scattering parameter, according to this Rational Matrix, produce equivalent electrical circuit; (C) Rational Matrix of the non-common extreme value form of scattering parameter is converted to state space matrices; (F) the state space matrices substitution Hamiltonian Matrices being converted to is analyzed, according to the eigenwert of Hamiltonian Matrices, whether comprise pure imaginary number value judges whether the Rational Matrix of the non-common extreme value form of scattering parameter meets passivity, thereby judge whether equivalent electrical circuit meets passivity requirement.

Compared to prior art, scattering parameter passivity analysis system and method provided by the present invention, can carry out passivity analysis to the equivalent electrical circuit producing according to scattering parameter, to judge whether equivalent electrical circuit meets the requirements.

Accompanying drawing explanation

Fig. 1 is the applied environment figure of scattering parameter passivity analysis system of the present invention preferred embodiment.

Fig. 2 is the functional block diagram of scattering parameter passivity analysis system preferred embodiment in Fig. 1.

Fig. 3 is the process flow diagram of scattering parameter passivity analysis method of the present invention preferred embodiment.

Fig. 4 is the schematic diagram of four port difference transmission lines in circuit.

Fig. 5 is that the scattering parameter of four ports in Fig. 4 is through the vectorial matching schematic diagram of 8 pairs of extreme value-residual values.

Fig. 6 is that four port scattering parameters shown in Fig. 4 utilize the present invention to analyze the schematic diagram of the pure imaginary number eigenwert of resulting Hamiltonian matrix.

Main element symbol description

Circuit	10
		Measuring instrument	20
Calculation element	30
		S parameter passivity analysis system	31
S Parameter File	32
		Processor	33
Storer	34
		Parameter read module	311
Vector Quasi compound module	312
		Matrix conversion module	313
Passivity analysis module	314
		Equivalent electrical circuit	315
Difference transmission lines	L 1、L2
		Port	1、2、3、4

Embodiment

Consulting shown in Fig. 1, is the applied environment figure of scattering parameter of the present invention (scattering parameters is called for short S parameter) passivity analysis system 31 (hereinafter to be referred as system 31) preferred embodiment.This system 31 is applied to calculation element 30.This calculation element 30 also comprises processor 33 and storer 34.As shown in Figure 1, this calculation element 30 is connected with measuring instrument 20, and this measuring instrument 20, for the signal of transmission in measurement circuit 10, obtains scattering parameter (scattering parameters, be called for short S parameter) file 32, and this S Parameter File 32 is stored in to storer 34.

This system 31 is for producing the Rational Matrix of S parameter to this S Parameter File 32, the Rational Matrix of S parameter is converted to state space matrices, and bring this state space matrices into Hamiltonian matrix and analyze, according to the eigenwert of Hamiltonian matrix, whether comprise pure imaginary number value judges whether the Rational Matrix of S parameter meets passivity, thereby whether the design of decision circuitry 10 meets passivity requirement.

Described storer 34 is also for the sequencing code of storage system 31.Processor 33 is carried out described sequencing code, and the above-mentioned functions of system 31 is provided.

In the present embodiment, the passivity circuit that this circuit 10 is design.For example, this circuit 10 may comprise difference transmission lines L1 and L2 as shown in Figure 4.This measuring instrument 20 is network analyzer.Described calculation element 30 can be PC, notebook, and server, workstation, or other has the electronic installation of data processing function.

Consulting shown in Fig. 2, is the functional block diagram of system 31.This system 31 comprises parameter read module 311, Vector Quasi compound module 312, matrix conversion module 313 and passivity analysis module 314.Function below in conjunction with the functional module of the method flow illustrative system 31 shown in Fig. 3.

As shown in Figure 3, be the process flow diagram of S parameter passivity analysis method of the present invention preferred embodiment.

Step S301, parameter read module 311 reads S Parameter File 32, and this original S Parameter File 32 is included in the S parameter value measuring from circuit 10 each ports under different frequency.For example, suppose that circuit 10 comprises difference transmission lines L1 and L2 as shown in Figure 4, this difference transmission lines L1 and L2 comprise four port ones, 2,3,4.When from a signal of port one transmitting, can reflect a signal from port one simultaneously, port 2,3,4 all can receive a signal, and the signal that port one reflects is denoted as S11 with the ratio of the signal that port one is launched, and is called again reflection parameters; The signal that port 2 receives is denoted as S12 with the ratio of the signal that port one is launched, and is called again insertion loss parameter; The signal that port 3 receives is denoted as S13 with the ratio of the signal that port one is launched, and is called again near-end cross parameter; The signal that port 4 receives is denoted as S14 with the ratio of the signal that port one is launched, and is called again far-end cross talk parameter.All these parameters include amplitude and phase place, and amplitude curve can change along with the variation of the frequency transmitting.Therefore, in the present embodiment, S parameter comprises: reflection parameters, insertion loss parameter, near-end cross parameter and far-end cross talk parameter etc.

In general, the rational function of the S parameter that vector matching produces has (1a) and (1b) two kinds of forms, wherein (1a) is common extreme value (common-pole) form, is (1b) non-common extreme value (non-common-pole) form:

$S\left(s\right)\≈\left(\underset{m=1}{\overset{M}{\mathrm{\Σ}}}\frac{\left[\begin{array}{cccc}{r}_m^{\mathrm{1,1}}& r_m^{\mathrm{1,2}}& ...& r_m^{1,N}\\ r_m^{\mathrm{2,1}}& r_m^{\mathrm{2,2}}& ...& r_m^{2,N}\\ ...& ...& ...& ...\\ r_m^{N,1}& r_m^{N,2}& ...& r_m^{N,N}\end{array}\right]}{s+p_m}\right)+\left[\begin{array}{cccc}{d}_m^{\mathrm{1,1}}& d_m^{\mathrm{1,2}}& ...& d_m^{1,N}\\ d_m^{\mathrm{2,1}}& d_m^{\mathrm{2,2}}& ...& d_m^{2,N}\\ ...& ...& ...& ...\\ d_m^{N,1}& d_m^{N,2}& ...& d_m^{N,N}\end{array}\right]---\left(1a\right)$

$S\left(s\right)\≈\left(\underset{m=1}{\overset{M}{\mathrm{\Σ}}}\left[\begin{array}{cccc}\frac{r_m^{\mathrm{1,1}}}{s+p_m^{1,1}}& \frac{r_m^{\mathrm{1,2}}}{s+p_m^{\mathrm{1,2}}}& ...& \frac{r_m^{1,N}}{s+p_m^{1,N}}\\ \frac{r_m^{\mathrm{2,1}}}{s+p_m^{\mathrm{2,1}}}& \frac{r_m^{\mathrm{2,2}}}{s+p_m^{\mathrm{2,2}}}& ...& \frac{r_m^{2,N}}{s+p_m^{2,N}}\\ ...& ...& ...& ...\\ \frac{r_m^{N,1}}{s+p_m^{N,1}}& \frac{r_m^{N,2}}{s+p_m^{N,2}}& ...& \frac{r_m^{N,N}}{s+p_m^{N,N}}\end{array}\right]\right)+\left[\begin{array}{cccc}{d}_m^{\mathrm{1,1}}& d_m^{\mathrm{1,2}}& ...& d_m^{1,N}\\ d_m^{\mathrm{2,1}}& d_m^{\mathrm{2,2}}& ...& d_m^{2,N}\\ ...& ...& ...& ...\\ d_m^{N,1}& d_m^{N,2}& ...& d_m^{N,N}\end{array}\right]---\left(1b\right)$

Wherein, M represents control accuracy, and N represents port number, r _mrepresent residual value, p _mrepresent extreme value, s=ω=2 π f represents angular frequency, d _mfor constant.

As shown in Figure 5, be the S parameter of four ports through 8 the vectorial fitting result to (M=8) extreme value-residual value.Transverse axis in Fig. 5 represents frequency, and the longitudinal axis represents amplitude.By Fig. 5, can be learnt, use the rational function curve of the common extreme value form (By Common pole) of the S parameter that vectorial matching produces to depart from the original amplitude curve of S parameter far away, and use the rational function of the non-common extreme value form (By Non-common pole) of the S parameter that vectorial matching produces more to approach the original amplitude curve (Original System) of S parameter.That is to say, use the rational function error of common extreme value form of the S parameter that vectorial matching produces larger, use the rational function precision of non-common extreme value form of the S parameter that vectorial matching produces higher.Therefore the rational function that, the present invention is directed to non-common extreme value form is analyzed.

Step S303, Vector Quasi compound module 312 carries out respectively to each S parameter the Rational Matrix that vectorial matching produces the non-common extreme value form of S parameter, produces the equivalent electrical circuit 315 of circuit 10 according to this Rational Matrix.

Shown in (2), be that S parameter obtains the Rational Matrix of M to the non-common extreme value form of extreme value-residual value through vectorial matching as follows:

$S\left(s\right)\≈\hat{S}\left(s\right)=\left[\begin{array}{cccc}{\hat{S}}_{11}\left(s\right)& {\hat{S}}_{12}\left(s\right)& ...& {\hat{S}}_{1N}\left(s\right)\\ {\hat{S}}_{21}\left(s\right)& {\hat{S}}_{22}\left(s\right)& ...& {\hat{S}}_{2N}\left(s\right)\\ ...& ...& ...& ...\\ {\hat{S}}_{N1}\left(s\right)& {\hat{S}}_{N2}\left(s\right)& ...& {\hat{S}}_{\mathrm{NN}}\left(s\right)\end{array}\right]---\left(2\right)$

Wherein, represent extreme value, its real number (real) form as shown in (3a), represent residual value, its plural number (complex) form is as shown in (3b):

$\hat{S}{r}_{p,q}\left(s\right)=\underset{u=1}{\overset{U}{\mathrm{\Σ}}}\frac{r_u^{p,q}}{s+p_u^{p,q}}---\left(3a\right)$

Wherein, and for real number.

$\hat{S}{c}_{p,q}\left(s\right)=\underset{v=1}{\overset{V}{\mathrm{\Σ}}}\frac{\mathrm{Re}\left(r_v^{p,q}\right)+\mathrm{Im}\left(r_v^{p,q}\right)j}{s+\mathrm{Re}\left(p_v^{p,q}\right)+\mathrm{Im}\left(p_v^{p,q}\right)j}+\frac{\mathrm{Re}\left(r_v^{p,q}\right)-\mathrm{Im}\left(r_v^{p,q}\right)j}{s+\mathrm{Re}\left(p_v^{p,q}\right)-\mathrm{Im}\left(p_v^{p,q}\right)j}---\left(3b\right)$

Wherein, u+2V=M, and

Step S305, matrix conversion module 313 is converted to state space matrices by the Rational Matrix of the non-common extreme value form of S parameter.For example, matrix conversion module 313 is converted to the state space matrices shown in following formula (4) by the Rational Matrix shown in formula (2) in conjunction with (3a) and expression way (3b):

jωX(jω)＝AX(jω)+BU(jω) (4)

Y(jω)＝CX(jω+DU(jω)，

Wherein, A, B, C and D matrix represent with formula (5) respectively:

$A=\left[\begin{array}{cc}{A}_r& 0\\ 0& A_c\end{array}\right],$ $B=\left[\begin{array}{c}{B}_r\\ B_c\end{array}\right],$ C＝[C _r C _c]，

$D=\begin{array}{c}\left[\begin{array}{cccc}{d}^{\mathrm{1,1}}& d^{\mathrm{1,2}}& ...& d^{1,N}\\ d^{\mathrm{2,1}}& d^{\mathrm{2,2}}& ...& d^{2,N}\\ ...& ...& ...& ...\\ d^{N,1}& d^{N,1}& ...& d^{N,N}\end{array}\right]---\left(5\right)\end{array}$

A in formula (5) _r, B _rand C _rspace matrix for the real number form of extreme value-residual value in formula (3a), can represent by formula (6):

$A_r(i,j)=\left\{\begin{array}{cc}{p}_u^{p,q},& i=\mathrm{j\; and\; i}=(U\·N)(p-1)+U(q-1)+u\\ 0& o.w.\end{array}\right.$

$B_r(i,j)=\left\{\begin{array}{cc}1,& i=(U\·N)(p-1)+U(q-1)+\mathrm{u\; and\; j}=q\\ 0& o.w.\end{array}\right.$

$C_r(i,j)=\left\{\begin{array}{cc}{r}_u^{p,q},& i=\mathrm{p\; and\; j}=(U\·N)(p-1)+U(q-1)+u\\ 0& o.w.\end{array}\right.---\left(6\right)$

A wherein _rfor (NNU) * (NNU), B _rfor (NU) * N, C _rfor the sparse matrix of N * (NU), o.w represents other condition (otherwise).

A in formula (5) _c, B _cand C _cspace matrix for the plural form of extreme value-residual value in formula (3b), can represent by formula (7):

$A_c(i,j)=\left\{\begin{array}{cc}\mathrm{Re}\left(p_v^{p,q}\right),& i=\mathrm{j\; and\; i}=\mathrm{\Ψ}(p,q,v),i=\mathrm{\Ψ}(p,q,v)-1\\ \mathrm{Im}\left(p_v^{p,q}\right),& i=\mathrm{\Ψ}(p,q,v)-1,j=\mathrm{\Ψ}(p,q,v)\\ -\mathrm{Im}\left(p_v^{p,q}\right),& i=\mathrm{\Ψ}(p,q,v),j=\mathrm{\Ψ}(p,q,v)-1\\ 0,& o.w.\end{array}\right.$

$B_c(i,j)=\left\{\begin{array}{cc}2,& i=\mathrm{\Ψ}(p,q,v)-1\mathrm{and\; j}=q\\ 0& o.w.\end{array}\right.$

$C_c(i,j)=\left\{\begin{array}{cc}\mathrm{Re}\left(r_v^{p,q}\right),& i=\mathrm{p\; and\; j}=\mathrm{\Ψ}-1\\ \mathrm{Im}\left(r_v^{p,q}\right),& i=\mathrm{p\; and\; j}=\mathrm{\Ψ}\\ 0,& o.w.\end{array}\right.---\left(7\right)$

Wherein, (p-1)+2V (the q-1)+2v of Ψ (p, q, v)=(2VN), A _cfor (NN2V) * (NN2V), B _rfor (N2V) * N, C _rfor the sparse matrix of N * (N2V).

Matrix conversion module 313, in conjunction with above-mentioned formula (5), (6), (7), can obtain the A in the state space matrices shown in formula (4), B, C, the expression formula of D.

Step S307, passivity analysis module 314 is state space matrices substitution Hamilton (Hamiltonian) matrix being converted to, the Hamiltonian matrix H as shown in (8):

$H=\left[\begin{array}{cc}A-{\mathrm{BR}}^{-1}{D}^{T}C& -{\mathrm{BR}}^{-1}{B}^T\\ C^T{Q}^{-1}C& -A^T+C^T{\mathrm{DR}}^{-1}{B}^T\end{array}\right]---\left(8\right)$

Wherein, R=D ^td-IQ=DD ^t-I, I is unit matrix.

Step S309,314 pairs of Hamiltonian matrixes of passivity analysis module are analyzed, and according to the eigenwert of Hamiltonian matrix, whether comprise that pure imaginary number value judges whether the Rational Matrix of the non-common extreme value form of S parameter meets passivity.

For example, if the Hamiltonian matrix H eigenwert of passivity analysis module 314 judgements as shown in (8) do not comprise pure imaginary number value, the Rational Matrix of the non-common extreme value form of the S parameter shown in passivity analysis module 314 judgment formulas (2) meets passivity.Otherwise, if the Hamiltonian matrix H eigenwert of passivity analysis module 314 judgements as shown in (8) comprises pure imaginary number value, the Rational Matrix of the non-common extreme value form of the S parameter shown in passivity analysis module 314 judgment formulas (2) does not meet passivity.

As shown in Figure 6, for four port S parameters shown in Fig. 4 obtain after the Rational Matrix of the non-common extreme value form as shown in formula (2) the vectorial matching of (M=8) extreme value-residual value through 8, in conjunction with substitution formula (8) after formula (3) to (7), the pure imaginary number eigenwert of resulting Hamiltonian matrix.Wherein,

Step S311, whether passivity analysis module 314 meets passivity according to the Rational Matrix of the non-common extreme value form of S parameter and judges whether equivalent electrical circuit 315 meets passivity requirement.

For example, if the Rational Matrix of the non-common extreme value form of S parameter does not meet passivity, passivity analysis module 314 judgement equivalent electrical circuit 315 do not meet passivity requirement, and the vectorial fitting result of Vector Quasi compound module 312 needs to adjust.If the Rational Matrix of the non-common extreme value form of S parameter meets passivity, equivalent electrical circuit 315 meets passivity requirement, can be applied to subsequent conditioning circuit sunykatuib analysis, for example, be applied to circuit simulation software Hspice and carry out time-domain analysis.

Claims (8)

1. a scattering parameter passivity analysis system, is applied to calculation element, and this calculation element is connected with measuring instrument, and the signal transmitting in this measuring instrument measurement circuit, obtains scattering parameter file, it is characterized in that, this system comprises:

Parameter read module, for reading scattering parameter file, this scattering parameter file is included in the scattering parameter measuring from each port of circuit under different frequency;

Vector Quasi compound module, for each scattering parameter being carried out respectively to the Rational Matrix of the non-common extreme value form of vectorial matching generation scattering parameter, produces equivalent electrical circuit according to this Rational Matrix, and described non-common extreme value form is:

$S\left(s\right)\≈\left(\begin{array}{c}\underset{m=1}{\overset{M}{\mathrm{\Σ}}}\left[\begin{array}{cccc}\frac{r_m^{\mathrm{1,1}}}{s+p_m^{\mathrm{1,1}}}& \frac{r_m^{\mathrm{1,2}}}{s+p_m^{\mathrm{1,2}}}& ...& \frac{r_m^{1,N}}{s+p_m^{1,N}}\\ \frac{r_m^{\mathrm{2,1}}}{s+p_m^{\mathrm{2,1}}}& \frac{r_m^{\mathrm{2,2}}}{s+p_m^{\mathrm{2,2}}}& ...& \frac{r_m^{2,N}}{s+p_m^{2,N}}\\ ...& ...& ...& ...\\ \frac{r_m^{N,1}}{s+p_m^{N1,}}& \frac{r_m^{N,2}}{s+p_m^{N,2}}& ...& \frac{r_m^{N,N}}{s+p_m^{N,N}}\end{array}\right]\end{array}\right)+\left[\begin{array}{cccc}{d}_m^{\mathrm{1,1}}& d_m^{\mathrm{1,2}}& ...& d_m^{1,N}\\ d_m^{\mathrm{2,1}}& d_m^{\mathrm{2,2}}& ...& d_m^{2,N}\\ ...& ...& ...& ...\\ d_m^{N,1}& d_m^{N,2}& ...& d_m^{N,N}\end{array}\right]$

Wherein, M represents control accuracy, and N represents port number, r _mrepresent residual value, p _mrepresent extreme value, s=ω=2 π f represents angular frequency, d _mfor constant;

Matrix conversion module, for being converted to state space matrices by the Rational Matrix of the non-common extreme value form of scattering parameter;

Passivity analysis module, for the state space matrices substitution Hamiltonian Matrices being converted to is analyzed, according to the eigenwert of Hamiltonian Matrices, whether comprise pure imaginary number value judges whether the Rational Matrix of the non-common extreme value form of scattering parameter meets passivity, to judge whether equivalent electrical circuit meets passivity requirement.

2. scattering parameter passivity analysis system as claimed in claim 1, it is characterized in that, if the eigenwert of Hamiltonian Matrices comprises pure imaginary number value, the Rational Matrix of the non-common extreme value form of passivity analysis module judgement scattering parameter does not meet passivity, and then judgement equivalent electrical circuit does not meet passivity requirement, if the eigenwert of Hamiltonian Matrices does not comprise pure imaginary number value, the Rational Matrix of the non-common extreme value form of passivity analysis module judgement scattering parameter meets passivity, and then judgement equivalent electrical circuit meets passivity requirement.

3. scattering parameter passivity analysis system as claimed in claim 1, is characterized in that, described calculation element is PC, server, or workstation.

4. scattering parameter passivity analysis system as claimed in claim 1, is characterized in that, described measuring instrument is network analyzer.

5. a scattering parameter passivity analysis method, is applied to calculation element, and this calculation element is connected with measuring instrument, and the signal transmitting in this measuring instrument measurement circuit, obtains scattering parameter file, it is characterized in that, the method comprises:

Read scattering parameter file, this scattering parameter file is included in the scattering parameter measuring from each port of circuit under different frequency;

Each scattering parameter is carried out respectively to the Rational Matrix that vectorial matching produces the non-common extreme value form of scattering parameter, according to this Rational Matrix, produce equivalent electrical circuit, described non-common extreme value form is:

$S\left(s\right)\≈\left(\begin{array}{c}\underset{m=1}{\overset{M}{\mathrm{\Σ}}}\left[\begin{array}{cccc}\frac{r_m^{\mathrm{1,1}}}{s+p_m^{\mathrm{1,1}}}& \frac{r_m^{\mathrm{1,2}}}{s+p_m^{\mathrm{1,2}}}& ...& \frac{r_m^{1,N}}{s+p_m^{1,N}}\\ \frac{r_m^{\mathrm{2,1}}}{s+p_m^{\mathrm{2,1}}}& \frac{r_m^{\mathrm{2,2}}}{s+p_m^{\mathrm{2,2}}}& ...& \frac{r_m^{2,N}}{s+p_m^{2,N}}\\ ...& ...& ...& ...\\ \frac{r_m^{N,1}}{s+p_m^{N1,}}& \frac{r_m^{N,2}}{s+p_m^{N,2}}& ...& \frac{r_m^{N,N}}{s+p_m^{N,N}}\end{array}\right]\end{array}\right)+\left[\begin{array}{cccc}{d}_m^{\mathrm{1,1}}& d_m^{\mathrm{1,2}}& ...& d_m^{1,N}\\ d_m^{\mathrm{2,1}}& d_m^{\mathrm{2,2}}& ...& d_m^{2,N}\\ ...& ...& ...& ...\\ d_m^{N,1}& d_m^{N,2}& ...& d_m^{N,N}\end{array}\right]$

Wherein, M represents control accuracy, and N represents port number, r _mrepresent residual value, p _mrepresent extreme value, s=ω=2 π f represents angular frequency, d _mfor constant;

The Rational Matrix of the non-common extreme value form of scattering parameter is converted to state space matrices;

The state space matrices substitution Hamiltonian Matrices being converted to is analyzed, according to the eigenwert of Hamiltonian Matrices, whether comprise pure imaginary number value judges whether the Rational Matrix of the non-common extreme value form of scattering parameter meets passivity, to judge whether equivalent electrical circuit meets passivity requirement.

6. scattering parameter passivity analysis method as claimed in claim 5, it is characterized in that, if the eigenwert of Hamiltonian Matrices comprises pure imaginary number value, the Rational Matrix that judges the non-common extreme value form of scattering parameter does not meet passivity, and then judgement equivalent electrical circuit does not meet passivity requirement, if the eigenwert of Hamiltonian Matrices does not comprise pure imaginary number value, the Rational Matrix that judges the non-common extreme value form of scattering parameter meets passivity, and then judgement equivalent electrical circuit meets passivity requirement.

7. scattering parameter passivity analysis method as claimed in claim 5, is characterized in that, described calculation element is PC, server, or workstation.

8. scattering parameter passivity analysis method as claimed in claim 5, is characterized in that, described measuring instrument is network analyzer.