EEA算法的Matlab实现

最新推荐文章于 2021-09-25 16:55:32 发布

鲁小志

最新推荐文章于 2021-09-25 16:55:32 发布

阅读量1.5k

点赞数 1

分类专栏：其它算法_除卡尔曼和机器学习 matlab相关文章标签：算法 ARMA模型

本文链接：https://blog.csdn.net/lu_mingzhi/article/details/53413987

版权

matlab相关同时被 2 个专栏收录

14 篇文章 2 订阅

订阅专栏

其它算法_除卡尔曼和机器学习

3 篇文章 0 订阅

订阅专栏

clear all;
close all;
N=8000;%运算循环的次数
x=zeros(3,N);
x(1,:)=randn(1,N);%产生N个正态分布高斯白噪声
x(2,:)=rand(1,N);%产生N个均匀分布的随机信号
x(3,:)=sin(0.01*pi*N);%产生正弦信号
d=zeros(1,N);%真实信号
e=zeros(1,N);%误差值
u=[0.1,0.15,0.2;0.1,0.15,0.2;0.1,0.15,0.2];%设定三种信号的步长
coefficient=[1,-0.5,0.3,-0.2,1,-0.8,0.6,-0.3];%真实系数
w=zeros(N,8);%系数估计值
for n0=1:3%for1开始
for n1=1:length(u(n0,:))%for2开始
    figure((n0-1)*3+n1);
for n=5:N%for3开始
    D=[fliplr(d(1,n-4:n-1)),fliplr(x(n0,n-3:n))];
    d(1,n)=D*coefficient.';
    e(1,n)=d(1,n)-w(n-1,:)*D.';
    w(n,:)=w(n-1,:)+u(n0,n1)*conj(e(1,n))*D;
end%for3结束
for n2=1:8%for4开始
    subplot(2,4,n2);
    plot(1:N,w(:,n2));
    xlabel('n');
    ylabel({['信号',num2str(n0),'系数估计值','(u=',num2str(u(n0,n1)),')']});
end%for4结束
end%for2结束
end%for1结束