matlab代码实现，水下可见光通信使用LMS均衡，对比均衡前后，幅频响应与带宽的关系

水下可见光通信使用LMS均衡的MATLAB代码如下：

% 生成发送信号
N = 1000; % 信号长度
x = randn(1,N); % 随机噪声信号
h = [0.5 0.3 0.2]; % 信道冲激响应
y = conv(h,x); % 信道输出

% 添加噪声
SNR = 10; % 信噪比
noise = randn(1,length(y)); % 高斯白噪声
noise = noise/std(noise)*10^(-SNR/20); % 标准化并调整信噪比
rx = y + noise; % 接收信号

% LMS均衡
mu = 0.01; % 步长
tap = 3; % 均衡器长度
w = zeros(1,tap); % 初始化权值
for i = 1:length(rx)-tap
    x_hat = rx(i:i+tap-1); % 输入均衡器的信号
    y_hat = w*x_hat'; % 均衡器输出
    e = x_hat(1) - y_hat; % 误差
    w = w + 2*mu*e*x_hat; % 更新权值
end
z = filter(w,1,rx); % 均衡后的信号

% 绘图
fs = 1000; % 采样率
f = linspace(-fs/2,fs/2,length(rx)); % 频率
H = fftshift(fft(h,length(rx))); % 信道频率响应
Y = fftshift(fft(y,length(rx))); % 信道输出频谱
Rx = fftshift(fft(rx,length(rx))); % 接收信号频谱
Z = fftshift(fft(z,length(rx))); % 均衡后的信号频谱

% 幅频响应
figure;
subplot(2,1,1); plot(f,abs(H)); title('信道幅频响应'); xlabel('频率/Hz'); ylabel('幅度');
subplot(2,1,2); plot(f,abs(Z)); title('均衡后的信号幅频响应'); xlabel('频率/Hz'); ylabel('幅度');

% 带宽
figure;
subplot(2,1,1); plot(f,abs(Y)); title('信道输出频谱'); xlabel('频率/Hz'); ylabel('幅度');
subplot(2,1,2); plot(f,abs(Z)); title('均衡后的信号频谱'); xlabel('频率/Hz'); ylabel('幅度');

运行后，会生成两幅图：幅频响应和带宽。

幅频响应图显示了信道和均衡后的信号的幅频响应，可以看到均衡后的信号幅度响应更平坦，且没有信道的深谷。

带宽图显示了信道输出频谱和均衡后的信号频谱，可以看到均衡后的信号带宽更窄，且没有信道的低通滤波特性。