matlab在电路中的全响应程序
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
A=[1 0 0 1 1 0;0 0 1 0 -1 1;0 1 0 -1 0 -1] %输入相关矩阵A
Is=[2;0;-5;0;0;0] %输入电流源列向量
Us=[0;3;0;-2;0;0] %输入电压源列向量
Y=[1 0 0 0 0 0 ;-2 1/3 0 0 0 0 ;0 0 1/4 0 0 0;0 0 0 1/5 0 1/3 ;0 0 0 0 1/2 0;0 0 0 0 0 1/8] %输入导纳矩阵Y
Z=inv(Y)%求解阻抗矩阵Z
Yn=A*Y*A' %求解节点导纳矩阵Yn
Un=inv(Yn)*(A*Is-A*Y*Us) %求解节点电压Un
U=A'*Un %求解之路电压U
I=Y*(U+Us)-Is %求解支路电流I
结果如下：
Z =
1.0000 0 0 0 0 0
6.0000 3.0000 0 0 0 0
0 0 4.0000 0 0 0
0 0 0 5.0000 0 -13.3333
0 0 0 0 2.0000 0
0 0 0 0 0 8.0000
Yn =
1.7000 -0.1667 -0.5333
-0.5000 0.8750 -0.1250
-2.2000 -0.4583 0.9917
Un =
-7.8339
-13.7855
-25.1626
U =
-7.8339
-25.1626
-13.7855
17.3287
5.9516
11.3772
I =
-9.8339
8.2803
1.5536
6.8581
2.9758
1.4221
A=input('相关矩阵A=:')
Is=input('电流源Is= :')
Us=input('电压源Us=:' )
display ('是否含有受控源？是（a=1)否（a=0)')
a=input('输入a=:' )
if(a==0) %判断语句
Z=input('阻抗矩阵Z=:')
Y=inv(Z) %求解导纳矩阵Y
else (a==1)
Y=input('导纳矩阵Y=:')
Z=inv(Y) %求解阻抗矩阵Z
end
Yn=A*Y*A' %求解节点导纳矩阵Yn
Un=inv(Yn)*(A*Is-A*Y*Us) %求解节点电压Un
U=A'*Un %求解之路电压U
I=Y*(U+Us)-Is %求解支路电流I
含受控源的运行结果如下
相关矩阵A=:[1 0 0 1 1 0;0 0 1 0 -1 1;0 1 0 -1 0 -1]
电流源Is= :[2;0;-5;0;0;0]
电压源Us=:[0;3;0;-2;0;0]
是否含有受控源？是（a=1)否（a=0)
输入a=:1
导纳矩阵Y=:[1 0 0 0 0 0 ;-2 1/3 0 0 0 0 ;0 0 1/4 0 0 0;0 0 0 1/5 0 1/3;0 0 0 0 1/2 0;0 0 0 0 0 1/8]
Z =
1.0000 0 0 0 0 0
6.0000 3.0000 0 0 0 0
0 0 4.0000 0 0 0
0 0 0 5.0000 0 -13.3333
0 0 0 0 2.0000 0
0 0 0 0 0 8.0000
Yn =
1.7000 -0.1667 -0.5333
-0.5000 0.8750 -0.1250
-2.2000 -0.4583 0.9917
Un =
-7.8339
-13.7855
-25.1626
U =
-7.8339
-25.1626
-13.7855
17.3287
5.9516
11.3772
I =
-9.8339
8.2803
1.5536
6.8581
2.9758
1.4221
A=[1 -1 0;-1 0 1] %输入相关矩阵A
Is=[-2;-1;0] %输入电流源列向量
Us=[2;0;-3] %输入电压源列向量
Z=[2 0 0;0 0.8 0;0 0 3] %输入阻抗矩阵Z
Y=inv(Z) %求解导纳矩阵Y
Yn=A*Y*A'%求解节点导纳矩阵Yn
Un=inv(Yn)*(A*Is-A*Y*Us) %求解节点电压Un U=A'*Un%求解支路电压U
I=Y*(U+Us)-Is%求解支路电流I
运行结果如下：
Yn =
1.7500 -0.5000
-0.5000 0.8333
Un =
0.2759
4.9655
U =
-4.6897
-0.2759
4.9655
I =
0.6552
0.6552
0.6552
(要有表达式RU194页)
A=input('相关矩阵A=:')
Is=input('电流源Is= :')
Us=input('电压源Us=:' )
display ('是否含有受控源？是（a=1)否（a=0)')