function?[NKkstns0xssta0flamdadt]=signal_model（）;

N=12;%陣元數(shù)
f=5.926*10^9;%信號(hào)中心頻率Hz
lamda=（3*10^8）/f;%信號(hào)波長(zhǎng)
d=0.5*lamda;%陣元間距
K=100*N;%采樣快拍數(shù)
sta0=0;%信號(hào)入射角
SNR=20;%輸入SNR:dB

sgm2s=1*10^（SNR/10）;%輸入的信號(hào)功率
%***************信號(hào)產(chǎn)生***************
fs=2*f;%采樣頻率
Ts=1/fs;%采樣間隔
k=[1:K];%快拍數(shù)向量
t=Ts*k;%時(shí)間向量
st=sqrt（sgm2s）*exp（j*2*pi*f*t）;%信號(hào)
n=[1:N]‘;
s0=exp（j*2*pi*（n-1）*（d/lamda）*sin（sta0*pi/180））;%信號(hào)的導(dǎo)向矢量
Signal=s0*st;%信號(hào)陣列快拍矩陣
%*************************噪聲*************************
%這里的噪聲為高斯白噪聲，服從標(biāo)準(zhǔn)正態(tài)分布
randn（‘state‘0）;%為了確保每次所產(chǎn)生的隨機(jī)矩陣，均相同。
Noise=randn（NK）;
%******************干擾*******************
%在這里假設(shè)干擾與信號(hào)的頻率相同，只是入射角有所不同。
JNR=5;
thta_J=30;%干擾信號(hào)的入射角
Jam_P=1*10^（JNR/10）;%干擾信號(hào)的功率
Jamming_0=sqrt（Jam_P）*exp（j*2*pi*f*t+pi/6）;%參考陣元接收到的干擾信號(hào)（初始相位設(shè)定為一固定值pi/6）
a_thta_J=exp（j*2*pi*（n-1）*d*sin（thta_J*pi/180）/lamda）;%干擾信號(hào)的導(dǎo)向矢量
Jamming=a_thta_J*Jamming_0;%各陣元接收到的干擾信號(hào)的向量。
x=Signal+Jamming+Noise;
sta=[-90:0.5:90];%角度搜索范圍
A=exp（i*2*pi*（n-1）*（d/lamda）*sin（sta*pi/180））;%搜索角度范