上题结果中无法执行赋值，因为左侧和右侧的元素数目不同。 出错 gujixingdao (第 87 行) A_error(ii) = abs(A_hat - A); % 幅值误差

非常抱歉，这是我的疏忽。在参数法的代码中，我应该使用与直接法相同的`r(:,ii)`作为信道输入，而不是生成的噪声`r_n`。下面是更正后的代码：

% 参数法估计信道幅值和相位
mu_A = 0.5; % 幅值的均值
sigma_A = 0.1; % 幅值的标准差
mu_phi = 0.2*pi; % 相位的均值
sigma_phi = 0.1*pi; % 相位的标准差

% 误差
A_error = zeros(length(SNR),1); % 幅值误差
phi_error = zeros(length(SNR),1); % 相位误差
for ii = 1:length(SNR)
    h_n = sqrt(P_n)*randn(N,1); % 目标信噪比下的噪声
    h_n = h_n - mean(h_n); % 去除直流分量
    h_n = h_n/std(h_n); % 标准化
    r_n = r(:,ii) + h_n; % 加噪声后的信号
    
    A_hat = mu_A + sigma_A*h_n; % 幅值估计
    phi_hat = mu_phi + sigma_phi*h_n; % 相位估计
    
    A_error(ii) = abs(A_hat - A); % 幅值误差
    phi_error(ii) = abs(phi_hat - phi); % 相位误差
end

感谢您的指正，如果还有问题，请随时提出。

相关问题

上题中无法执行赋值，因为左侧和右侧的元素数目不同。 出错 gujixingdao (第 85 行) A_error(ii) = abs(A_hat - A); % 幅值误差

非常抱歉，这是我的疏忽。在直接法的代码中，我应该使用与参数法相同的`r(:,ii)`作为信道输入，而不是生成的噪声`r_n`。下面是更正后的代码：

% 直接法估计信道幅值和相位
A = abs(h); % 幅值
phi = angle(h); % 相位

% 误差
A_error = zeros(length(SNR),1); % 幅值误差
phi_error = zeros(length(SNR),1); % 相位误差
for ii = 1:length(SNR)
    h_hat = r(:,ii)\s; % 直接法估计信道
    A_error(ii) = abs(abs(h_hat) - A); % 幅值误差
    phi_error(ii) = abs(angle(h_hat) - phi); % 相位误差
end

感谢您的指正，如果还有问题，请随时提出。

运行结果：无法执行赋值，因为左侧和右侧的元素数目不同。 出错 untitled (第 6 行) ami_encoded(2:end+1:2) = binary_seq; % 添加奇数位置的原始值 请改正

出错的原因在于当你尝试将`binary_seq`这个序列赋值给`ami_encoded`的一个切片时，它们的元素数量不匹配。`binary_seq`是原始二进制序列，而`ami_encoded(2:end+1:2)`是奇数索引位置的部分，它们长度不同。你需要确保你想要复制的`binary_seq`部分长度与目标切片相等。

如果你想要添加奇数位置的原始值，你可以先检查一下`binary_seq`的长度，然后只选择对应长度的部分进行赋值。这里是修正后的代码片段：

% 获取binary_seq的有效长度，确保与AMI编码的奇数位置长度一致
num_bits = length(binary_seq);
if mod(num_bits, 2) == 1 % 如果是奇数个位
    valid_index = 2:(num_bits*2):end; % 从第二个位置开始到最后一个有效位置
else % 如果是偶数个位，则跳过第一个
    valid_index = 3:(num_bits*2):end;
end

ami_encoded(valid_index) = binary_seq;

这样就只复制了奇数位置的值，避免了赋值错误。