%Samp5_19
N=5;  %濾波器階數
Wn=1000;  %邊界頻率
[b,a]=besself(N,Wn); %設計Bessel濾波器
figure(1)
[H,w]=freqs(b,a,512);   %繪出濾波器的復數頻率特性
magH=abs(H);phaH=unwrap(angle(H));  %求幅頻響應和相頻響應
subplot(2,1,1),plot(w/(2*pi),20*log10(magH));  %以頻率為橫軸繪幅頻響應
grid on;xlabel('頻率/Hz');ylabel('振幅/dB');
subplot(2,1,2),plot(w/2/pi,angle(H)*180/pi);%繪制相頻響應
grid on;xlabel('頻率/Hz');ylabel('相位/^o')
figure(2)    %新的繪圖框
dt=1/10000;   %信號采樣間隔
f1=100;f2=1000;  %輸入兩個頻率成分
t=0:dt:0.1;
x=sin(2*pi*f1*t)+0.5*cos(2*pi*f2*t); %輸入信號
H=[tf(b,a)];   %濾波器在MATLAB系統中的表示
[y,t1]=lsim(H,x,t);  %模擬濾波器的輸出
subplot(2,1,1),plot(t,x),title('輸入信號')  %繪出輸入信號
subplot(2,1,2),plot(t1,y)  %繪出輸出信號
title('輸出信號'),xlabel('時間/s')