资源描述：
《峰值指标提取程序.doc》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、一、峰值指标提取程序%计算信号的峰值指标，输入信号为1024维。functiony=fengzhizhibiao(x)maxs=max(abs(x));jfg=sqrt(sum(x.^2)/1024);y=maxs/jfg;end二．newrbe函数的程序1.参数为峰值指标和峭度时%%初始化，清零closeall;clear;clc;formatcompact;%%载入训练样本数据文件[filename,pathname]=uigetfile(...{'*.ast','AST-files(*.AST)';...'*.*','AllFiles(
2、*.*)'},...'选择所有训练样本数据',...'MultiSelect','on');num_test=length(filename);fork=1:num_testdata1=load([pathname,filename{k}]);fzzb(k)=fengzhizhibiao(data1);qd(k)=kurtosis(data1);a=strncmp('nq1.ast',filename{k},1);b=strncmp('wq1.ast',filename{k},1);ifa==1t(:,k)=[001];elseifb==1t
3、(:,k)=[010];elset(:,k)=[100];endendx=[fzzb;qd];x2=x';%%%%载入测试样本数据文件[filename,pathname]=uigetfile(...{'*.ast','AST-files(*.AST)';...'*.*','AllFiles(*.*)'},...'选择所有测试样本数据',...'MultiSelect','on');num_test=length(filename);fork=1:num_testdata1=load([pathname,filename{k}]);fzzb1
4、(k)=fengzhizhibiao(data1);qd1(k)=kurtosis(data1);endx1=[fzzb1;qd1];x3=x1';%%spread=1;rbfnet=newrbe(x,t,spread);%%ty=sim(rbfnet,x1);t1=t't2=ty'2.参数为峰值指标时%%初始化，清零closeall;clear;clc;formatcompact;%%载入训练样本数据文件[filename,pathname]=uigetfile(...{'*.ast','AST-files(*.AST)';...'*.*'
5、,'AllFiles(*.*)'},...'选择所有训练样本数据',...'MultiSelect','on');num_test=length(filename);fork=1:num_testdata1=load([pathname,filename{k}]);fzzb(k)=fengzhizhibiao(data1);a=strncmp('nq1.ast',filename{k},1);b=strncmp('wq1.ast',filename{k},1);ifa==1t(:,k)=[001];elseifb==1t(:,k)=[010]
6、;elset(:,k)=[100];endendx=[fzzb];x2=x';%%%%载入测试样本数据文件[filename,pathname]=uigetfile(...{'*.ast','AST-files(*.AST)';...'*.*','AllFiles(*.*)'},...'选择所有测试样本数据',...'MultiSelect','on');num_test=length(filename);fork=1:num_testdata1=load([pathname,filename{k}]);fzzb1(k)=fengzhizhi
7、biao(data1);endx1=[fzzb1];x3=x1';%%spread=1;rbfnet=newrbe(x,t,spread);%%ty=sim(rbfnet,x1);t1=t't2=ty'3.参数为峭度时%%初始化，清零closeall;clear;clc;formatcompact;%%载入训练样本数据文件[filename,pathname]=uigetfile(...{'*.ast','AST-files(*.AST)';...'*.*','AllFiles(*.*)'},...'选择所有训练样本数据',...'Multi
8、Select','on');num_test=length(filename);fork=1:num_testdata1=load([pathname,filena
