文章目錄

• 0.符號說明
• 1.如何根據(jù)連續(xù)系統(tǒng)建立差分方程
• • 1.1.獲取連續(xù)系統(tǒng)的傳遞函數(shù)
  • 1.2.獲取離散系統(tǒng)的傳遞函數(shù)
  • 1.3.轉(zhuǎn)換為差分方程
• 2.基本PID控制原理
• 3.比較PID輸出，分析參數(shù)產(chǎn)生的影響
• 4.改進PID算法（遇限削弱積分法）
• 5.simulink仿真

0.符號說明

y(k)——系統(tǒng)響應輸出的離散值

u(k)——數(shù)字PID控制輸出的離散值

r(k)——期望輸出的離散值（事先已知），在本例中為常數(shù)（即階躍輸入）

e(k)——e(k)=r(k)-y(k)，為期望值-實際值，是單位負反饋的誤差比較信號

注：圖片來源于百度百科

1.如何根據(jù)連續(xù)系統(tǒng)建立差分方程

1.1.獲取連續(xù)系統(tǒng)的傳遞函數(shù)

線性定常系統(tǒng)的控制中，PID是個非常常見的控制方式，如果可以通過Matlab仿真出PID的控制效果圖，那么對系統(tǒng)設計時的實時調(diào)試將會容易得多。在這里我們將會以一個利用系統(tǒng)辨識參數(shù)的PID設計為為例展示Matlab仿真PID的過程。
首先需要對一個未知的系統(tǒng)的參數(shù)進行辨識，以延遲環(huán)節(jié)可以忽略不計的電機調(diào)速系統(tǒng)為例。將時間戳導入xdata向量，對應的時刻轉(zhuǎn)速導入ydata向量，進行系統(tǒng)辨識

鏈接：Matlab的系統(tǒng)辨識

我們就以上文鏈接中辨識的系統(tǒng)傳遞函數(shù)為例：
$$G(s)=\frac{0.998}{0.021s+1}$$因此通過tf函數(shù)建立系統(tǒng)結(jié)構(gòu)體如下：

sys=tf(0.998,[0.021,1]); %建立被控對象傳遞函數(shù)，即式4.1

1.2.獲取離散系統(tǒng)的傳遞函數(shù)

由于是數(shù)字PID仿真，我們需要選取一個采樣時間，本案例選用的是0.005s（注意，采樣周期應該小于系統(tǒng)純滯后時間的0.1倍）。在對其進行數(shù)字PID控制前，我們需要將這個系統(tǒng)離散化：

ts=0.005; %采樣時間=0.005s dsys=c2d(sys,ts,'z'); %離散化

dsys即我們根據(jù)采樣周期離散化的Z變換系統(tǒng)。首先我們需要提取這個Z變化d那系統(tǒng)的參數(shù)方便后面的計算：

[num,den]=tfdata(dsys,'v');%'v'代表強制以向量的格式（默認為元胞數(shù)組）輸出num和den

1.3.轉(zhuǎn)換為差分方程

求解出的Z變換表達式為$dsys=\frac{num(1)?z+num(2)}{den(1)?z+den(2)}=\frac{0.2114}{z?0.7881}$
在PID仿真的過程中我們需要求解出時域表達式 ，因此需要借助差分方程解決，對于以下的Z變換：

\begin{equation}
Y(z)=dsys\cdot U(z)=\frac{num(2)}{den(1)\cdot z+den(2)}\cdot U(z)
\label{eq:Sample1}
\end{equation}

\begin{equation}
zY(z)+den(2)Y(z)=num(1)zU(z)+num(2)U(z)
\label{eq:Sample2}
\end{equation}
對上式進行反Z變換，可以得到以下的差分方程：

\begin{equation}
y(k+1)+den(2)y(k)=num(1)u(k+1)+num(2)u(k)
\label{eq:Sample3}
\end{equation}

