簡介這份資源面向控制系統(tǒng)學(xué)習(xí)者和MATLAB/Simulink入門用戶聚焦如何在MATLAB環(huán)境中繪制S平面上的三維傳遞函數(shù)圖幫助理解動(dòng)態(tài)系統(tǒng)的頻率響應(yīng)與穩(wěn)定性。內(nèi)容圍繞傳遞函數(shù)建模、S平面極點(diǎn)零點(diǎn)分布、tf函數(shù)用法以及3D曲面可視化展開適合正在學(xué)習(xí)自動(dòng)控制原理或需要直觀觀察系統(tǒng)頻域特性的讀者。壓縮包共2個(gè)文件包含1個(gè)m腳本和1個(gè)txt說明整體約2KB其中m腳本用于定義傳遞函數(shù)、設(shè)定頻率范圍并計(jì)算響應(yīng)txt文件提供許可相關(guān)信息。已有1819人學(xué)習(xí)下載說明該主題在基礎(chǔ)控制仿真中具有一定關(guān)注度。通過這份資源讀者可以掌握從傳遞函數(shù)到三維曲面圖的完整實(shí)現(xiàn)思路理解bode、freqresp與surf等函數(shù)的配合方式并借助可視化結(jié)果分析系統(tǒng)在不同頻率下的行為為后續(xù)系統(tǒng)設(shè)計(jì)與性能優(yōu)化提供參考。1. 傳遞函數(shù)圖從“能畫出來”到“能講清楚”的那道坎很多人第一次在 MATLAB 里敲完bode(tf(num, den))看到一條漂亮的幅頻曲線彈出來就以為傳遞函數(shù)圖這件事已經(jīng)翻篇了。真正上手做項(xiàng)目才發(fā)現(xiàn)畫出來只是起點(diǎn)伯德圖要看截止頻率和相位裕度奈奎斯特圖要判斷包圍圈數(shù)根軌跡要讀增益臨界點(diǎn)零極點(diǎn)圖要區(qū)分左右半平面——同一套傳遞函數(shù)換一張圖就是換一個(gè)診斷視角。標(biāo)題里的“3傳遞函數(shù)圖”我理解成三層遞進(jìn)先把傳遞函數(shù)建模對再把三類常用圖伯德、奈奎斯特、根軌跡畫準(zhǔn)最后把圖上的關(guān)鍵參數(shù)讀出來用于控制器整定。這套流程適合做控制系統(tǒng)、電力電子、電機(jī)驅(qū)動(dòng)、過程控制的工程師也適合剛學(xué)自動(dòng)控制原理、想用 MATLAB 把課本結(jié)論跑一遍的學(xué)生。下面按“建?！媹D—讀圖—避坑—進(jìn)階”的順序把每一步落到可復(fù)現(xiàn)的命令和參數(shù)上。2. 傳遞函數(shù)建模先把 num/den 寫對再談畫圖2.1 三種建模方式與適用邊界MATLAB 里表示傳遞函數(shù)常見做法有三種系數(shù)向量、tf對象、zpk零極點(diǎn)增益形式。系數(shù)向量最原始num [1 2]、den [1 3 2]分別對應(yīng)分子 $s2$ 和分母 $s^23s2$但向量本身不帶采樣時(shí)間、不帶輸入輸出名一旦系統(tǒng)變復(fù)雜就容易搞混。tf對象是主力sys tf(num, den)之后可以直接sys.InputName、sys.TimeUnit還能做series、parallel、feedback運(yùn)算。zpk形式在已知零極點(diǎn)時(shí)最省事zpk([-1], [-2 -3], 5)表示零點(diǎn)在 -1、極點(diǎn)在 -2 和 -3、增益為 5。選哪種取決于你手里有什么。從微分方程出發(fā)先寫系數(shù)向量再轉(zhuǎn)tf從物理參數(shù)出發(fā)比如 RLC 電路先算零極點(diǎn)再轉(zhuǎn)zpk從實(shí)驗(yàn)數(shù)據(jù)出發(fā)先tfest辨識(shí)再轉(zhuǎn)tf。我一般會(huì)統(tǒng)一轉(zhuǎn)成tf對象再畫圖因?yàn)閎ode、nyquist、rlocus對tf對象的支持最完整返回的句柄也方便后續(xù)改線寬、加標(biāo)注。% 方式一系數(shù)向量直接構(gòu)造 num [1 2]; % 分子 s 2 den [1 3 2]; % 分母 s^2 3s 2 sys1 tf(num, den); % 方式二零極點(diǎn)增益形式 z [-1]; % 零點(diǎn) p [-2 -3]; % 極點(diǎn) k 5; % 增益 sys2 zpk(z, p, k); % 方式三由兩個(gè)簡單環(huán)節(jié)串聯(lián) s tf(s); % 先聲明 s 變量 sys3 (s 2) / (s^2 3*s 2); % 統(tǒng)一轉(zhuǎn)成 tf 對象方便后續(xù)畫圖 sys tf(sys3); disp(sys);這段代碼的關(guān)鍵在s tf(s)這一行。很多人不知道 MATLAB 支持這種“聲明式”寫法結(jié)果每次都要手算多項(xiàng)式乘法。用s變量寫出來的表達(dá)式可讀性接近數(shù)學(xué)公式(s2)/(s^23*s2)一眼就能看出結(jié)構(gòu)。tf(sys3)做一次類型歸一化避免后面bode和nyquist因?yàn)閷ο箢愋筒灰恢聢?bào)錯(cuò)。參數(shù)上注意den的系數(shù)順序是從高階到低階漏掉某一項(xiàng)要補(bǔ)零比如 $s^21$ 寫成[1 0 1]寫成[1 1]就變成 $s1$ 了這是新手最常翻車的地方。2.2 離散系統(tǒng)與采樣時(shí)間別讓 z 域和 s 域混著畫如果系統(tǒng)是數(shù)字控制器或者采樣數(shù)據(jù)傳遞函數(shù)要寫成 z 域形式tf(num, den, Ts)第三個(gè)參數(shù)是采樣時(shí)間。Ts -1表示未指定采樣時(shí)間畫出來的圖會(huì)按連續(xù)系統(tǒng)處理頻率軸單位是 rad/sTs 0.001表示 1 ms 采樣頻率軸會(huì)變成 Hz 或者歸一化頻率取決于你用的繪圖函數(shù)。Ts 0.001; % 采樣周期 1 ms num_d [0.5 0.3]; % 離散分子 den_d [1 -0.8 0.15]; % 離散分母 sys_d tf(num_d, den_d, Ts); % 離散傳遞函數(shù) % 畫伯德圖時(shí)指定頻率單位 figure; bode(sys_d, {2*pi*10, 2*pi*10000}); % 頻率范圍 10 Hz 到 10 kHz grid on;這里{2*pi*10, 2*pi*10000}是頻率范圍單位是 rad/s。很多人直接寫{10, 10000}結(jié)果圖只畫到 10 rad/s看不到高頻段。離散系統(tǒng)畫伯德圖時(shí)MATLAB 默認(rèn)會(huì)做雙線性變換頻率軸在奈奎斯特頻率附近會(huì)翹起來這是正常現(xiàn)象不是 bug。如果要做控制器設(shè)計(jì)建議先用c2d把連續(xù)控制器離散化再和離散被控對象串聯(lián)不要混著畫。2.3 從物理系統(tǒng)到傳遞函數(shù)一個(gè) RLC 例子假設(shè)有一個(gè)串聯(lián) RLC 電路輸入電壓 $u$輸出電容電壓 $y$根據(jù)基爾霍夫定律可以寫出 $\frac{Y(s)}{U(s)} \frac{1}{LCs^2 RCs 1}$。取 $L1$ mH、$C100$ μF、$R10$ Ω代入得 $LC 10^{-7}$、$RC 10^{-3}$。L 1e-3; % 1 mH C 100e-6; % 100 uF R 10; % 10 ohm num_rlc [1]; den_rlc [L*C, R*C, 1]; sys_rlc tf(num_rlc, den_rlc); % 查看自然頻率和阻尼比 [wn, zeta] damp(sys_rlc); fprintf(自然頻率 wn %.2f rad/s\n, wn(1)); fprintf(阻尼比 zeta %.4f\n, zeta(1));damp返回的自然頻率和阻尼比是讀圖前的重要參考。阻尼比小于 1 說明系統(tǒng)有諧振峰伯德圖上會(huì)在 $wn$ 附近鼓一個(gè)包阻尼比大于 1 說明過阻尼幅頻曲線單調(diào)下降。知道這兩個(gè)數(shù)后面看伯德圖時(shí)心里就有底不會(huì)把諧振峰誤當(dāng)成噪聲。3. 三類傳遞函數(shù)圖的畫法與讀圖要點(diǎn)3.1 伯德圖幅頻、相頻與裕度標(biāo)注伯德圖是最常用的頻域圖bode(sys)默認(rèn)畫幅頻和相頻兩張子圖。要讀截止頻率、相位裕度、幅值裕度直接用margin(sys)更省事它會(huì)在圖上標(biāo)出穿越頻率和裕度。figure; margin(sys_rlc); % 自動(dòng)標(biāo)注幅值裕度和相位裕度 grid on; % 手動(dòng)獲取裕度數(shù)值 [Gm, Pm, Wcg, Wcp] margin(sys_rlc); fprintf(幅值裕度 Gm %.2f dB\n, 20*log10(Gm)); fprintf(相位裕度 Pm %.2f deg\n, Pm); fprintf(相位穿越頻率 Wcg %.2f rad/s\n, Wcg); fprintf(幅值穿越頻率 Wcp %.2f rad/s\n, Wcp);margin返回的Gm是線性值轉(zhuǎn)成 dB 要取20*log10。Pm單位是度一般要求 45° 到 60° 之間太小超調(diào)大太大響應(yīng)慢。Wcp是幅頻曲線穿過 0 dB 的頻率也就是閉環(huán)帶寬的近似值。如果Pm是負(fù)數(shù)說明系統(tǒng)不穩(wěn)定這時(shí)候別急著調(diào)控制器先檢查建模有沒有把右半平面極點(diǎn)漏掉。畫伯德圖時(shí)頻率范圍很關(guān)鍵。默認(rèn)范圍可能太窄看不到高頻衰減也可能太寬低頻段擠成一團(tuán)。我一般用bode(sys, {0.1, 1e6})手動(dòng)指定先粗看一遍再縮范圍。另外bode支持多系統(tǒng)對比bode(sys1, sys2, sys3)會(huì)把三條曲線畫在同一張圖上用不同顏色區(qū)分適合比較不同參數(shù)下的頻響。3.2 奈奎斯特圖包圍圈數(shù)與穩(wěn)定性判斷奈奎斯特圖把頻率從 $-\infty$ 到 $\infty$ 的 $G(j\omega)$ 畫在復(fù)平面上穩(wěn)定性判據(jù)是 $Z N P$其中 $P$ 是開環(huán)右半平面極點(diǎn)數(shù)$N$ 是包圍 $(-1, j0)$ 點(diǎn)的圈數(shù)$Z$ 是閉環(huán)右半平面極點(diǎn)數(shù)。nyquist(sys)畫圖nyquistplot可以加選項(xiàng)。figure; nyquist(sys_rlc); grid on; title(Nyquist Plot of RLC System); % 判斷閉環(huán)穩(wěn)定性 P sum(real(pole(sys_rlc)) 0); % 開環(huán)右半平面極點(diǎn)數(shù) disp([開環(huán)右半平面極點(diǎn)數(shù) P , num2str(P)]);對于這個(gè) RLC 例子開環(huán)極點(diǎn)都在左半平面$P0$只要曲線不包圍 $(-1, j0)$閉環(huán)就穩(wěn)定。實(shí)際讀圖時(shí)要注意MATLAB 默認(rèn)只畫正頻率部分負(fù)頻率部分是對稱的判斷包圍圈數(shù)時(shí)要腦補(bǔ)完整曲線。如果曲線正好穿過 $(-1, j0)$說明系統(tǒng)臨界穩(wěn)定這時(shí)候margin返回的裕度接近零需要加阻尼或者降增益。奈奎斯特圖容易看花眼尤其是曲線繞來繞去的時(shí)候。我的習(xí)慣是先用margin看裕度再用nyquist確認(rèn)包圍關(guān)系兩者對不上就回去檢查傳遞函數(shù)有沒有寫錯(cuò)。另外nyquist支持ShowFullContour選項(xiàng)可以把負(fù)頻率也畫出來但圖會(huì)更亂一般不需要。3.3 根軌跡圖增益變化與極點(diǎn)走向根軌跡畫的是閉環(huán)極點(diǎn)隨開環(huán)增益 $K$ 從 0 到 $\infty$ 變化的軌跡。rlocus(sys)直接畫rlocfind可以交互式選點(diǎn)讀增益。figure; rlocus(sys_rlc); grid on; title(Root Locus of RLC System); % 交互式選取增益 % [K, poles] rlocfind(sys_rlc); % 取消注釋后鼠標(biāo)點(diǎn)選rlocus畫出的每條分支從一個(gè)開環(huán)極點(diǎn)出發(fā)終止于一個(gè)開環(huán)零點(diǎn)或無窮遠(yuǎn)。分支數(shù)等于極點(diǎn)數(shù)。讀圖時(shí)關(guān)注三點(diǎn)軌跡是否穿過虛軸決定穩(wěn)定增益范圍、分離點(diǎn)在哪里決定阻尼變化、漸近線角度決定高頻走向。rlocfind會(huì)返回你點(diǎn)選位置對應(yīng)的增益和閉環(huán)極點(diǎn)適合做“給定阻尼比反推增益”的活。如果系統(tǒng)有多個(gè)參數(shù)變化比如同時(shí)調(diào) $K_p$ 和 $K_i$根軌跡就不夠用了得用sisotool或者自己寫循環(huán)掃參數(shù)。根軌跡的局限性在于它只適合單參數(shù)變化多參數(shù)場景下容易誤導(dǎo)。3.4 零極點(diǎn)圖一眼看清左右半平面pzmap(sys)把零點(diǎn)和極點(diǎn)畫在復(fù)平面上零點(diǎn)用圓圈極點(diǎn)用叉號(hào)。這張圖信息密度最高但最容易被忽略。figure; pzmap(sys_rlc); grid on; title(Pole-Zero Map of RLC System); % 獲取零極點(diǎn)數(shù)值 z zero(sys_rlc); p pole(sys_rlc); disp(零點(diǎn)); disp(z); disp(極點(diǎn)); disp(p);讀零極點(diǎn)圖時(shí)右半平面的極點(diǎn)直接判不穩(wěn)定右半平面的零點(diǎn)會(huì)導(dǎo)致非最小相位特性比如初始反向響應(yīng)。虛軸上的極點(diǎn)說明系統(tǒng)有振蕩模態(tài)共軛極點(diǎn)對對應(yīng)一個(gè)二階振蕩環(huán)節(jié)。零極點(diǎn)距離虛軸越近對動(dòng)態(tài)響應(yīng)影響越大。如果零點(diǎn)和極點(diǎn)靠得很近可能互相抵消實(shí)際系統(tǒng)階數(shù)比傳遞函數(shù)看起來低這時(shí)候要考慮模型降階。4. 避坑與排查傳遞函數(shù)圖里最容易翻車的五個(gè)地方4.1 現(xiàn)象伯德圖幅頻曲線在低頻段翹起來原因傳遞函數(shù)里有積分環(huán)節(jié)$1/s$低頻增益趨于無窮MATLAB 默認(rèn)頻率范圍沒覆蓋到足夠低的頻段或者den最后一項(xiàng)寫成了非零值把積分環(huán)節(jié)變成了慣性環(huán)節(jié)。解決檢查den最后一項(xiàng)是否為零。如果是積分環(huán)節(jié)den [1 0]對應(yīng) $1/s$den [1 1]對應(yīng) $1/(s1)$。畫圖時(shí)用bode(sys, {1e-3, 1e3})把低頻范圍拉寬確認(rèn)斜率是不是 -20 dB/dec。4.2 現(xiàn)象奈奎斯特圖包圍圈數(shù)數(shù)不對原因MATLAB 默認(rèn)只畫正頻率負(fù)頻率部分需要手動(dòng)補(bǔ)全或者傳遞函數(shù)有虛軸上的極點(diǎn)nyquist會(huì)畫一個(gè)無窮大半徑的圓弧容易看漏。解決用nyquist(sys, {wmin, wmax})指定頻率范圍避開虛軸極點(diǎn)。如果極點(diǎn)就在虛軸上先做小擾動(dòng)移到左半平面再畫。判斷包圍圈數(shù)時(shí)從 $(-1, j0)$ 點(diǎn)向曲線引射線數(shù)穿越次數(shù)順時(shí)針減逆時(shí)針。4.3 現(xiàn)象根軌跡分支數(shù)對不上極點(diǎn)數(shù)原因傳遞函數(shù)分子分母有公因子minreal沒做MATLAB 把零極點(diǎn)對消的環(huán)節(jié)也算進(jìn)去了。解決畫圖前先sys_min minreal(sys)把可對消的零極點(diǎn)去掉。minreal的容差默認(rèn)是sqrt(eps)如果零極點(diǎn)靠得很近但不完全相等可以手動(dòng)指定容差minreal(sys, 1e-6)。4.4 現(xiàn)象離散系統(tǒng)伯德圖頻率軸單位混亂原因tf對象的Ts沒設(shè)對或者畫圖時(shí)用了bode但頻率范圍給的是 Hz 而不是 rad/s。解決確認(rèn)sys.Ts的值Ts -1是連續(xù)系統(tǒng)Ts 0是離散系統(tǒng)。畫圖時(shí)統(tǒng)一用 rad/s或者用bode(sys, {2*pi*fmin, 2*pi*fmax})把 Hz 轉(zhuǎn)成 rad/s。離散系統(tǒng)在奈奎斯特頻率附近幅頻曲線會(huì)翹這是雙線性變換的正常現(xiàn)象不要試圖“修平”。4.5 現(xiàn)象margin返回的裕度是負(fù)數(shù)但系統(tǒng)仿真穩(wěn)定原因margin算的是開環(huán)頻域裕度如果系統(tǒng)是非最小相位或者有延遲頻域裕度和時(shí)域穩(wěn)定性可能不一致。另外margin對多輸入多輸出系統(tǒng)只取第一個(gè)通道容易誤判。解決先用step(feedback(sys, 1))看階躍響應(yīng)確認(rèn)時(shí)域穩(wěn)定。如果是非最小相位系統(tǒng)頻域裕度只能做參考最終以時(shí)域仿真為準(zhǔn)。MIMO 系統(tǒng)用margin要指定通道或者改用sigma畫奇異值圖。5. 進(jìn)階技巧把傳遞函數(shù)圖用成設(shè)計(jì)工具5.1 用sisotool做交互式控制器設(shè)計(jì)sisotool(sys)打開控制系統(tǒng)設(shè)計(jì)器左邊是根軌跡或伯德圖右邊是階躍響應(yīng)可以拖拽零極點(diǎn)實(shí)時(shí)看效果。這個(gè)工具適合做 PID 整定和超前滯后校正。我一般先用sisotool粗調(diào)拿到零極點(diǎn)位置后再用pidtune或者手動(dòng)寫tf復(fù)現(xiàn)。% 打開交互式設(shè)計(jì)器 % sisotool(sys_rlc); % 取消注釋后運(yùn)行 % 用 pidtune 自動(dòng)整定 PID [C, info] pidtune(sys_rlc, PID); disp(C); disp(info); % 閉環(huán)階躍響應(yīng) sys_cl feedback(C * sys_rlc, 1); figure; step(sys_cl); grid on;pidtune返回的info里有相位裕度和截止頻率可以和你手動(dòng)設(shè)計(jì)的結(jié)果對比。pidtune默認(rèn)目標(biāo)是 60° 相位裕度如果系統(tǒng)本身阻尼很小整定出來的增益可能偏保守需要手動(dòng)加大。5.2 批量掃參數(shù)畫一族伯德圖做魯棒性分析時(shí)經(jīng)常要看好幾個(gè)參數(shù)下的頻響。用循環(huán)生成傳遞函數(shù)數(shù)組再一次性畫出來。R_values [5, 10, 20, 50]; % 不同電阻值 figure; hold on; for i 1:length(R_values) R_i R_values(i); den_i [L*C, R_i*C, 1]; sys_i tf(1, den_i); bode(sys_i, {10, 1e5}); end legend(R5, R10, R20, R50); grid on; hold off;hold on之后連續(xù)bode會(huì)把曲線疊在同一張圖上。注意每次bode都會(huì)刷新頻率范圍如果范圍不一致后面的圖會(huì)覆蓋前面的。解決辦法是第一次調(diào)用時(shí)指定范圍后續(xù)用bode(sys_i)不帶范圍或者統(tǒng)一用bode(sys_i, w)傳入固定的頻率向量w logspace(1, 5, 500)。5.3 從圖反推傳遞函數(shù)tfest與頻響數(shù)據(jù)如果手里只有實(shí)測頻響數(shù)據(jù)可以用tfest辨識(shí)傳遞函數(shù)。先把數(shù)據(jù)整理成iddata對象再指定極點(diǎn)數(shù)。% 假設(shè) freq 是頻率向量mag 和 phase 是幅值和相位 % resp mag .* exp(1j * deg2rad(phase)); % data idfrd(resp, freq, Ts); % sys_est tfest(data, 2); % 辨識(shí) 2 階傳遞函數(shù) % bode(sys_est, data); % 對比辨識(shí)結(jié)果和實(shí)測數(shù)據(jù)tfest的階數(shù)選擇很關(guān)鍵階數(shù)太低擬合不準(zhǔn)太高會(huì)過擬合。我一般從 2 階開始試看擬合百分比sys_est.Report.Fit.FitPercent到 85% 以上就夠用。辨識(shí)出來的傳遞函數(shù)要檢查零極點(diǎn)是否在左半平面右半平面的零極點(diǎn)可能是噪聲引起的需要手動(dòng)剔除。5.4 一個(gè)我常犯的錯(cuò)誤早年做電機(jī)控制時(shí)我把電流環(huán)的傳遞函數(shù)寫成 $1/(LsR)$畫伯德圖看截止頻率挺合適結(jié)果一上電就振蕩。后來才發(fā)現(xiàn)漏了 PWM 延遲和采樣保持環(huán)節(jié)實(shí)際傳遞函數(shù)多了一個(gè) $e^{-T_d s}$相位在截止頻率附近掉了 30 多度裕度根本不夠。從那以后我養(yǎng)成了一個(gè)習(xí)慣畫完理論伯德圖一定再疊一條含延遲的近似曲線用pade把延遲展開成有理函數(shù)對比兩條曲線的相位差。如果差超過 20°設(shè)計(jì)就得留更多裕度。這個(gè)習(xí)慣幫我省了不少后悔藥。希望幫到你。本文還有配套的精品資源點(diǎn)擊獲取