matlab-(含教程)基于随机共振发的低信噪比信号去噪matlab仿真

function [noiseFinal] = createNoise(duration,samplingRate, amplitude, lf, hf, plotWav) npoints = duration * samplingRate; noise = (amplitude/3)* randn( 1, npoints ); time = (0:npoints-1) / samplingRate; X = [1; exp(i*2*pi*rand(npoints/2-1,1)); 1]; freqbins = (0:npoints/2)'/npoints*samplingRate; X(find((freqbins < lf) | (freqbins > hf))) = 0; X = [X; conj(flipud(X(2:end-1)))]; noise = real(ifft(X)); noiseFinal = amplitude*noise/prctile(noise, 99.999); timePlot = 1/samplingRate:1/samplingRate:duration; binwidth = 1/2*std(noiseFinal); numbins = round(range(noiseFinal)/binwidth); if isequal(plotWav,'on') h = figure; subplot(1,2,1); plot(timePlot,noiseFinal); xlabel('Duration in sec'); ylabel('Amplitude in mA'); title(['Noise for Frequencies between ' num2str(lf) 'Hz to ' num2str(hf) 'Hz (Sampling Rate: ' num2str(samplingRate) 'Hz)']) subplot(1,2,2); histfit(noiseFinal,100) xlabel('Amplitudes in mA'); ylabel('Frequency of Amplitudes'); title(['Distribution of Amplitudes over whole signal']) if ~exist('signals','dir') mkdir('signals') end set(h,'PaperPositionMode','auto'); set(h,'PaperOrientation','landscape'); set(h,'Position',[50 50 1200 800]); print(h,['signals/noise_' num2str(lf) '_' num2str(hf) '_' num2str(samplingRate) '_' num2str(amplitude) '_' num2str(duration) '.pdf'],'-dpdf','-r0') audiowrite(['signals/noise_' num2str(lf) '_' num2str(hf) '_' num2str(samplingRate) '_' num2str(amplitude) '_' num2str(duration) '.wav'], noiseFinal, samplingRate) end end