生成 ls mmse信道估计matlab代码

以下是在Matlab中实现LS和MMSE信道估计的代码：

LS信道估计：

function [h_hat_ls] = ls_channel_est(y,x,p)
% y: 接收到的信号
% x: 发送的信号
% p: 信道的长度
% h_hat_ls: LS估计的信道

% 构建Toeplitz矩阵
X = toeplitz(x,[x(1) zeros(1,p-1)]);
Y = y(p:length(y));
Y = Y';

% LS估计信道
h_hat_ls = inv(X'*X)*X'*Y;
end

MMSE信道估计：

function [h_hat_mmse] = mmse_channel_est(y,x,p,sigma_n)
% y: 接收到的信号
% x: 发送的信号
% p: 信道的长度
% sigma_n: 噪声方差
% h_hat_mmse: MMSE估计的信道

% 构建Toeplitz矩阵
X = toeplitz(x,[x(1) zeros(1,p-1)]);
Y = y(p:length(y));
Y = Y';

% MMSE估计信道
R = X'*X;
R_inv = inv(R);
h_hat_mmse = R_inv*X'*(Y)/(sigma_n + Y'*Y - Y'*X*R_inv*X'*Y);
end

这些函数的输入参数`y`和`x`是接收和发送信号。`p`是信道的长度，`sigma_n`是噪声方差。这些函数的输出是估计的信道系数`h_hat_ls`和`h_hat_mmse`。