




版權說明:本文檔由用戶提供并上傳,收益歸屬內容提供方,若內容存在侵權,請進行舉報或認領
文檔簡介
1、參考教材中例子設計一包含了模糊技術與PID技術的混合智能控制器,其被控對象為:采樣時間為1ms,編寫matlab仿真程序,確定其在階躍輸入的響應結果,并與經典PID控制仿真結果相比較。要求詳細描述控制系統(tǒng)的設計,控制系統(tǒng)工作流程,模糊系統(tǒng)中的輸入輸出的隸屬函數(shù)設計及其采用的模糊規(guī)則,分析仿真結果并進行總結。表1 kp的模糊規(guī)則表eckpeNBNMNSZOPSPMPBNBPBPBPMPMPSZOZONMPBPBPMPSPSZONSNSPMPMPMPSZONSNSZOPMPMPSZONSNMNMPSPSPSZONSNSNMNMPMPSZONSNMNMNMNBPBZOZONMNMNMNBNB表2 k
2、i的模糊規(guī)則表kieceNBNMNSZOPSPMPBNBNBNBNMNMNSZOZONMNBNBNMNSNSZOZONSNBNMNSNSZOPSPSZONMNMNSZOPSPMPMPSMNNSZOPSPSPMPBPMZOZOPSPSPMPBPBPBZOZOPSPMPMPBPB表3 kd的模糊規(guī)則表kdeceNBNMNSZOPSPMPBNBPSNSNBNBNBNMPSNMPSNSNBNMNMNSZONSZONSNMNMNSNSZOZOZONSNSNSNSNSZOPSZOZOZOZOZOZOZOPMPBNSPSPSPSPSPBPBPBPMPMPSPSPSPBKp,ki,kd的模糊控制規(guī)則表建立好以
3、后,可根據(jù)以下方法進行kp,ki,kd的自適應校正。將系統(tǒng)誤差e和誤差變化ec變化范圍定義為模糊集上的論域,即e,ec = -3,-2,-1,0,1,2,3,其模糊子集為e,ec = NB,NM,NS,ZO,PS,PM,PB,子集中元素分別代表負大,負中,負小,零,正小,正中,正大。應用模糊合成推理設計PID參數(shù)的整定算法。第k個采樣時間的整定為Kp(k)=kp0+kp(k)Ki(k)=ki0+ki(k)Kd(k)=kd0+kd(k)在線運行過程中,控制系統(tǒng)通過對模糊邏輯規(guī)則的結果處理、查表和運算,完成對PID參數(shù)的在線自校正。其工作流程圖如下圖所示。入口e(k),ec(k)模糊化e(k-1)
4、=e(k)PID控制器輸出計算當前kp,ki,kd模糊整定kp,ki,kdec(k)=e(k)-e(k-1)e(k)=r(k)-y(k)取當前采樣值返回圖1 誤差的隸屬函數(shù)圖2 誤差變化率的隸屬函數(shù)圖3 kp的隸屬函數(shù)圖4 ki的隸屬函數(shù)圖5 kd的隸屬函數(shù)圖6 模糊系統(tǒng)fuzzpid.fis的結構圖7 模糊推理系統(tǒng)的動態(tài)仿真環(huán)境在程序PID_b.m中,利用所設計的模糊系統(tǒng)fuzzpid.fis進行PID控制參數(shù)的整定,并利用模糊PID控制進行階躍響應,在第300個采樣時間時控制器輸出端加上1.0的干擾,響應結果及PID控制參數(shù)的自適應變化如圖8到13所示。圖8 模糊PID控制階躍響應圖9 模
5、糊PID控制誤差響應圖10 控制器輸入u圖11 kp的自適應調整圖12 ki的自適應調整圖13 kd的自適應調整在對三階線性系統(tǒng)的控制中,利用穩(wěn)定邊界法進行參數(shù)整定的經典PID控制的超調量比模糊PID控制的超調量要大,但模糊PID控制存在一定的穩(wěn)態(tài)誤差。模糊控制用模糊集合和模糊概念描述過程系統(tǒng)的動態(tài)特性,根據(jù)模糊集和模糊邏輯來做出控制決策,它在解決復雜控制問題方面有很大的潛力,可以動態(tài)地適應外界環(huán)境的變化。附錄1 模糊系統(tǒng)設計程序PID_a%Fuzzy Turnning PID Controlclear all;close all;a=newfis('fuzzpid');a=a
6、ddvar(a,'input','e',-3,3); a=addmf(a,'input',1,'NB','zmf',-3,-1);a=addmf(a,'input',1,'NM','trimf',-3,-2,0);a=addmf(a,'input',1,'NS','trimf',-3,-1,0);a=addmf(a,'input',1,'Z','trimf',-2,0,2)
7、;a=addmf(a,'input',1,'PS','trimf',-1,1,3);a=addmf(a,'input',1,'PM','trimf',0,2,3);a=addmf(a,'input',1,'PB','smf',1,3);a=addvar(a,'input','ec',-3,3);a=addmf(a,'input',2,'NB','zmf',-3,-1);a=
8、addmf(a,'input',2,'NM','trimf',-3,-2,0);a=addmf(a,'input',2,'NS','trimf',-3,-1,0);a=addmf(a,'input',2,'Z','trimf',-2,0,2);a=addmf(a,'input',2,'PS','trimf',-1,1,3);a=addmf(a,'input',2,'PM',
9、'trimf',0,2,3);a=addmf(a,'input',2,'PB','smf',1,3);a=addvar(a,'output','kp',-0.3,0.3);a=addmf(a,'output',1,'NB','zmf',-0.3,-0.1);a=addmf(a,'output',1,'NM','trimf',-0.3,-0.2,0);a=addmf(a,'output',1
10、,'NS','trimf',-0.3,-0.1,0);a=addmf(a,'output',1,'Z','trimf',-0.2,0,0.2);a=addmf(a,'output',1,'PS','trimf',-0.1,0.1,0.3);a=addmf(a,'output',1,'PM','trimf',0,0.2,0.3);a=addmf(a,'output',1,'PB','
11、smf',0.1,0.3);a=addvar(a,'output','ki',-0.06,0.06);a=addmf(a,'output',2,'NB','zmf',-0.06,-0.02);a=addmf(a,'output',2,'NM','trimf',-0.06,-0.04,0);a=addmf(a,'output',2,'NS','trimf',-0.06,-0.02,0.02);a=addmf(a,&
12、#39;output',2,'Z','trimf',-0.04,0,0.04);a=addmf(a,'output',2,'PS','trimf',-0.02,0.02,0.06);a=addmf(a,'output',2,'PM','trimf',0,0.04,0.06);a=addmf(a,'output',2,'PB','smf',0.02,0.06);a=addvar(a,'output'
13、,'kd',-3,3);a=addmf(a,'output',3,'NB','zmf',-3,-1);a=addmf(a,'output',3,'NM','trimf',-3,-2,0);a=addmf(a,'output',3,'NS','trimf',-3,-1,0);a=addmf(a,'output',3,'Z','trimf',-2,0,2);a=addmf(a,'out
14、put',3,'PS','trimf',-1,1,3);a=addmf(a,'output',3,'PM','trimf',0,2,3);a=addmf(a,'output',3,'PB','smf',1,3);rulelist = 1 1 7 1 5 1 1;1 2 3 1 3 1 1;1 3 6 2 1 1 1;1 4 6 2 1 1 1;1 5 5 3 1 1 1;1 6 4 4 2 1 1;1 7 4 4 5 1 1;2 1 7 1 5 1 1;2 2
15、 7 1 3 1 1;2 3 6 2 1 1 1;2 4 5 3 2 1 1;2 5 5 3 2 1 1;2 6 4 4 3 1 1;2 7 3 4 4 1 1;3 1 6 1 4 1 1;3 2 6 2 3 1 1;3 3 6 3 2 1 1;3 4 5 3 2 1 1;3 5 4 4 3 1 1;3 6 3 5 3 1 1;3 7 3 5 4 1 1;4 1 6 2 4 1 1;4 2 6 2 4 1 1;4 3 5 3 3 1 1;4 4 4 4 3 1 1;4 5 3 5 3 1 1;4 6 2 6 3 1 1;4 7 2 6 4 1 1;5 1 5 2 4 1 1;5 2 5 3 4
16、 1 1;5 3 4 4 4 1 1;5 4 3 5 4 1 1;5 5 3 5 4 1 1;5 6 2 6 4 1 1;5 7 2 7 4 1 1;6 1 5 4 7 1 1;6 2 4 4 5 1 1;6 3 3 5 5 1 1;6 4 2 5 5 1 1;6 5 2 6 5 1 1;6 6 2 7 5 1 1;6 7 1 7 7 1 1;7 1 4 4 7 1 1;7 2 4 4 6 1 1;7 3 2 5 6 1 1;7 4 2 6 6 1 1;7 5 2 6 5 1 1;7 6 1 7 5 1 1;7 7 1 7 7 1 1;a = addrule(a,rulelist);a = s
17、etfis(a,'DefuzzMethod','centroid');writefis(a,'fuzzpid');a = readfis('fuzzpid');figure(1);plotmf(a,'input',1);figure(2);plotmf(a,'input',2);figure(3);plotmf(a,'output',1);figure(4);plotmf(a,'output',2);figure(5);plotmf(a,'output'
18、;,3);figure(6);plotfis(a);fuzzy fuzzpid;showrule(a);ruleview fuzzpid;附錄2 模糊控制程序PID_b%Fuzzy PID Controlclose all;clear all;a=readfis('fuzzpid');ts=0.001;sys=tf(4.23,1,1.64,8.46);days=c2d(sys,ts,'tustin');num,den=tfdata(days,'v');u_1=0.0;u_2=0.0;u_3=0.0;y_1=0;y_2=0;y_3=0;x=0,0,
19、0'e_1=0;ec_1=0;kp0=0.40;kd0=1.0;ki0=0.0;for k=1:1:3000;time(k)=k*ts;r(k)=1.0;%Using fuzzy inference to tuning PIDk_pid=evalfis(e_1,ec_1,a);kp(k)=kp0+k_pid(1);ki(k)=ki0+k_pid(2);kd(k)=kd0+k_pid(3);u(k)=kp(k)*x(1)+kd(k)*x(2)+ki(k)*x(3);if k=300 %adding disturbance(1.0v at time 0.3s) u(k)=u(k)+1.0;endy(k)=-den(2)*y_1-den(3)*y_2-den(4)*y_3+num(1)*u(k)+num(2)*u_1+num(3)*u_2+num(4)*u_3;e(k)=r(k)-y(k);%Return of PID parameters%u_3=u_2;u_2=u_1;u_1=u(k);y_3=y_2;y_2-
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網頁內容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
- 4. 未經權益所有人同意不得將文件中的內容挪作商業(yè)或盈利用途。
- 5. 人人文庫網僅提供信息存儲空間,僅對用戶上傳內容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內容本身不做任何修改或編輯,并不能對任何下載內容負責。
- 6. 下載文件中如有侵權或不適當內容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 法院辦公室集體活動方案
- 消防展會活動方案
- 活動線上推廣策劃方案
- 活動預熱美食節(jié)活動方案
- 消防救援站掃雪活動方案
- 2024-2025學年青海省海北市九年級化學第一學期期末達標測試試題含解析
- 2024年山東省棗莊薛城區(qū)五校聯(lián)考九年級化學第一學期期末經典模擬試題含解析
- 山東省濰坊市諸城龍源學校2024年七上數(shù)學期末預測試題含解析
- 福建省泉州德化縣聯(lián)考2024年七年級數(shù)學第一學期期末綜合測試試題含解析
- 企業(yè)智慧化進程中零售解決方案全案規(guī)劃
- 臨時工的免責協(xié)議書
- 《急救知識普及》課件
- 《員工的七個習慣》課件
- 防御性駕駛全例培訓
- 《應急救援知識》課件
- 智慧物業(yè)管理方案
- 2024年注塑工作總結與計劃
- 現(xiàn)實表現(xiàn)材料模板
- (新插圖)人教版五年級下冊數(shù)學 6-3-1 分數(shù)加減混合運算 知識點梳理課件
- 人教版初中生物知識點匯總
- 2022年宜賓市敘州區(qū)區(qū)內外考試選調在編在職教師考試真題
評論
0/150
提交評論