%《MATLAB在電子信息課程中的應(yīng)用》第五章例5.17程序q517
% 交流諧振電路
% 電子工業(yè)出版社出版  陳懷琛 吳大正 高西全合著  2001年10月

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; 		% 折算內(nèi)阻
f0=1/(2*pi*sqrt(C*L));
Q0=sqrt(L/C)/r,r0=L/C/r;			% 空載Q0值
re=r0*rse/(r0+rse); 				% 折算內(nèi)阻與回路電阻的并聯(lián)
Q=Q0*re/r0,B=f0/Q,				% 實際Q值和通帶
s=log10(f0);
f=logspace(s-.1,s+.1,501);w=2*pi*f;		% 設(shè)定計算的頻率范圍及數(shù)組
z1e=r1+j*w*L; z2e=r2+1./(j*w*C); % 等效單回路中兩個電抗支路的阻抗
ze = 1./(1./z1e+1./z2e+1./rse);	 % 等效單回路中三個支路的并聯(lián)阻抗
subplot(2,1,1),loglog(w,abs(ze)),grid	% 畫對數(shù)幅頻特性
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),