【射频噪声干扰】射频噪声信号调频+调幅+调相干扰（幅度谱和频谱图）【含Matlab源码 2897期】.zip

% 清空变量，关闭所有窗口，情况控制台窗口 clear;close all;clc; %% 射频噪声干扰 %信号长度 N=1024; n=1:1:N; %系统采样率 f_s=10e6; T_s=1/f_s; A=1; B_n=2e6;%射频信号要产生的带宽 f_0=2e6; % 调制噪声参数 delta_F=B_n/2;%射频信号要产生的带宽B_n时，噪声要产生的带宽 N_n=round(N*(B_n/(2*f_s)));%频谱上采样点数 S=random('Normal',0,1,1,N)+1i*random('Normal',0,1,1,N);%高斯白噪声频谱,中值0，标准差1，维度(1，N) S_n=S; for i=N_n:1:N-N_n-1 S_n(i)=0; end s_n=ifft(S_n);%对做逆傅里叶变换就可以得到时域的有一定带宽的噪声信号 s_n=s_n/std(s_n,0,2);%归一化,方差为1 %产生需要调制的射频噪声干扰s_0 s_0=exp(1j*2*pi*f_0*n*T_s); s=s_0.*s_n; subplot(3,1,1); plot(n/f_s,s_0); xlabel('us');title('原始信号'); %时域波形 figure(1); 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('射频噪声干扰功率谱');