怎么破解PSPICE里面器件模型，例如传输线模型，有高手支招吗？

1945007625

这个是理想传输线的Pspice等效电路，但是我看NGSPICE里面的理想传输线模型不是节点矩阵法实现的，下面是传输线源程序
) V& M+ A3 U: h9 j1 ]5 s7 C8 F' d#include "cktdefs.h"
#include "tradefs.h"
0 W% K' v) Z4 Z2 r#include "sperror.h"
& P+ P$ l# o+ G#include "suffix.h"
1 R4 Q7 ~' g7 E7 r3 A- m
3 a& \3 Q# q' B* u$ u
& Y1 }5 _8 Y4 W6 u" N, Kint
9 O# {7 P7 D1 r9 v; Y% yTRAaccept(ckt,inModel)
    register CKTcircuit *ckt;
    GENmodel *inModel;
8 \+ o! Y5 N/ E1 ^1 c{
1 j8 ]# C* R) g6 \    register TRAmodel *model = (TRAmodel *)inModel;
    register TRAinstance *here;
    register int i=0,j;
    double v1,v2,v3,v4;
    double v5,v6,d1,d2,d3,d4;
8 n' g6 @$ p. V( ]    double *from,*to;
6 d, Q6 D/ k( Q7 d    int error;
, B0 C8 f: P5 Y+ J6 A* [9 g. C4 \! B. h

! T0 {0 w; F& A& K7 d& N5 T' I& t    /*  loop through all the transmission line models */
! z' d% G; k- K, D! \    for( ; model != NULL; model = model->TRAnextModel ) {
! I: H6 ~7 _8 S+ J* D8 Z# J
2 \2 q$ x6 q  V% |9 \; ?        /* loop through all the instances of the model */
        for (here = model->TRAinstances; here != NULL ;
' D2 ^; [6 T4 J/ h" y5 h; s) |                here=here->TRAnextInstance) {
2 d6 _( Y" D( M% Z3 f5 W6 [            if( (ckt->CKTtime - here->TRAtd) > *(here->TRAdelays+6)) {
% s3 x: y' _, n- S% c. f                /* shift! */
                for(i=2;i<here->TRAsizeDelay &&
. b" `* N& X0 n  @6 G9 r( ~$ D                    (ckt->CKTtime - here->TRAtd > *(here->TRAdelays+3*i));i++)
3 r# O5 O8 j- {* ]; O$ {5 J                        { /* loop does it all */ ; }
                i -= 2;
                for(j=i;j<=here->TRAsizeDelay;j++) {
                    from = here->TRAdelays + 3*j;
! G5 C. U3 h' q) J8 t, V( R                    to = here->TRAdelays + 3*(j-i);
                    *(to) = *(from);
                    *(to+1) = *(from+1);
" M. x% W. ]3 y- e                    *(to+2) = *(from+2);
                }
3 r6 M: M; u' o+ F: F9 C: N) y                here->TRAsizeDelay -= i;
            }
            if(ckt->CKTtime - *(here->TRAdelays+3*here->TRAsizeDelay) >
                    ckt->CKTminBreak) {
                if(here->TRAallocDelay <= here->TRAsizeDelay) {
/ Z" z/ v  |3 _7 h) ?" f2 g8 \  z2 I( X                    /* need to grab some more space */
                    here->TRAallocDelay += 5;
                    here->TRAdelays = (double *)REALLOC((char *)here->TRAdelays,
                        (here->TRAallocDelay+1)*3*sizeof(double));
                }
2 I& k! W9 N4 m" S) g2 T; C                here->TRAsizeDelay ++;
                to = (here->TRAdelays  +3*here->TRAsizeDelay);
                *to = ckt->CKTtime;
                to = (here->TRAdelays+1+3*here->TRAsizeDelay);
                *to = ( *(ckt->CKTrhsOld + here->TRAposNode2)
2 p8 t" B6 `2 l: d3 ~- {/ V                     -*(ckt->CKTrhsOld + here->TRAnegNode2))
: j* q9 J' N& Y; R' k. O                    + *(ckt->CKTrhsOld + here->TRAbrEq2)*
                        here->TRAimped;
                *(here->TRAdelays+2+3*here->TRAsizeDelay) =
( T: D. ^3 ~6 e" @5 G7 o* T                    ( *(ckt->CKTrhsOld + here->TRAposNode1)
) @% |5 k! ~4 k( m                     -*(ckt->CKTrhsOld + here->TRAnegNode1))
                    + *(ckt->CKTrhsOld + here->TRAbrEq1)*
                        here->TRAimped;
#ifdef NOTDEF
                v1 = *(here->TRAdelays+1+3*here->TRAsizeDelay);
                v2 = *(here->TRAdelays+1+3*(here->TRAsizeDelay-1));
6 F$ v  P% x% S7 D                v3 = *(here->TRAdelays+2+3*here->TRAsizeDelay);
  A& k( N+ ?; c# _! U$ F                v4 = *(here->TRAdelays+2+3*(here->TRAsizeDelay-1));
$ p! V) V$ H% X7 y& Y& X- s$ b" {% V6 L' t                if( (FABS(v1-v2) >= 50*ckt->CKTreltol*
$ g4 P$ x: m0 z, x  R                        MAX(FABS(v1),FABS(v2))+50*ckt->CKTvoltTol) ||
! A/ g* i) F- r  o                    (FABS(v3-v4) >= 50*ckt->CKTreltol*
                        MAX(FABS(v3),FABS(v4))+50*ckt->CKTvoltTol)  ) {
$ L: t, x+ Y  G0 {2 a/ N                    /* changing - need to schedule after delay */
                    /*printf("%s:  at %g set for %g and %g\n",here->TRAname,
: S3 A% p4 ?2 ]- O3 V0 ^                        ckt->CKTtime,
                        ckt->CKTtime+here->TRAtd,
4 B8 Z+ ?  E& F$ x# ?9 K5 {0 y                        *(here->TRAdelays+3*here->TRAsizeDelay-3)+
6 x* L5 D1 a) I# h                                here->TRAtd);*/
                    error = CKTsetBreak(ckt,ckt->CKTtime+here->TRAtd);
                    if(error) return(error);
                    /* also set for break after PREVIOUS point */
# O8 e+ i; m$ M9 S& z+ U                    error = CKTsetBreak(ckt,
8 G: `/ s) S3 O- r) j                            *(here->TRAdelays+3*here->TRAsizeDelay -3) +
+ A% m9 y# A( G% N                            here->TRAtd);
8 \0 O9 K9 l4 F8 F                    CKTbreakDump(ckt);
* H0 c  r& E4 @" L" J                    if(error) return(error);
                }
3 [& M3 Q8 D% w0 ]3 y( F, G- v#else
$ Z! o7 \9 v5 e& e. j7 b' j4 _                v1 = *(here->TRAdelays+1+3*here->TRAsizeDelay);
( U* Y% w! U+ U# u! {                v2 = *(here->TRAdelays+1+3*(here->TRAsizeDelay-1));
5 [( K* U7 u1 ~: X( T; I                v3 = *(here->TRAdelays+1+3*(here->TRAsizeDelay-2));
* N6 R8 e/ E. s- [8 Z( S8 i4 |: \                v4 = *(here->TRAdelays+2+3*here->TRAsizeDelay);
; k; p; t7 a5 S& M                v5 = *(here->TRAdelays+2+3*(here->TRAsizeDelay-1));
                v6 = *(here->TRAdelays+2+3*(here->TRAsizeDelay-2));
                d1 = (v1-v2)/ckt->CKTdeltaOld[0];
                d2 = (v2-v3)/ckt->CKTdeltaOld[1];
5 \. W9 v9 \$ t/ b( b                d3 = (v4-v5)/ckt->CKTdeltaOld[0];
                d4 = (v5-v6)/ckt->CKTdeltaOld[1];
. ~1 c. }) Q7 Q0 ~                /*printf("%s: at %g derivs are %g, %g and %g, %g\n",
                        here->TRAname,ckt->CKTtime,d1,d2,d3,d4);*/
                if( (FABS(d1-d2) >= here->TRAreltol*MAX(FABS(d1),FABS(d2))+
5 k) w2 {  k; \5 A7 M% u                        here->TRAabstol) ||
/ L5 ]( S% S% O  i  j; v1 {                        (FABS(d3-d4) >= here->TRAreltol*MAX(FABS(d3),FABS(d4))+
                        here->TRAabstol) ) {
+ l6 d+ j. [! B1 g+ o                    /* derivitive changing - need to schedule after delay */
                    /*printf("%s:  at %g set for %g\n",here->TRAname,
5 z2 @( e8 l' E0 z, T                        ckt->CKTtime,
) f5 R+ X* h, }: v" h                        *(here->TRAdelays+3*here->TRAsizeDelay-3)+here->TRAtd
                        );*/
                    /*printf("%g, %g, %g -> %g, %g \n",v1,v2,v3,d1,d2);*/
  s# c2 T0 o' p+ @8 U                    /*printf("%g, %g, %g -> %g, %g \n",v4,v5,v6,d3,d4);*/
, p  c& ^6 U+ L; ^                    /* also set for break after PREVIOUS point */
                    /*printf("setting break\n");*/
1 s  `7 N5 U1 p$ q7 `* N                    error = CKTsetBreak(ckt,
& \" _* |) O- y8 i                            *(here->TRAdelays+3*here->TRAsizeDelay -3) +
                            here->TRAtd);
                    if(error) return(error);
                }
6 f) K5 F7 H% D" U3 R6 B7 q#endif /*NOTDEF*/
( ^5 |, ^8 Y" F, c            }
        }
! y0 R: S! S; l. p/ d+ K    }
    return(OK);
}

* A* s) C7 o6 T; F1 T
$ e8 H# p" K4 P$ K

芦根苏木

我也想知道，看看别人怎么说
2 {* o8 V8 t6 v8 u! f