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

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

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

%【天線陣】：
X=[1234];Y=[1234];Z=[1234];??%設置天線陣單元的坐標，單位：m；

%【系統白噪聲】：
RND=50;??%設定系統白噪聲的信噪比，單位：dB；

%?【FFT相關設定】：
fc=40*fs;%采樣頻率設定；
N=2048;%采樣點數設定；

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

SIG=As*sin（2*pi*fs*t）;??%存儲有用信號時域函數；
DIS=A（1）*sin（W（1）*t）;
for?i=2:m
????dis=A（i）*sin（W（i）*t）;
????DIS=cat（1DISdis）;
end???%此循環是將m個干擾源的時域函數存入m行矩陣DIS中；


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