节点电压法matlab

%利用matlab编写的节点电压法解电路电压NUM=5;%thenumberofthenodesR=ones(NUM,NUM);%存储电阻的矩阵I=zeros(NUM,1);%存储电流源的矩阵fora=1:NUMforb=1:NUMR(a,b)=realmax;%令矩阵中的值等于浮点数最大值endendpara=1;whilepara==1%选择输入type=menu('要输入的选项','电阻','电流源','结束');switchtypecase1node1=input('元件的第一个节点:');node2=input('元件的第二个节点:');parameter=input('输入电阻/欧姆:');R(node1,node2)=parameter;R(node2,node1)=parameter;case2node1=input('元件的第一个节点:');node2=input('元件的第二个节点:');parameter=input('电流源/毫安:');I(node1,1)=parameter;I(node2,1)=-parameter;case3para=0;%退出endendA=zeros(NUM,NUM);%电导矩阵B=zeros(NUM,1);%电流源矩阵tracer=1;fora=1:NUMforb=1:NUMifa~=bA(a,a)=A(a,a)+1/R(a,b);%节点的总跨导endifb~=aA(a,b)=-1/R(a,b);%互导endendendfora=1:NUMifI(a,1)~=0B(a,1)=I(a,1);endendA(:,1)=0;%第一个节点为0combine=[A,B];answer=rref(combine);fora=1:NUM-1%displaythenodevoltagescreen=strcat('Voltageatnode',num2str(a+1));disp(screen);disp(answer(a,NUM+1));end选择节点1为地。结果截图如下：