% 生成线性调频信号和主瓣噪声干扰信号 fs = 1000; % 采样率 T = 1; % 信号时长 t = 0:1/fs:T-1/fs; % 时间序列 f0 = 50; % 起始频率 f1 = 200; % 终止频率 x = chirp(t,f0,T,f1); % 线性调频信号 n = 0.5*sin(2*pi*100*t); % 主瓣噪声干扰信号 s = x + n; % 混合信号 % 盲源分离 s_hat = fastica(s); % 使用fastICA算法 % 计算方向图 fc = 100; % 中心频率 N = 1024; % FFT长度 d = 0.5; % 元件间距 theta = -90:0.1:90; % 方向角 w = ones(size(s_hat,1),1); % 权重 S = fft(s_hat,N,2); % FFT P = abs(sum(repmat(w,1,N).*S.*exp(-1j*2*pi*fc*(0:N-1)*d*sin(theta*pi/180)),1)).^2; % beam pattern MATLAB报错显示矩阵维度不一致

根据你提供的代码，我看到主要问题出现在计算方向图的部分。在计算权重w时，你使用了s_hat的行数来创建一个列向量，但是在计算FFT时，你使用了N来计算每个信号的FFT，这会导致FFT结果的维度与权重w不匹配，从而导致矩阵维度不一致的错误。

为了解决这个问题，你可以在计算FFT时使用s_hat的列数来保持维度匹配，即将这行代码：

S = fft(s_hat,N,2); % FFT

替换为：

S = fft(s_hat,size(s_hat,2),2); % FFT

这样做可以确保FFT结果的维度与权重w相同，从而避免矩阵维度不一致的错误。