如何用matlab实现频域块LMS算法

先对t离散化，然后FFT，再画个幅频和相频就可以了。

不过貌似你这个式子有点问题，randn产生1*1024的高斯噪声已经指定了s的维数了，所以t必须去1024个离散值，限定了t的取值了。

你可以参考下面的程序，建议你看懂之后自己编一下。

这是对sin(x)求幅频相频

%用fft求幅频相频

clc%清空

clearall%清除所有变量

closeall%关闭所有窗口

A=10%振幅

fw=50%固有频率

phi=pi/3%相位

step=1000

t=0:1/step:10*pi%时间t

y=A*sin(2*pi*fw*t+phi)%正弦函数y

f=step*(0:256)/512%频率

subplot(3,1,1)%三行一列第一幅图

plot(t,y)%绘制图形

xlabel('t/s','fontsize',13)%横坐标显示t/s，字号13

ylabel('y','fontsize',13)%纵坐标显示y，字号13

title('正弦函数曲线','fontsize',13)%显示标题

Y=fft(y,512)%对y进行傅里叶变换

subplot(3,1,2)%三行一列第二幅图

plot(f,abs(Y(1:257)))%绘制图形

xlabel('f/Hz','fontsize',13)%横坐标显示f/Hz，字号13

ylabel('幅值','fontsize',13)%纵坐标显示幅值，字号13

title('幅频特性曲线','fontsize',13)%显示标题

[value,index]=max(abs(Y))%将abs（Y)最大值点的横坐标、纵坐标分别赋给

value和index

text(f(index),value,sprintf('maxpiont=(%f,%f)',f(index),value),'fonts

ize',13)%显示最大值点坐标

subplot(3,1,3)%三行一列第三幅图

plot(f,angle(Y(1:257))*180/pi)%绘制图形

xlabel('f/Hz','fontsize',13)%横坐标显示f/Hz，字号13

ylabel('相位/°','fontsize',13)%纵坐标显示相位/°，字号13

title('相位特性曲线','fontsize',13)%显示标题

LMS算法步骤：

1,、设置变量和参量：

X(n)为输入向量，或称为训练样本

W(n)为权值向量

e(n)为偏差

d(n)为期望输出

y(n)为实际输出

η为学习速率

n为迭代次数

2、初始化，赋给w(0)各一个较小的随机非零值，令n=0

3、对于一组输入样本x(n)和对应的期望输出d，计算

e(n)=d(n)-X（n）

W(n+1)=W(n)+ηX(n)e(n)

4、判断是否满足条件，若满足算法结束，若否n增加1，转入第3步继续执行。

---