用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)
for m=1:(a-1)       %循环输入支路关系：关联且离开节点为+1、关联且进入节点为-1、无关联为0  
    for n=1:b
        A(m,n)=input('请输入支路关系：') %将输入的数字循环代替A矩阵中的0元素
    end
end
for k=1:b  %循环输入各支路中的阻值，有电感应先消互感
    fprintf('请输入第%g支路电阻:',k)
    Z(k,k)=input('R=')  %将输入的数字循环代替Z矩阵中的0元素
end
    for k=1:b %循环输入各支路中的电压值
    fprintf('请输入第%g支路电压:',k)
    Us(k,1)=input('Us=')   %将输入的数字循环代替Us矩阵中的0元素
    end
    for k=1:b  %循环输入各支路中的电压值
    fprintf('请输入第%g支路电流源数值:',k)
    Is(k,1)=input('Is=')  %将输入的数字循环代替Is矩阵中的0元素
    end
r=input('请输入受控源个数：')   %无受控源输入0
Y=Z'    %导纳矩阵Y
if r==0  %无受控源情况
J=A*Is-A*Y*Us
Yn=A*Y*A'
Un=inv(Yn)*J  %节点电压
U=A'*Un
I=Y*(U+Us)-Is  %节点电流
else       %多受控源情况
    for t=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





注：
1111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111
输入节点数：3
输入支路数：5
输入支路关系：1 0 -1 1 0 -1 -1 0 0 1（每输入一个按一下回车，下同）
输入各支路中的电压值：0 0 5 0 0
输入各支路中的电流值：0 0 0 0 0（或直接按回车）
输入各支路中的电阻：1 1 1 1 1
输入受控源个数：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轴范围
if R>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(['电容电压'],['电感电流'])
elseif R==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('时间t/s')
ylabel('电压 U/V  电流 I/A')
legend(['电容电压'],['电感电流'])
else   %欠阻尼情况
a=M   %衰减系数
w0=sqrt(N) %固有振荡角频率
wd=sqrt(w0*w0-a*a)   %衰减振荡角频率    
k1=u0
k2=1/w0*(a*u0+i0/C)
k=sqrt(k1*k1+k2*k2)
r=atan(k1/k2)
uc=k*exp(-a*t).*sin(wd*t+r)   %电容电压
il=-C*k*a*exp(-a*t).*sin(wd*t+r)+C*k*wd*exp(-a*t).*cos(wd*t+r)  %电容电流
plot(t,uc,'k-',t,il,'b-')  %输出并设置图像黑色实线为Uc图像、蓝色实线为il图像
xlabel('时间t/s')
ylabel('电压 U/V  电流 I/A')
legend(['电容电压'],['电感电流'])
end

  注：
  222222222222222222222222222222222222222222222222222222222222222222222222222222222222
  输入电阻值R=？？？（过阻尼、临界阻尼、欠阻尼依次为：30、20、10）
  输入电容值C=：0.004
  输入电感值L=：0.4
  输入初始状态电压值u0：6
  输入初试状态电流值i0：
  输入时间常数T：好像是0.6



第三道题：

   P=input('请输入功率:')
  U=input('请输入电压:')
  w=input('请输入功率因数:')
  C=input('请输入补偿电容：')
    r1=acos(w)  %初始相位差
    b=tan(r1)-C*U*U*314/P
    r=atan(b)   %接入补偿电容后的相位差
    w1=cos(r)   %接入补偿电容后的功率因数
    I=P/U/w     %初始电流
    I2=P/U/w1   %接入补偿电容后的电流
    t=0:0.04:2  %设置步长与x轴范围
    U0=U*cos(314*t*pi/180)  %初始电压表达式
    I0=I*cos(314*t*pi/180-r1) %初始电流表达式
    U1=U0  %接入补偿电容后电压表达式
    I1=I2*cos(314*t*pi/180+r) %接入补偿电容后电流表达式
    subplot(2,1,1)  %将当前图形窗口分为2行1列2个子图，指定一号图
    plot(t,U0,'r-',t,I0,'b-') %输出图像，红色实线表示U0、蓝色实线表示I0
    subplot(2,1,2)  %将当前图形窗口分为2行1列2个子图，指定二号图
    plot(t,U1,'r-',t,I1,'b-') %输出图像，红色表示U1、蓝色表示I1
   
   
   
   
   
   
    注：
    3333333333333333333333333333333333333333333333333333333333333333333333333333333333333333
输入功率:10000

输入电压:220

输入功率因数:0.6

输入补偿电容：？？？（8.7733e-004是W1为1的情况，输入的应是小于、等于、大于该数三种情况）


第四道题：

R=2000  %初始电阻
C=0.00000005  %初始电感
U=10   %初始电压
f=2000  %初始频率
w=2*pi*f  %角频率
t=0:0.001:0.1  %设置步长与x轴范围
U1=1/(j*w*C)/(R+1/(j*w*C))*40/pi   %输出电压U中相应谐波的复数U1、U2、U3（取展开的谐波序列前三项）
U2=1/(3*j*w*C)/(R+1/(3*j*w*C))*40/pi
U3=1/(5*j*w*C)/(R+1/(5*j*w*C))*40/pi
c1=abs(U1) %取输出电压U中相应谐波的复数U1、U2、U3的绝对值
c2=abs(U2)
c3=abs(U3)
a1=real(U1) %取输出电压U中相应谐波的复数U1、U2、U3的实部
a2=real(U2)
a3=real(U3)
b1=imag(U1) %取输出电压U中相应谐波的复数U1、U2、U3的虚部
b2=imag(U2)
b3=imag(U3)
r1=asin(b1/c1) %各相位差
r2=asin(b2/c2)
r3=asin(b3/c3)
u1=c1*sin(w*t*pi/180+r1) %在展开的谐波序列前三项电压作用下的电感电压
u2=c2*sin(3*w*t*pi/180+r2)
u3=c3*sin(5*w*t*pi/180+r3)
U=u1+u2+u3 %电感电压
plot(t,U'b-') %输出图像，以蓝色实线输出
xlabel('wt')
ylabel('电压 U')