這其實是我的dsp課程的課程設計,老實說,乙個學期的dsp課程並沒有使我對dsp晶元增加了多少了解,更多的收穫是終於知道之前看如天書的數字訊號處理是用來幹嘛的了,然而,這也並沒有什麼卵用,也絲毫沒有減少我對傅利葉的恨意,哈哈,說說笑。
設計乙個低通濾波器,這其實是乙個很簡單很簡單的課程設計,甚至,還稱不上課程設計,叫實驗也不為過。我之所以想把它寫出來,一是臨近畢業,想把這之前做的一些有意思的東西分享一下,否則他們只能靜靜的躺在硬碟裡慢慢被遺忘,被刪掉。二,當然是造福各位還「掙扎」在天堂裡的「莘莘學子」啦。插一句,最近在實習,才發現,學校才真的是最接近天堂的地方,還在天堂裡的請好好珍惜。其實我只是想讓後來的師弟師妹們初步了解一下像傅利葉變換是怎麼在實際中應用的,而不只是背幾個不知何解的公式,想當初我就是背一下頻域相乘等同於時域卷積就去考訊號與系統了。。。。,第一次嘗試去寫東西,寫不好莫怪。
首先,奉上大神傑作好了,看完還不懂就去打死作者吧,與我無關,拜拜。。。。
老實說,要是早一點看到大神的傑作,可能就不用雲裡霧裡的就去考試了,總的來說,我們之所以要用傅利葉變換,一般訊號在時域裡面都是由幾個不同頻率和幅度的正弦波合成的,波形雜亂無章,我們並不能對它進行直接的處理,而經過傅利葉變換轉到頻域裡面就能直觀地看到每個正弦波的頻率,在頻域裡面對各個頻率成分進行處理,如刪減某個頻率,在經過反傅利葉變換就會得到我們想要的訊號。要說一句的是,時域到頻域,再由頻域到時域只是人理解的乙個過程,並不是說在dsp晶元中執行的程式就是先進行傅利葉變換,處理完再反變換的過程。頻域和時域是兩個不同的空間,反映的是同乙個物體,在乙個空間對這個物體進行處理自然在另乙個空間也會看到它的變化,舉個不恰當的例子,時域和頻域就如乙個人的身前和身後,在這個人身後胡一巴掌,你在前面看到 的不是胡一巴掌,對,你看到的是這個人吐了一口血。。。。。。而傅利葉變換就是這一巴掌和這一口血之間的聯絡。這樣就可以很好的解釋「頻域相乘等同於時域卷積」這句話了,時域裡面進行卷積運算要經過翻轉,移位,相乘等,然而這一切在頻域中就相當只進行了相乘運算。
廢話不多說,上圖,上圖之前還是先奉上高手作品, ,這裡面提到了數字訊號處理的第一步,抽樣,抽樣的具體原理推導原諒我也說不清,只能說抽樣在實際中就是常說的模/數轉換,就是通過一塊a/d晶元將模擬訊號轉化為一系列的數字量讓dsp去處理,下圖中的黃圈圈住的就是本實驗用的a/d晶元tlc320ad50,
前面說過,傅利葉表明的是任何訊號都是由不同頻率和幅度的正弦波合成的,下面就來合成兩個不同頻率的正弦波訊號,原諒我matlab臨時學的,手癢用用。。。。
可以看到,頻域圖明確的表明我們的訊號只有1000hz頻率分量,同樣,乙個4khz的正弦波訊號如下,
可以預見,如果把這兩個頻率分量的正弦波訊號合成,它的頻域是什麼樣子的,如下
就是這樣,時域和頻域之間還是有很多很有趣的聯絡的,有興趣的可以關心一下,基本來說,傅利葉變換是很強大很有用處的,寫這個更多是由於看了前面第一篇博文有感而寫,至今從網上博文受益良多,有時候自己動手做做,雖然簡單,但看到預期結果是還是很開心。
matlab設計低通濾波器
輸入 被白雜訊汙染的正弦訊號,fs 100khz,訊號頻率為10khz,雜訊訊號為20khz,現在要濾掉20khz的正弦訊號 clear all fs 100000 t 0 1 fs 0.003 f1 10000 f2 20000 signal1 sin 2 pi f1 t signal2 sin ...
IIR和FIR濾波器設計低通濾波器
設原始訊號為 採用iir濾波器和fir濾波器設計低通濾波器,比較兩類濾波器的濾波結果。要求 採用matlab語言實現,並分析結果。首先對原始訊號進行頻譜分析,確定濾波器引數。通過分析加噪訊號的頻譜,雜訊訊號為 原始訊號為 iir濾波器設計 按照數字濾波器技術指標 通帶邊界頻率wp 通帶最大衰減 阻帶...
用MATLAB設計低通濾波器
濾波器有很多種設計方式。最簡單的,就是用已有的buttord計算出最符合條件的巴特沃思濾波器的階數n和截止頻率,再用butter計算n階巴特沃斯數字濾波器系統函式分子 分母多項式的係數向量b a。用freqz函式畫出濾波器的幅頻 相頻圖。用filter濾波。例 規定 wp為通帶截止頻率 ws為阻帶截...