bilinaer.m

%impulse response
freq=5;
phi=pi/3;
sample=20*freq;
dt=1/sample;
t=0:dt:10;
t=t';
X_ps=3*sin(2*pi*freq*t+phi);%primary signal noise
plot(t,X_ps);
%secondary path impulse response
A=2;
zeta=0.05;
Wn=20;
wd=sqrt(1-zeta^2);
phil=pi/3;
Y_sec_Ir=A*exp(-zeta*Wn*t);
figure(6)
plot(t,Y_sec_Ir);
title('Impulse Response of secondary source');
xlabel('t');
ylabel('pressure');
Xr=sin(2*pi*freq*t);
N=60;
mu1=7*power(10,-6.1);
mu2=7*power(10,-6.5);
mu3=7*power(10,-6.65);
M=length(t);
Aw=zeros(1,N);
Bw=zeros(1,N);
Cw=zeros(N,N);
x_buf=zeros(N,1);
y_buf=zeros(M,1);
y_buf_last=zeros(N,1);
y_buf_vector_last=zeros(1,N);
r_buf=zeros(1,N);
r_buf1=zeros(1,N);
r_1=zeros(M,1);
r_2=zeros(M,1);
E_t=zeros();
y_signal=zeros(M,1);
v_n=zeros(M,1);
v_n1=zeros(M,1);
v_buf=zeros(1,N);
x_buf1=zeros(N,1);
control_signal=zeros();
T2=zeros();
for n=1:M
         x_buf=[Xr(n);x_buf(1:end-1)];%storing refernce signal
   
         r_1=[Xr(n);r_1(1:M-1)];
         r=sum(r_1.*Y_sec_Ir);%convultion with secondary impulse response
         r_buf=[r,r_buf(1:end-1)];
         y_sig=sum(Aw*x_buf(1:N))+sum(Bw*y_buf_last)+sum(y_buf_last.*(Cw*x_buf(1:N)));
         %disp(y_sig);
         y_signal=[y_sig;y_signal(1:end-1)];
         for i=0:N-1
             for j=1:N
                 if(n>N)
                     v_n(n)=Xr(n-i).*y_buf(n-j);
                 end
             end
         end
         v_n1=[v_n(n);v_n1(1:M-1)];
         v=sum(r_1.*Y_sec_Ir);
         v_buf=[v,v_buf(1:end-1)];


         %sz1=size(y_sig);
         %disp(sz1);
         y_buf=[y_sig;y_buf(1:M-1)];
         r1=sum(y_buf.*Y_sec_Ir);
    %disp(r1);
         r_buf1=[r1,r_buf1(1:end-1)];
         e=X_ps(n)+r1;


    %sz=size(Xr(n));
    %disp(sz);
   
    %sz1=size(r_1);
    %disp(sz1);
    
    %sz=size(u);
    %disp(sz);
    
    %y_buf=[y_signal;y_buf(1:N)];
    %sz1=size(y_buf);
    %disp(sz1);

    Aw=Aw-mu1*e*r_buf;
    Bw=Bw-mu2*e*y_buf_vector_last;
    Cw=Cw-mu3*e*v_buf;
    y_buf_last=y_buf(1:N);
    y_buf_vector_last=r_buf1(1:N);
    E_t(n)=e;
    %sz1=size(r_buf1);
    %disp(r_buf1);
    %r1=sum(r_2*Y_sec_Ir);
    figure(3)
    plot(y_signal,'*--');
    n_iter=n;
    control_signal(n)=y_sig;
    T2(n)=t(n);
    


end
figure(2)
%subplot(2,1,1)
plot(t,X_ps,'b',T2,E_t,'r');
title(  'Waveform Obtained by Bilinear' ) ;
xlabel ( ' time [ second ] ' ) ;
ylabel ( 'Pressure ' ) ;