南京理工大学建模辨识与基础实验报告.doc

课程实验报告 姓名：_ 学号： — 花名册编号： 成绩 实验一 一、 实验目标： ?掌握Mat lab编程的基础知识以及一般操作； ?利用Mat lab编程实现M序列。 二、 实验过程： 1、参照例2. 4,编程产生一组M序列和一组连续型伪随机信号，画出信号的 图形。要求它们的幅值和长度可以随意调整，并产生的M序列信号存入文件。 实验程序： %用于生成序列M clear; % Register二[1 111]; % C=[0 10 1]; Register =[111111]; C=[0 0 0 0 1 1]; Ampl i = l ;改变幅度 Bmpli=5; Delta=l; step=0.01; n= 1ength (C); N=2b-1; LenthM=10*N;改变长度 for i=l:LenthM sum=0; for j二1:n surrFsum+Regi ster (j)*C(j); end Regi sterN二mod(sum, 2); M(i)=Register(n); for j二n:-1:2 Register (j) =Register (j-1); end; Regis ter ⑴二 Regi sterN; end for i=l:LenthM Y(l, i)=Delta*i; if (M(i)==0) Y(2, i)=Ampl i ; else Y(2, i)=-Ampli; end end NewY(l, 1)=0. 000000000; NewY(2, 1)=Y(2, 1); CurrentTime=step; j=2； for i=l:LenthM wh i1e CurrentTimeY(1, i)-0. 00001 NewY(1, j)=CurrentTime; NewY(2, j)二Y(2, i); j二j+1； CurrentTime二CurrentTime+step; end end stairs(M); hold on plot(0, Ampli+0. 0001);hold on; plot(0, -Ampli-0. 0001); save MQ Y 实验结果: 生成的连续型伪随机信号: Figure 1 Figure 1 文件 文件(F) 蜩(E) M(V)価入⑴ 1M(T) AQD(D) SD(W) 助(H) % %◎◎电 M?辽□ So100200300400500600 % %◎◎电 M?辽□ S o 100 200 300 400 500 600 700 改变幅值、长度后: Q Figure 1 文件(F)啊(E) M(V) ?A(l)工M(T) A3|(D)窗口 (W) m 助(H) 、 oao^i^i^^o??x*iaiDo Ta 2、计算该伪随机序列的自相关函数。具体地，若采样时间间隔为Is,计算 Rx(-0. 5), Rx(0), R(0. 5), Rx(1), Rx(1.5)。 Rx(-0. 5)=1.5079 Rx(0)=l R(0. 5)=1.5079 Rx(l)=-0. 0159 Rx(1.5)=-0. 0159 实验二 一、 实验冃标： ?掌握利用相关分析法辨识系统的脉冲响应函数。 二、 实验过程： 用计算机模拟仿真如下算例，辨识系统的脉冲响应函数被识系统的传递函数 假设为 251?46以+ 251?46以+ 14?5庁+】 输入信号采用6阶M序列，幅值a=l, △二0。输出采样周期为1 So 实验程序： r=8; for m=l:N temp=0; for k二N+l:N+r*N temp=temp+Y(2, k~m)*z(k); end Ruz (m)二temp/(r*N); end C=-Ruz(N); for i二1:N GM(1, i) = (i-l)*Delta; GM(2, i)=N*(Ruz(i)+C)/((N+l)*Amp1i 2); end plot (GM(1 , :), GM(2 , :), r- ) ;hold on; den二[2];num二[51. 46 14. 5 1]; sys二tf (den, num); impulse(sys); 实验结果： 有噪声情况： Figure 1 — 口 文件(F) 彌8(E) ?I(V)描入(I) IMCD 蠶 U(D)童口 (W)群助(H) MBA A \O ? ?復□目 0 12 O 8pnu-du/v ? 020 10 20 30 40 50 60 Time (seconds) Impulse Response 08 0604 02 O 70 无噪声情况: Figure」 w盘m) 82?(_) H?3 ?S(H) ffi Amplitude b 8 rf -mpu-se