WSJT-X/lib/fchisq0.f90

real function fchisq0(y,npts,a)

  real y(npts),a(4)
  
!  rewind 51
  chisq = 0.
  do i=1,npts
     x=i
     z=(x-a(3))/(0.5*a(4))
     yfit=a(1)
     if(abs(z).lt.3.0) then
        d=1.0 + z*z
        yfit=a(1) + a(2) * (1.0/d - 0.1)
     endif
     chisq=chisq + (y(i) - yfit)**2
!     write(51,3001) i,y(i),yfit,y(i)-yfit
!3001 format(i5,3f10.4)
  enddo
  fchisq0=chisq

  return
end function fchisq0
Various changes to JT65 decoding, all potentially temporary. 1. Measure Doppler width by fitting a (modified) Lorentzian. 2. Don't call "slope" in sync65(). 3. New definition of "sync1". 4. Get snr from sync1. git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/wsjtx@6536 ab8295b8-cf94-4d9e-aec4-7959e3be5d79 2016-03-17 14:52:06 -04:00			`real function fchisq0(y,npts,a)`

			`real y(npts),a(4)`

Add submode to argument list of JT9 subroutine decode(). Add test program fersum; several tweaks to fer65. Code cleanup in lorentzian, and return chisqr in a(5). git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/wsjtx@6567 ab8295b8-cf94-4d9e-aec4-7959e3be5d79 2016-04-01 13:57:08 -04:00			`! rewind 51`
Various changes to JT65 decoding, all potentially temporary. 1. Measure Doppler width by fitting a (modified) Lorentzian. 2. Don't call "slope" in sync65(). 3. New definition of "sync1". 4. Get snr from sync1. git-svn-id: svn+ssh://svn.code.sf.net/p/wsjt/wsjt/branches/wsjtx@6536 ab8295b8-cf94-4d9e-aec4-7959e3be5d79 2016-03-17 14:52:06 -04:00			`chisq = 0.`
			`do i=1,npts`
			`x=i`
			`z=(x-a(3))/(0.5*a(4))`
			`yfit=a(1)`
			`if(abs(z).lt.3.0) then`
			`d=1.0 + z*z`
			`yfit=a(1) + a(2) * (1.0/d - 0.1)`
			`endif`
			`chisq=chisq + (y(i) - yfit)**2`
			`! write(51,3001) i,y(i),yfit,y(i)-yfit`
			`!3001 format(i5,3f10.4)`
			`enddo`
			`fchisq0=chisq`

			`return`
			`end function fchisq0`