crosstalk_f.m

该程序用于计算n根平行导线间的串扰,程序结构及说明见程序内。源可以是单频的或脉冲的（这里给出的是频域的）。频率、波形、线长、电源及负载阻抗均可在用MATLAB编写的语句中方便地修改。但必须输入导线系统

% File Name: Crosstalk_f.m
% Calculate Cross talk of 3 parallel conductors at specified frequncy
%Part 1: Initializing data;
  clear
  Ncond=3;                          % Number of conductors;
  Lc=0.1;                           % Length of condutors(m);
  C=[53.232   -3.6051  -23.423;
     -3.6051  45.246   -10.996;
     -23.423  -10.996   55.679]*1E-12;      % Capacitance coeffiecient matrix(F);
       
  L=1.0e-006 *[0.7808    0.2049    0.4120;
               0.2049    0.7862    0.2735;
               0.4120    0.2735    0.7800]; % Inductance coeffiecient matrix (H);
  
  Z0=[50 0  0;      
      0  50 0
      0  0  50];                            % Inner impendance;
  Zl1=50;                                   % Load impendance;
  Zl2=50;
  Zl3=50;
  Zl=[Zl1   0    0;
      0     Zl2  0;
      0     0    Zl3];           
  
   Us=1;                                    % Applied voltage;
   f=1E9;                                   % Frequency;
    
   
 % Part 2: Calculate propagation constsnt ［Γ］and velocity;
  
  [T,D]=eig(C*L)  
  for i=1:Ncond
     Lamuda(i)=D(i,i);    
  end;
  
  w=2*pi*f;
  for k=1:Ncond
    Gama(k,k)=j*w*sqrt(Lamuda(k));
 end;
 for i=1:3
    v(i)=w/imag(Gama(i,i))     
 end;
 
 % Part 3: Calculating crosstalk voltage at start and end points of No.1,3 ;
  
  Y=j*w*C;
  E1p=exp(Lc*Gama);
  E2p=exp(-Lc*Gama);
  for k=1:Ncond
     E1(k,k)=E1p(k,k);
     E2(k,k)=E2p(k,k);
  end;
  
  Ai11=0.5*inv(Y)*T*(E1+E2)*inv(T)*Y;
  Ai12=-0.5*inv(Y)*T*Gama*(E1-E2)*inv(T);
  Ai21=-0.5*T*(E1-E2)*inv(Gama)*inv(T)*Y;
  Ai22=0.5*T*(E1+E2)*inv(T);
  Ai=[Ai11 Ai12; 
      Ai21 Ai22];
  V0=Us*[0 1 0]';       % Source applyied at No.2 conductor;
  I_start=inv(Zl*Ai22-Zl*Ai21*Z0-Ai12+Ai11*Z0)*(Ai11-Zl*Ai21)*V0;
  I_end=(Ai22-Ai21*Z0)*I_start+Ai21*V0;
  V_start=V0-Z0*I_start;
  V_end=Zl*I_end;

  U1_start=V_start(1);        % No.1 conductor;
  U1_end=V_end(1);
  U1s_Amp=abs(U1_start)       % Amplitude (v);
  U1s_Ang=angle(U1_start);    % Angle (rad/s);
  U1e_Amp=abs(U1_end)         % Amplitude (v);
  U1e_Ang=angle(U1_end);
   
  U3_start=V_start(3);        % No.3 conductor;
  U3_end=V_end(3);      
  U3s_Amp=abs(U3_start)       % Amplitude (v);
  U3s_Ang=angle(U3_start);    % Angle (rad/s);
  U3e_Amp=abs(U3_end)         % Amplitude (v);
  U3e_Ang=angle(U3_end);

  % The end;