PID神經(jīng)元網(wǎng)絡(luò)解耦控制算法-多變量系統(tǒng)控制.doc

《PID神經(jīng)元網(wǎng)絡(luò)解耦控制算法-多變量系統(tǒng)控制.doc》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在教育資源-天天文庫。

1、%%清空環(huán)境變量clcclear%%網(wǎng)絡(luò)結(jié)構(gòu)初始化rate1=0.006;rate2=0.001;%學(xué)習(xí)率k=0.3;K=3;y_1=zeros(3,1);y_2=y_1;y_3=y_2;%輸出值u_1=zeros(3,1);u_2=u_1;u_3=u_2;%控制率h1i=zeros(3,1);h1i_1=h1i;%第一個控制量h2i=zeros(3,1);h2i_1=h2i;%第二控制量h3i=zeros(3,1);h3i_1=h3i;%第三個空置量x1i=zeros(3,1);x2i=x1i;x3i=

2、x2i;x1i_1=x1i;x2i_1=x2i;x3i_1=x3i;%隱含層輸出%權(quán)值初始化k0=0.03;%第一層權(quán)值w11=k0*rand(3,2);w12=k0*rand(3,2);w13=k0*rand(3,2);%第二層權(quán)值w21=k0*rand(1,9);w22=k0*rand(1,9);w23=k0*rand(1,9);%值限定ynmax=1;ynmin=-1;%系統(tǒng)輸出值限定xpmax=1;xpmin=-1;%P節(jié)點輸出限定qimax=1;qimin=-1;%I節(jié)點輸出限定qdmax=1;

3、qdmin=-1;%D節(jié)點輸出限定uhmax=1;uhmin=-1;%輸出結(jié)果限定%%網(wǎng)絡(luò)迭代優(yōu)化fork=1:1:200%%控制量輸出計算%--------------------------------網(wǎng)絡(luò)前向計算--------------------------%系統(tǒng)輸出y1(k)=(0.4*y_1(1)+u_1(1)/(1+u_1(1)^2)+0.2*u_1(1)^3+0.5*u_1(2))+0.3*y_1(2);y2(k)=(0.2*y_1(2)+u_1(2)/(1+u_1(2)^2)+0.4

4、*u_1(2)^3+0.2*u_1(1))+0.3*y_1(3);y3(k)=(0.3*y_1(3)+u_1(3)/(1+u_1(3)^2)+0.4*u_1(3)^3+0.4*u_1(2))+0.3*y_1(1);r1(k)=0.7;r2(k)=0.4;r3(k)=0.6;%控制目標(biāo)%系統(tǒng)輸出限制yn=[y1(k),y2(k),y3(k)];yn(find(yn>ynmax))=ynmax;yn(find(yn

5、);yn(2)];x3o=[r3(k);yn(3)];%隱含層x1i=w11*x1o;x2i=w12*x2o;x3i=w13*x3o;%比例神經(jīng)元P計算xp=[x1i(1),x2i(1),x3i(1)];xp(find(xp>xpmax))=xpmax;xp(find(xp

6、,h2i(2),h3i(2)];qi=qi_1+xi;qi(find(qi>qimax))=qimax;qi(find(qiqdmax))=qdmax;qd(find(qd

7、h1i(3)=qd(1);h2i(3)=qd(2);h3i(3)=qd(3);%輸出層計算wo=[w21;w22;w23];qo=[h1i',h2i',h3i'];qo=qo';uh=wo*qo;uh(find(uh>uhmax))=uhmax;uh(find(uh

8、r=[r1(k)-y1(k);r2(k)-y2(k);r3(k)-y3(k)];error1(k)=error(1);error2(k)=error(2);error3(k)=error(3);J(k)=0.5*(error(1)^2+error(2)^2+error(3)^2);%調(diào)整大小ypc=[y1(k)-y_1(1);y2(k)-y_1(2);y3(k)-y_1(3)];uhc=[u_1(1)-u_2(1);u_1(2