小小文 发表于 2011-1-2 11:51

1/2对信号进行1/2谱切片处理 下面是程序 麻烦高手帮我改一下谢谢

N=512
fs=100000
x=x0311170536(:,2);%信号
%估计三阶累积量
K=64;
M=fix(N/K);
%分段去均值
for i=0:k-1
    y=x(i*M+1:(i+1)*M);
    y=y-mean(y);
    xk(i+1,:)=y;
end
nfft=256;
for t=-(nfft/2-1):nfft/2
    c(t+nfft/2)=third_slice_cumulant(xk,t);
end
Y=fft(c,nfft);
P=Y.*conj(Y)/nfft;
plot(f,p(1:nfft/2);
xlabel('频率')
ylabel('幅值')

gogspe 发表于 2012-7-8 22:45

回复 1 # 小小文 的帖子

还在研究不
f=fs*(0:nfft/2-1)/nfft;少了这个,别外是报什么错呢

gogspe 发表于 2012-7-23 18:12

代吗不全,N呢
页: [1]
查看完整版本: 1/2对信号进行1/2谱切片处理 下面是程序 麻烦高手帮我改一下谢谢