GPC.rar_大作业_matlab GPC_matlab例程_数值算法/人工智能下载-pudn.com

<html xmlns="http://www.w3.org/1999/xhtml"> <head> <meta charset="utf-8"> <meta name="generator" content="pdf2htmlEX"> <meta http-equiv="X-UA-Compatible" content="IE=edge,chrome=1"> <link rel="stylesheet" href="https://static.pudn.com/base/css/base.min.css"> <link rel="stylesheet" href="https://static.pudn.com/base/css/fancy.min.css"> <link rel="stylesheet" href="https://static.pudn.com/prod/directory_preview_static/622b630215da9b288b48602d/raw.css"> <script src="https://static.pudn.com/base/js/compatibility.min.js"></script> <script src="https://static.pudn.com/base/js/pdf2htmlEX.min.js"></script> <script> try{ pdf2htmlEX.defaultViewer = new pdf2htmlEX.Viewer({}); }catch(e){} </script> <title></title> </head> <body> <div id="sidebar" style="display: none"> <div id="outline"> </div> </div> <div id="pf1" class="pf w0 h0" data-page-no="1"><div class="pc pc1 w0 h0"><img class="bi x0 y0 w1 h1" alt="" src="https://static.pudn.com/prod/directory_preview_static/622b630215da9b288b48602d/bg1.jpg"><div class="t m0 x1 h2 y1 ff1 fs0 fc0 sc0 ls0 ws0">1</div><div class="c x0 y2 w0 h3"><div class="t m0 x2 h4 y3 ff2 fs1 fc0 sc1 ls0 ws0">被控对象与作<span class="_ _0"></span>业要求</div><div class="t m0 x3 h5 y4 ff2 fs2 fc0 sc0 ls0 ws0">一个二阶非线<span class="_ _1"></span>性对象：</div><div class="t m0 x2 h5 y5 ff2 fs2 fc0 sc0 ls0 ws0">要求：</div><div class="t m0 x2 h5 y6 ff2 fs2 fc0 sc0 ls0 ws0">⒈<span class="_"> </span>建立一个<span class="_ _2"> </span>T-S<span class="_ _2"> </span>模型充分逼近它。<span class="_ _3"> </span>给初值，<span class="_ _4"> </span>用随机数或<span class="_ _1"></span>变化数生成，</div><div class="t m0 x3 h5 y7 ff2 fs2 fc0 sc0 ls0 ws0">获得一批样本<span class="_ _1"></span>。用某种方法建立<span class="_ _1"></span>规则前件。正交最<span class="_ _1"></span>小二乘法辨识参数<span class="_ _1"></span>得<span class="_ _5"> </span>T-S<span class="_ _5"> </span>模型。</div><div class="t m0 x2 h5 y8 ff2 fs2 fc0 sc0 ls0 ws0">⒉<span class="_"> </span>基于建立的<span class="_ _5"> </span>T-<span class="_ _1"></span>S<span class="_ _5"> </span>模型，设计<span class="_ _5"> </span>GPC<span class="_ _6"> </span>来控制原系统。</div><div class="t m0 x3 h4 y9 ff2 fs1 fc0 sc1 ls0 ws0">T-<span class="_ _0"></span>S<span class="_ _5"> </span>模糊<span class="_ _0"></span>辨识</div><div class="t m0 x3 h5 ya ff1 fs2 fc0 sc0 ls0 ws0">T-S<span class="_ _5"> </span><span class="ff2">模糊规则<span class="_ _1"></span>表示如下</span></div><div class="t m0 x4 h5 yb ff2 fs2 fc0 sc0 ls0 ws0">：IF<span class="_ _7"> </span>is<span class="_ _8"> </span>and<span class="_ _9"> </span>is<span class="_ _a"> </span>,......，and<span class="_ _8"> </span>is</div><div class="t m0 x5 h5 yc ff2 fs2 fc0 sc0 ls0 ws0">Then</div></div><div class="c x0 yd w0 h6"><div class="t m0 x2 h5 ye ff2 fs2 fc0 sc0 ls0 ws0">式中<span class="_ _b"> </span>表示第<span class="_ _6"> </span><span class="ff1">i<span class="_ _6"> </span></span>条规则<span class="_ _1"></span>；<span class="_ _b"> </span>表示模糊<span class="_ _1"></span>子集；<span class="_ _c"> </span>为模型的<span class="_ _6"> </span><span class="ff1">j<span class="_ _6"> </span></span>个输入变量<span class="_ _1"></span>，<span class="_ _d"></span>其中<span class="_ _6"> </span><span class="ff1">j=1</span>，<span class="_ _e"></span><span class="ff1">2<span class="ff2">，<span class="_ _e"></span><span class="ff1">…<span class="_ _1"></span>m<span class="_ _1"></span><span class="ff2">，</span></span></span></span></div></div><div class="c x0 yf w0 h6"><div class="t m0 x2 h5 y10 ff1 fs2 fc0 sc0 ls0 ws0">m<span class="_ _6"> </span><span class="ff2">为输入变量的个<span class="_ _1"></span>数；<span class="_ _d"></span><span class="ff1">i<span class="_ _d"></span><span class="ff2">（<span class="ff1">i<span class="_ _5"> </span>=<span class="_ _5"> </span>1</span>，<span class="_ _e"></span><span class="ff1">2<span class="ff2">，<span class="_ _d"></span><span class="ff1">…<span class="_ _1"></span><span class="ff2">，<span class="_ _d"></span><span class="ff1">n<span class="_ _1"></span><span class="ff2">）<span class="_ _d"></span>为规则数；<span class="_ _f"> </span>为第<span class="_ _6"> </span><span class="ff1">i<span class="_ _5"> </span></span>条模<span class="_ _1"></span>糊规则的输出值<span class="_ _1"></span>；<span class="_ _10"> </span>为</span></span></span></span></span></span></span></span></span></div></div><div class="c x0 y2 w0 h3"><div class="t m0 x2 h5 y11 ff2 fs2 fc0 sc0 ls0 ws0">结论参数。</div><div class="t m0 x2 h5 y12 ff2 fs2 fc0 sc1 ls0 ws0">1.<span class="_ _11"> </span>输入空间的划分</div><div class="t m0 x3 h5 y13 ff2 fs2 fc0 sc0 ls0 ws0">系统输入<span class="_ _6"> </span>x∈[-15<span class="_"> </span>15],其中，<span class="_ _1"></span>在这种划分中，每<span class="_ _1"></span>个模糊集的曲线为<span class="_ _1"></span>：</div></div></div><div class="pi" data-data='{"ctm":[1.611811,0.000000,0.000000,1.611811,0.000000,0.000000]}'></div></div> </body> </html>

