1、精選優(yōu)質(zhì)文檔-傾情為你奉上成績(jī)廣州航海學(xué)院計(jì)算機(jī)控制的理論基礎(chǔ) 實(shí)驗(yàn)報(bào)告專(zhuān)心-專(zhuān)注-專(zhuān)業(yè)專(zhuān)業(yè)班級(jí) 電氣16 姓名學(xué)號(hào) 馮 實(shí)驗(yàn)名稱(chēng) PID控制及仿真 實(shí)驗(yàn)日期 2019.12.20 實(shí)驗(yàn)組號(hào) 指導(dǎo)教師 李老師 （報(bào)告內(nèi)容包括：實(shí)驗(yàn)?zāi)康摹?shí)驗(yàn)設(shè)備及器材、實(shí)驗(yàn)步驟、實(shí)驗(yàn)數(shù)據(jù)、圖標(biāo)及曲線處理、實(shí)驗(yàn)小結(jié)等）一、實(shí)驗(yàn)?zāi)康?了解常用的計(jì)算機(jī)控制算法；2熟悉模擬PID控制的基本原理及控制參數(shù)的作用；3掌握數(shù)字PID控制算法及其仿真。二、實(shí)驗(yàn)設(shè)備及器材PC機(jī)一臺(tái)三、實(shí)驗(yàn)步驟及內(nèi)容1. 位置式PID控制仿真1.1位置式PID控制程序如下：1.1.1創(chuàng)建PlantModel的M文件并保存function dy=P

2、lantModel(t,y,flag,para) %定義PlantModel函數(shù)u=para;J=0.0067;B=0.1;dy=zeros(2,1);dy(1)=y(2);dy(2)=-(B/J)*y(2)+(1/J)*u;1.1.2位置式PID程序如下clear all; %清除所有close all; %關(guān)閉所有ts=0.001; %設(shè)置采樣時(shí)間xk=zeros(2,1); %設(shè)置矩陣e_1=0;u_1=0;for k=1:1:2000 %設(shè)置循環(huán)函數(shù)以及采樣時(shí)刻 time(k)=k*ts; rin(k)=0.50*sin(1*2*pi*k*ts); %設(shè)置輸入為正弦函數(shù) para=u_

3、1; tSpan=0 ts; tt,xx=ode45('PlantModel',tSpan,xk, ,para); %調(diào)用PlantModel函數(shù) xk=xx(length(xx),:); %輸出xx的長(zhǎng)度 yout(k)=xk(1); e(k)=rin(k)-yout(k); %求出e(k)的表達(dá)式 de(k)=(e(k)-e_1)/ts; %求出de(k)的表達(dá)式 u(k)=20.0*e(k)+0.50*de(k); %求出u(k)的表達(dá)式 if u(k)>10.0 %通過(guò)if函數(shù)取u（k）上限值或者下限值 u(k)=10.0; end u_1=u(k); e_1=e

4、(k);endfigure(1); %第一漲圖片的命名plot(time,yout,'k'); %輸出橫坐標(biāo)為時(shí)間t，縱坐標(biāo)為輸出響應(yīng)，并且曲線為黑色xlabel('time(s)'),ylabel('yout');figure(2); %第一漲圖片的命名plot(time,rin-yout,'k'); %輸出橫坐標(biāo)為時(shí)間t，縱坐標(biāo)為偏差變化，并且曲線為黑色xlabel('time(s)'),ylabel('error'); %設(shè)置橫坐標(biāo)為時(shí)間t，縱坐標(biāo)為偏差變化1.1.3仿真結(jié)果如下輸出相應(yīng)曲線偏

5、差曲線1.2位置式PID控制小結(jié)通過(guò)仿真結(jié)果可以看出只要PID參數(shù)設(shè)置合適，參考輸入為正弦函數(shù)信號(hào)，系統(tǒng)輸出能夠較好地跟蹤參考輸入。2. 增量式PID控制仿真。2.1增量式PID控制程序如下clear all;%清除所有close all; %關(guān)閉所有ts=0.001; %設(shè)置采樣時(shí)間sys=tf(400,1,50,0); %設(shè)置連續(xù)對(duì)象的傳遞函數(shù)dsys=c2d(sys,ts,'z'); %傳遞函數(shù)進(jìn)行z變換num,den=tfdata(dsys,'v'); %設(shè)置傳遞函數(shù)的分子分母變量u_1=0.0;u_2=0.0;u_3=0.0;y_1=0;y_2=0;y

6、_3=0;x=0,0,0'error_1=0;error_2=0;for k=1:1:1000 %設(shè)置仿真時(shí)間 time(k)=k*ts; rin(k)=1.0; %設(shè)置輸入為單位階躍信號(hào) kp=8; %分別設(shè)置比例、積分、微分參數(shù) ki=0.10; kd=10; du(k)=kp*x(1)+kd*x(2)+ki*x(3); %求出du（k）的表達(dá)式 u(k)=u_1+du(k); %求出u(k)的表達(dá)式 if u(k)>=10 %通過(guò)if函數(shù)取u（k）上限值或者下限值 u(k)=10; endif u(k)<=-10 u(k)=-10;endyout(k)=-den(2)

7、*y_1-den(3)*y_2+num(2)*u_1+num(3)*u_2;error=rin(k)-yout(k); %計(jì)算偏差u_3=u_2;u_2=u_1;u_1=u(k); %數(shù)據(jù)儲(chǔ)存y_3=y_2;y_2=y_1;y_1=yout(k);x(1)=error-error_1; %計(jì)算比例項(xiàng)x(2)=error-2*error_1+error_2; %計(jì)算積分項(xiàng)x(3)=error; %計(jì)算微分項(xiàng)error_2=error_1;error_1=error;endplot(time,rin,'b',time,yout,'r'); %輸出橫坐標(biāo)為時(shí)間t，縱坐標(biāo)為階躍信號(hào)，并且曲線為紅色xlabel('time(s)');ylabel('rin,yout'); %設(shè)置橫坐標(biāo)為時(shí)間t，縱坐標(biāo)為階躍信號(hào)，2.2仿真結(jié)果如下2.3增量式PID控制小結(jié)由仿真結(jié)果可以看出參考輸入為階躍信號(hào)，只要增量式PID參數(shù)設(shè)置合適，系統(tǒng)輸出也能夠較好地跟蹤參考輸入。四、實(shí)驗(yàn)總結(jié)通過(guò)本次實(shí)驗(yàn)，了解常用的計(jì)算機(jī)控制算法；熟悉模擬PID控制的基本原理及控制參數(shù)的作用；掌握數(shù)字PID控制算法及其仿真。同時(shí)