function[y]=twt_model(x,twtdata,ibo)
pwrin=twtdata(:,1);pwrout=twtdata(:,2);phaseout=twtdata(:,3);
n=length(x);ibo_real=10^(ibo/10);
instpwr=0.5*(abs(x).^2);
avepwr=mean(instpwr);
%instpwr=insper*(ibo_real/avepwr);
instpwr=insper*(ibo_real/avepwr);
outpwr=interp1(pwrin,pwrout,instpwr);
outmag=sqrt(2*outpwr);
outphase=(pi/180)*interp1(pwrin,phaseout,instpwr);
y=outmag.*exp(sqrt(-1)*(outphase+angle(x)));