<div id="pf2" class="pf w0 h0" data-page-no="2"><div class="pc pc2 w0 h0"><img class="bi x0 y0 w1 h1" alt="" src="https://static.pudn.com/prod/directory_preview_static/622b630215da9b288b48602d/bg2.jpg"><div class="t m0 x1 h2 y1 ff1 fs0 fc0 sc0 ls0 ws0">2</div><div class="c x0 y2 w0 h3"><div class="t m0 x3 h5 y14 ff2 fs2 fc0 sc0 ls0 ws0">如图所示：</div></div><div class="t m0 x2 h5 y15 ff1 fs2 fc0 sc1 ls0 ws0">2<span class="_ _0"></span>.<span class="_ _c"> </span><span class="ff2">输入输出模型</span></div><div class="t m0 x3 h5 y16 ff2 fs2 fc0 sc0 ls0 ws0">据所定划分，<span class="_ _1"></span>每个变量划分都<span class="_ _1"></span>为<span class="_ _5"> </span><span class="ff1">4<span class="_ _5"> </span></span>个模糊集，则共<span class="_ _1"></span>有<span class="_ _5"> </span><span class="ff1">64<span class="_ _5"> </span></span>个规则。选<span class="_ _1"></span>取如下<span class="_ _5"> </span><span class="ff1">4<span class="_ _6"> </span></span>个规则起</div><div class="t m0 x2 h5 y17 ff2 fs2 fc0 sc0 ls0 ws0">作用。</div><div class="t m0 x3 h7 y18 ff1 fs2 fc0 sc0 ls0 ws0">R1<span class="_ _5"> </span>:<span class="_ _5"> </span>if<span class="_ _5"> </span>x1<span class="_ _5"> </span>is<span class="_ _5"> </span>E1<span class="_ _5"> </span>and<span class="_ _6"> </span>x2<span class="_ _5"> </span>is<span class="_ _5"> </span>F1<span class="_ _5"> </span>and<span class="_ _5"> </span>x3<span class="_ _6"> </span>is<span class="_ _12"> </span>G1</div><div class="t m0 x5 h7 y19 ff1 fs2 fc0 sc0 ls0 ws0">then<span class="_ _5"> </span>y1<span class="_ _5"> </span>=<span class="_ _5"> </span>a1<span class="_ _1"></span>1<span class="_ _5"> </span>+<span class="_ _5"> </span>a1<span class="_ _1"></span>2*x1<span class="_ _5"> </span>+<span class="_ _5"> </span>a13*x2<span class="_ _5"> </span>+<span class="_ _6"> </span>a<span class="_ _0"></span>14*x3</div><div class="t m0 x3 h7 y1a ff1 fs2 fc0 sc0 ls0 ws0">R2<span class="_ _5"> </span>:<span class="_ _5"> </span>if<span class="_ _5"> </span>x1<span class="_ _5"> </span>is<span class="_ _5"> </span>E2<span class="_ _5"> </span>and<span class="_ _6"> </span>x2<span class="_ _5"> </span>is<span class="_ _5"> </span>F2<span class="_ _5"> </span>and<span class="_ _5"> </span>x3<span class="_ _6"> </span>is<span class="_ _12"> </span>G2</div><div class="t m0 x5 h7 y1b ff1 fs2 fc0 sc0 ls0 ws0">then<span class="_ _5"> </span>y2<span class="_ _5"> </span>=<span class="_ _5"> </span>a21<span class="_ _5"> </span>+<span class="_ _6"> </span>a22*x1<span class="_ _5"> </span>+<span class="_ _5"> </span>a23*x2<span class="_ _5"> </span>+<span class="_ _6"> </span>a<span class="_ _0"></span>24*<span class="_ _1"></span>x3</div><div class="t m0 x3 h7 y1c ff1 fs2 fc0 sc0 ls0 ws0">R3<span class="_ _5"> </span>:<span class="_ _5"> </span>if<span class="_ _5"> </span>x1<span class="_ _5"> </span>is<span class="_ _5"> </span>E3<span class="_ _5"> </span>and<span class="_ _6"> </span>x2<span class="_ _5"> </span>is<span class="_ _5"> </span>F3<span class="_ _5"> </span>and<span class="_ _5"> </span>x3<span class="_ _6"> </span>is<span class="_ _12"> </span>G3</div><div class="t m0 x5 h7 y1d ff1 fs2 fc0 sc0 ls0 ws0">then<span class="_ _5"> </span>y3<span class="_ _5"> </span>=<span class="_ _5"> </span>a31<span class="_ _5"> </span>+<span class="_ _6"> </span>a32*x1<span class="_ _5"> </span>+<span class="_ _5"> </span>a33*x2<span class="_ _5"> </span>+<span class="_ _6"> </span>a<span class="_ _0"></span>34*<span class="_ _1"></span>x3</div><div class="t m0 x3 h7 y1e ff1 fs2 fc0 sc0 ls0 ws0">R4<span class="_ _5"> </span>:<span class="_ _5"> </span>if<span class="_ _5"> </span>x1<span class="_ _5"> </span>is<span class="_ _5"> </span>E4<span class="_ _5"> </span>and<span class="_ _6"> </span>x2<span class="_ _5"> </span>is<span class="_ _5"> </span>F4<span class="_ _5"> </span>and<span class="_ _5"> </span>x3<span class="_ _6"> </span>is<span class="_ _12"> </span>G4</div><div class="t m0 x5 h7 y1f ff1 fs2 fc0 sc0 ls0 ws0">then<span class="_ _5"> </span>y4<span class="_ _5"> </span>=<span class="_ _5"> </span>a41<span class="_ _5"> </span>+<span class="_ _6"> </span>a42*x1<span class="_ _5"> </span>+<span class="_ _5"> </span>a43*x2<span class="_ _5"> </span>+<span class="_ _6"> </span>a<span class="_ _0"></span>44*<span class="_ _1"></span>x3</div><div class="t m0 x3 h5 y20 ff2 fs2 fc0 sc0 ls0 ws0">则模型输出为<span class="_ _1"></span>：</div><div class="t m0 x2 h5 y21 ff2 fs2 fc0 sc0 ls0 ws0">其中，</div><div class="t m0 x2 h5 y22 ff1 fs2 fc0 sc1 ls0 ws0">3<span class="_ _0"></span>.<span class="_ _c"> </span><span class="ff2">参数辨识程序及结果</span></div><div class="t m0 x2 h5 y23 ff2 fs2 fc0 sc0 ls0 ws0">yf(1)=y(1);yf(2)<span class="_ _1"></span>=y(2);</div><div class="t m0 x2 h5 y24 ff2 fs2 fc0 sc0 ls0 ws0">y(1)=0;y(2)=0;u(<span class="_ _1"></span>1)=0;u(2)=0;ee=0<span class="_ _1"></span>.01;</div><div class="t m0 x2 h5 y25 ff2 fs2 fc0 sc0 ls0 ws0">N=8;L=3;p=1;q=0.<span class="_ _1"></span>7;P=8;<span class="_ _11"> </span>na=2;<span class="_ _13"> </span>nb=0;</div><div class="t m0 x2 h5 y26 ff2 fs2 fc1 sc0 ls0 ws0">for<span class="_"> </span><span class="fc0">step=1:N0;</span></div><div class="t m0 x6 h5 y27 ff2 fs2 fc0 sc0 ls0 ws0">A(1)=1;</div><div class="t m0 x6 h5 y28 ff2 fs2 fc0 sc0 ls0 ws0">A(2)=-theta(5<span class="_ _1"></span>:8)*Beta(:,k-2);</div><div class="t m0 x6 h5 y29 ff2 fs2 fc0 sc0 ls0 ws0">A(3)=-theta(9<span class="_ _1"></span>:12)*Beta(:,k-2)<span class="_ _1"></span>;</div><div class="t m0 x6 h5 y2a ff2 fs2 fc0 sc0 ls0 ws0">b(1)=theta(13<span class="_ _1"></span>:16)*Beta(:,k-2)<span class="_ _1"></span>;</div><div class="t m0 x5 h5 y2b ff2 fs2 fc0 sc0 ls0 ws0">e(1)=1;</div><div class="t m0 x5 h5 y2c ff2 fs2 fc1 sc0 ls0 ws0">for<span class="_"> </span><span class="fc0">i=1:na</span></div><div class="t m0 x7 h5 y2d ff2 fs2 fc0 sc0 ls0 ws0">F(1,i)=A(i)-<span class="_ _1"></span>A(i+1);</div></div><div class="pi" data-data='{"ctm":[1.611811,0.000000,0.000000,1.611811,0.000000,0.000000]}'></div></div>