【数字信号去噪】基于matlab奇异值分解（SVD）数字信号降噪【含Matlab源码 1020期】.zip

% 声波信号降噪处理 function [TempData,Sample] =ReadData(FileName) %数据读取，数据格式为特检院检测仪存储格式 % FileName='SR10-100021-SPL-5-4.dat'; fid=fopen(FileName,'rb'); %读取文件头 % 采样时间结构 Sample_Time = struct('year',{2009},'month',{1}, 'day',{22},'hour',{12},'minute',{12},'second',{12}) ;%定义采样时间格式 SensorAddr=fread(fid,1,'int32');%传感器地址 fseek(fid,4,'bof');%数据文件头长度为44字节 Sample_Time.year=fread(fid,1,'int32'); Sample_Time.month=fread(fid,1,'int32'); Sample_Time.day=fread(fid,1,'int32'); Sample_Time.hour=fread(fid,1,'int32'); Sample_Time.minute=fread(fid,1,'int32'); Sample_Time.second=fread(fid,1,'int32'); GroupOrder=fread(fid,1,'int8'); %读取的批次，0xff表示实时采样数据 startPack=fread(fid,1,'int8'); %起始包号，0～255 totalPack=fread(fid,1,'int8'); %总包数，0～255，0表示256包 location=fread(fid,1,'int8'); %/GPS定位信息 %设备配置信息结构 Device_Config Device_Config.group=fread(fid,1,'int8'); %批号 Device_Config.year=fread(fid,1,'int8');%年 Device_Config.month=fread(fid,1,'int8');%月 Device_Config.day=fread(fid,1,'int8');%日 Device_Config.hour=fread(fid,1,'int8');%时 Device_Config.minute=fread(fid,1,'int8');%分 Device_Config.second=fread(fid,1,'int8');%秒 Device_Config.repeat=fread(fid,1,'int8');%次数 Device_Config.interval=fread(fid,1,'int8');%间隔 Device_Config.sampleLen=fread(fid,1,'int8');%采样长度 Device_Config.speed=fread(fid,1,'int8');%采样速度 Device_Config.gain=fread(fid,1,'int8');%增益 fseek(fid,44,'bof');%数据文件头长度为44字节 Count=261120;%数据长度 [TempData,Count]=fread(fid,'int16'); TempData=TempData*2.5/2048; %得到数据长度 switch Device_Config.sampleLen case {8} Sample.Len='512KB'; case {9} Sample.Len='512KB'; case {10} Sample.Len='512B'; case {11} Sample.Len='1024B'; end % //采样速度 switch Device_Config.speed case {1} Sample.Frequency = '500k'; case {2} Sample.Frequency = '250k'; case {3} Sample.Frequency = '125k'; case {4} Sample.Frequency = '62.5k'; end % //增益 switch Device_Config.gain case {0} Sample.Gain = 1; case {1} Sample.Gain = 4; case {2} Sample.Gain = 16; case {3} Sample.Gain = 64; case {4} Sample.Gain = 10; case {5} Sample.Gain = 40; case {6} Sample.Gain = 160; case {7} Sample.Gain = 640; case {8} Sample.Gain = 100; case {9} Sample.Gain = 400; case {10} Sample.Gain = 1600; case {11} Sample.Gain = 6400; end