gsshpzsh.m

高斯噪声产生的代码,含有射频噪声干扰的代码

close all;
clear all;
clc;
t=1e-3;%产生高斯分布的噪声%
fs=1e7;
t1=0:1/fs:t-1/fs;
n=length(t1);
randn('state',0);
u=randn(1,n);

wp=[2.5e6 3.5e6];%加滤波器
ws=[2e6 4e6];
rp=1;
rs=60;
[N,wn]=buttord(wp/(fs/2),ws/(fs/2),rp,rs);
[b,a]=butter(N,wn);
u1=filter(b,a,u);

J=abs(fft(u1,512));%噪声信号频谱分析
f=fs*(0:256)/512;

figure,
subplot(2,2,1),plot(u),title('高斯分布噪声'),xlabel('t(单位：s)');
subplot(2,2,2),hist(u,-4:0.02:4),title('高斯分布信号直方图');
subplot(2,2,3),plot(u1),title('滤波后噪声信号'),xlabel('t(单位：s)');
subplot(2,2,4),plot(f,J(1:257)),title('射频噪声干扰信号频谱图');