bispeci_bai.m

可以改变信噪比 分析在不同信噪比下 高斯白噪声中正弦信号的高阶谱情况 用于检测信号

close all;clc;clear;
fs=5000; %采样频率
t=0:1/fs:1.0-1/fs;N=1.0*fs; %采样点数 
Am=0.8; %输入信号幅度
s=Am*sin(2*pi*68*t);%输入信号,频率68Hz
A=1;
noi=A*randn(1,N);%高斯白噪声
%noi=wgn(1,N,10);
d=s+noi;
%d=awgn(s,12,'measured'); %信号的混合
SNR=10*log10(Am*Am); %信噪比
SNR

% subplot(311);plot(s);
% subplot(312);plot(noi);
% subplot(313);plot(d);

%间接法双谱估计
%bispeci返回用间接法从有限个观测信号中估计出的双谱bispec，并且用等高线图显示
figure;
subplot(221);bspec1=bispeci(noi,21,64,0,'unbiased',128,1);title('bispectrum of noi');

subplot(222);[bspec2,waxis2]=bispeci(d,21,64,0,'unbiased',128,1);title('bispectrum of d');

% subplot(223);mesh(abs(bspec1));title('noi');  %网格图
% subplot(224);mesh(abs(bspec2));title('d');

figure;
subplot(121);mesh(abs(bspec1));title('noi');  %网格图

subplot(122);mesh(abs(bspec2));title('d');
figure;
subplot(121);mesh(abs(bspec1));title('noi');
view(-37.5,0);
subplot(122);mesh(abs(bspec2));title('d');
view(-37.5,0);


%直接法双谱估计 
%bispecd返回估计双谱矩阵，原点在中心，轴的方向是向下与向右
figure;
subplot(221);bspec3=bispecd(noi,128,1,64,0);title('bispectrum of noi');
figure;
subplot(222);[bspec4,waxis4]=bispecd(d,128,1,64,0);title('bispectrum of d');
figure;
subplot(223);mesh(abs(bspec3));title('noi');  %网格图
subplot(224);mesh(abs(bspec4));title('d');