今天犯了乙個很低階的錯誤，但是不容易發現，在此做記錄提醒自己。
低通濾波器的理想衝激響應是：
hlp=sin(wc*n)/π/n
於是，我這樣寫了matlab程式：
m=-256:256;
wc=pi/4;
hn=sin(wc*m)./(pi*m);
win=hamming(513);
b=hn.*win;
這樣做會出現錯誤：如果你對b做fft，那麼fft(b)的結果全是nan的複數。
為什麼會這樣呢？
注意到m可能等於0，這是hn是的分母為0，這是造成錯誤的根源所在。
所以，此時當m=0時，hn的對應值應該去當m→0時的hn極限，這樣就不會出現問題了。
用MATLAB設計低通濾波器
濾波器有很多種設計方式。最簡單的，就是用已有的buttord計算出最符合條件的巴特沃思濾波器的階數n和截止頻率，再用butter計算n階巴特沃斯數字濾波器系統函式分子 分母多項式的係數向量b a。用freqz函式畫出濾波器的幅頻 相頻圖。用filter濾波。例 規定 wp為通帶截止頻率 ws為阻帶截...
FIR濾波器設計中的窗函式
使用窗函式的作用在於減少頻譜洩漏，由於在實際應用中，通常都是fir濾波器，這就不能使濾波器的係數達到sinc函式上的無限大的抽樣，從而實現在頻域的方波形狀的低通濾波，人們想到使用窗函式的方法對輸入的取樣點進行加窗處理。常用的窗函式有漢寧窗，漢明窗，布萊克慢窗，凱澤窗等。它們的評價指標就在於過渡帶寬度...
matlab中的高斯低通濾波器 加性高斯色雜訊通道
對於高斯白雜訊，不同時刻的樣本值是不相關的。為方便，在實際應用中通常將在訊號上疊加高斯白雜訊，對一些演算法效能進行評估。依據評估結果，指導工程設計實現。在實際工程中，高斯白雜訊的不相關性並一定能夠很好近似反映實際使用的雜訊環境。這會使得利用高斯白雜訊獲得的效能評估結果通常是實際工程問題的過估計。本節...