复杂电磁环境下的抗干扰建模

%-------------------------------用户设定内容-------------------------------
%【有用信号】：
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