%《MATLAB在電子信息課程中的應用》第五章例5.9程序q509
% 交流電路加多頻信號：戴維南定理
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著  2001年10月

clear, format compact
w=[eps,1,2];us=[10,10,0];Is=[5,0,5];	% 置頻率及相應的電壓，電流數組
z1=1./(0.5*w*j);z2=[2,2,2];z3=[2,2,2];z4=1*w*j;	%各阻抗數組
uoc=(z2./(z1+z2)-z4./(z3+z4)).*us;		% 求開路電壓
zeq=z3.*z4./(z3+z4)+z1.*z2./(z1+z2);	% 求等效內阻
u=Is.*zeq+uoc;				% 按等效電路求u
disp('	 w       um       phi  ')		% 顯示各分量的值
disp([w',abs(u'),angle(u')*180/pi])