系统辨识及其matlab仿真(一些噪声与辨识算法))

《系统辨识及其matlab仿真(一些噪声与辨识算法))》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、beecq工程师精华0积分180帖子86水位180技术分0系统辨识及其matlab仿真（一些噪声和辨识算法)）参考文献：[1]系统辨识及其matlab仿真。侯媛彬，汪梅，王立琦。科学出版社2004。[2]过程辨识。方崇智。清华大学出版社，1987。[3]系统辨识讲义。伍清河，北京理工大学，2005。（见http://www.baisi.net/post/view?bid=51&id=204235&sty=1&tpg=2&age=0主要内容：【1】随机序列产生程序【2】白噪声产生程序【3】M序列产生程序【4】二阶系统一次性完成最小二乘辨识程序【

2、5】实际压力系统的最小二乘辨识程序【6】递推的最小二乘辨识程序【7】增广的最小二乘辨识程序【8】梯度校正的最小二乘辨识程序【9】递推的极大似然辨识程序【10】Bayes辨识程序【11】改进的神经网络MBP算法对噪声系统辨识程序【12】多维非线性函数辨识程序的Matlab程序【13】模糊神经网络解耦Matlab程序【14】F-检验法部分程序……待更新注：帖子上出现为“：）”！2007-5-3013:04#1  popoi工程师精华0积分189帖子91水位189技术分0【1】随机序列产生程序A=6;x0=1;M=255;fork=1:100x2=

3、A*x0;x1=mod(x2,M);v1=x1/256;v(:,k)=v1;x0=x1;v0=v1;endv2=vk1=k;%grapherk=1:k1;plot(k,v,k,v,'r');xlabel('k'),ylabel('v');title('(0,1)均匀分布的随机序列')2007-5-3014:09#2  xlc73工程师精华0积分173帖子80水位161技术分0【2】白噪声产生程序A=6;x0=1;M=255;f=2;N=100;fork=1:Nx2=A*x0;x1=mod(x2,M);v1=x1/256;v(:,k)=(v1-

4、0.5)*f;x0=x1;v0=v1;endv2=vk1=k;%grapherk=1:k1;plot(k,v,k,v,'r');xlabel('k'),ylabel('v');title('(-1,+1)均匀分布的白噪声')2007-5-3016:05#3  fbqclq工程师精华0积分195帖子98水位195技术分0【3】M序列产生程序X1=1;X2=0;X3=1;X4=0;%移位寄存器输入Xi初T态（0101），Yi为移位寄存器各级输出m=60;%置M序列总长度fori=1:m%1#Y4=X4;Y3=X3;Y2=X2;Y1=X1;X4=Y

5、3;X3=Y2;X2=Y1;X1=xor(Y3,Y4);%异或运算ifY4==0U(i)=-1;elseU(i)=Y4;endendM=U%绘图i1=ik=1:1:i1;plot(k,U,k,U,'rx')xlabel('k')ylabel('M序列')title('移位寄存器产生的M序列')2007-5-3017:28#4  lywsd2000工程师精华0积分210帖子103水位210技术分0【4】二阶系统一次性完成最小二乘辨识程序%FLch3LSeg1u=[-1,1,-1,1,1,1,1,-1,-1,-1,1,-1,-1,1,1];%系统

6、辨识的输入信号为一个周期的M序列z=zeros(1,16);%定义输出观测值的长度fork=3:16z(k)=1.5*z(k-1)-0.7*z(k-2)+u(k-1)+0.5*u(k-2);%用理想输出值作为观测值endsubplot(3,1,1)%画三行一列图形窗口中的第一个图形stem(u)%画出输入信号u的经线图形subplot(3,1,2)%画三行一列图形窗口中的第二个图形i=1:1:16;%横坐标范围是1到16，步长为1plot(i,z)%图形的横坐标是采样时刻i,纵坐标是输出观测值z,图形格式为连续曲线subplot(3,1,3)

7、%画三行一列图形窗口中的第三个图形stem(z),gridon%画出输出观测值z的经线图形，并显示坐标网格u,z%显示输入信号和输出观测信号%L=14%数据长度HL=[-z(2)-z(1)u(2)u(1);-z(3)-z(2)u(3)u(2);-z(4)-z(3)u(4)u(3);-z(5)-z(4)u(5)u(4);-z(6)-z(5)u(6)u(5);-z(7)-z(6)u(7)u(6);-z(8)-z(7)u(8)u(7);-z(9)-z(8)u(9)u(8);-z(10)-z(9)u(10)u(9);-z(11)-z(10)u(11)

8、u(10);-z(12)-z(11)u(12)u(11);-z(13)-z(12)u(13)u(12);-z(14)-z(13)u(14)u(13);-z(15)-