Matlab?freqz函數(shù)使用

(2012-10-22 13:42:03) 轉(zhuǎn)載▼

標(biāo)簽： 雜談

freqz函數(shù)計算線性系統(tǒng)的頻率響應(yīng),包括幅頻響應(yīng)和相頻響應(yīng),基本輸入為線性系統(tǒng)的AMMA模型系數(shù)向量,一個典型的AMMA模型為

%??????????????jw??????????????-jw?????????????-jmw
%???????jw?B(e)??? b(1) +b(2)e + .... + b(m+1)e
%????H(e) = ---- = ------------------------------------
%??????????????jw??????????????-jw?????????????-jnw
%???????????A(e)??? a(1) +a(2)e + .... + a(n+1)e

其中,向量A為IIR系統(tǒng)分母向量,對FIR系統(tǒng)為1;向量B為分子向量

函數(shù)使用零散傅立葉變換(fft)計算系統(tǒng)響應(yīng),計算范圍為0~pi,即單位圓的上半部分,因為默認(rèn)情況下對于實(shí)數(shù)系數(shù)系統(tǒng)響應(yīng)關(guān)于DC對稱.計算公式為上式的頻域離散化形式.計算完成后在對目標(biāo)解析度使用算術(shù)插值方法進(jìn)行向量填充,填充形式取決于附加輸入?yún)?shù)信息.

?

freqz的基本調(diào)用格式為freqz(b,a),此時函數(shù)內(nèi)部默認(rèn)進(jìn)行512點(diǎn)算術(shù)插值,所以輸出是512點(diǎn)的幅度值和對應(yīng)的0~pi自變量向量,如果需要改變向量長度便于繪圖,則可以附加維度信息,例如freqz(b,a,N)則進(jìn)行N點(diǎn)插值

?

如果系統(tǒng)不是實(shí)系數(shù)的而是復(fù)數(shù)系統(tǒng),則系統(tǒng)響應(yīng)不再對稱,0~pi的范圍不足以反映整體特性,此時有兩種選擇:

1. 繪制整個主值區(qū)(-pi~pi),調(diào)用格式是增加字串型參數(shù)'whole',例如freqz(b,a,n,'whole');

2. 繪制定制區(qū)域,則直接給定范圍,該范圍為給定自變量向量,例如 freqz(b,a,[-pi:2*pi/n:pi])

使用上述調(diào)用格式,輸出為-pi~pi范圍內(nèi)的離散傅立葉變換計算結(jié)果,如果為繪圖方便要標(biāo)識頻率信息,則調(diào)用時需要附加采樣頻率信息,例如freqz(b,a,[-fs/2:(fs/n):fs/2],fs);

?

freqz的基本使用方法是利用函數(shù)返回值進(jìn)行后續(xù)分析處理,此時函數(shù)不會自動繪制曲線,如果直接調(diào)用而不指定返回接收變量,則函數(shù)自動按照通常標(biāo)準(zhǔn)繪制幅頻/相頻曲線,幅頻曲線為對數(shù)格式

?

可以使用type freqz分析源代碼

分享：

27

喜歡

閱讀(16746)┊ 評論 (7)┊ 收藏(2)┊轉(zhuǎn)載(8) ┊喜歡▼┊打印┊舉報 已投稿到：

排行榜

轉(zhuǎn)載列表：

轉(zhuǎn)載

轉(zhuǎn)載是分享博文的一種常用方式...

前一篇：小寶的近照 后一篇：Perl?遞歸遍歷目錄，找出目錄下以及所有子目錄中的文件名 評論 重要提示：警惕虛假中獎信息 [發(fā)評論]

• jiewangcong

  總喜歡長篇大論，不知所云

  2012-10-23 ?10:42舉報刪除分享回復(fù)(0)

• 愛珊瑚鼠

  講的真的很好?贊一個！

  2013-3-30 ?23:52舉報刪除分享回復(fù)(0)

• 揚(yáng)州白言

  freqz(b,a,[-fs/2:(fs/n):fs/2],fs);這種語法好像查不到啊？好像只有freqz(b,a,N,fs)這種用法啊??

  2013-7-29 ?08:44舉報刪除分享回復(fù)(0)

• mits

  講的很好啊

  2013-10-30 ?10:03舉報刪除分享回復(fù)(0)

• mits

  但是為什么我設(shè)B=[1?0];
  
  
  A?=[1?-1];
  
  
  [H,?w]?=?freqz(B,?A,?[-4000:1:4000],?8000);
  得到的只有之流分量...

  2013-10-30 ?10:09舉報刪除分享回復(fù)(0)

• 吳胡韡

  沒源程序?沒真相

  2014-1-8 ?18:21舉報刪除分享回復(fù)(0)

• mykilla

  很有幫助?感謝！

  2015-11-3 ?09:36舉報刪除分享回復(fù)(0)

發(fā)評論 id="postCommentIframe" frameborder="0" style="border:1px solid #C7C7C7;height:158px;width:448px;maring-top:1px;background-color:white;" src="http://blog.sina.com.cn/main_v5/ria/blank2.html"> 更多>>

登錄名： 密碼： 找回密碼 注冊 記住登錄狀態(tài)

昵???稱：

分享到微博 ???評論并轉(zhuǎn)載此博文

按住左邊滑塊，拖動完成上方拼圖 點(diǎn)擊上圖按鈕并沿道路拖動到終點(diǎn)處 移動到此開始驗證

匿名評論 發(fā)評論

以上網(wǎng)友發(fā)言只代表其個人觀點(diǎn)，不代表新浪網(wǎng)的觀點(diǎn)或立場。

<?前一篇小寶的近照 后一篇?>Perl?遞歸遍歷目錄，找出目錄下以及所有子目錄中的文件名 《新程序員》：云原生和全面數(shù)字化實(shí)踐50位技術(shù)專家共同創(chuàng)作，文字、視頻、音頻交互閱讀

總結(jié)

以上是生活随笔為你收集整理的matlab freqz函数使用的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。