从IEPE传声器获取和分析声压数据
本示例演示如何从130E20型PCB®IEPE阵列麦克风获取和显示声压数据。该传感器正在记录音叉在中C (261.626 Hz)频率下产生的声压。
发现支持麦克风的设备万博1manbetx
方法返回的表中访问该设备,可以发现支持麦克风测量的设万博1manbetx备daqlist
命令。在本例中,麦克风连接在CompactDAQ Chassis NI cDAQ-9178上的National Instruments®设备NI 9234的0通道上,并带有IDcDAQ1Mod3
.
D = daqlist(“倪”)
d = 12×4表的DeviceID DeviceInfo描述模型 ___________ __________________________________ _____________ ____________________ " cDAQ1Mod1”“国家仪器NI 9205”“NI 9205”[1×1采集。DeviceInfo] "cDAQ1Mod2" "美国国家仪器NI 9263" "NI 9263" [1×1 daq。DeviceInfo] "cDAQ1Mod3" "美国国家仪器NI 9234" "NI 9234" [1×1 daq。DeviceInfo] "cDAQ1Mod4" "美国国家仪器NI 9201" "NI 9201" [1×1 daq。DeviceInfo] "cDAQ1Mod5" "国家仪器NI 9402" "NI 9402" [1×1 daq。DeviceInfo] "cDAQ1Mod6" "美国国家仪器NI 9213" "NI 9213" [1×1 daq。DeviceInfo] "cDAQ1Mod7" "美国国家仪器NI 9219" "NI 9219" [1×1 daq。DeviceInfo] "cDAQ1Mod8" "美国国家仪器NI 9265" "NI 9265" [1×1 daq。DeviceInfo] "Dev1" "美国国家仪器公司pci -6363" " pci -6363" [1×1 daq。DeviceInfo]"Dev2" "National Instruments NI ELVIS II" "NI ELVIS II" [1×1 daq.DeviceInfo] "Dev3" "National Instruments PCIe-6363" "PCIe-6363" [1×1 daq.DeviceInfo] "Dev4" "National Instruments PCIe-6363" "PCIe-6363" [1×1 daq.DeviceInfo]
deviceInfo = d{3,“DeviceInfo”}
deviceInfo = ni: National Instruments ni 9234(设备ID: 'cDAQ1Mod3')模拟输入支持:-5.0到+5.0伏范围速万博1manbetx率从1000.0到5200.0扫描/秒4通道('ai0','ai1','ai2','ai3')“电压”,“加速度计”,“麦克风”,“IEPE”测量类型该模块位于“cDAQ-9178”机箱的插槽3,名称为“cDAQ1”。
创建数据采集并添加麦克风通道
创建数据采集并添加通道麦克风
测量类型。
Dq = daq(“倪”);Ch = addinput(dq,“cDAQ1Mod3”,“ai0”,“麦克风”);
设置传感器属性
设置麦克风通道灵敏度
属性设置为传感器数据表中指定的值。对于该传感器,灵敏度
值为0.037伏特/帕斯卡。检查通道属性以查看设备配置中的变化。
ch.灵敏度= 0.037;ch
ch =指数类型设备通道测量类型范围名称 _____ ____ ___________ _______ ___________________ ______________________ _______________ 1“人工智能”“cDAQ1Mod3”“ai0”“麦克风(Diff)”"-200到+200帕斯卡" "cDAQ1Mod3_ai0"
配置并启动采集
将采集扫描速率设置为每秒51200次扫描,然后使用读
来获取四秒钟的数据。
dq。费率= 51200;Tt = read(dq, seconds(4));t = tt.时间;data = tt.cDAQ1Mod3_ai0;
分析数据
使用audioplayer
播放获取的麦克风信号
p = audioplayer(data, dq.Rate);玩(p);
在时间域内检查数据
情节(t)数据);ylabel (声压(帕斯卡));
频域数据检验
使用fft
并通过以下参数计算输入数据的单侧振幅谱:
计算信号长度(表中样本或条目的数量)
计算nfft
计算振幅和频率
Len =高度(tt);NFFT = 2^nextpow2(len);Y = fft(data,nfft)/len;f = dq.Rate/2*linspace(0,1,nfft/2+1);A = 2*abs(y(1:nfft/2+1));
绘制单面振幅谱
情节(f);xlim (1000 [0]);包含(的频率(赫兹));ylabel (“振幅”);
图中显示了261.626赫兹的峰值。这与音叉的频率相匹配。