最小二乘辨识.rar_最小二乘辨识_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/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">班级：信研<span class="ff2 sc1">0701 </span>姓名：李海霞<span class="ff2 sc1"> </span>学号：</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">作业<span class="ff2 sc1">2</span>：产生白噪声和<span class="ff2 sc1">M</span>序列</div><div class="t m0 x1 h3 y5 ff2 fs0 fc0 sc1 ls0 ws0">1.<span class="_ _0"> </span><span class="ff1 sc0">程序说明：</span></div><div class="t m0 x3 h3 y6 ff1 fs0 fc0 sc1 ls0 ws0">本程序用<span class="ff3">MATLAB</span>语言编写<span class="ff3">.</span>程序运行完成可以显示白噪声图形和白噪声信</div><div class="t m0 x3 h3 y7 ff1 fs0 fc0 sc1 ls0 ws0">号数据，数据保存在’白噪声序列<span class="ff3">.txt’</span>文件中。</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">程序代码：</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">初值</span>x<span class="ff1">取正奇数</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">伪随机序列是</span>e(i)=x(i)/M<span class="ff1">，其中</span>x(i)=A*x(i)+k*M<span class="_ _1"></span><span class="ff1">。两边同除一个</span>M<span class="ff1">。</span></span></div><div class="t m0 x4 h3 y13 ff1 fs0 fc2 sc1 ls0 ws0">根据同余式定理得到<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">白噪声的均值</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">白噪声的方差</span></span></div><div class="t m0 x1 h3 y16 ff5 fs0 fc0 sc1 ls0 ws0">n=600; <span class="fc2">%<span class="ff1">白噪声的序列长度：</span>( <span class="ff1">在</span>1<span class="ff1">到</span>1000<span class="ff1">之间</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">白噪声序列</span>.txt'</span>,<span class="fc3">'wt'</span>); <span class="fc2">%<span class="ff1">产生一个空白文件，用于存放数据</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">是四舍五入到最近的整数</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">均匀分布的白噪声，用变换抽样法转换为正态分布的白噪声</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">写入白噪声</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">输出白噪声图形</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">白噪声信号</span>'</span>) </div><div class="t m0 x1 h3 y29 ff5 fs0 fc0 sc1 ls0 ws0">title(<span class="fc3">'<span class="ff1">正态分布白噪声</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">程序运行结果：</span></div><div class="t m0 x1 h8 y2c ff3 fs1 fc0 sc1 ls0 ws0">3.1 <span class="ff1">产生白噪声图形：</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/6268b4f84f8811599e1782ad/bg2.jpg"><div class="c x0 y1 w2 h2"><div class="t m0 x1 h8 y2d ff3 fs1 fc0 sc1 ls0 ws0">3.2 <span class="ff1">产生白噪声信号数据</span></div><div class="t m0 x1 h7 y2e ff3 fs0 fc0 sc1 ls0 ws0">6.612451e-001 5.774923e-001 4.792902e-001 -1.683940e-001 </div><div class="t m0 x1 h7 y2f ff3 fs0 fc0 sc1 ls0 ws0">7.288249e-002 1.575171e-002 -1.460384e-002 6.120653e-002 -</div><div class="t m0 x1 h7 y30 ff3 fs0 fc0 sc1 ls0 ws0">1.113070e-001 5.508005e-002 -1.950977e-002 -1.364850e-001 </div><div class="t m0 x1 h7 y31 ff3 fs0 fc0 sc1 ls0 ws0">8.878680e-002 3.012075e-002 -2.957884e-002 5.373468e-002 -</div><div class="t m0 x1 h7 y32 ff3 fs0 fc0 sc1 ls0 ws0">3.597785e-002 5.740693e-002 4.419843e-003 1.119452e-001 </div><div class="t m0 x1 h7 y33 ff3 fs0 fc0 sc1 ls0 ws0">1.403342e-001 -1.762979e-002 -1.302758e-001 -9.629335e-002 </div><div class="t m0 x1 h7 y34 ff3 fs0 fc0 sc1 ls0 ws0">-5.101620e-002 -7.917252e-003 -3.297326e-002 -1.521888e-001 </div><div class="t m0 x1 h7 y35 ff3 fs0 fc0 sc1 ls0 ws0">8.243325e-002 -4.356131e-002 -1.308579e-001 8.989907e-002 -</div><div class="t m0 x1 h7 y36 ff3 fs0 fc0 sc1 ls0 ws0">1.867038e-002 3.254982e-003 4.693020e-002 -5.067337e-002 -</div><div class="t m0 x1 h7 y37 ff3 fs0 fc0 sc1 ls0 ws0">1.292082e-001 1.155992e-001 4.529656e-003 5.654237e-002 -</div><div class="t m0 x1 h7 y38 ff3 fs0 fc0 sc1 ls0 ws0">1.556838e-001 -5.583714e-002 -8.669206e-002 4.079805e-002 -</div><div class="t m0 x1 h7 y39 ff3 fs0 fc0 sc1 ls0 ws0">1.673424e-001 7.605478e-002 -5.124867e-003 -7.658467e-002 -</div><div class="t m0 x1 h7 y3a ff3 fs0 fc0 sc1 ls0 ws0">3.278001e-002 2.651477e-002 -1.684770e-001 -5.747737e-002 -</div><div class="t m0 x1 h7 y3b ff3 fs0 fc0 sc1 ls0 ws0">2.768529e-002 -1.459871e-001 -1.010220e-001 1.081117e-001 </div><div class="t m0 x1 h7 y3c ff3 fs0 fc0 sc1 ls0 ws0">3.651037e-002 -7.220661e-003 -7.997269e-002 -9.720626e-002 </div><div class="t m0 x1 h7 y3d ff3 fs0 fc0 sc1 ls0 ws0">-1.026337e-001 -2.993016e-002 8.088102e-002 1.235801e-001 -</div><div class="t m0 x1 h7 y3e ff3 fs0 fc0 sc1 ls0 ws0">1.443281e-001 -3.308421e-002 -1.372234e-001 1.280606e-001 </div><div class="t m0 x1 h7 y3f ff3 fs0 fc0 sc1 ls0 ws0">2.565054e-002 2.674790e-001 -9.064962e-002 1.535033e-001 </div></div></div><div class="pi" data-data='{"ctm":[1.611850,0.000000,0.000000,1.611850,0.000000,0.000000]}'></div></div>