exn812.m

MATLAB及在电子信息课程中的应用.感觉很有用,本来也不花钱,有想用的来拿吧

%《MATLAB及其在理工课程中的应用指南》(十一五规划版)第八章例8－1－2程序
% 直流电路的暂态计算
% 西安电子科技大学出版社出版  陈怀琛编著 2007年5月
%
% 电路暂态
r1=3;us=18;is=3;r2=12;r3=6;C=1;	% 给出原始数据
uc0=-12;ir20=uc0/r2;ir30=uc0/r3; 	% 算出初值
ic0=is-ir20-ir30;
ir2f=is*r3/(r2+r3); 			% 算出终值
ir3f=is*r2/(r2+r3);
ucf=ir2f*r2; icf=0; 					
t=[[-2:0]-eps,0:15];			% 注意时间数组的设置，在t=0附近设两个点
uc(1:3)=-12;ir2(1:3)=is;			% t<0时的值
T = r2*r3/(r2+r3)*C;			% 求充电时常数	
uc(4:19)=ucf+(uc0-ucf)*exp(-t(4:19)/T);	%
ir2(4:19)=ir2f+(ir20-ir2f)*exp(-t(4:19)/T);	% 用三要素法求输出
subplot(2,1,1);h1=plot(t,uc),		% 绘图
grid,set(h1,'linewidth',2)					% 加大线宽
subplot(2,1,2),h2=plot(t,ir2);
grid,set(h2,'linewidth',2)