从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赫兹的峰值。这与音叉的频率相匹配。