小弟初学matlab  今日师兄让偶识别一个信号的衰减。

clear

clc

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 一些参数 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

st=1024;                               % 采样频率 1024Hz

TT=300;                               % 采样周期

nn=st*TT;                            % 数据的采集长度=频率*时间

%注意：数据长度一般是1024的整数倍，当采样时间*频率不满足条件时，会延长采样时间自动满足1024的整数倍.

%例如采样时间为60s，数据长度应为3840，最后自动延长至4096,故采样时间为64S

nyquist=0.5*st;                      % 最高Nyquist频率=1/2采样频率

df=1/nn*st;                          % 最低起始频率

duan=40;                             % 分段数

kk=20;                                %采用随机减量技术时的周期数

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 测量数据读取 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

fid1=fopen('G:\matlab\work\try.uff','rt');

%fid2=fopen('G:\matlab\work\try.uff','rt');

for i=1:13

c=fgets(fid1); %行读取指令。

end

for j=1:nn/1024;       % 共300秒，200次1024采样，

b=fscanf(fid1,'%f',[1,1024]);  % 读取文件fid

Y(1024*(j-1)+1:1024*j)=b;

for i=1:15

c=fgets(fid1);   %  进行行读取，循环次，抵消说明文字

end

end

Y=Y-mean(Y);

//该段程序读取的是一个1024HZ*300秒的信号(一个样本)，uff文件除了前面13行是一些注释性的参数外，后面的就是信号的加速度数据。

用MATLAB不能直接算Y的平方和，所以小弟不知道如何计算Y的均方根值，请教各位前辈。指点下小弟，不甚感激！

(Y的时域图如下)

[本帖最后由 mp2625 于 2009-11-6 22:58 编辑]

1111.jpg

(32.63 KB, 下载次数: 0)

2009-11-6 22:58 上传