本帖最后由 wdhd 于 2016-9-19 13:08 编辑
二次相位耦合有人研究过吗?
我的程序不知道毛病在哪,请高手指教.二次相位一种是耦合,一种是非耦合, 我只给我了.耦合的情况,发现两种情况结果相同,不知道为什么?
clear;
%产生数据
N=128*8;
fs=100;
f=0:fs/N:fs/N*(N-1);
n=0:N-1;
f1=10;
fai1=2*pi*randn;
f2=30;
fai2=2*pi*randn;
f3=40;
fai3=fai1+fai2;
e=randn(1,N);
t=[0:1/fs:(N-1)/fs];
s1=cos(2*pi*f1*t+fai1);
s2=cos(2*pi*f2*t+fai2);
s3=cos(2*pi*f3*t+fai3);
%x=s1+s2+0.5*s3+s1.*s2+e;
x=s1+s2+s3+e;
X=fft(x,N);
X=abs(X);
X=20*log10(X);
figure(1);
subplot(211);
plot(f(1:N/2),X(1:N/2));
Y=zeros;
u=rand(1,64*N);u=u-mean(u);
for i=0:63
x=s1+s2+s3+10*u(1+i*N:i*N+N);
X=fft(x,N);
X=abs(X);
Y=Y+X;
end
Y=Y/64;
Y=20*log10(Y);
subplot(212);plot(f(1:N/2),Y(1:N/2));
xlabel('Hz')
%双谱估计平画图
figure(2);
x=s1+s2+s3+e;
[bspec,waxis]=bispecd(x,16);
title('双谱估计平画图')
figure(3);
bspec_abs=abs(bspec);
mesh(bspec_abs);
title('双谱估计三维图')
figure(4);
[bic,waxis] = bicoher (x,16);
mesh(bic)
title('双谱双相干图')
|