%《MATLAB及其在理工課程中的應用指南》(十一五規劃版)第八章例8－1－6程序
% 交流諧振電路
% 西安電子科技大學出版社出版  陳懷琛編著 2007年5月
%
r1=2; r2=3; L1=0.75e-3; L2=0.25e-3; C=1000e-12;rs=28200;
L=L1+L2;r=r1+r2;rse=rs*(L/L1)^2; 		% 折算內阻
f0=1/(2*pi*sqrt(C*L))
Q0=sqrt(L/C)/r,r0=L/C/r;			% 空載Q0值
re=r0*rse/(r0+rse); 				% 折算內阻與回路電阻的并聯
Q=Q0*re/r0,B=f0/Q,				% 實際Q值和通帶
s=log10(f0);
f=logspace(s-.1,s+.1,501);w=2*pi*f;		% 設定計算的頻率范圍及數組
z1e=r1+j*w*L; z2e=r2+1./(j*w*C); % 等效單回路中兩個電抗支路的阻抗
ze = 1./(1./z1e+1./z2e+1./rse);	 % 等效單回路中三個支路的并聯阻抗
subplot(2,1,1),loglog(w,abs(ze)),grid	% 畫對數幅頻特性
axis([min(w),max(w),0.9*min(abs(ze)),1.1*max(abs(ze))])
subplot(2,1,2),semilogx(w,angle(ze)*180/pi),grid % 畫相頻特性
axis([min(w),max(w),-100,100])
fh=w(find(abs(1./(1./z1e+1./z2e))>50000))/2/pi; % % 求幅特性大于50k的頻帶
fhmin=min(fh),fhmax=max(fh),
set(gcf,'color','w')