zishiying.m

来自「这是我师兄做的空间谱估计的内容。这只是一小部分程序」· M 代码 · 共 19 行

clear;clc
%set up parameters
N=100;%输入序列个数
delta=0.0000001;
x=rand(1,N,[-10,10]);
b=1:N;
a=-1:N-1;
d=filter(b,a,x);
[h,y]=lms(x,d,delta,N);
e=d-y;
[h1,w]=freqz(b,a)
h1f=abs(h1);
clf
figure(1);
plot(w,h1f);title('未知系统幅频特性曲线');
[h2,w]=freqz(h,1);
h2f=abs(h2);
figure(2);
plot(w,h2f);title('FIR系统幅频特性曲线');