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

clear, format compact
R=2; C=0.5; T=R*C;uc0=4;			% 輸入元件參數
Um=10; w=2; Zc=1/j/w/C;		% 輸入給定參數
t=0:0.1:10;						% 設定時間數組
us=Um*cos(w*t);				% 設定激勵信號
Ha=abs(Zc/(Zc+R));Hph=angle(Zc/(Zc+R));
ucst=Um*Ha*cos(w*t+Hph);			% 穩態分量計算
ucp0=ucst(1);				% 穩態分量的初始值
uctr=(uc0-ucp0)*exp(-t/T);	% 暫態分量
uc=uctr+ucst;					% 總的uc為兩項之和
%把三種數據畫在一張圖上
plot(t,uc,'-',t,uctr,':',t,ucst,'-.'),grid
%set(h,'Linewidth',2)			% 線條加粗
legend('uc','uctr','ucst')		% 用圖例標注