《系统辨识及其matlab仿真一些噪声和辨识算法》由会员分享，可在线阅读，更多相关《系统辨识及其matlab仿真一些噪声和辨识算法（9页珍藏版）》请在金锄头文库上搜索。
1、1】随机序列产生程序2】白噪声产生程序【3】M序列产生程序4】二阶系统一次性完成最小二乘辨识程序5】实际压力系统的最小二乘辨识程序6】递推的最小二乘辨识程序7】增广的最小二乘辨识程序8】梯度校正的最小二乘辨识程序9】递推的极大似然辨识程序10】Bayes辨识程序11】改进的神经网络MBP算法对噪声系统辨识程序12】多维非线性函数辨识程序的Matlab程序【13】模糊神经网络解耦Matlab程序【14】F-检验法部分程序【1】随机序列产生程序A=6;x0=1;M=255;fork=1:100x2=A*x0;x1=mod(x2,M);v1=x1/256;v(:,k)=v1;x0=x1;v0=v1;
2、endv2=vk1=k;%grapherk=1:k1;plot(k,v,k,v,r);xlabel(k),ylabel(v);title(O,l)均匀分布的随机序列)【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-0.5)*f;x0=x1;v0=v1;endv2=vk1=k;%grapherk=1:k1;plot(k,v,k,v,r);xlabel(k),ylabel(v);title(-l,+l)均匀分布的白噪声)【3】M序列产生程序Xl=l;X2=0;X3=l;X4=0
3、;%移位寄存器输入Xi初T态(0101),Yi为移位寄存器各级输出m=60;%置M序列总长度fori=1:m%1#Y4=X4;Y3=X3;Y2=X2;Y1=X1;X4=Y3;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序列)【4】二阶系统一次性完成最小二乘辨识程序%FLch3LSeg1u=-1,1,-1,1,1,1,1,-1,-1,-1,1,-1,-1,1,1;%系统辨识的输
4、入信号为一个周期的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)%画三行一列图形窗口中的第三个图形stem(z),grido
5、n%画出输出观测值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)u(10);-z(12)-z(11)u(12)u(11);-z(13)-z(12)u(13)u(12);-z(14
6、)-z(13)u(14)u(13);-z(15)-z(14)u(15)u(14)%给样本矩阵HL赋值ZL=z(3);z(4);z(5);z(6);z(7);z(8);z(9);z(10);z(11);z(12);z(13);z(14);z(15);z(16)%给样本矩阵zL赋值%calculatingparameters%计算参数c1=HL*HL;c2=inv(c1);c3=HL*ZL;c=c2*c3%计算并显示%DISPLAYPARAMETERSa1=c(1),a2=c(2),bl=c(3),b2=c(4)%从中分离出并显示al、a2、bl、b2%End【5】实际压力系统的最小二乘辨识程序%
7、FLch3LSeg2clear%工作间清零V=54.3,61.8,72.4,88.7,118.6,194.0%赋初值V,并显示P=61.2,49.5,37.6,28.4,19.2,10.1%赋初值P,并显示%logP=-alpha*logV+logbeita=-logV,1alpha,log(beita)=HL*sita%注释P、V之间的关系fori=1:6;%循环变量的取值为从1到6Z(i)=log(P(i);%赋系统的输出采样值end%循环结束ZL=Z%给zL赋值HL=-log(V(1),1;-log(V(2),1;-log(V(3),1;-log(V(4),1;-log(V(5),1;-
8、log(V(6),1%给HL赋值%calculatingparameters%计算参数c1=HL*HL;c2=inv(c1);c3=HL*ZL;c4=c2*c3%计算%SeparationofParameters%分离变量alpha=c4(1)%为c4的第1个元素beita=exp(c4(2)%为以自然数为底的c4的第2个元素的指数【6】递推的最小二乘辨识程序%FLch3RLSeg3clear%清理工作间变量L=15;%M序列的周期y1=1;y2=1;y3=1;y4=0;%四个移位积存器的输出初始值fori=1;%开始循环，长度为Lx1=xor(y3,y4);%第一个移位积存器的输入是第3个与
9、第4个移位积存器的输出的或x2=y1;%第二个移位积存器的输入是第3个移位积存器的输出x3=y2;%第三个移位积存器的输入是第2个移位积存器的输出x4=y3;%第四个移位积存器的输入是第3个移位积存器的输出y(i)=y4;%取出第四个移位积存器幅值为0和1的输出信号，ify(i)0.5,u(i)=-0.03;%如果M序列的值为1时，辨识的输入信号取-0.03elseu(i)=0.03;%当M序列的值为0时，辨识的输入信号取0.03end%小循环结束y1=x1;y2=x2;y3=x3;y4=x4;%为下一次的输入信号做准备end%大循环结束，产生输入信号ufigure(1);%第1个图形stem
10、(u),gridon%以径的形式显示出输入信号并给图形加上网格z(2)=0;z(1)=0;%取z的前两个初始值为零fork=3:15;%循环变量从3到15z(k)=1.5*z(k-1)-0.7*z(k-2)+u(k-1)+0.5*u(k-2);%给出理想的辨识输出采样信号end%RLS递推最小二乘辨识c0=0.0010.0010.0010.001;%直接给出被辨识参数的初始值,即一个充分小的实向量p0=10人6*eye(4,4);%直接给出初始状态P0，即一个充分大的实数单位矩阵E=0.000000005;%相对误差E=0.000000005c=c0,zeros(4,14);%被辨识参数矩阵的
11、初始值及大小e=zeros(4,15);%相对误差的初始值及大小fork=3:15;%开始求Kh1=-z(k-1),-z(k-2),u(k-1),u(k-2);x=h1*p0*h1+1;x1=inv(x);%开始求K(k)k1=p0*h1*x1;%求出K的值d1=z(k)-h1*c0;c1=c0+k1*d1;%求被辨识参数ce1=c1-c0;%求参数当前值与上一次的值的差值e2=e1./c0;%求参数的相对变化e(:,k)=e2;%把当前相对变化的列向量加入误差矩阵的最后一列c0=c1;%新获得的参数作为下一次递推的旧参数c(:,k)=c1;%把辨识参数c列向量加入辨识参数矩阵的最后一列p1=
12、p0-k1*k1*h1*p0*h1+1;%求出p(k)的值p0=p1;%给下次用ife20.5,u(i)=-1;%M序列的值为1时，辨识的输入信号取-1elseu(i)=1;%M序列的值为0时，辨识的输入信号取1endy1=x1;y2=x2;y3=x3;y4=x4;%为下一次的输入信号作准备endfigure(1);%画第一个图形subplot(2,1,1);%画第一个图形的第一个子图stem(u),gridon%画出M序列输入信号v=randn(1,60);%产生一组60个正态分布的随机噪声subplot(2,1,2);%画第一个图形的第二个子图plot(v),gridon;%画出随机噪声信
13、号R=corrcoef(u,v);%计算输入信号与随机噪声信号的相关系数r=R(1,2);%取出互相关系数u%显示输入型号V%显示噪声型号z=zeros(7,60);zs=zeros(7,60);zm=zeros(7,60);zmd=zeros(7,60);%输出采样、不考虑噪声时系统输出、不考虑噪声时模型输出、模型输出矩阵的大小z(2)=0;z(1)=0;zs(2)=0;zs(1)=0;zm(2)=0;zm(1)=0;zmd(2)=0;zmd(1)=0;%给输出采样、不考虑噪声时系统输出、不考虑噪声时模型输出、模型输出赋初值%增广递推最小二乘辨识c0=0.0010.0010.0010.0010.0010.0010.001;%直
