%《MATLAB及其在理工課程中的應用指南》(十一五規劃版)第八章例8－2－3程序
% 運算放大器有限增益和頻率響應對電路閉環特性的影響
% 西安電子科技大學出版社出版  陳懷琛編著 2007年7月
%
Z2=[20,100,500]*1000; Z1=2000;	 % 設定元件參數
A0=input('A0= (書上取2e6) ');
w1=input('w1= (書上先取500，后取50) ');
w2=input('w2= (書上取2e6) ');
w3=input('w3= (書上取5e7) ');
w=logspace(2,8);		 % 設定頻率數組
b=A0*w1*w2*w3; a=poly([-w1,-w2,-w3]);	 % 列出運算放大器分子分母系數向量
A=polyval(b,j*w)./polyval(a,j*w);	 % 求放大器開環頻率響應
for i=1:3			 % 循環計算三種Z2的閉環響應
   Z12(i)=Z2(i)/Z1;
   H(i,:)=-Z12(i)./(1+(1+Z12(i))./A);	% 放大器閉環響應
   % Hdb=20*log10(abs(H));		% 化為分貝
   semilogx(w,abs(H(i,:)),'linewidth',2),hold on	% 畫出頻率-增益曲線
end
v=axis;axis(v);				% 保持w坐標,
semilogx(w,abs(A),'r--','linewidth',2)			% 畫出開環頻率-增益響應
grid,hold off
set(gcf,'color','w')