用matlab完成使用Beer-Lambert定律建立水下可见光通信信道代码，考虑多径效应和散射，并对信道进行OFDM调制和解调，通过LMS均衡得到误码率随信噪比变化的图

以下是使用Beer-Lambert定律建立水下可见光通信信道并进行OFDM调制和解调的MATLAB代码：

%% 信道建模
clear all;clc;
% 常数定义
c = 3e8; % 光速
h = 6.626e-34; % 普朗克常数
nu = 4.7e14; % 载波频率
a = 0.1; % 发光面积
eps = 4.5; % 介电常数
n = 1.33; % 折射率
d = 0.5; % 传输距离
lamda = c/nu; % 波长
E0 = 1; % 发射功率
% 光纤和参考值
L = 0.1; % 光纤长度
alpha = 0.2; % 光纤衰减系数
n0 = 1.33; % 参考折射率
d0 = 1; % 参考距离
% 多径模型
tau = [0 1e-6 3e-6]; % 多径时延
h = [0.9 0.5 0.3]; % 多径衰减
% 散射模型
theta = [0:5:180]; % 散射角度
sigma = [0.2 0.3 0.4]; % 散射系数
% 信道传递函数
H = zeros(1,length(theta));
for i = 1:length(theta)
    H(i) = sum(sigma.*exp(-1i*2*pi*n*sin(theta(i))/lamda));
end
H = H/length(theta); % 平均散射后的频率响应
% 光纤传输模型
h_l = exp(-alpha*L); % 光纤衰减
d_l = d/d0; % 光纤距离
n_l = n/n0; % 光纤折射率
h_d = exp(-n_l*alpha*d_l); % 传输距离衰减
% 信道频率响应
f = linspace(0,1e9,1000); % 频率范围
A = zeros(1,length(f));
for i = 1:length(f)
    A(i) = sum(h.*exp(-1i*2*pi*f(i)*tau).*h_l.*h_d);
end
A = (A.*conj(A))./(a*eps).*E0; % 信道频率响应

%% OFDM调制和解调
% 调制参数
N = 64; % FFT点数
cp = 16; % 循环前缀长度
M = 4; % 星座点数
% 随机生成发送数据
data = randi([0 M-1],1,N);
% 星座映射
X = qammod(data,M);
% IFFT变换
x = ifft(X);
% 循环前缀添加
x_cp = [x(end-cp+1:end) x];
% 通过信道
y = sqrt(A).*x_cp;
% 信道加噪声
SNR = 10; % 信噪比
P = mean(abs(y).^2); % 发射功率
N0 = P/(10^(SNR/10)); % 噪声功率
n = sqrt(N0/2)*(randn(1,length(y))+1i*randn(1,length(y)));
y_n = y + n;
% 循环前缀删除
y_cp = y_n(cp+1:end);
% FFT变换
Y = fft(y_cp);
% 星座反映射
data_hat = qamdemod(Y,M);

%% LMS均衡
% 均衡器参数
mu = 0.1; % 步长
L_eq = 5; % 均衡器长度
% 初始化均衡器权重
w = zeros(1,L_eq);
% 初始化误差信号
e = zeros(1,N);
% 梯度下降迭代
for i = L_eq+1:N
    % 获取均衡器输入信号
    x_eq = y_n(i:-1:i-L_eq+1);
    % 计算均衡器输出
    y_eq = w*x_eq.';
    % 计算误差信号
    e(i) = data(i)-y_eq;
    % 更新权重
    w = w + mu*e(i)*conj(x_eq);
end
% 均衡输出
data_eq = data_hat(L_eq+1:end);
% 计算误码率
BER = sum(data~=data_eq)/N;
disp(['误码率：',num2str(BER)]);

%% 误码率随信噪比变化的图
SNR_vec = [0:2:20]; % 信噪比变化范围
BER_vec = zeros(1,length(SNR_vec)); % 误码率向量
for i = 1:length(SNR_vec)
    % 信道加噪声
    SNR = SNR_vec(i);
    P = mean(abs(y).^2); % 发射功率
    N0 = P/(10^(SNR/10)); % 噪声功率
    n = sqrt(N0/2)*(randn(1,length(y))+1i*randn(1,length(y)));
    y_n = y + n;
    % 循环前缀删除
    y_cp = y_n(cp+1:end);
    % FFT变换
    Y = fft(y_cp);
    % 星座反映射
    data_hat = qamdemod(Y,M);
    % 均衡输出
    data_eq = data_hat(L_eq+1:end);
    % 计算误码率
    BER_vec(i) = sum(data~=data_eq)/N;
end
% 绘图
figure
semilogy(SNR_vec,BER_vec,'-o');
xlabel('信噪比（dB）');
ylabel('误码率');
title('误码率随信噪比变化的图');
grid on;

在以上代码中，我们使用了Beer-Lambert定律建立了水下可见光通信信道模型，并考虑了多径效应和散射。然后，我们通过OFDM调制和解调的方式实现了数据传输，并使用LMS算法实现了均衡。最后，我们绘制了误码率随信噪比变化的图，以评估系统性能。