clear all, close all, clc
load ..\yamashita\stiction_model\all_FIC_loops_nominal.mat
x = OP_coll{5}(1:4096);y = MV_coll{5}(1:4096);
x = x - mean(x); y = y - mean(y);
e = SP_coll{5}(1:4096) - MV_coll{5}(1:4096);
e = e - mean(e);
bicoherence(e)
xlabel('frequency/sampling frequency')
ylabel('frequency/sampling frequency')
zlabel('bic^2');
pause
subplot(4,2,[1 3])
f = [0.0594 0.1750]
xf = wfilter(x,f);yf = wfilter(y,f);
plot(xf,yf,'o-')
axis square
xlabel('OP_f'),ylabel('MV_f'),title('filter = [0.0594 0.1750]')
subplot(4,2,5)
plot(xf),ylabel('OP_f')
subplot(4,2,7)
plot(yf),ylabel('MV_f')
subplot(4,2,[2 4])
% f = [0.004 0.1750]
f = [0.004 0.03]
xf = wfilter(x,f);yf = wfilter(y,f);
plot(xf,yf,'o-')
axis square
xlabel('OP_f'),ylabel('MV_f'),title(['filter = ' num2str(f)])
subplot(4,2,6)
plot(xf),ylabel('OP_f')
subplot(4,2,8)
plot(yf),ylabel('MV_f')
