利用相关分析法辨识脉冲响应.docVIP

利用相关分析法辨识脉冲响应 1 实验方案设计 1.1 生成输入数据和噪声 用M序列作为辨识的输入信号用M序列作为辨识的输入信号 3.1 均匀分布随机数生成函数 function sita=U(N) %生成N个[0 1]均匀分布随机数 A=179; x0=11; M=2^15; for k=1:N x2=A*x0; x1=mod(x2,M); v1=x1/(M+1); v(:,k)=v1; x0=x1; end sita=v; end 3.2 正态分布白噪声生成函数 function v=noise(aipi) %生成正态分布N(0,sigma) sigma=1; %标准差 for k=1:length(aipi) ksai=0; for i=1:12 temp=mod(i+k,length(aipi))+1; ksai=ksai+aipi(temp); end v(k)=sigma*(ksai-6); end end 3.3 M序列生成函数 function [Np r M]=createM(n,a) %生成长度为n的M序列,周期为Np，周期数为r x=[1 1 1 1 1 1]; %初始化初态 for i=1:n y=x; x(2:6)=y(1:5);