function?[NKkstns0xssta0flamdadt]=signal_model（）;

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

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