因为某些原因把

I2((mm-1)*M2+nn)=exp(j*2*pi*(mm-1)*d*sin(ThetaJ2)*cos(PhiJ2)/lamda)*exp(j*2*pi*(nn-1)*d*sin(ThetaJ2)*sin(PhiJ2)/lamda);

改写成了

I2((mm-1)*M2+nn)=exp(j*(2*pi*(mm-1)*d*sin(ThetaJ2)*cos(PhiJ2)/lamda+2*pi*(nn-1)*d*sin(ThetaJ2)*sin(PhiJ2)/lamda));

基本就是exp(ja)*exp(jb)=exp(j*(a+b))这么简单的改动，但是报错了，想不通哪里错了(本想这样改写成C语言的时候简单些)

全部代码是这样的：

clc

clear all

close all

% close all

freqs=1.268e9;                       %工作频率

c=3e8;                               %光速

lamda=c/freqs;                       %波长

d=0.5*lamda;                           %单元间距

M1=2;

M2=2;

M=M1*M2;                             %天线阵元数

N=9;                                 %时延步长

pd=32;                             %快拍数

pdd = 1024;

%% digital parameter

fs=62e6;                              %采样频率

fc=15.48e6;

h = fir1(32,20e6/fs);

h = conv(h,h);

%%

ThetaJ2=30*pi/180;                  %干扰方向

PhiJ2=40*pi/180;

freqJ2=5e5;                        %第1个干扰的频率

J2NR=60;

AJ2=sqrt(10^(J2NR/10));

Jam2one=AJ2*exp(j*(2*pi*freqJ2*(1:1:pdd)/fs));

% %------------------------------------------第一个干扰模型的建立--------------------

ThetaJ1=60*pi/180;                  %干扰方向

PhiJ1=70*pi/180;

freqJ1=15e5;                        %第2个干扰的频率

J1NR=60;

AJ1=sqrt(10^(J1NR/10));

Jam1one=AJ1*exp(j*(2*pi*freqJ1*(1:1:pdd)/fs));

I2=zeros(M,1);

I1=zeros(M,1);

%--------------------------------------方阵---------------------------------

for mm=1:M1

for nn=1:M2

I2((mm-1)*M2+nn)=exp(j*2*pi*(mm-1)*d*sin(ThetaJ2)*cos(PhiJ2)/lamda)*exp(j*2*pi*(nn-1)*d*sin(ThetaJ2)*sin(PhiJ2)/lamda);

I1((mm-1)*M2+nn)=exp(j*2*pi*(mm-1)*d*sin(ThetaJ1)*cos(PhiJ1)/lamda)*exp(j*2*pi*(nn-1)*d*sin(ThetaJ1)*sin(PhiJ1)/lamda);

end

end

Jam2=I2*Jam2one;

Jam1=I1*Jam1one;

JAM2=Jam2.';

JAM1=Jam1.';

% JAM4=Jam4.';

%

% %------------------------------------------产生噪声-------------------------

noise=sqrt(1/2)*randn(pdd,M)+j*sqrt(1/2)*randn(pdd,M);

% %------------------------------------------产生接收信号-------------------------

%%

Receive=JAM1+JAM2+noise;

% Receive=JAM1+noise;

Rbak = Receive(1:32,:);

ddcen = 1;

if(ddcen==1)

%

Receive = conj(Receive');

R1 = Receive(1,:);

R2 = Receive(2,:);

R3 = Receive(3,:);

R4 = Receive(4,:);

NN = length(R1);

nn = 0:1:NN-1;

coho = exp(1i*2*pi*fc*nn/fs);

%%

IfR1 = real(R1.*coho);

IfR2 = real(R2.*coho);

IfR3 = real(R3.*coho);

IfR4 = real(R4.*coho);

coho = exp(-1i*2*pi*fc*nn/fs);

% IfR1 = R1.*coho;

% IfR2 = R2.*coho;

% IfR3 = R3.*coho;

% IfR4 = R4.*coho;

%% ddc

IfR1 = coho.*IfR1;

IfR2 = coho.*IfR2;

IfR3 = coho.*IfR3;

IfR4 = coho.*IfR4;

%% 滤波

S1 = conv(IfR1,h);

S2 = conv(IfR2,h);

S3 = conv(IfR3,h);

S4 = conv(IfR4,h);

%%

S1 = S1(257:320);

S2 = S2(257:320);

S3 = S3(257:320);

S4 = S4(257:320);

Receive = [S1;S2;S3;S4];

Receive = conj(Receive');

% Receive = reshape(Receive,1,size(Receive,1)*size(Receive,2));

end

% 功率倒置 %

%------------------------------------------常规求最优权值方法-------------------------

% tic;

RR=zeros(M,M);

RR=((Receive.')*(Receive.')')/pd;

R_R=Receive(:,1)+Receive(:,2)+Receive(:,3)+Receive(:,4);

y_w0=Receive.'*(Receive.')'/pd ;                         %未加权值前的输出

deta=[1,zeros(1,(M-1))].';

wopt=inv(RR)*deta*inv(deta'*inv(RR)*deta)   %功率倒置

% wopt_fi = floor(wopt/max(abs(wopt))*(2^(BW-1)));

%------------------------------------------画二维方向图----------------------

cida=-pi/2:pi/180:pi/2;

phi=-pi/2:pi/180:pi/2;

ff=10e6:0.05e6:20e6;

for jj=1:1:length(cida)                                         %空间角度扫描

for gg=1:1:length(phi)

w_s_scan=zeros(M,1);

sss=cida(jj);

ttt=phi(gg);

for mm=1:1:M1

for nn=1:1:M2

w_s_scan((mm-1)*M2+nn)=exp(j*2*pi*(mm-1)*d*sin(sss)*cos(ttt)/lamda)*exp(j*2*pi*(nn-1)*d*sin(sss)*sin(ttt)/lamda);

end                                                      %空间导向矢量 1*M

end

F2(jj,gg)=abs(wopt'*w_s_scan);                               %阵列方向图

end

end

Fmax2=max(max(F2));

Fa2=F2./Fmax2;

Fa2=20*log10(Fa2);

figure(1);

mesh(180*cida/pi,180*phi/pi,Fa2.');

% title('空域功率倒置');

ylabel('Phi方向');

xlabel('Theta方向');

zlabel('增益(dB)');

grid on;