\begin{equation}
y(k+1)=-den(2)y(k)+num(1)u(k+1)+num(2)u(k)
\label{eq:Sample4}
\end{equation}
位置型PID仿真時實際上可以不需要保存前一個數(shù)據(jù)(u(k)和y(k))，增量型PID必須要保存前一個數(shù)據(jù)。這里我們使用了位置型PID，但仍然利用$u_1$和$y_1$保存了上一個數(shù)據(jù)，僅僅是為了演示這一過程。\begin{equation}
y(k+1)=-den(2)y(k)+num(1)u(k+1)+num(2)u(k)
\end{equation}
可以轉(zhuǎn)換為下面的式子：
\begin{equation}
y(k)=-den(2)y_1+num(1)u(k)+num(2)u_1
\label{eq:Sample5}
\end{equation}
我們的差分方程就這樣建立完畢。注意，此差分方程僅僅是描述系統(tǒng)模型的運算規(guī)律的，和我們的控制無關。因此是y(k)和u(k)的映射關系。我們下面的控制則是利用負反饋信號e(k)導出u(k)的輸出，求解的是控制器u(k)的序列值。

2.基本PID控制原理

以位置型PID控制為例。將連續(xù)的PID控制轉(zhuǎn)換為數(shù)字式時，微分環(huán)節(jié)被用差分代替，積分環(huán)節(jié)被累加和代替，比例環(huán)節(jié)則保持不變。差分的實現(xiàn)非常簡單，只需要用$e(k+1)?e(k)$即$e(k)?e_1$等效即可。積分的實現(xiàn)在每一次運算的后面都累加原來的誤差，即Ee=Ee+e_1;即可。PID的控制器輸出$$u(k)=Kp?e(k)+Kd?(e(k)?e_1)+Ki?Ee$$
PID控制器構(gòu)造完畢，我們需要通過r(k)和y(k)得到e(k)，再通過e(k)得出u(k)，進而再求解出y(k)，再結(jié)合r(k)求解出e(k),…以此循環(huán)，求解出離散的響應點。
詳細的代碼如下：

ts=0.005; %采樣時間=0.005s sys=tf(0.998,[0.021,1]); %建立被控對象傳遞函數(shù)，即式4.1 dsys=c2d(sys,ts,'z'); %離散化 [num,den]=tfdata(dsys,'v'); % e_1=0; %前一時刻的偏差 Ee=0; %累積偏差 u_1=0.0; %前一時刻的控制量 y_1=0; %前一時刻的輸出 %PID參數(shù) kp=0.22; ki=0.13; kd=0; u=zeros(1,1000);%預先分配內(nèi)存 time=zeros(1,1000);%時刻點（設定1000個） for k=1:1:1000time(k)=k*ts; %時間參數(shù)r(k)=1500; %期望值y(k)=-1*den(2)*y_1+num(2)*u_1+num(1)*u(k);%系統(tǒng)響應輸出序列e(k)=r(k)-y(k); %誤差信號u(k)=kp*e(k)+ki*Ee+kd*(e(k)-e_1); %系統(tǒng)PID控制器輸出序列Ee=Ee+e(k); %誤差的累加和u_1=u(k); %前一個的控制器輸出值y_1=y(k); %前一個的系統(tǒng)響應輸出值e_1=e(k); %前一個誤差信號的值 end %（僅繪制過渡過程的曲線，x坐標限制為[0,1]） p1=plot(time,r,'-.');xlim([0,1]);hold on;%指令信號的曲線（即期望輸入） p2=plot(time,y,'--');xlim([0,1]);%不含積分分離的PID曲線 hold on;

輸出的PID控制曲線如下：

3.比較PID輸出，分析參數(shù)產(chǎn)生的影響

