%-------------------------------用戶設(shè)定內(nèi)容-------------------------------
%【有用信號(hào)】：
As=10;??%設(shè)定有用信號(hào)的振幅，單位：?V；
fs=3*10^8;??%設(shè)定有用信號(hào)的頻率單位：Hz；
phi=unifrnd（0pi/2）;
theta=unifrnd（02*pi）;??%設(shè)定有用信號(hào)的來波方向，服從均勻分布，單位：rad；

%【干擾信號(hào)】：
SNR=20;??%設(shè)定信噪比，單位：dB；
m=100;??%設(shè)定干擾源的數(shù)量m；
A=unifrnd（0.5*As*10^（-SNR/20）2*As*10^（-SNR/20）1m）;??%設(shè)定干擾源的振幅A，服從均勻分布，單位：V；
W=normrnd（6*pi*10^82*10^91m）;??%設(shè)定干擾源的角頻率W，以有用信號(hào)頻率為均值，服從正態(tài)分布，單位：rad/s；
PHI=unifrnd（0pi/21m）;
THETA=unifrnd（02*pi1m）;??%設(shè)定干擾源的來波方向，服從均勻分布，單位：rad；

%【信道】：
PHIh=unifrnd（-pipi1m+1）;??%設(shè)定信道產(chǎn)生的相位，服從均勻分布，單位：rad；
Th=unifrnd（00.051m+1）;??%設(shè)定信道產(chǎn)生的時(shí)延，服從均勻分布，單位：s；

%【天線陣】：
X=[1234];Y=[1234];Z=[1234];??%設(shè)置天線陣單元的坐標(biāo)，單位：m；

%【系統(tǒng)白噪聲】：
RND=50;??%設(shè)定系統(tǒng)白噪聲的信噪比，單位：dB；

%?【FFT相關(guān)設(shè)定】：
fc=40*fs;%采樣頻率設(shè)定；
N=2048;%采樣點(diǎn)數(shù)設(shè)定；

%---------------------------------程序部分---------------------------------
n=0:N-1;
t=n/fc;??%設(shè)定觀察圖像的時(shí)間序列；
f=（n*fc/N）*10^-6;??%設(shè)定觀察圖像的頻率序列；

SIG=As*sin（2*pi*fs*t）;??%存儲(chǔ)有用信號(hào)時(shí)域函數(shù)；
DIS=A（1）*sin（W（1）*t）;
for?i=2:m
????dis=A（i）*sin（W（i）*t）;
????DIS=cat（1DISdis）;
end???%此循環(huán)是將m個(gè)干擾源的時(shí)域函數(shù)存入m行矩陣DIS中；


Hs=SIG*exp（1i*PHIh（m+1））*exp（-1i*W（1）*Th（m+1））;??%有用信號(hào)通過信道響應(yīng)；
dis=DIS（1:）;
H=dis*exp（1i*PHIh（1））*exp（-1i*W（1）*Th（1））;
for?i=2:m
????dis=DIS（i:）;
????h=dis