|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
求助:求取时域下的特征值,出现错误,有请诸位i帮忙
如题:我用一下一段程序,求取振动信号(用加速度传感器采集)时域下的六个特征值,出现错误,还请大家帮我看看这段程序
帮忙修改下,小弟感激不尽!!!
就此先道谢一声!:handshake :handshake
程序:
clear all;
load zc01.dat;
x01=zc01;
load zc02.dat;
x02=zc02;
load zc03.dat;
x03=zc03;
load zc04.dat;
x04=zc04;
load zc05.dat;
x05=zc05;
load zc06.dat;
x06=zc06;
load zc07.dat;
x07=zc07;
load zc08.dat;
x08=zc08;
load zc09.dat;
x09=zc09;
load zc10.dat;
x10=zc10;
load zc11.dat;
x11=zc11;
load zc12.dat;
x12=zc12;
X=[x01 x02 x03 x04 x05 x06 x07 x08 x09 x10 x11 x12]';
%求样本时域特征
function Y=get_sample(X,N,4096);%X样本矩阵;N欲求样本数;数据点数4096
i=1;
for i=1:N
sample(:,i)=X(N*(i-1)+1+4096):N*i+4096,:);%获得样本
max_eg=max(sample(:,i));%求出峰值
min_eg=abs(min(sample(:,i)));
if(min_eg>max_eg)
max_eg=min_eg;
end
Y(1,i)=max_eg;
Y(2,i)=var(sample(:,i))^0.5;%球均方根
Y(3,i)=skewness(sample(:,i));%计算偏度
Y(4,i)=kurtosis(sample(:,i))-3;%计算峭度
Y(5,i)=mean(sample(:,i));%均值
Y(6,i)=var(sample(:,i),1);%均方差
end
运行显示错误:
??? function Y=get_sample(X,N,4096);%X样本矩阵;N欲求样本数;数据点数4096
|
Error: Function definitions are not permitted at the prompt or in scripts. |
|