matlab freqz函数使用
Matlab?freqz函數(shù)使用
(2012-10-22 13:42:03) 轉(zhuǎn)載▼標(biāo)簽： 雜談 |
freqz函數(shù)計(jì)算線性系統(tǒng)的頻率響應(yīng),包括幅頻響應(yīng)和相頻響應(yīng),基本輸入為線性系統(tǒng)的AMMA模型系數(shù)向量,一個(gè)典型的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)分母向量,對(duì)FIR系統(tǒng)為1;向量B為分子向量
函數(shù)使用零散傅立葉變換(fft)計(jì)算系統(tǒng)響應(yīng),計(jì)算范圍為0~pi,即單位圓的上半部分,因?yàn)槟J(rèn)情況下對(duì)于實(shí)數(shù)系數(shù)系統(tǒng)響應(yīng)關(guān)于DC對(duì)稱.計(jì)算公式為上式的頻域離散化形式.計(jì)算完成后在對(duì)目標(biāo)解析度使用算術(shù)插值方法進(jìn)行向量填充,填充形式取決于附加輸入?yún)?shù)信息.
?
freqz的基本調(diào)用格式為freqz(b,a),此時(shí)函數(shù)內(nèi)部默認(rèn)進(jìn)行512點(diǎn)算術(shù)插值,所以輸出是512點(diǎn)的幅度值和對(duì)應(yīng)的0~pi自變量向量,如果需要改變向量長(zhǎng)度便于繪圖,則可以附加維度信息,例如freqz(b,a,N)則進(jìn)行N點(diǎn)插值
?
如果系統(tǒng)不是實(shí)系數(shù)的而是復(fù)數(shù)系統(tǒng),則系統(tǒng)響應(yīng)不再對(duì)稱,0~pi的范圍不足以反映整體特性,此時(shí)有兩種選擇:
1. 繪制整個(gè)主值區(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)的離散傅立葉變換計(jì)算結(jié)果,如果為繪圖方便要標(biāo)識(shí)頻率信息,則調(diào)用時(shí)需要附加采樣頻率信息,例如freqz(b,a,[-fs/2:(fs/n):fs/2],fs);
?
freqz的基本使用方法是利用函數(shù)返回值進(jìn)行后續(xù)分析處理,此時(shí)函數(shù)不會(huì)自動(dòng)繪制曲線,如果直接調(diào)用而不指定返回接收變量,則函數(shù)自動(dòng)按照通常標(biāo)準(zhǔn)繪制幅頻/相頻曲線,幅頻曲線為對(duì)數(shù)格式
?
可以使用type freqz分析源代碼
分享：27
喜歡
閱讀(16746)┊ 評(píng)論 (7)┊ 收藏(2)┊轉(zhuǎn)載(8) ┊喜歡▼┊打印┊舉報(bào)| 排行榜 |
轉(zhuǎn)載列表：
轉(zhuǎn)載是分享博文的一種常用方式...
前一篇：小寶的近照 后一篇：Perl?遞歸遍歷目錄，找出目錄下以及所有子目錄中的文件名 評(píng)論 重要提示：警惕虛假中獎(jiǎng)信息 [發(fā)評(píng)論]-
jiewangcong
總喜歡長(zhǎng)篇大論，不知所云2012-10-23 ?10:42舉報(bào)刪除分享回復(fù)(0)
-
愛珊瑚鼠
講的真的很好?贊一個(gè)！2013-3-30 ?23:52舉報(bào)刪除分享回復(fù)(0)
-
揚(yáng)州白言
freqz(b,a,[-fs/2:(fs/n):fs/2],fs);這種語(yǔ)法好像查不到啊？好像只有freqz(b,a,N,fs)這種用法啊??2013-7-29 ?08:44舉報(bào)刪除分享回復(fù)(0)
-
mits
講的很好啊2013-10-30 ?10:03舉報(bào)刪除分享回復(fù)(0)
-
mits
但是為什么我設(shè)B=[1?0];
A?=[1?-1];
[H,?w]?=?freqz(B,?A,?[-4000:1:4000],?8000);
得到的只有之流分量...2013-10-30 ?10:09舉報(bào)刪除分享回復(fù)(0)
-
吳胡韡
沒源程序?沒真相2014-1-8 ?18:21舉報(bào)刪除分享回復(fù)(0)
-
mykilla
很有幫助?感謝！2015-11-3 ?09:36舉報(bào)刪除分享回復(fù)(0)
登錄名： 密碼： 找回密碼 注冊(cè) 記住登錄狀態(tài)
昵???稱：
分享到微博 ???評(píng)論并轉(zhuǎn)載此博文
按住左邊滑塊，拖動(dòng)完成上方拼圖 點(diǎn)擊上圖按鈕并沿道路拖動(dòng)到終點(diǎn)處 移動(dòng)到此開始驗(yàn)證 匿名評(píng)論 發(fā)評(píng)論以上網(wǎng)友發(fā)言只代表其個(gè)人觀點(diǎn)，不代表新浪網(wǎng)的觀點(diǎn)或立場(chǎng)。
<?前一篇小寶的近照 后一篇?>Perl?遞歸遍歷目錄，找出目錄下以及所有子目錄中的文件名 《新程序員》：云原生和全面數(shù)字化實(shí)踐50位技術(shù)專家共同創(chuàng)作，文字、視頻、音頻交互閱讀總結(jié)
以上是生活随笔為你收集整理的matlab freqz函数使用的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。
- 上一篇： 关于2016年房价的思考
- 下一篇： matlab添加多个legend