红外光谱信噪比
翁老爷子的新书《傅里叶变换红外光谱分析》(第2版)中,有一段对红外仪器信噪比的无奈描述:
“红外仪器的信噪比是衡量一台仪器性能好坏的一项非常重要的技术指标。但是信噪比的测量方法目前没有统一的、公认的标准,因此,各个红外仪器公司所给定的仪器信噪比没有可比性。每个红外仪器公司都有信噪比的测量方法,因此,信噪比指标的验收只能按照仪器公司的验收方法进行验收。”
看来这个“红外信噪比”真个是乱花渐欲迷人眼,让人雾里看花隔一层啊!但是zwyu我充分发挥超人的大打特打、死缠烂打、穷追猛打的精神,欲对“红外信噪比”进行一次非官方、全方位的刨根问底,追踪探秘。各位好奇同学请跟进!
正文
信噪比(signal-to-noise ratio,简记为SNR ),顾名思义,就是信号值与噪声值的比,这一比值当然是越高越好。可是,翻遍《GB/T 21186-2007 傅立叶变换红外光谱仪》,《GB/T 6040-傅里叶变换公式证明
2002 红外光谱分析方法通则》(见红外光谱相关标准与检定规程大合集)以及其他的一些行业性、地方性的检定规程(国家级的傅里叶变换红外光谱仪检定规程至今还未出台),甚至中国药典,愣是不到关于信噪比的只言片语的定义。信噪比指标对红外仪器性能的评判很重要,怎么会不呢?且慢,注意标准中屡屡提到的“基线噪声”(100%T线噪声)XXXX:1或1:XXXX,还往往标了P-P或RMS,这不就是我们熟悉的信噪比的表示方法吗?哈哈,总算到你了。
艰难的看过标准上的描述(没办法,中国国标写的水平就是高!?),为了各位同学能够顺利读懂,我将它写为白话现代汉语版:
红外信噪比,是通过基线(100%T线)噪声来表征。也就是,在样品室中不放样品的情况下(空光路),测得一条假定理想的100%T透射光谱。信号,当然就是100%T了,如果没有噪声,那么这条光谱将是一条严格的纵坐标为100%T的直线,但是,实际情况是噪声总是存在的,这就使得这条光谱的各个波数点上的值不见得一定是100%T,可能高一些(比如100.1%T),也可能低一些(比如99.9%T)。P-P(峰-峰值)噪声的意思就是说刚才测得的那条光谱在某一段波数区间内(比如2200~2100cm-1)的最大值与最小值之差,比如说是
100.1%T-99.9%T=0.2%T。前面说了,信号是假定为100%T,那么,根据信噪比的定义,信号值/噪声值,比如100%T/0.2%T=500(注意此处单位相消,也就是说,信噪比用信号噪声比值表示的话,是一个无量纲的数)。此时,我们可以说,这台红外光谱仪的信噪比是500:1。换句话说,我们知道了P-P(峰-峰值)噪声,我们也就自然知道了P-P值信噪比;同理,我们知道了P-P值信噪比,比如500:1,那么我们很自然的也能利用噪声=信号/信噪比,即100%T/500=0.2%T,得到P-P噪声值的大小为0.2%T。
有人说,为了避免小概率事件的发生(此君是迷,鉴定完毕!),噪声值应该用更具代表性和统计性的RMS(均方根值)噪声来表示。那啥是RMS呢?我不得不祭出万恶的数学公式(霍金一部《时间简史》,只用了一个公式。我这个小小的原创这么早就出公式了。我不如霍金。。。)
设{Y1, Y2, Y3, …YN}为给定波数区间内N个连续波数点对应的纵坐标值(按照前述条件下,为一系列%T透过率值),则这些值的均值为:
均方根(root mean square,简记为RMS)偏差为:
如果不用公式,通俗地讲,均方根值就是一组数的平方的平均值的平方根;均方根偏差就是一组数与这组数均值之差的平方的平均值的平方根。所以,你瞧,我早早放弃了只用文字叙述,还是看看万恶的公式吧。显然,用上式求得的一条光谱在某波数(横坐标)区间内全部N个数据点纵坐标值的均方根偏差就作为了RMS噪声的度量。
一般对红外光谱来讲,P-P(峰-峰值)噪声会比RMS(均方根值)噪声大5倍左右,换句话说,RMS噪声的绝对数值更小,换算成信噪比时就更大,所以你发现用RMS值表示的信噪比往往看起来都很漂亮也就不奇怪了,因为它比P-P值表示的信噪比大了5倍(而且,显然参与运算的波数点越多,这一倍数还会增加)。
上面的“基线噪声”都是用了100%T基线,对应的是透射光谱的透过率表示形式;国际上越来越多的地方采用透射光谱的吸光度表示形式,此时的“基线”自然变成了0A基线。该“零基线”上的噪声单位,显然也就变成了A(吸光度单位,有时写做AU)。此时,计算P-P噪声和RMS噪声的方法与前面完全一样。但是,因为这些基线都是在样品室中不放样品的情况下(空光路)测得的,所以此时的信号应该是0A,如果直接计算信噪比的话,0/噪声=0,显然得不到明确的有意义数值。所以有很多同学这个地方就会糊涂了,由吸光度表示的基线噪声值,怎么得到信噪比?在此,zwyu独家奉献推导过程(呵呵,反正市面上所有的资料里都没写,可能觉得太简单了吧。):
版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198,我们将在24小时内删除。
发表评论