1/2对信号进行1/2谱切片处理 下面是程序 麻烦高手帮我改一下谢谢
N=512fs=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('幅值') 回复 1 # 小小文 的帖子
还在研究不
f=fs*(0:nfft/2-1)/nfft;少了这个,别外是报什么错呢 代吗不全,N呢
页:
[1]