【噪声干扰】噪声调幅干扰（幅度谱和频谱图）【含Matlab源码 2899期】

⛄一、获取代码方式

获取代码方式1：
完整代码已上传我的资源：【噪声干扰】基于matlab噪声调幅干扰（幅度谱和频谱图）【含Matlab源码 2899期】
点击上面蓝色字体，直接付费下载，即可。

获取代码方式2：
付费专栏Matlab信号处理（初级版）

备注：
点击上面蓝色字体付费专栏Matlab信号处理（初级版），扫描上面二维码，付费29.9元订阅海神之光博客付费专栏Matlab信号处理（初级版），凭支付凭证，私信博主，可免费获得1份本博客上传CSDN资源代码（有效期为订阅日起，三天内有效）；
点击CSDN资源下载链接：1份本博客上传CSDN资源代码

⛄二、简介

有源干扰是现代电子战中的一种主要干扰形式。它在干扰制导系统的同时, 又干扰无线电引信, 会得到更显著的干扰效果。尤其是对制导系统的干扰难以奏效时, 干扰武器系统的最后一个工作环节——无线电引信, 可起到事半功倍的效果。对无线电引信的人为有源干扰与一般雷达的干扰基本相同, 大体分为压制式干扰和回答式欺骗干扰两大类。

对引信威胁最大的干扰是压制式干扰中的阻塞式干扰, 阻塞式干扰发射宽频带的干扰信号, 可对频带内的引信同时进行干扰, 为此要求干扰机发射宽频谱的大功率干扰信号, 而噪声调频干扰具有宽的干扰带宽和较大的噪声功率, 是目前对雷达、引信、通信进行阻塞式干扰中最常用的干扰形式。因此, 本文研究噪声调频干扰信号, 利用计算机产生其时域波形、频谱、功率谱, 分析噪声调频信号的特点;同时, 计算脉冲多谱勒引信天线接收到的噪声调频干扰有效噪声功率;利用计算机产生混频器及解码输出端干扰信号的时域信号;最后, 将计算机产生的噪声调频信号与干扰模拟器产生的噪声调频干扰信号比较, 验证本文所产生的噪声调频干扰仿真的正确性。

1 噪声调频干扰信号时域表达式、频谱和功率谱仿真分析
噪声调频干扰信号最常见的是射频振荡的频率与调制噪声电压ξ (t) 呈线性关系, 噪声调频干扰电压的时域表达式为[1]:

uj (t) =Ecos[ω0t+∫t0kfξ (t) dt] (1)

式中:E为电压振幅, 为常数 (忽略其寄生调幅) ;kf为调频指数;ξ (t) 为高斯噪声, 其均值为0, 方差为σ2ξ, 在频带Ω0内为均匀频谱。

设调制噪声电压ξ (t) 是高斯噪声, 其幅度概率密度分布为高斯函数。

由于噪声调频干扰的角频率与ξ (t) 呈线性关系, 故瞬时角频率或角频偏的概率密度亦应为高斯分布, 其均方根值为:

其中fe为有效频偏。由此得噪声调频干扰的频谱函数为:

噪声调频干扰的功率谱为:

噪声调频干扰信号的频带宽Δfj, 可按半功率电平对应的两个频率之间的差值来定义, 即:

利用计算机产生式 (1) 的噪声调频干扰时域信号uj (t) , 并应用快速傅立叶变换在计算机上产生噪声调频干扰的频谱信号。图1是中心频率为4.3 GHz、调频指数为0.2 GHz/V、调制噪声功率为2.0 dBW的噪声调频时域信号和同步频谱信号。图2是计算机利用式 (3) 产生的调制噪声功率分别为2.0 dBW和4.0 dBW, 调频指数为0.2 GHz/V, 中心频率为4.35 GHz的噪声调频频谱信号曲线。图3是计算机利用式 (3) 产生的中心频率分别为4.35 GHz和2.0 GHz, 调制噪声功率2.0 dBW, 调频指数为0.2 GHz/V的噪声调频频谱信号曲线。图4为我所研制的7.5 cm频端噪声调频引信干扰模拟器产生的噪声调频干扰信号频谱。

⛄三、部分源代码

%--------------------噪声调幅干扰------------------
%信号长度
N=10000;
n=1:1:N;
%系统采样率
f_s=100e6;
T_s=1/f_s;
B_n=2e6;%射频信号要产生的带宽
f_0=10e6;
m_A=1;%有效调制系数
% 调制噪声参数
delta_F=B_n/2;%射频信号要产生的带宽B_n时，噪声要产生的带宽
N_n=round(N*(B_n/(2f_s)));%频谱上采样点数
S_n=random(‘Normal’,0,1,1,N)+1irandom(‘Normal’,0,1,1,N);%高斯白噪声频谱,中值0，标准差1，维度(1，N)
for i=N_n:1:N-N_n-1
s=(A+real(s_n)).cos(2pif_0n*T_s);

figure(1);
subplot(3,1,1);
plot(n/f_s,s_0);
xlabel(‘us’);title(‘原始信号’);
subplot(3,1,2);
plot(n/f_s,real(s));
xlabel(‘us’);title(‘噪声调幅干扰时域波形’);
%幅度谱
subplot(3,1,3);
k = -(N-1)/2:(N-1)/2;
f = k/N*f_s;
plot(f,abs(fftshift(fft(s)))/N);
xlabel(‘hz’);title(‘噪声调幅干扰功率谱’);

⛄四、运行结果

⛄五、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1]夏红娟,陈潜.噪声调频干扰信号仿真及应用[J].上海铁道大学学报. 2000,(06)

3 备注
简介此部分摘自互联网，仅供参考，若侵权，联系删除