sel.m

EDA程序： 将spice网表转化为用于计算电路各种参数所需要的相应稀疏矩阵。

function [matall] = SEL(n,matall,num1,num2,value)
 %n=1:R       2:M 3:C  4:I
 if (n==1)  
         if num1~=0
         matall.G(num1,num1)=matall.G(num1,num1)+1/value;
         end
         if num2~=0
         matall.G(num2,num2)=matall.G(num2,num2)+1/value;
         end
         if (num1~=0) &(num2~=0)
         matall.G(num1,num2)=matall.G(num1,num2)-1/value;
         matall.G(num2,num1)=matall.G(num2,num1)-1/value;
         end   
 elseif (n==2)
         if (num1~=0) & (num2~=0)
           matall.C(num1,num2)=matall.C(num1,num2)-value;
           matall.C(num2,num1)=matall.C(num2,num1)+value;
         end
         
elseif (n==3)
         
         if num1~=0
         matall.C(num1,num1)=matall.C(num1,num1)+value;
         end
         if num2~=0
         matall.C(num2,num2)=matall.C(num2,num2)+value;
         end
         if (num1~=0) &(num2~=0)
         matall.C(num1,num2)=matall.C(num1,num2)-value;
         matall.C(num2,num1)=matall.C(num2,num1)-value;
         end
 elseif (n==4)
      if num1~=0
          matall.B(num1,1)=matall.B(num1,1)-value;
      end

      if num2~=0
            matall.B(num2,1)=matall.B(num2,1)+value;
      end 
end