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)')
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' %求解节点电压Un
Un=inv(Yn)*(A*Is-A*Y*Us) %求解节点电压Un %求解支路电流I
U=A'*Un %求解支路电压U
I=Y*(U+Us)-Is %求解支路电流I
不含受控源的运行结果如下：
是否含有受控源？是（a=1)否（a=0)
输入a=:0
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页)
封面设计题目：电路信号处理综合课程设计课程设计任务书
一、解决问题
1．大规模电路的计算机辅助分析
2．动态电路的计算机辅助分析
3．谐振电路频率响应的分析
设计说明书写作：
一、课程设计任务书（第1页）四设计总结五参考文献
二、目录（第3页）
三、设计第一题
A程序B例子C最后运行结果D对曲线进行注解。