%《MATLAB在電子信息課程中的應用》第五章例5.4程序q504
% 一階電路開關的暫態計算
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著  2001年10月

% 電路暫態
clear,format compact
%figure(1),fg521,pause
r1=3;us=18;is=3;r2=12;r3=6;C=1;	% 給出原始數據
uc0=-12;ir20=uc0/r2;ir30=uc0/r3; 	% 算出初值ir20及uc0
ic0=is-ir20-ir30;
ir2f=is*r3/(r2+r3); 			% 算出終值ir2f及ucf
ir3f=is*r2/(r2+r3);
ucf=ir2f*r2; icf=0; 					
% 注意時間數組的設置，在t=0及10附近設兩個點
t=[-2,-1,0-eps,0+eps,1:9,10-eps,10+eps,11:20];  % 這語句與書上不同,已修改.
% 原語句t=[-2-eps:0-eps,0:9,10-eps,10+eps,11:20];只適用于MATLAB5.3以下版本.
figure(2),plot(t),grid,pause 	% 找出時間與數組下標的關系，t=10+eps對應下標15
uc(1:3)=-12;ir2(1:3)=3;			% t<0時的值
T = r2*r3/(r2+r3)*C;			% 求充電時常數	
uc(4:14)=ucf+(uc0-ucf)*exp(-t(4:14)/T);	%
ir2(4:14)=ir2f+(ir20-ir2f)*exp(-t(4:14)/T);	% 用三要素法求輸出
uc(15)=uc(14);ir2(15)=is;		% 求t=10+eps時的各初值
ucf2=-12;ir2f=is;					% 求uc和ir2在新區間終值ucf2和ir2f
T2=r1*r3/(r1+r3)*C;						% t=10+eps到t=20區間的時常數
uc(15:25)=ucf2+(uc(15)-ucf2)*exp(-(t(15:25)-t(15))/T2);	% 再用三要素法求輸出
ir2(15:25)=is;
figure(3),subplot(2,1,1);h1=plot(t,uc);		% 繪uc圖
grid,set(h1,'linewidth',3)					% 加大線寬
subplot(2,1,2),h2=plot(t,ir2);		% 繪ir2圖
grid,set(h2,'linewidth',3)