卷积噪声的matlab仿真

clear all close all clc fs=100e6; Ts=1/fs; Tp=26e-6; t=0:Ts:Tp-Ts; f=10e6; signal=sin(2*pi*f*t); signalfft=fft(signal); figure; subplot(211) plot(t,signal) subplot(212) f=[0:length(signalfft)-1]*fs/length(signalfft); plot(f,abs(signalfft)) noise=randn(1,length(signal)); %噪声经过带通滤波器 load lvbo3040.mat b=lvbo3040; a=1; [h,f1]=freqz(b,a,101,'whole',fs); mag=abs(h); figure; plot(f1,mag) u=filter(b,a,noise);%生成带限高斯白噪声p=2.7917w Nfft=2^nextpow2(length(u)); j2=fft(u); figure; subplot(2,1,2),plot(f,abs(j2)) convnoise=conv(u,signal); convnoisefft=fft(convnoise); f=[0:length(convnoisefft)-1]*fs/length(convnoisefft); figure; plot(f,abs(convnoisefft))