matlab實(shí)現(xiàn)波特圖

1、用matlab實(shí)現(xiàn)傳遞函數(shù)的波特圖波特圖是線性非時(shí)變系統(tǒng)的傳遞函數(shù)對頻率的半對數(shù)坐標(biāo)圖，其橫軸頻率以對數(shù)尺度（log scale）表示，利用波特圖可以看出系統(tǒng)的頻率響應(yīng)。波特圖一般是由二張圖組合而成，一張幅頻圖表示頻率響應(yīng)增益的分貝值對頻率的變化，另一張相頻圖則是頻率響應(yīng)的相位對頻率的變化。波特圖可以用電腦軟件（matlab）或儀器繪制，也可以自行繪制。利用波特圖可以看出在不同頻率下，系統(tǒng)增益的大小及相位，也可以看出大小及相位隨頻率變化的趨勢。波特圖的圖形和系統(tǒng)的增益，極點(diǎn)、零點(diǎn)的個(gè)數(shù)及位置有關(guān)，只要知道相關(guān)的資料，配合簡單的計(jì)算就可以畫出近似的波特圖，這是使用波特圖的好處波特圖的分析波特圖表

2、明了一個(gè)電路網(wǎng)絡(luò)對不同頻率信號的放大能力。但是在電子電路中，這種圖有可能比較麻煩，一方面，要表示一個(gè)網(wǎng)絡(luò)在低頻和高頻下的所有情況，那么橫軸（頻率軸會很長）。此外，一般放大電路的放大倍數(shù)可能達(dá)到幾百，使得縱軸也很長。第三，這樣畫出的圖形往往是很不規(guī)則曲線。把曲線做直線化處理。畫圖所依據(jù)的式子中會得到fL fH的數(shù)值。得出的波特圖也應(yīng)該在fL和fH處出現(xiàn)拐角（此點(diǎn)所在的頻率稱為截?cái)囝l率），不過這樣處理會產(chǎn)生一定的誤差。理論計(jì)算可知：在截?cái)囝l率處真實(shí) 值與估計(jì)值有3dB的誤差。 在斜率不為0的直線處要標(biāo)明斜率。標(biāo)明出每十倍頻程放大倍數(shù)的變化情況。 經(jīng)過這三種簡化，波特圖的曲線就是由一條折線組成看起來

3、非常舒服。雖然經(jīng)過處理造成了誤差，但已經(jīng)成為一種標(biāo)準(zhǔn)。已知傳遞函數(shù)G（s）=30（1+0.2s）s（s2+16s+100)要求w在10-2和103之間的波特圖。用matlab實(shí)現(xiàn)：Matlab實(shí)現(xiàn)的源程序：num=6,30; den=1 16 100 0; w=logspace(-2,3,500) mag,pha=bode(num,den,w); magdB=20*log10(mag); subplot(211); semilogx(w,magdB); grid on; title('Bode Diagram'); xlabel('Frequency(red/sec)'); ylabel('Gain dB'); subplot(212); semilogx(w,pha); grid on; xlabel('Frequency