xiaji99
2012-03-27, 14:34
function[cfp]=chafunP(ud)
%P波相捡拾问题
%编写2006.11.27 马强
%函数[cfp]=chafun(ud,ew,ns)
%ud 垂直向记录
%ew 东西向记录
%ns 南北向记录
%cfp 特征函数
ntud=max(size(ud));
%一、只用垂直分量,目的为捡拾P波;
%特征函数为CFp=(X(k))^2+(X(k)-X(k-1))^2 (吴逸民,1998);
cfp=zeros(1,ntud+1);
cfp(1)=ud(1)^2;
ud(ntud+1)=ud(ntud); %增加一个递归点
for k=2:ntud+1
cfp(k)=ud(k)^2+(ud(k)-ud(k-1))^2;
end
cfp=cfp(1:ntud);
这是chafunP.m文件
我在命令窗输入结果是
>> chafunP
??? Input argument "ud" is undefined.
Error in ==> chafunP at 13
ntud=max(size(ud));
那怎么定义ud啊 ud是一个数据记录文件 里面有很多数据
%P波相捡拾问题
%编写2006.11.27 马强
%函数[cfp]=chafun(ud,ew,ns)
%ud 垂直向记录
%ew 东西向记录
%ns 南北向记录
%cfp 特征函数
ntud=max(size(ud));
%一、只用垂直分量,目的为捡拾P波;
%特征函数为CFp=(X(k))^2+(X(k)-X(k-1))^2 (吴逸民,1998);
cfp=zeros(1,ntud+1);
cfp(1)=ud(1)^2;
ud(ntud+1)=ud(ntud); %增加一个递归点
for k=2:ntud+1
cfp(k)=ud(k)^2+(ud(k)-ud(k-1))^2;
end
cfp=cfp(1:ntud);
这是chafunP.m文件
我在命令窗输入结果是
>> chafunP
??? Input argument "ud" is undefined.
Error in ==> chafunP at 13
ntud=max(size(ud));
那怎么定义ud啊 ud是一个数据记录文件 里面有很多数据