%this program calculates the iir from the assembly program


clear all

a4=  0.6272;
a3= -0.6108;
a2=  0.1412;
a1= -0.0255;
b4=  4.0;
b3=  6.0;
b2=  4.0;
b1=  1.0;
scale=0.05078125;


x(1:30)=[
5000.0; 4333.0; 5465.0; 13556.0; 7423.0; -5000.0; -4333.0; -5465.0; -13556.0; -7423.0;
5000.0; 4333.0; 5465.0; 13556.0; 7423.0; -5000.0; -4333.0; -5465.0; -13556.0; -7423.0;
5000.0; 4333.0; 5465.0; 13556.0; 7423.0; -5000.0; -4333.0; -5465.0; -13556.0; -7423.0];
w(1:4)=0;

for i=1:30

  w(i+4)=x(i)*scale+w(i+3)*a4+w(i+2)*a3+w(i+1)*a2+w(i)*a1;
  y(i)=w(i+4)+      w(i+3)*b4+w(i+2)*b3+w(i+1)*b2+w(i)*b1;
end;

    fin = fopen('output_matlab.dat', 'wt');
        for i=1:30
            fprintf(fin, '%f\n', y(i));
        end;
    fclose(fin);