gdmu.m

在matlab下实现系统辨识设计

x  = randn(1,700);                     % 作为系统的输入
       b  = fir1(63,0.5)                       % 确定未知系统
       n  = 0.1*randn(1,700);                  % 加入噪声
       d  = filter(b,1,x)+n;                     % 期望信号
mu =0.003;                             % 固定步长值
        h = adaptfilt.lms(64,mu);                % 模拟建模
     [y,e] = filter(h,x,d);                       % 模拟系统输出
        S =abs(h.Coefficients.')-abs(b.');           % 误差取绝对值作差
        subplot(3,2,1); plot(1:700,d);axis([0,700,-3,3]);
        title(‘未知系统输出’);xlabel(‘采样点’); ylabel(‘输出’); 
        subplot(3,2,3); plot(1:700,y);axis([0,700,-3,3]);
        title(‘模拟系统输出’);xlabel(‘采样点’); ylabel(‘输出’); 
        subplot(3,2,5); plot(1:700,e);axis([0,700,-3,3]);
        title(‘误差值’);xlabel(‘采样点’); ylabel(‘误差’); 
        subplot(3,2,2);stem([b.']);axis([0,63,-0.1,0.6]);
        title(‘未知系统采样点’);xlabel(‘采样点’); ylabel(‘数值’); 
        subplot(3,2,4);stem([h.Coefficients.']);axis([0,63,-0.1,0.6]);
        title(‘模拟系统采样点’);xlabel(‘采样点’); ylabel(‘数值’);
        subplot(3,2,6);stem([s]);axis([0,63,-0.03,0.03]);
        title(‘误差比对’);xlabel(‘采样点’); ylabel('误差大小');