用MatLAB做电路课程设计.pdf

用MatLAB做电路课程设计

06年年底做的电路课程设计，当时第一次接触MatLAB，破老师一点都不给
讲，只好自己一点点的学习总结。

题目现在楼忘记了，只剩下程序了。

第一道题：

a=input('请输入节点数：')%输入节点个数
b=input('请输入支路数：')%输入支路个数
A=zeros(a-1,b)%创建元素全为零关联矩阵A、阻抗矩阵Z、电压
源矩阵Us、电流源矩阵Is
Z=zeros(b,b)
Us=zeros(a,1)
Is=zeros(a,1)
form=1:(a-1)%循环输入支路关系：关联且离开节点为+1、关联且进入
节点为-1、无关联为0
forn=1:b
A(m,n)=input('请输入支路关系：')%将输入的数字循环代替A矩阵
中的0元素
end
end
fork=1:b%循环输入各支路中的阻值，有电感应先消互感
fprintf('请输入第%g支路电阻:',k)
Z(k,k)=input('R=')%将输入的数字循环代替Z矩阵中的0元素
end
fork=1:b%循环输入各支路中的电压值
fprintf('请输入第%g支路电压:',k)
Us(k,1)=input('Us=')%将输入的数字循环代替Us矩阵中的0元素
end
fork=1:b%循环输入各支路中的电压值
fprintf('请输入第%g支路电流源数值:',k)
Is(k,1)=input('Is=')%将输入的数字循环代替Is矩阵中的0元素
end
r=input('请输入受控源个数：')%无受控源输入0
Y=Z'%导纳矩阵Y
ifr==0%无受控源情况
J=A*Is-A*Y*Us
Yn=A*Y*A'
Un=inv(Yn)*J%节点电压
U=A'*Un
I=Y*(U+Us)-Is%节点电流
else%多受控源情况
fort=1:r
control=input('请输入被控制支路：')
uncontrol=input('请输入控制支路：')
math=input('请输入受控系数：')
Y(control,uncontrol)=math%将输入的数字循环代替Y矩阵中的相应0元素
end
J=A*Is-A*Y*Us
Yn=A*Y*A'
Un=inv(Yn)*J%节点电压
U=A'*Un
I=Y*(U+Us)-Is%节点电流
end

注：
111111111111111111111111111111111111111111111111111111111111111111111
1111111111111111111111111111111
输入节点数：3
输入支路数：5
输入支路关系：10-110-1-1001（每输入一个按一下回车，下同）
输入各支路中的电压值：00500
输入各支路中的电流值：00000（或直接按回车）
输入各支路中的电阻：11111
输入受控源个数：1
输入被控制支路:2
输入控制支路：4
输入受控系数：2

第二道题：

clear
R=input('请输入电阻值R=：')
C=input('请输入电容值C=：')
L=input('请输入电感值L=：')
u0=input('请输入初始状态电压值u0：')
i0=input('请输入初试状态电流值i0：')
T=input('请输入时间常数T：')
M=R/2/L
N=1/L/C
t=0:T/100:T%设置步长与x轴范围
ifR>2*sqrt(L/C)%过阻尼情况
O=abs(M*M-N)
s1=-M+sqrt(O)%两个特征值
s2=-M-sqrt(O)
uc=u0/(s1-s2)*(s1*exp(s2*t)-s2*exp(s1*t))+i0/C/(s1-
s2)*(exp(s1*t)-exp(s2*t))%电容电压
il=u0*s1*s2*C/(s1-s2)*(exp(s2*t)-exp(s1*t))+i0/(s1-
s2)*(s1*exp(s1*t)-s2*exp(s2*t))%电容电流
plot(t,uc,'k-',t,il,'b-')%输出并设置图像,黑色实线为Uc图像、蓝
色实线为il图像
xlabel('时间t/s')
ylabel('电压U/V电流I/A')
legend(['电容电压'],['电感电流'])
elseifR==2*sqrt(L/C)%临界阻尼情况
a=-M%特征值
uc=u0*(1+a*t).*exp(a*t)+i0/C*t.*exp(a*t)%电容电压
il=-u0*a*a*C*t*diag(exp(a*t))+i0*(1+a*t)*diag(exp(a*t))%电容电
流
plot(t,uc,'k-',t,il,'b-')%输出并设置图像,黑色实线为Uc图像、蓝
色实线为il图像
xlabel(