一個基本的PID就完成了。下面如果我們想要知道修改PID的三個參數(shù)kp,ki,kd會帶來什么效果，只需要在程序中修改即可。為了方便起見，我們建立一個PID的數(shù)組，kp,ki,kd每次都取數(shù)組的一個值，然后設定一個大循環(huán)開始循環(huán)仿真。再利用subplot輸出子圖的方式將所有的PID效果都輸出到一個圖進行對比。該代碼根據(jù)上述代碼修改已經(jīng)很容易，PID比較圖的代碼如下：

close all PID=[0.22,0.13,0;0.4,0.13,0;0.4,0.25,0;0.8,0.23,0.4;0.8,0.2,1;0.7,0.2,0.9];%初始化PID參數(shù) for pid=1:1:6 ts=0.005; %采樣時間=0.005s sys=tf(0.998,[0.021,1]); %建立被控對象傳遞函數(shù)，即式4.1 dsys=c2d(sys,ts,'z'); %離散化 [num,den]=tfdata(dsys,'v'); % e_1=0; %前一時刻的偏差 Ee=0; %累積偏差 u_1=0.0; %前一時刻的控制量 y_1=0; %前一時刻的輸出 %PID參數(shù) kp=PID(pid,1); ki=PID(pid,2); kd=PID(pid,3); u=zeros(1,1000); time=zeros(1,1000); for k=1:1:1000time(k)=k*ts; %時間參數(shù)r(k)=1500; %給定量y(k)=-1*den(2)*y_1+num(2)*u_1+num(1)*u(k);e(k)=r(k)-y(k); %偏差u(k)=kp*e(k)+ki*Ee+kd*(e(k)-e_1); Ee=Ee+e(k); u_1=u(k); y_1=y(k); e_1=e(k); end subplot(2,3,pid); p1=plot(time,r,'-.');xlim([0,1]);hold on; p2=plot(time,y,'--');xlim([0,1]); title(['Kp=',num2str(kp),' Ki=',num2str(ki),' Kd= ',num2str(kd)]); hold on; end

輸出的子圖矩陣如下：

可以發(fā)現(xiàn)，修改Kp會造成上升時間的縮短，但是有可能也會帶來較大的超調(diào)。積分的增加是一個嚴重的滯后環(huán)節(jié)，會減小相位裕度，也會帶來超調(diào)（超調(diào)量并不是絕對的，相對于較小的Kp可能會產(chǎn)生較大的超調(diào)，而Kp較大時超調(diào)會減小（例如第一行的1圖和2圖的對比））。然而積分的引入也是必要的，否則將會很長時間無法削弱誤差e(k)（例如第二行第二個圖）。微分的引入相當于一個超前校正，會減少超調(diào)，但是過渡的微分很可能會造成尾部振蕩，系統(tǒng)逐漸變得不穩(wěn)定。因此微分和積分之間需要一個平衡，當滿足這個平衡的時候，系統(tǒng)幾乎沒有振蕩，同時響應速度也較快。（第一行的圖3是積分過多，產(chǎn)生超調(diào)，第二行的圖1和圖3就比較理想）
綜合上述，PID的調(diào)節(jié)經(jīng)驗可以歸結(jié)為以下幾點：

• Kp較小時，系統(tǒng)對微分和積分環(huán)節(jié)的引入較為敏感，積分會引起超調(diào)，微分可能會引起振蕩，而振蕩劇烈的時候超鐵也會增加。
• Kp增大時，積分環(huán)節(jié)由于滯后產(chǎn)生的超調(diào)逐漸減小，此時如果想要繼續(xù)減少超調(diào)可以適當引入微分環(huán)節(jié)。繼續(xù)增大Kp系統(tǒng)可能會不太穩(wěn)定，因此在增加Kp的同時引入Kd減小超調(diào)，可以保證在Kp不是很大的情況下也能取得較好的穩(wěn)態(tài)特性和動態(tài)性能。
• Kp較小時，積分環(huán)節(jié)不宜過大，Kp較大時積分環(huán)節(jié)也不宜過小（否則調(diào)節(jié)時間會非常地長），在下面這個例子中我們還會介紹到，當使用分段PID，在恰當?shù)臈l件下分離積分，可以取得更好的控制效果。原因在于在穩(wěn)態(tài)誤差即將滿足要求時，消除了系統(tǒng)的滯后。因此系統(tǒng)超調(diào)會明顯減少。本例中采樣的抗積分飽和的方法是遇限削弱積分法。

