% Matlab 程式設計 基礎篇，葉倍宏 編著，全華
% Matlab 程式設計 應用篇，葉倍宏 編著，全華
% 以下程式摘自以上書籍，敬請多多選購。 有關頻率響應，亦可參考 bode 指令 (繪製波德圖)
%  	OPA : 低通濾波器
	clear;
%	constant
	R1 = 1.0e3;
   	R2 = 100e3;
    C2 = 2*pi*1.59e-9;
%  求臨界頻率 WH     
	WH = 1/(R2*C2);			
% cal
	 W = logspace(-1,6);
	 amp=1./(1+j*W./WH);      
	   	at = -(R2 / R1) * amp;
	   	atp = 20*log10(abs(at));
% plotting
subplot(2,1,1);
    semilogx(W, atp); 
    grid;
	axis([1.0e-1 1.0e6 0 50]);
    ylabel('|T|');
subplot(2,1,2);
    semilogx(W, angle(at)*57.3);        
    grid;
	axis([1.0e-1 1.0e6 90 180]);
    xlabel('f (Hz)');
    ylabel('angle');