vl.m

来自「EDA程序： 将spice网表转化为用于计算电路各种参数所需要的相应稀疏矩阵。」· M 代码 · 共 26 行

function [matall] = VL(type,matall,num1,num2,num3,value)
if type==1
        if num1~=0
         matall.G(num1,num3)=matall.G(num1,num3)+1;  %bu shi -1?
         matall.G(num3,num1)=matall.G(num3,num1)+1;
         end
         if num2~=0
         matall.G(num2,num3)=matall.G(num2,num3)+1;  
         matall.G(num3,num2)=matall.G(num3,num2)+1;  %bu shi -1?
         end
     matall.B(num3,1)=matall.B(num3,1)+value;
       
elseif type==2
        if (num1~=0) & (num3~=0)
         matall.G(num1,num3)=matall.G(num1,num3)-1;     % -1 +1?  
         matall.G(num3,num1)=matall.G(num3,num1)+1;
         end
         if (num2~=0) &(num3~=0)
         matall.G(num2,num3)=matall.G(num2,num3)+1;
         matall.G(num3,num2)=matall.G(num3,num2)-1;
         end
     if num3~=0
    matall.C(num3,num3)=matall.C(num3,num3)+value;     
     end
end