filtering - NaN values in Android filter -

-

i use following function filter signal. code works order 4, when use higher order, 5, nan values in output. highlly appreciated

void filt(int ord, double a[], double b[], int np,  double x[], arraylist<double> array)     {         int i,j;        // y[0]=b[0]*x[0];         array.add(0, b[0]*x[0]);         (i=1;i<ord+1;i++)         {              array.add(i, 0.0);             (j=0;j<i+1;j++)                array.add(i, array.get(i)+b[j]*x[i-j]);             (j=0;j<i;j++)                 array.add(i, array.get(i)-a[j+1]*array.get(i-j-1));         }         (i=ord+1;i<np;i++)         {            array.add(i, 0.0);                 (j=0;j<ord+1;j++)                 array.add(i, array.get(i)+b[j]*x[i-j]);                  (j=0;j<ord;j++)                 array.add(i, array.get(i)-a[j+1]*array.get(i-j-1));         }      }

i tried same functionality of matlab, , got nan values. dig in matlab documentation , found following explains nan values:

"for higher order filters (possibly starting low order 8), numerical problems due roundoff errors may occur when forming transfer function using [b,a] syntax." therefor butterworth filter unstable higher orders.

to on problem in matlab, try following:

n = 6; wn = [2.5e6 29e6]/500e6; ftype = 'bandpass';  % transfer function design [b,a] = butter(n,wn,ftype); h1=dfilt.df2(b,a);      % unstable filter.  % zero-pole-gain design [z, p, k] = butter(n,wn,ftype); [sos,g]=zp2sos(z,p,k); h2=dfilt.df2sos(sos,g);  % plot , compare results hfvt=fvtool(h1,h2,'frequencyscale','log'); legend(hfvt,'tf design','zpk design')

Comments

Post a Comment

Popular posts from this blog

c# - Unity IoC Lifetime per HttpRequest for UserStore -

-
i'm trying clean default implementation of accountcontroller.cs comes out of box in new mvc5/owin security implementation. have modified constructor this: private usermanager<applicationuser> usermanager; public accountcontroller(usermanager<applicationuser> usermanager) { this.usermanager = usermanager; } also, have created lifetime manager unity looks this: public class httpcontextlifetimemanager<t> : lifetimemanager, idisposable { private httpcontextbase _context = null; public httpcontextlifetimemanager() { _context = new httpcontextwrapper(httpcontext.current); } public httpcontextlifetimemanager(httpcontextbase context) { if (context == null) throw new argumentnullexception("context"); _context = context; } public void dispose() { this.removevalue(); } public override object get...
Read more

Change the color of an oval at click in Java AWT -

-
i have draw ovals in java, , @ click change color. beginning tried change color after 20 ms, doesn't work. my code is: public class mycomponentnew extends frame { public graphics2d g2d; public mycomponentnew(string title) { super(title); setsize(400, 550); } @override public void paint(graphics g) { this.g2d = (graphics2d) g; this.g2d.setcolor(color.red); this.g2d.filloval(10, 55, 50, 100); } public void changecolor () { this.g2d.setcolor(color.blue); this.g2d.filloval(10, 55, 50, 100); } } and in class main method have: mycomponentnew m; m = new mycomponentnew("fereastra cu baloane"); m.setvisible(true); m.addwindowlistener(new windowadapter() { @override public void windowclosing(windowevent we) { system.exit(0); } }); try { thread.sleep(20); } catch(interruptedexception e) {} m.changecolor(); the color of oval remains red. ...
Read more

I am trying to solve the error message 'incompatible ranks 0 and 1 in assignment' in a fortran 95 program. -

-
i writing program in fortran find velocity of parachuting person in relation time. keep getting error can't fix. new programming , appreciated. the error is v(i+1)=v(i)+[32-((c*v(i)*v(i))/m)]*(h) 1 error: incompatible ranks 0 , 1 in assignment @ (1) and program is program para integer :: real :: v(11) !velocity real :: q !initial velocity real :: h !time step real :: c !drag coefficient real :: m !mass ! gravity equal 32 ft/s^2 write (*,*)'enter time step' read(*,*)h write(*,*)'enter initial velocity' read(*,*)q write(*,*)'enter drag coefficient' read(*,*)c write(*,*)'enter mass' read(*,*)m i=1,10 ! 1 10, 1 being interval. end v(i+1)=v(i)+[32-((c*v(i)*v(i))/m)]*(h) q=v(1) end program you cannot use [] normal parenthesis in expressions. array constructor, [ items ] means array items elements. end do should after line.
Read more