【噪声干扰】基于matlab噪声调相干扰（幅度谱和频谱图）【含Matlab源码 2898期】.zip

% 清空变量，关闭所有窗口，情况控制台窗口 clear;close all;clc; %% %--------------------噪声调相干扰------------------ %信号长度 N=10000; n=1:1:N; %系统采样率 f_s=100e6; T_s=1/f_s; B_n=2e6;%射频信号要产生的带宽 f_0=10e6; K_PM=5; A=1; % 调制噪声参数 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);%归一化 s_n=real(s_n); %调相干扰 s_0=A*cos(2*pi*f_0*n*T_s); s=A*cos(2*pi*f_0*n*T_s+K_PM*s_n); %一直要产生B_n带宽的噪声 D=K_PM*std(s_n,0,2); if D>5 delta_F=2*sqrt(2*log(2))*sqrt((D^2*B_n^2)/3); %fprintf('%6.2f\n',D^2); else D<0.1 delta_F=2*B_n; end %图像 figure(1); 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('噪声调相干扰功率谱');