# 最小二乘辨识.rar

• zui556
了解作者
• matlab
开发工具
• 83KB
文件大小
• rar
文件格式
• 0
收藏次数
• 10 积分
下载积分
• 4
下载次数
• 2020-03-29 15:40
上传日期

• 递推最小二乘程序
• whitenoise.txt
9KB
• WLS.m
2.4KB
• WLS.asv
2.2KB
• Data.txt
1.6KB
• work3.asv
660B
• 白噪声序列.txt
9KB
• RFF.asv
2.4KB
• M序列.txt
6.9KB
• 最小二乘辨识.doc
66.5KB
• 产生白噪声和M序列.doc
81KB
• Matlab中文论坛--助努力的人完成毕业设计.url
183B
• 使用帮助：新手必看.htm
3.5KB

<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/6268b4f84f8811599e1782ad/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/6268b4f84f8811599e1782ad/bg1.jpg"><div class="c x0 y1 w2 h2"><div class="t m0 x1 h3 y2 ff1 fs0 fc0 sc0 ls0 ws0">&#29677;&#32423;&#65306;&#20449;&#30740;<span class="ff2 sc1">0701 </span>&#22995;&#21517;&#65306;&#26446;&#28023;&#38686;<span class="ff2 sc1"> </span>&#23398;&#21495;&#65306;</div><div class="t m0 x1 h4 y3 ff2 fs0 fc0 sc1 ls0 ws0">2007000669 </div><div class="t m0 x2 h3 y4 ff1 fs0 fc0 sc0 ls0 ws0">&#20316;&#19994;<span class="ff2 sc1">2</span>&#65306;&#20135;&#29983;&#30333;&#22122;&#22768;&#21644;<span class="ff2 sc1">M</span>&#24207;&#21015;</div><div class="t m0 x1 h3 y5 ff2 fs0 fc0 sc1 ls0 ws0">1.<span class="_ _0"> </span><span class="ff1 sc0">&#31243;&#24207;&#35828;&#26126;&#65306;</span></div><div class="t m0 x3 h3 y6 ff1 fs0 fc0 sc1 ls0 ws0">&#26412;&#31243;&#24207;&#29992;<span class="ff3">MATLAB</span>&#35821;&#35328;&#32534;&#20889;<span class="ff3">.</span>&#31243;&#24207;&#36816;&#34892;&#23436;&#25104;&#21487;&#20197;&#26174;&#31034;&#30333;&#22122;&#22768;&#22270;&#24418;&#21644;&#30333;&#22122;&#22768;&#20449;</div><div class="t m0 x3 h3 y7 ff1 fs0 fc0 sc1 ls0 ws0">&#21495;&#25968;&#25454;&#65292;&#25968;&#25454;&#20445;&#23384;&#22312;&#8217;&#30333;&#22122;&#22768;&#24207;&#21015;<span class="ff3">.txt&#8217;</span>&#25991;&#20214;&#20013;&#12290;</div><div class="t m0 x1 h5 y8 ff3 fs1 fc0 sc1 ls0 ws0"> </div><div class="t m0 x1 h3 y9 ff4 fs0 fc0 sc1 ls0 ws0">2. <span class="ff1 sc0">&#31243;&#24207;&#20195;&#30721;&#65306;</span></div><div class="t m0 x1 h6 ya ff5 fs0 fc1 sc1 ls0 ws0">global<span class="fc0"> M f d n P</span></div><div class="t m0 x1 h6 yb ff5 fs0 fc0 sc1 ls0 ws0">e=zeros(1,1001);</div><div class="t m0 x1 h6 yc ff5 fs0 fc0 sc1 ls0 ws0">u=zeros(1,1000);</div><div class="t m0 x1 h6 yd ff5 fs0 fc0 sc1 ls0 ws0">c=zeros(1,1000);</div><div class="t m0 x1 h6 ye ff5 fs0 fc0 sc1 ls0 ws0">A=179;</div><div class="t m0 x1 h6 yf ff5 fs0 fc0 sc1 ls0 ws0">P=3.1415926;</div><div class="t m0 x1 h6 y10 ff5 fs0 fc0 sc1 ls0 ws0">M=2^35;</div><div class="t m0 x1 h3 y11 ff5 fs0 fc0 sc1 ls0 ws0">x=11; <span class="fc2">%<span class="ff1">&#21021;&#20540;</span>x<span class="ff1">&#21462;&#27491;&#22855;&#25968;</span></span></div><div class="t m0 x1 h3 y12 ff5 fs0 fc0 sc1 ls0 ws0">e(1)=x/M;<span class="fc2">%<span class="ff1">&#20266;&#38543;&#26426;&#24207;&#21015;&#26159;</span>e(i)=x(i)/M<span class="ff1">&#65292;&#20854;&#20013;</span>x(i)=A*x(i)+k*M<span class="_ _1"></span><span class="ff1">&#12290;&#20004;&#36793;&#21516;&#38500;&#19968;&#20010;</span>M<span class="ff1">&#12290;</span></span></div><div class="t m0 x4 h3 y13 ff1 fs0 fc2 sc1 ls0 ws0">&#26681;&#25454;&#21516;&#20313;&#24335;&#23450;&#29702;&#24471;&#21040;<span class="ff5">e(i)=A*e(i)+k*e(i)</span></div><div class="t m0 x1 h3 y14 ff5 fs0 fc0 sc1 ls0 ws0">f=0; <span class="fc2">%<span class="ff1">&#30333;&#22122;&#22768;&#30340;&#22343;&#20540;</span></span></div><div class="t m0 x1 h3 y15 ff5 fs0 fc0 sc1 ls0 ws0">d=0.1; <span class="fc2">%<span class="ff1">&#30333;&#22122;&#22768;&#30340;&#26041;&#24046;</span></span></div><div class="t m0 x1 h3 y16 ff5 fs0 fc0 sc1 ls0 ws0">n=600; <span class="fc2">%<span class="ff1">&#30333;&#22122;&#22768;&#30340;&#24207;&#21015;&#38271;&#24230;&#65306;</span>( <span class="ff1">&#22312;</span>1<span class="ff1">&#21040;</span>1000<span class="ff1">&#20043;&#38388;</span>)</span></div><div class="t m0 x1 h6 y17 ff5 fs0 fc2 sc1 ls0 ws0"> </div><div class="t m0 x1 h3 y18 ff5 fs0 fc0 sc1 ls0 ws0">fid = fopen(<span class="fc3">'<span class="ff1">&#30333;&#22122;&#22768;&#24207;&#21015;</span>.txt'</span>,<span class="fc3">'wt'</span>); <span class="fc2">%<span class="ff1">&#20135;&#29983;&#19968;&#20010;&#31354;&#30333;&#25991;&#20214;&#65292;&#29992;&#20110;&#23384;&#25918;&#25968;&#25454;</span>.</span></div><div class="t m0 x1 h6 y19 ff5 fs0 fc0 sc1 ls0 ws0"> </div><div class="t m0 x1 h6 y1a ff5 fs0 fc1 sc1 ls0 ws0">for<span class="fc0"> i=2:n+1</span></div><div class="t m0 x1 h3 y1b ff5 fs0 fc0 sc1 ls0 ws0"> e(i)=A*e(i-1)-fix(A*e(i-1)); <span class="fc2">%round<span class="ff1">&#26159;&#22235;&#33293;&#20116;&#20837;&#21040;&#26368;&#36817;&#30340;&#25972;&#25968;</span>,</span></div><div class="t m0 x1 h6 y1c ff5 fs0 fc1 sc1 ls0 ws0">end<span class="fc0"> </span></div><div class="t m0 x1 h6 y1d ff5 fs0 fc0 sc1 ls0 ws0"> </div><div class="t m0 x1 h3 y1e ff5 fs0 fc1 sc1 ls0 ws0">for<span class="fc0"> a=1:n <span class="fc2">%<span class="ff1">&#22343;&#21248;&#20998;&#24067;&#30340;&#30333;&#22122;&#22768;&#65292;&#29992;&#21464;&#25442;&#25277;&#26679;&#27861;&#36716;&#25442;&#20026;&#27491;&#24577;&#20998;&#24067;&#30340;&#30333;&#22122;&#22768;</span></span></span></div><div class="t m0 x1 h6 y1f ff5 fs0 fc0 sc1 ls0 ws0"> u(a)=sqrt((-2)*(log(e(a))))*cos(2*P*e(a+1)); </div><div class="t m0 x1 h6 y20 ff5 fs0 fc0 sc1 ls0 ws0"> c(a)=f+d*u(a);</div><div class="t m0 x1 h3 y21 ff5 fs0 fc0 sc1 ls0 ws0"> fprintf(fid,<span class="fc3">' %d'</span>,c(a)); <span class="fc2">%<span class="ff1">&#20889;&#20837;&#30333;&#22122;&#22768;</span></span></div><div class="t m0 x1 h6 y22 ff5 fs0 fc1 sc1 ls0 ws0">end</div><div class="t m0 x1 h6 y23 ff5 fs0 fc1 sc1 ls0 ws0"> </div><div class="t m0 x1 h6 y24 ff5 fs0 fc0 sc1 ls0 ws0">fclose(fid);</div><div class="t m0 x1 h6 y25 ff5 fs0 fc0 sc1 ls0 ws0">a=1:n</div><div class="t m0 x1 h3 y26 ff5 fs0 fc0 sc1 ls0 ws0">plot(a,c(a)); <span class="fc2">%<span class="ff1">&#36755;&#20986;&#30333;&#22122;&#22768;&#22270;&#24418;</span></span></div><div class="t m0 x1 h6 y27 ff5 fs0 fc0 sc1 ls0 ws0">xlabel(<span class="fc3">'a'</span>) </div><div class="t m0 x1 h3 y28 ff5 fs0 fc0 sc1 ls0 ws0">ylabel(<span class="fc3">'<span class="ff1">&#30333;&#22122;&#22768;&#20449;&#21495;</span>'</span>) </div><div class="t m0 x1 h3 y29 ff5 fs0 fc0 sc1 ls0 ws0">title(<span class="fc3">'<span class="ff1">&#27491;&#24577;&#20998;&#24067;&#30333;&#22122;&#22768;</span>'</span>)</div><div class="t m0 x1 h7 y2a ff5 fs0 fc0 sc1 ls0 ws0"> <span class="ff3"> </span></div><div class="t m0 x1 h3 y2b ff4 fs0 fc0 sc1 ls0 ws0">3. <span class="ff1 sc0">&#31243;&#24207;&#36816;&#34892;&#32467;&#26524;&#65306;</span></div><div class="t m0 x1 h8 y2c ff3 fs1 fc0 sc1 ls0 ws0">3.1 <span class="ff1">&#20135;&#29983;&#30333;&#22122;&#22768;&#22270;&#24418;&#65306;</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>

相关推荐