最小二乘参数辨识matlab程序_matlab_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/62798b24517cd20ea4d44206/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/62798b24517cd20ea4d44206/bg1.jpg"><div class="c x0 y1 w2 h2"><div class="t m0 x1 h3 y2 ff1 fs0 fc0 sc0 ls0 ws0">SI3.1_remindianzu.m</div><div class="t m0 x1 h3 y3 ff1 fs0 fc0 sc0 ls0 ws0">clear all</div><div class="t m0 x1 h3 y4 ff1 fs0 fc0 sc0 ls0 ws0">close all</div><div class="t m0 x1 h3 y5 ff1 fs0 fc0 sc0 ls0 ws0">clc</div><div class="t m0 x1 h4 y6 ff1 fs0 fc0 sc0 ls0 ws0">T=[20.5 26 32.7 40 51 61 73 80 88 95.7];%<span class="ff2">温度</span></div><div class="t m0 x1 h4 y7 ff1 fs0 fc0 sc0 ls0 ws0">R=[765 790 826 850 873 910 942 980 1010 1032];%<span class="ff2">阻值</span></div><div class="t m0 x1 h3 y8 ff1 fs0 fc0 sc0 ls0 ws0">[m,n]=size(T);</div><div class="t m0 x1 h3 y9 ff1 fs0 fc0 sc0 ls0 ws0">figure </div><div class="t m0 x1 h3 ya ff1 fs0 fc0 sc0 ls0 ws0">plot(T<span class="_ _0"></span>,R,'b+')</div><div class="t m0 x1 h3 yb ff1 fs0 fc0 sc0 ls0 ws0">t=0;</div><div class="t m0 x1 h3 yc ff1 fs0 fc0 sc0 ls0 ws0">z=0;</div><div class="t m0 x1 h3 yd ff1 fs0 fc0 sc0 ls0 ws0">tz=0;</div><div class="t m0 x1 h3 ye ff1 fs0 fc0 sc0 ls0 ws0">tt=0;</div><div class="t m0 x1 h3 yf ff1 fs0 fc0 sc0 ls0 ws0">for i=1:n</div><div class="t m0 x1 h3 y10 ff1 fs0 fc0 sc0 ls0 ws0"> t=t+T(i);</div><div class="t m0 x1 h3 y11 ff1 fs0 fc0 sc0 ls0 ws0"> tt=tt+T(i)*T(i);</div><div class="t m0 x1 h3 y12 ff1 fs0 fc0 sc0 ls0 ws0"> z=z+R(i);</div><div class="t m0 x1 h3 y13 ff1 fs0 fc0 sc0 ls0 ws0"> tz=tz+T(i)*R(i);</div><div class="t m0 x1 h3 y14 ff1 fs0 fc0 sc0 ls0 ws0">end</div><div class="t m0 x1 h3 y15 ff1 fs0 fc0 sc0 ls0 ws0">a=(tt*z-t*tz)/(n*tt-t*t);</div><div class="t m0 x1 h3 y16 ff1 fs0 fc0 sc0 ls0 ws0">b=(n*tz-t*z)/(n*tt-t*t);</div><div class="t m0 x1 h3 y17 ff1 fs0 fc0 sc0 ls0 ws0">R1=a+70*b;</div><div class="t m0 x1 h4 y18 ff1 fs0 fc0 sc0 ls0 ws0">%<span class="ff2">最小二乘拟合</span></div><div class="t m0 x1 h3 y19 ff1 fs0 fc0 sc0 ls0 ws0">A=polyfit(T<span class="_ _0"></span>,R,1);</div><div class="t m0 x1 h3 y1a ff1 fs0 fc0 sc0 ls0 ws0">z=polyval(A,T);</div><div class="t m0 x1 h4 y1b ff1 fs0 fc0 sc0 ls0 ws0">%<span class="ff2">画图</span></div><div class="t m0 x1 h3 y1c ff1 fs0 fc0 sc0 ls0 ws0">figure</div><div class="t m0 x1 h3 y1d ff1 fs0 fc0 sc0 ls0 ws0">plot(T<span class="_ _0"></span>,z);</div><div class="t m0 x1 h3 y1e ff1 fs0 fc0 sc0 ls0 ws0">figure</div><div class="t m0 x1 h3 y1f ff1 fs0 fc0 sc0 ls0 ws0">plot(T<span class="_ _0"></span>,R,'b+')</div><div class="t m0 x1 h3 y20 ff1 fs0 fc0 sc0 ls0 ws0">hold on</div><div class="t m0 x1 h3 y21 ff1 fs0 fc0 sc0 ls0 ws0">plot(T<span class="_ _0"></span>,z,'r');</div><div class="t m0 x1 h3 y22 ff1 fs0 fc0 sc0 ls0 ws0">hold off</div><div class="t m0 x1 h3 y23 ff1 fs0 fc0 sc0 ls0 ws0">SI3.2_LS.m</div><div class="t m0 x1 h3 y24 ff1 fs0 fc0 sc0 ls0 ws0">clear all</div><div class="t m0 x1 h3 y25 ff1 fs0 fc0 sc0 ls0 ws0">close all</div><div class="t m0 x1 h3 y26 ff1 fs0 fc0 sc0 ls0 ws0">clc</div><div class="t m0 x1 h3 y27 ff1 fs0 fc0 sc0 ls0 ws0">randn('seed',100)</div><div class="t m0 x1 h4 y28 ff1 fs0 fc0 sc0 ls0 ws0">v=randn(1,16); %<span class="ff2">产生一组<span class="_ _1"> </span></span>16<span class="_ _1"> </span><span class="ff2">个<span class="_ _1"> </span></span>N(<span class="_ _2"></span>0,1)<span class="ff2">的高斯分布的随机噪声</span></div><div class="t m0 x1 h4 y29 ff1 fs0 fc0 sc0 ls0 ws0">% M<span class="_ _1"> </span><span class="ff2">序列产生程序</span></div><div class="t m0 x1 h4 y2a ff1 fs0 fc0 sc0 ls0 ws0">L=15;% M<span class="_ _1"> </span><span class="ff2">序列的周期</span></div><div class="t m0 x1 h4 y2b ff1 fs0 fc0 sc0 ls0 ws0">y1=1;y2=1;y3=1;y4=0;%<span class="ff2">四个移位积存器的输出初始值</span></div></div></div><div class="pi" data-data='{"ctm":[1.611850,0.000000,0.000000,1.611850,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/62798b24517cd20ea4d44206/bg2.jpg"><div class="c x0 y1 w2 h2"><div class="t m0 x1 h3 y2 ff1 fs0 fc0 sc0 ls0 ws0">for i=1:L;</div><div class="t m0 x1 h3 y3 ff1 fs0 fc0 sc0 ls0 ws0"> x1=xor(y3,y4);</div><div class="t m0 x2 h3 y4 ff1 fs0 fc0 sc0 ls0 ws0">x2=y1;</div><div class="t m0 x2 h3 y5 ff1 fs0 fc0 sc0 ls0 ws0">x3=y2;</div><div class="t m0 x1 h3 y2c ff1 fs0 fc0 sc0 ls0 ws0"> x4=y3;</div><div class="t m0 x1 h3 y2d ff1 fs0 fc0 sc0 ls0 ws0"> y(i)=y4;</div><div class="t m0 x1 h3 y8 ff1 fs0 fc0 sc0 ls0 ws0"> if y(i)>0.5,u(i)=-5;</div><div class="t m0 x1 h3 y9 ff1 fs0 fc0 sc0 ls0 ws0"> else u(i)=5;</div><div class="t m0 x1 h3 ya ff1 fs0 fc0 sc0 ls0 ws0"> end</div><div class="t m0 x1 h3 yb ff1 fs0 fc0 sc0 ls0 ws0"> y1=x1;y2=x2;y3=x3;y4=x4;</div><div class="t m0 x1 h3 yc ff1 fs0 fc0 sc0 ls0 ws0">end</div><div class="t m0 x1 h3 yd ff1 fs0 fc0 sc0 ls0 ws0">figure</div><div class="t m0 x1 h3 ye ff1 fs0 fc0 sc0 ls0 ws0">stem(u),grid on</div><div class="t m0 x1 h4 y2e ff1 fs0 fc0 sc0 ls0 ws0">title('<span class="ff2">输入信号<span class="_ _1"> </span></span>M<span class="_ _1"> </span><span class="ff2">序列</span>')</div><div class="t m0 x1 h4 y2f ff1 fs0 fc0 sc0 ls0 ws0">% <span class="ff2">最小二乘辨识程序</span> </div><div class="t m0 x1 h4 y30 ff1 fs0 fc0 sc0 ls0 ws0">z=zeros(1,16); %<span class="ff2">定义输出观测值的长度</span></div><div class="t m0 x1 h3 y12 ff1 fs0 fc0 sc0 ls0 ws0">for k=3:16 </div><div class="t m0 x1 h4 y31 ff1 fs0 fc0 sc0 ls0 ws0"> z(k)=-1.5*z(k-1)-0.7*z(k-2)+u(k-1)+0.5*u(k-2)+1*v(k); %<span class="ff2">观测值</span></div><div class="t m0 x1 h3 y14 ff1 fs0 fc0 sc0 ls0 ws0">end</div><div class="t m0 x1 h3 y15 ff1 fs0 fc0 sc0 ls0 ws0">figure(2)</div><div class="t m0 x1 h3 y16 ff1 fs0 fc0 sc0 ls0 ws0">plot([1:16],z)</div><div class="t m0 x1 h4 y32 ff1 fs0 fc0 sc0 ls0 ws0">title('<span class="ff2">输出观测值</span>')</div><div class="t m0 x1 h3 y33 ff1 fs0 fc0 sc0 ls0 ws0">figure(3)</div><div class="t m0 x1 h3 y19 ff1 fs0 fc0 sc0 ls0 ws0">stem(z),grid on</div><div class="t m0 x1 h4 y34 ff1 fs0 fc0 sc0 ls0 ws0">title('<span class="ff2">输出观测值<span class="_ _1"> </span></span>z<span class="_ _1"> </span><span class="ff2">的经线图形</span>')</div><div class="t m0 x1 h4 y1b ff1 fs0 fc0 sc0 ls0 ws0">%<span class="ff2">给样本系数矩阵</span></div><div class="t m0 x1 h3 y1c ff1 fs0 fc0 sc0 ls0 ws0">H=[-z(2) -z(1) u(2)<span class="_ _3"></span> u(1);-z(3) -z(2) u(3) u<span class="_ _3"></span>(2);-z(4) -z(3) u(4) u(3);-<span class="_ _3"></span>z(5) -z(4) u(5) u(4);-z(6) <span class="_ _3"></span>-z(5) u(6) u(5);-z(7)</div><div class="t m0 x1 h3 y1d ff1 fs0 fc0 sc0 ls0 ws0">-z(6) <span class="_ _3"></span>u(7)<span class="_ _3"></span> u(6);-z(8)<span class="_ _3"></span> <span class="_ _3"></span>-z(7) <span class="_ _3"></span>u(8)<span class="_ _3"></span> u(7);-z(9)<span class="_ _3"></span> -z<span class="_ _3"></span>(8) <span class="_ _3"></span>u(9)<span class="_ _3"></span> u(8);-z(10)<span class="_ _3"></span> <span class="_ _3"></span>-z(9) <span class="_ _3"></span>u(10)<span class="_ _3"></span> u(9);-z(1<span class="_ _2"></span>1) <span class="_ _3"></span>-z(10)<span class="_ _3"></span> u<span class="_ _3"></span>(1<span class="_ _2"></span>1) u(10);-z(12)</div><div class="t m0 x1 h3 y1e ff1 fs0 fc0 sc0 ls0 ws0">-z(1<span class="_ _2"></span>1) u(12) u(1<span class="_ _0"></span>1<span class="_ _3"></span>);-z(13) -z(12) u(13) u(12);-z(14) -z(13) u(14) u(13);-z(15) -z(14) u(15) u(<span class="_ _2"></span>14)]<span class="_ _3"></span>;</div><div class="t m0 x1 h4 y35 ff1 fs0 fc0 sc0 ls0 ws0">%<span class="ff2">给出样本观测矩阵</span></div><div class="t m0 x1 h3 y20 ff1 fs0 fc0 sc0 ls0 ws0">Z=[z(3);z(4);z(5);z(6);z(7);z(8);z(9);z(10);z(1<span class="_ _2"></span>1);z(12);z(13);z(14);z(15);z(16)] </div><div class="t m0 x1 h4 y36 ff1 fs0 fc0 sc0 ls0 ws0">%<span class="ff2">计算参数</span></div><div class="t m0 x1 h3 y22 ff1 fs0 fc0 sc0 ls0 ws0">c=inv(H'*H)*H'*Z;</div><div class="t m0 x1 h4 y37 ff1 fs0 fc0 sc0 ls0 ws0">%<span class="ff2">分离参数</span></div><div class="t m0 x1 h3 y38 ff1 fs0 fc0 sc0 ls0 ws0">a1=c(1), a2=c(2), b1=c(3), b2=c(4)</div><div class="t m0 x1 h3 y24 ff1 fs0 fc0 sc0 ls0 ws0">SI3.1_RLS.m</div><div class="t m0 x1 h3 y25 ff1 fs0 fc0 sc0 ls0 ws0">clear all</div><div class="t m0 x1 h3 y26 ff1 fs0 fc0 sc0 ls0 ws0">close all</div><div class="t m0 x1 h3 y27 ff1 fs0 fc0 sc0 ls0 ws0">clc</div><div class="t m0 x1 h4 y28 ff1 fs0 fc0 sc0 ls0 ws0">%<span class="ff2">产生<span class="_ _1"> </span></span>N(0,1)<span class="ff2">正态分布的随机噪声</span></div><div class="t m0 x1 h3 y39 ff1 fs0 fc0 sc0 ls0 ws0">randn('seed',100);</div><div class="t m0 x1 h3 y3a ff1 fs0 fc0 sc0 ls0 ws0">v=randn(1,60); </div></div></div><div class="pi" data-data='{"ctm":[1.611850,0.000000,0.000000,1.611850,0.000000,0.000000]}'></div></div>