粉红噪声matlab产生

粉红噪声是一种无规则的噪声，具有连续的噪声谱，其功率与频率呈反比，频率成分比白噪声多，在固定的倍频程带宽内，声能量相同的。如低频125~250Hz带宽与高频2K~4K带宽内能量是相同的。

白噪声属于全频噪声，主要用来测量电声设备，粉红噪声是根据人耳的听觉特性与习惯
设计出的测量噪声，主要用来模拟人的听觉，测量厅堂的声学特性

function out = create_pink_noise(Fs, Sec, Amp)% Creates a pink noise signal and saves it as a wav file
%
% Usage: create_noise(Fs, Sec, Amp);
%
%        Fs is the desired sampling rate
%        Sec is the duration of the signal in seconds
%        Amp is the amplitude in dB of the signal (0dB to -144dB)
%
% Author: sparafucile17 06/14/02%error trapping
if((Amp > 0) || (Amp < -144))error('Amplitude is not within the range of 0dB to -144dB');
end%Create Whitenoise
white_noise = randn((Fs*Sec)+1,1);%Apply weighted sum of first order filters to approximate a -10dB/decade
%filter.  This is Paul Kellet's "refined" method (a.k.a instrumentation
%grade)  It is accurate to within +/-0.05dB above 9.2Hz
b=zeros(7,1);
for i=1:((Fs*Sec)+1)b(1) = 0.99886 * b(1) + white_noise(i) * 0.0555179;b(2) = 0.99332 * b(2) + white_noise(i) * 0.0750759;b(3) = 0.96900 * b(3) + white_noise(i) * 0.1538520;b(4) = 0.86650 * b(4) + white_noise(i) * 0.3104856;b(5) = 0.55000 * b(5) + white_noise(i) * 0.5329522;b(6) = -0.7616 * b(6) - white_noise(i) * 0.0168980;pink_noise(i) = b(1) + b(2) + b(3) + b(4) + b(5) + b(6) + b(7) + white_noise(i) * 0.5362;b(7) = white_noise(i) * 0.115926;
end%Normalize to +/- 1
if(abs(min(pink_noise)) > max(pink_noise))  pink_noise = pink_noise / abs(min(pink_noise));
elsepink_noise = pink_noise / max(pink_noise);
end%Normalize to prevent positive saturation (We can't represent +1.0)
pink_noise = pink_noise /abs(((2^31)-1)/(2^31));%Scale signal to match desired level
pink_noise = pink_noise * 10^(Amp/20);%Output noise signal
out = pink_noise(1:end-1);

这里使用了两种 Paul Kellet 的滤波方法降低误差，滤波系数是

Here are some new lower-order pink noise filter coefficients.These have approximately equiripple error in decibels from 20hz to 20khz at a 44.1khz
sampling rate.1st order, ~ +/- 3 dB error (not recommended!)
num = [0.05338071119116  -0.03752455712906]
den = [1.00000000000000  -0.97712493947102]2nd order, ~ +/- 0.9 dB error
num = [ 0.04957526213389  -0.06305581334498   0.01483220320740 ]
den = [ 1.00000000000000  -1.80116083982126   0.80257737639225 ]

最上面的使用更高阶次的实现滤波相关减小误差

粉红噪声matlab产生相关推荐

Matlab生成白噪和粉噪
音频测试中常用音源包括:白噪.粉噪(粉噪由于模拟人耳,常用来测量扬声器响度) 生成音源文件的方式:CoolEdit.Matlab function WhiteNoise = generateWhite ...
白噪声，有色噪声的定义、特性及其MATLAB仿真
一.白噪声白噪声(white noise)是指功率谱密度在整个频域内是常数的噪声. 所有频率具有相同能量密度的随机噪声称为白噪声.白噪声是指在较宽的频率范围内,各等带宽的频带所含的噪声功率谱密度相等 ...
Matlab应变片仿真,一种基于Matlab/Adams联合仿真的真实路谱再现系统和方法与流程...
本发明属于汽车系统动力学仿真技术领域,特别是一种基于Matlab/Adams联合仿真的真实路谱再现系统和方法. 背景技术: 汽车系统动力学仿真技术是汽车设计制造中一项不可或缺的技术,尤其是在汽车操纵稳 ...
matlab如何代码实现原理,色噪声原理及matlab代码实现
色噪声原理及matlab实现 1.实验目的: ⑴了解随机信号自身的特性,包括均值(数学期望).均方值.方差.相关函数.概率密度.频谱及功率谱密度等. (2)了解色噪声的基本概念和分析方法,掌握用mat ...
色噪声原理及matlab代码实现,色噪声原理及matlab代码实现
色噪声原理及matlab代码实现色噪声原理及 matlab 实现1.实验目的:⑴ 了解随机信号自身的特性,包括均值(数学期望) .均方值.方差.相关函数.概率密度.频谱及功率谱密度等.(2)了解色噪 ...
Audition频率分析的Matlab实现代码：对同一个音源，matlab计算结果与Audition相同
Audition是音频工程师广泛使用的音频分析处理软件.Audition的频率分析模块能够得到音频序列的频谱图,通过编写matlab,实现类似Audition的结果. 大致思路是:对序列进行分帧,加窗 ...
matlab 倍频,matlab倍频程分析
关于倍频程的概念声音信号的频率范围为 20Hz 到 20KHz,在声音信号频谱分析一般不需要对每个频率成分进行具体分析.为了方便起见,人们把 20Hz 到 20KHz 的声频...... 第四章 ...
matlab中daubechie小波,[转载]小波滤波器
滤波器组完美重构与小波快速算法. 前面的分析可以知道Vj相当于在j分辨率的逼近,Vj-1相当于j-1分辨率的逼近,这样Wj-1相当于两个分辨率逼近的差.在高分辨率下,我们可以用f在(2^j*t)的采样 ...
matlab加有色噪声,这么讲解高斯白噪声，还有色噪声，谁都会懂
原标题:这么讲解高斯白噪声,还有色噪声,谁都会懂百度百科上解释为 "高斯白噪声,幅度分布服从高斯分布,功率谱密度服从均匀分布",听起来有些晦涩难懂,下面结合例子通俗而详细地介绍一 ...

粉红噪声matlab产生

粉红噪声matlab产生相关推荐

最新文章

热门文章