[xiaji99]

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是一个数据记录文件 里面有很多数据