matlab圖形繪制ppt課件-資料下載頁
2026-01-05 03:50本頁面
【正文】 頻率特性的 MATLAB函數(shù) 。 所以這里有必要 給大家介紹一下如何運用 MATLAB對頻率響應(yīng) 函數(shù)的計算方法。 (1)polyval()因為 G(jw)的分子 、 分母均為有理 多項式 ， 可用多項式計算指令 ployval()計算系 統(tǒng)的頻率響應(yīng)，其調(diào)用格式為 Y=polyval(P,X) 2022/2/11 哈爾濱理工大學(xué)機(jī)械動力工程學(xué)院 藺勇智制作 38 【 說明 】 ?P是多項式系數(shù)向量 (降冪排列 )， X是自變量 (應(yīng)設(shè)為 jw)， Y是返回的計算結(jié)果 (復(fù)數(shù)數(shù)組 )。 ?對 Y利用 abs()、 angle()即可求出系統(tǒng)的 幅頻特性和相頻特性。 2022/2/11 哈爾濱理工大學(xué)機(jī)械動力工程學(xué)院 藺勇智制作 39 例如：繪制系統(tǒng) 的幅頻 特性曲線和相頻特性曲線。 程序如下： num=[11 11]。den=[1 15 4 0]。 w=::*pi。%產(chǎn)生頻率向量 Gw=polyval(num,j*w)./polyval(den,j*w)。%產(chǎn)生頻率特性 mag=abs(Gw)。%計算幅頻特性 theta=angle(Gw)。%計算相頻特性 subplot(2,1,1),plot(w,mag) grid,title(39。幅頻特性 39。) ylabel(39。|G|39。) 21 1 ( 1 )()( 1 5 4 )sGss s s????2022/2/11 哈爾濱理工大學(xué)機(jī)械動力工程學(xué)院 藺勇智制作 40 subplot(2,1,2),plot(w,theta) grid,title(39。相頻特性 39。) xlabel(39。\omega(rad/s)39。),ylabel(39。deg39。) 2022/2/11 哈爾濱理工大學(xué)機(jī)械動力工程學(xué)院 藺勇智制作 41 (2) freqs()若已知系統(tǒng)傳遞函數(shù) ， 還可以用 freqs求系統(tǒng)的頻率相應(yīng) 。 其有以下幾種調(diào) 用格式 h=freqs(b,a,w) 指定實式角頻率 ,返回 響應(yīng)值。 [h,w]=freqs(b,a)自動確定 200個頻率 點，返回響應(yīng)值和對應(yīng)的角頻率向量。 [h,w]=freqs(b,a,f)指定頻率 (HZ)向量 ， 返 回響應(yīng)值和對應(yīng)的角頻率向量。 freqs(b,a,w)繪制對指定正實角頻率向量 的幅頻和相頻特性曲線。 2022/2/11 哈爾濱理工大學(xué)機(jī)械動力工程學(xué)院 藺勇智制作 42 【 說明 】 ?b,a均為系統(tǒng)傳遞函數(shù)的分子 、 分母的系 數(shù)向量。 ?在返回指令值的指令中 ， 需調(diào)用 abs()和 angle()求取幅頻和相頻特性。 ?第四種調(diào)用可直接繪制系統(tǒng)的幅頻和相頻 特性曲線 ， 其中幅頻特性曲線為全對數(shù)坐 標(biāo) ， 而相頻特性曲線為半對數(shù)坐標(biāo) ， 并且 可以不指定頻率向量。 2022/2/11 哈爾濱理工大學(xué)機(jī)械動力工程學(xué)院 藺勇智制作 43 例子：用 freqs()指令重繪上面例子所示系統(tǒng) 的頻率特性曲線。 程序如下： num=[11 11]。den=[1 15 4 0]。 w=::*pi。%產(chǎn)生頻率向量 freqs(num,den,w) 2022/2/11 哈爾濱理工大學(xué)機(jī)械動力工程學(xué)院 藺勇智制作 44 101100 1 3 0 1 2 0 1 1 0 1 0 0 9 0F r e q u e n c y ( r a d / s )Phase (degrees)101100101100101102F r e q u e n c y ( r a d / s )Magnitude2022/2/11 哈爾濱理工大學(xué)機(jī)械動力工程學(xué)院 藺勇智制作 45 若不指定頻率向量 w， 直接執(zhí) freqs(num,den), 結(jié)果如下所示 。 顯然該指令能夠自動確定系 統(tǒng)頻率響應(yīng)何時得頻率范圍。 程序如下： num=[11 11]。den=[1 15 4 0]。 freqs(num,den) 2022/2/11 哈爾濱理工大學(xué)機(jī)械動力工程學(xué)院 藺勇智制作 46 102101100101102 1 8 0 1 6 0 1 4 0 1 2 0 1 0 0 8 0F r e q u e n c y ( r a d / s )Phase (degrees)102101100101102105100105F r e q u e n c y ( r a d / s )Magnitude2022/2/11 哈爾濱理工大學(xué)機(jī)械動力工程學(xué)院 藺勇智制作 47 請大家指正 ,謝謝 !!!!!!
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號-1