最大信号法测角的matlab代码

以下是使用最大信号法测角的MATLAB代码示例：

% 生成信号
fs = 1000; % 采样率
t = 0:1/fs:1-1/fs; % 时间向量
f = 10; % 信号频率
theta = 30; % 信号入射角
s = exp(1j*2*pi*f*cosd(theta)*t); % 生成信号

% 添加噪声
SNR = 10; % 信噪比
noise = randn(size(s)); % 高斯白噪声
power_signal = norm(s)^2/length(s); % 信号功率
power_noise = norm(noise)^2/length(noise); % 噪声功率
scale_factor = sqrt(power_signal/(power_noise*10^(SNR/10))); % 计算噪声缩放因子
noise = noise*scale_factor; % 缩放噪声
x = s + noise; % 信号加噪声

% 最大信号法估计角度
N = length(x); % 信号长度
M = 8; % 阵元数
d = 0.5; % 阵元间距
theta_range = -90:0.1:90; % 角度搜索范围
P = zeros(size(theta_range)); % P(theta)函数
for ii = 1:length(theta_range)
    a = exp(1j*2*pi*f*d*cosd(theta_range(ii))*(0:M-1)'); % 阵列流形
    P(ii) = abs(a'*x/N)^2/abs(a'*a/N)^2; % 计算P(theta)函数
end
[~, index] = max(P); % 最大值索引
angle_est = theta_range(index); % 估计角度
disp(['实际角度：', num2str(theta), ' 估计角度：', num2str(angle_est)]);

此代码使用最大信号法估计了一个10 Hz的入射角为30度的信号，加入了10 dB的高斯白噪声，并使用8个阵元的均匀线阵进行测量。代码输出了实际角度和估计角度。