|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
%%%%%关于用互相关函数进行互谱分析
clear
load 'wind.txt'; %导入数据文件
M=2;
N=3;
%变换此处的M和N的值可以求不同点的互谱
X=wind(:,M);
Y=wind(:,N);
[Pxy,F] = cpsd(X,Y,[],[],512,6.25);
其中cpsd函数是进行互谱分析的,但其所用的方法是经典谱估计的welch方法,其谱线起伏度还是很大,我查了很多资料,好像现在进行的互谱分析用的都是cpsd,我想能不能用现代谱估计(例如burg)方法代替,程序如下:
M=2;
N=3;
%变换此处的M和N的值可以求不同点的互谱
X=wind(:,M);
Y=wind(:,N);
cxy = cohere(X,Y,256,6.25);
[Pxx1,F]=pburg(X,6,256,6.25);
[Pxx2,F]=pburg(Y,6,256,6.25);
for n=1:129
[Pxy,F]=(Pxx1(n)*Pxx2(n))*cxy(n);
end
其中红色部分就是改后的,先进行自谱估计,再根据相关函数与互谱、自谱的关系(用相关函数乘以两个信号的自谱)。
不知这种方法可行?望高手指教。
[ 本帖最后由 xulinzero 于 2008-10-10 17:06 编辑 ] |
|