久久久久久久999_99精品久久精品一区二区爱城_成人欧美一区二区三区在线播放_国产精品日本一区二区不卡视频_国产午夜视频_欧美精品在线观看免费

標(biāo)題: MATLAB Pid仿真實(shí)驗(yàn) [打印本頁]

作者: 51黑ss 時(shí)間: 2016-3-30 23:27
標(biāo)題: MATLAB Pid仿真實(shí)驗(yàn)

一實(shí)驗(yàn)要求:

學(xué)習(xí)PID控制器的基本原理，掌握PID參數(shù)的物理調(diào)節(jié)規(guī)律，通過仿真驗(yàn)證PID調(diào)節(jié)過程。選擇例題中的任意兩個(gè)實(shí)驗(yàn)進(jìn)行仿真研究，改變參數(shù)后觀察系統(tǒng)的動(dòng)態(tài)性能和穩(wěn)態(tài)性能的變化，確認(rèn)是否和理論分析一致。

二實(shí)驗(yàn)內(nèi)容

實(shí)驗(yàn)一：選擇chap-9.m作為pid調(diào)節(jié)。

1代碼：

%PID Controller with changing integration rate
clear all;
close all;
%Big time delay Plant
ts=20;
sys=tf([1],[60,1],'inputdelay',80);
dsys=c2d(sys,ts,'zoh');
[num,den]=tfdata(dsys,'v');

u_1=0;u_2=0;u_3=0;u_4=0;u_5=0;
y_1=0;y_2=0;y_3=0;
error_1=0;error_2=0;
ei=0;

for k=1:1:200
time(k)=k*ts;
rin(k)=1.0;  %Step Signal
%Linear model
yout(k)=-den(2)*y_1+num(2)*u_5;
error(k)=rin(k)-yout(k);
kp=0.45;kd=12;ki=0.0048;
A=0.4;B=0.6;
%T type integration
ei=ei+(error(k)+error_1)/2*ts;
M=1;
if M==1    %Changing integration rate
if abs(error(k))<=B
f(k)=1;
elseif abs(error(k))>B&abs(error(k))<=A+B
f(k)=(A-abs(error(k))+B)/A;
else
f(k)=0;
end
elseif M==2  %Not changing integration rate
f(k)=1;
end
u(k)=kp*error(k)+kd*(error(k)-error_1)/ts+ki*f(k)*ei;
if u(k)>=10
u(k)=10;
end
if u(k)<=-10
u(k)=-10;
end
%Return of PID parameters
u_5=u_4;u_4=u_3;u_3=u_2;u_2=u_1;u_1=u(k);
y_3=y_2;y_2=y_1;y_1=yout(k);
error_2=error_1;
error_1=error(k);
end
figure(1);
plot(time,rin,'b',time,yout,'r');
xlabel('time(s)');ylabel('rin,yout');
figure(2);
plot(time,f,'r');
xlabel('time(s)');ylabel('Integration rate f');

2實(shí)驗(yàn)結(jié)果

在沒有更改代碼的情況下，實(shí)驗(yàn)的效果如圖一所示。

圖一
改變kp

(1)增大kp,理論上可以看到pid調(diào)節(jié)的更快更陡。實(shí)際效果如圖2。Kp=3.

圖2

實(shí)際上產(chǎn)生了自激震蕩，說明kp不能調(diào)的過大，否則系統(tǒng)會不穩(wěn)定。

（2）減小kp，理論上系統(tǒng)調(diào)節(jié)緩慢，會產(chǎn)生超調(diào)量，調(diào)節(jié)曲線不平滑，有毛刺。Kp=0.1

圖3

實(shí)際上也證明比例系數(shù)太小，系統(tǒng)變化時(shí)調(diào)節(jié)相對緩慢。

改變kd

(1)增大kd，理論上系統(tǒng)穩(wěn)定時(shí)波動(dòng)很大，精度不高，且有很多毛刺與抖動(dòng)。Kd=30.

圖4

實(shí)際上系統(tǒng)開始調(diào)節(jié)時(shí)會有些尖峰，調(diào)節(jié)曲線不平滑。

（2）減小kd，理論上系統(tǒng)調(diào)節(jié)穩(wěn)定時(shí)變化緩慢，會有些超調(diào)量。Kd=3.

圖5

實(shí)際上也可以看出減小kd對系統(tǒng)影響不大，只是在系統(tǒng)剛進(jìn)入穩(wěn)定狀態(tài)時(shí)有一個(gè)小小的超調(diào)量。

改變ki

(1)增大ki,積分變大，系統(tǒng)系統(tǒng)變得不穩(wěn)定，抖動(dòng)波動(dòng)很厲害。

圖6

實(shí)際上可以其變化幅度大，超調(diào)量很高，系統(tǒng)很不穩(wěn)定。

（2）減小ki,系統(tǒng)調(diào)節(jié)變慢，系統(tǒng)可能無法調(diào)節(jié)。

