r(i)-h(i-l,:)*x(i-n+1:i),;h(i,:)=h(i-l,:)+e(l,i)*g(i,:);Cxx=Cxx-g(i,:)’*x(i-n+1:i)本Cxx;y(1,i)=h(i,:)*x(i-n+1:i)’;tr(1,i)=trace(Cxx);endfigure(l)plot(0:N~n,x(l,n:N)),gridtitleCx(k)inputsingnalinV’)xlabelCSamples')Figure(2)plot(0:N-n,xr(1,n:N),’r’),gridaxis([0800-1.21.2])title(,xr(k)referencesingnalinV’)xlabelCSamples')figure(3)plot(0:N-n,e(l,n:N)),holdonplot(0:N-n,y(l,n:N),’r’),holdongridtitle(5e(k)errorandy(k)outputinV’)xlabel(’Samples’)figure(4)plot(0:N-n,h(n:N,l)),holdonplot(0:N-n,h(n:N,2),’r’),holdoffgridtitleCa(n-l)anda(n~2)coeffcientsevolution')xlabel(’Samples’)figure(5)numl=fliplr(h(N,:));sysl=tf(numl,1,Ts);bode(sysl),holdofftitle(’Synthesizedfilter’)xlabelCFrequencyinrad/s’)ylabel(’Phaseindegree;ModuleindB’)figure(6)semi1ogy(0:N~n,tr(n:N)),gridtitle(’Cxxmatrixtrace’)xlabelCSamples')
全文预览