4.改進PID算法（遇限削弱積分法）

遇限削弱積分法的原理是
當$u(k)>u_{max}$時，若e(k)>0即輸出值還未到達指定值，則認為積分會帶來滯后，不再積分。
當$u(k)<0$時，若e(k)<0即輸出值超過了指定值，則認為積分會帶來滯后，不再積分。
在本案例中認為$u_{max}=r(k)$
改進PID算法如下（需要些兩個循環(huán)，當然也可以用一個循環(huán)，將其中的PID設為一個子過程調(diào)用）：

close all ts=0.005; %采樣時間=0.005s sys=tf(0.998,[0.021,1]); %建立被控對象傳遞函數(shù)，即式4.1 dsys=c2d(sys,ts,'z'); %離散化 [num,den]=tfdata(dsys,'v'); % e_1=0; %前一時刻的偏差 Ee=0; %累積偏差 u_1=0.0; %前一時刻的控制量 y_1=0; %前一時刻的輸出 %PID參數(shù) kp=0.22; ki=0.13; kd=0; u=zeros(1,1000); time=zeros(1,1000); for k=1:1:1000time(k)=k*ts; %時間參數(shù)r(k)=1500; %給定量y(k)=-1*den(2)*y_1+num(2)*u_1+num(1)*u(k);e(k)=r(k)-y(k); %偏差u(k)=kp*e(k)+ki*Ee+kd*(e(k)-e_1); Ee=Ee+e(k); u_1=u(k); y_1=y(k); e_1=e(k); end p1=plot(time,r,'-.');xlim([0,1]);hold on; p2=plot(time,y,'--');xlim([0,1]); hold on; a=1;%控制積分分離的二值數(shù) e_1=0;Ee=0;u_1=0.0;y_1=0;%重新初始化 for k=1:1:1000time(k)=k*ts; %時間參數(shù)r(k)=1500; %給定量y(k)=-1*den(2)*y_1+num(2)*u_1;e(k)=r(k)-y(k); %偏差u(k)=kp*e(k)+ki*Ee+kd*(e(k)-e_1); if ((u(k)>r(k)) && (e(k)>0))||((u(k)<0) && (e(k)<0))a=0;else a=1;end Ee=Ee+a*e(k); u_1=u(k); y_1=y(k); e_1=e(k); end p3=plot(time,y,'-');xlim([0,1]); title('含積分分離與不含積分分離的對比'); legend([p1,p2,p3],'指令信號','不含積分分離','含積分分離');

輸出的曲線對比圖如下：

可以發(fā)現(xiàn)，系統(tǒng)的超調(diào)量明顯減少了，調(diào)節(jié)時間也減少了一點。原因在于我們采用了分段PID的手段，既消除了穩(wěn)態(tài)誤差還削弱了積分環(huán)節(jié)帶來的滯后影響。

5.simulink仿真

需要的模塊名稱（不區(qū)分大小寫）如下：

• gain（參數(shù)分別為0.22和0.13/0.005）
• sum（參數(shù)分別為"|±"和"|++"）
• integrator
• scope
  注意：本文使用的是離散PID仿真，而simulink使用的是連續(xù)系統(tǒng)仿真，轉(zhuǎn)換PID參數(shù)時P參數(shù)不變，I參數(shù)應該除以仿真間隔Ts=0.005，D參數(shù)應該乘Ts。

以表中第一組PI參數(shù)為例：

得到的示波器曲線如下：

希望本文對您有幫助，謝謝閱讀。

總結(jié)

以上是生活随笔為你收集整理的Matlab仿真PID控制（带M文件、simulink截图和参数分析）的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

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