节点电压法matlab-豆柴文库

您所在位置：网站首页 / 节点电压法matlab.doc / 文档详情

在线预览结束，喜欢就下载吧，查找使用更方便

10 金币

下载文档

/ 2

下载提示

如果您无法下载资料，请参考说明：

1、部分资料下载需要金币，请确保您的账户上有足够的金币

2、已购买过的文档，再次下载不重复扣费

3、资料包下载后请先用软件解压，在使用对应软件打开

文本预览

%利用matlab编写的节点电压法解电路电压
NUM=5;%thenumberofthenodes
R=ones(NUM,NUM);%存储电阻的矩阵
I=zeros(NUM,1);%存储电流源的矩阵
fora=1:NUM
forb=1:NUM
R(a,b)=realmax;%令矩阵中的值等于浮点数最大值
end
end
para=1;
whilepara==1%选择输入
type=menu('要输入的选项','电阻','电流源','结束');
switchtype
case1
node1=input('元件的第一个节点:');
node2=input('元件的第二个节点:');
parameter=input('输入电阻/欧姆:');
R(node1,node2)=parameter;
R(node2,node1)=parameter;
case2
node1=input('元件的第一个节点:');
node2=input('元件的第二个节点:');
parameter=input('电流源/毫安:');
I(node1,1)=parameter;
I(node2,1)=-parameter;
case3
para=0;%退出
end
end
A=zeros(NUM,NUM);%电导矩阵
B=zeros(NUM,1);%电流源矩阵
tracer=1;
fora=1:NUM
forb=1:NUM
ifa~=b
A(a,a)=A(a,a)+1/R(a,b);%节点的总跨导
end
ifb~=a
A(a,b)=-1/R(a,b);%互导
end
end
end
fora=1:NUM
ifI(a,1)~=0
B(a,1)=I(a,1);
end
end

A(:,1)=0;%第一个节点为0
combine=[A,B];
answer=rref(combine);
fora=1:NUM-1%displaythenodevoltage
screen=strcat('Voltageatnode',num2str(a+1));
disp(screen);
disp(answer(a,NUM+1));
end
选择节点1为地。
结果截图如下：