q507.m

matlab在电子信息中的应用和实例中的源代码

%《MATLAB及在电子信息课程中的应用》第五章例5.7程序q507
% 直流电路的暂态计算：二阶欠阻尼电路
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月初版，2003年7月第二版
%
clear, format compact
L=0.5; C=0.02; R=12.5;		% 输入元件参数
uc0=1; iL0=0;
for R=1:10
    a=R/2/L;; w0=sqrt(1/(L*C)); % 输入给定参数
    a1=a-sqrt(a^2-w0^2);		% 方程的两个根
    a2=a+sqrt(a^2-w0^2);
    dt=0.01;t=0:dt:1;				% 设定时间数组
    % 方法1，用公式
    uc1=(a2*uc0+iL0/C)/(a2-a1)*exp(-a1*t);		% uc的第一个分量
    uc2=-(a1*uc0+iL0/C)/(a2-a1)*exp(-a2*t);		% uc的第二个分量
    iL1=-a1*C*(a2*uc0+iL0/C)/(a2-a1)*exp(-a1*t);		% iL的第一个分量
    iL2=a2*C*(a1*uc0+iL0/C)/(a2-a1)*exp(-a2*t);		% iL的第二个分量
    uc=uc1+uc2; iL=iL1+iL2;							% 把两个分量相加
    % 分别画出两种数据曲线
    figure(1),plot(t,uc),hold on,grid on
    set(gcf,'color','w')            % 设置图形背景色为白色
    figure(2),plot(t,iL),hold on,grid on
    set(gcf,'color','w')            % 设置图形背景色为白色
end
figure(1),grid,figure(2),grid
% 方法2，用拉普拉斯变换及留数法
for R=1:9
    num=[uc0,R/L*uc0+iL0/C];				% uc(s)的分子系数多项式
    den=[1,R/L,1/L/C];		% uc(s)的分母系数多项式
    [r,p,k]=residue(num,den);	% 求极点留数
    ucn=r(1)*exp(p(1)*t)+r(2)*exp(p(2)*t);	% 求时域函数
    iLn=C*diff(ucn)/dt;			% 对ucn求导得到电流iLn
    figure(3),plot(t,ucn),grid,hold on		% 绘曲线
    set(gcf,'color','w')            % 设置图形背景色为白色
    figure(4),plot(t(2:end),iLn),grid,hold on
    set(gcf,'color','w')            % 设置图形背景色为白色
end 
figure(3),grid,figure(4),grid