自适应波束形成matlab代码

自适应波束形成（Adaptive Beamforming）是一种信号处理技术，用于抑制干扰并增强感兴趣信号。下面是一个简单的自适应波束形成的 MATLAB 代码示例：

% 生成两个信号源
fs = 8000; % 采样率
t = 0:1/fs:1-1/fs; % 时间向量
f1 = 1000; % 第一个信号源频率
f2 = 2000; % 第二个信号源频率
s1 = cos(2*pi*f1*t); % 第一个信号源
s2 = cos(2*pi*f2*t); % 第二个信号源

% 生成两个干扰源
f3 = 3000; % 第一个干扰源频率
f4 = 4000; % 第二个干扰源频率
n1 = randn(size(t)); % 第一个白噪声干扰源
n2 = randn(size(t)); % 第二个白噪声干扰源
a1 = 0.5*sin(2*pi*f3*t); % 第一个正弦干扰源
a2 = 0.5*sin(2*pi*f4*t); % 第二个正弦干扰源

% 生成总信号
x = s1 + s2 + n1 + n2 + a1 + a2; % 总信号

% 自适应波束形成
N = 8; % 线阵元数
theta = [20 30]; % 信号源角度
theta = deg2rad(theta); % 角度转弧度
lambda = 1; % 波长
d = lambda/2; % 阵元间距
k = 2*pi/lambda; % 波数
W = ones(N,1); % 初始权值
for i = 1:100 % 迭代次数
    % 计算波束指向
    a = exp(-1i*k*d*(0:N-1)'*sin(theta));
    % 计算输出
    y = W'*a*x;
    % 计算误差
    e = a'*(y - x);
    % 更新权值
    W = W - 0.1*e.*conj(a);
end

% 绘制波束图
theta_range = -90:0.1:90;
theta_range = deg2rad(theta_range);
a_range = exp(-1i*k*d*(0:N-1)'*sin(theta_range));
P = abs(W'*a_range).^2;
plot(rad2deg(theta_range),P)
xlabel('Angle (degree)')
ylabel('Power')
title('Adaptive Beamforming')

上述代码生成了两个信号源和两个干扰源，然后进行了自适应波束形成。最后，绘制了波束图来显示自适应波束形成的效果。通过改变信号源角度，可以看到波束指向的变化。