圖7

系統(tǒng)無法達(dá)到設(shè)定的穩(wěn)態(tài)值，無法調(diào)節(jié)。

實(shí)驗(yàn)二：選擇char2_1.m作為實(shí)驗(yàn)

1代碼：

%Series System Control
clear all;
close all;
ts=2;
sys1=tf(1,[10,1]);
dsys1=c2d(sys1,ts,'z');
[num1,den1]=tfdata(dsys1,'v');
sys2=tf(1,[10,1]);
dsys2=c2d(sys2,ts,'z');
[num2,den2]=tfdata(dsys2,'v');
dph=1/zpk('z',ts);
Gc2=dph/(dsys2*(1-dph));
[nump,denp]=tfdata(Gc2,'v');
u1_1=0.0;u2_1=0.0;
y1_1=0;y2_1=0;
e2_1=0;ei=0;
for k=1:1:2000
time(k)=k*ts;
r1(k)=1;
%Linear model
y1(k)=-den1(2)*y1_1+num1(2)*y2_1;  %Main plant
y2(k)=-den2(2)*y2_1+num2(2)*u2_1;  %Assistant plant
error(k)=r1(k)-y1(k);
ei=ei+error(k);
u1(k)=1.2*error(k)+0.02*ei; %Main Controller
e2(k)=u1(k)-y2(k);          %Assistant Controller
u2(k)=-denp(2)*u2_1+nump(1)*e2(k)+nump(2)*e2_1;
d2(k)=0.01*rands(1);
u2(k)=u2(k)+d2(k);
%----------Return of PID parameters------------
u1_1=u1(k);
u2_1=u2(k);
e2_1=e2(k);
y1_1=y1(k);
y2_1=y2(k);
end
figure(1);    %Assistant Control
plot(time,u1,'b',time,y2,'r');
xlabel('time(s)');ylabel('u1,y2');
figure(2);    %Main Control
plot(time,r1,'b',time,y1,'r');
xlabel('time(s)');ylabel('r1,y1');
figure(3);
plot(time,d2,'r');
xlabel('time(s)');ylabel('disturbance');

2實(shí)驗(yàn)效果

在沒有更改遠(yuǎn)代碼的情況下，實(shí)驗(yàn)的結(jié)果如下所示

圖8

第3個(gè)圖是隨機(jī)數(shù)的隨時(shí)間的分布。

同理，改變kp，ki的值，系統(tǒng)會發(fā)生改變。

改變kp

（1）增大kp，系統(tǒng)變化會更快，會有超調(diào)量。Kp=3.

圖9

實(shí)際沒有發(fā)現(xiàn)系統(tǒng)有超調(diào)量，但在系統(tǒng)開始調(diào)節(jié)得時(shí)候會會產(chǎn)生尖峰脈沖，系統(tǒng)調(diào)節(jié)很快。

（2）減小kp，系統(tǒng)調(diào)節(jié)緩慢，但系統(tǒng)穩(wěn)定。

圖10

結(jié)果可以看出系統(tǒng)調(diào)節(jié)很穩(wěn)定，但調(diào)節(jié)時(shí)間有點(diǎn)長。

該變kd

（1）增大ki,系統(tǒng)很快達(dá)到穩(wěn)定，調(diào)節(jié)速度快.ki=0.2

圖11

結(jié)果顯示系統(tǒng)在轉(zhuǎn)換開始時(shí)系統(tǒng)產(chǎn)生了尖峰脈沖，然后很快地達(dá)到穩(wěn)定。

（2）減小ki，系統(tǒng)會很慢的達(dá)到穩(wěn)定.ki=0.005。

如圖12

結(jié)果顯示系統(tǒng)經(jīng)過相當(dāng)漫長的時(shí)間達(dá)到了穩(wěn)定狀態(tài)。

三實(shí)驗(yàn)感想

通過此在MATLAB上運(yùn)用pid，對pidd的實(shí)際作用有了基本了解。對于kp,要選取的合適，不能過大，否則會產(chǎn)生震蕩，過小可能會產(chǎn)生超常量，調(diào)節(jié)緩慢。對于kd,太大了產(chǎn)生毛刺抖動(dòng)，太小了會有些超常量。對于ki，太大了調(diào)節(jié)很快，但會產(chǎn)生過沖，尖峰脈沖，太小了調(diào)節(jié)很慢，可能無法調(diào)試達(dá)到穩(wěn)定狀態(tài)。

總結(jié)這三點(diǎn)，在調(diào)試pid時(shí)，kp應(yīng)從小到大調(diào)試，kd從小到大調(diào)試，ki從小到大調(diào)試。每次只調(diào)一個(gè)參數(shù)，直到最滿意的時(shí)候，再調(diào)下個(gè)參數(shù)。

歡迎光臨 (http://www.zg4o1577.cn/bbs/)