% exa080301_bispecd.m ,for example 8.3.1 
% to test bispecd.m 

clear;
% 信號的產生
N=256;
n=0:N-1;
f1=0.15;
fai1=pi/2;
f2=0.25;
fai2=pi/3;
f3=0.3;
fai3=pi/4;
en=randn(size(n));
s1=cos(2*pi*f1*n+fai1);
s2=cos(2*pi*f2*n+fai2);
s3=cos(2*pi*f3*n+fai3);
x=s1+s2+s3+en;

% 雙譜估計平面圖
[bspec,waxis]=bispecd(x,16);
title('雙譜平估計面圖')

% 雙譜估計三維圖
bspec_abs=abs(bspec);
figure
mesh(bspec_abs)
title('雙譜估計三維圖')