.MCAD 303010000 1 87 30 0 .CMD PLOTFORMAT 0 0 1 1 0 0 1 0 0 1 1 0 0 1 0 1 0 0 NO-TRACE-STRING 0 2 1 0 NO-TRACE-STRING 0 3 2 0 NO-TRACE-STRING 0 4 3 0 NO-TRACE-STRING 0 1 4 0 NO-TRACE-STRING 0 2 5 0 NO-TRACE-STRING 0 3 6 0 NO-TRACE-STRING 0 4 0 0 NO-TRACE-STRING 0 1 1 0 NO-TRACE-STRING 0 2 2 0 NO-TRACE-STRING 0 3 3 0 NO-TRACE-STRING 0 4 4 0 NO-TRACE-STRING 0 1 5 0 NO-TRACE-STRING 0 2 6 0 NO-TRACE-STRING 0 3 0 0 NO-TRACE-STRING 0 4 1 0 NO-TRACE-STRING 0 1 21 15 0 3 .CMD FORMAT rd=d ct=10 im=i et=6 zt=15 pr=3 mass length time charge temperature tr=0 vm=540024880 .CMD SET ORIGIN 0 .CMD SET TOL 0.001000000000000 .CMD SET PRNCOLWIDTH 8 .CMD SET PRNPRECISION 4 .CMD PRINT_SETUP 0.800000 0.500000 0.500000 0.500000 0 .CMD HEADER_FOOTER 1 1 *empty* *empty* *empty* 0 1 *empty* *empty* *empty* .CMD HEADER_FOOTER_FONT fontID=14 family=Arial points=10 bold=0 italic=0 underline=0 .CMD HEADER_FOOTER_FONT fontID=15 family=Arial points=10 bold=0 italic=0 underline=0 .CMD DEFINE_FONTSTYLE_NAME fontID=0 name=Variables .CMD DEFINE_FONTSTYLE_NAME fontID=1 name=Constants .CMD DEFINE_FONTSTYLE_NAME fontID=2 name=Text .CMD DEFINE_FONTSTYLE_NAME fontID=4 name=User^1 .CMD DEFINE_FONTSTYLE_NAME fontID=5 name=User^2 .CMD DEFINE_FONTSTYLE_NAME fontID=6 name=User^3 .CMD DEFINE_FONTSTYLE_NAME fontID=7 name=User^4 .CMD DEFINE_FONTSTYLE_NAME fontID=8 name=User^5 .CMD DEFINE_FONTSTYLE_NAME fontID=9 name=User^6 .CMD DEFINE_FONTSTYLE_NAME fontID=10 name=User^7 .CMD DEFINE_FONTSTYLE fontID=0 family=Times^New^Roman points=10 bold=0 italic=0 underline=0 .CMD DEFINE_FONTSTYLE fontID=1 family=Times^New^Roman points=10 bold=0 italic=0 underline=0 .CMD DEFINE_FONTSTYLE fontID=2 family=Arial points=10 bold=0 italic=0 underline=0 .CMD DEFINE_FONTSTYLE fontID=4 family=Arial points=10 bold=0 italic=0 underline=0 .CMD DEFINE_FONTSTYLE fontID=5 family=Courier^New points=10 bold=0 italic=0 underline=0 .CMD DEFINE_FONTSTYLE fontID=6 family=System points=10 bold=0 italic=0 underline=0 .CMD DEFINE_FONTSTYLE fontID=7 family=Script points=10 bold=0 italic=0 underline=0 .CMD DEFINE_FONTSTYLE fontID=8 family=Roman points=10 bold=0 italic=0 underline=0 .CMD DEFINE_FONTSTYLE fontID=9 family=Modern points=10 bold=0 italic=0 underline=0 .CMD DEFINE_FONTSTYLE fontID=10 family=Times^New^Roman points=10 bold=0 italic=0 underline=0 .CMD UNITS U=1 .CMD DIMENSIONS_ANALYSIS 0 0 .TXT 2 26 1 0 Cg a37.500000,48.000000,55 {\rtf1\ansi \deff0 {\fonttbl {\f0\fnil Arial;} {\f1\fnil Symbol;} } {\plain Signal{ Restoration with the Parametric Wiener Filter}} } .TXT 2 10 16 0 Cg a17.500000,54.000000,25 {\rtf1\ansi \deff0 {\fonttbl {\f0\fnil Arial;} {\f1\fnil Symbol;} } {\plain {Ken Castleman, 3/24/94}} } .TXT 0 32 2 0 Cg a10.750000,21.000000,13 {\rtf1\ansi \deff0 {\fonttbl {\f0\fnil Arial;} {\f1\fnil Symbol;} } {\plain PWF_1D{.MCD}} } .TXT 2 -30 17 0 Cg a10.500000,58.000000,17 {\rtf1\ansi \deff0 {\fonttbl {\f0\fnil Arial;} {\f1\fnil Symbol;} } {\plain {Section 16.2.4}} } .EQN 7 -35 3 0 {0:Np}NAME:20 .EQN 0 10 4 0 {0:u}NAME:0,(0.5)/({0:Np}NAME);0.5 .TXT 0 34 19 0 Cg a10.500000,47.000000,16 {\rtf1\ansi \deff0 {\fonttbl {\f0\fnil Arial;} {\f1\fnil Symbol;} } {\plain {The Gaussian:}} } .EQN 0 12 12 0 {0:G}NAME({0:\m}NAME,{0:\s}NAME,{0:x}NAME)~{0:if}NAME(|({0:x}NAME-{0:\m}NAME)<3*{0:\s}NAME,({0:e}NAME)^(-((((({0:x}NAME-{0:\m}NAME))^(2))/((2*{0:\s}NAME*{0:\s}NAME))))),0) .TXT 4 -58 23 0 Cg b86.500000,86.500000,94 {\rtf1\ansi \deff0 {\fonttbl {\f0\fnil Arial;} {\f1\fnil Symbol;} } {\plain The signal, f(x), has been smoothed by H(u) and corrupted by additive white noise of power N.} } .EQN 4 2 5 0 {0:\m}NAME~0 .EQN 0 7 10 0 {0:\s}NAME~0.2 .EQN 0 8 26 0 {0:N}NAME~0.1 .EQN 0 8 28 0 {0:k}NAME~100 .TXT 4 -24 24 0 Cg a14.875000,85.500000,19 {\rtf1\ansi \deff0 {\fonttbl {\f0\fnil Arial;} {\f1\fnil Symbol;} } {\plain The smoothing MTF:} } .EQN 0 17 6 0 {0:H}NAME({0:u}NAME):{0:G}NAME({0:\m}NAME,{0:\s}NAME,{0:u}NAME) .TXT 0 23 18 0 Cg a13.750000,87.000000,21 {\rtf1\ansi \deff0 {\fonttbl {\f0\fnil Arial;} {\f1\fnil Symbol;} } {\plain {The power spectra:}} } .EQN 0 16 7 0 {0:Pf}NAME({0:u}NAME):(2)/(1+{0:k}NAME*({0:u}NAME)^(2)) .EQN 0 18 8 0 {0:Pn}NAME({0:u}NAME):{0:N}NAME .TXT 5 -73 21 0 Cg a15.125000,87.000000,25 {\rtf1\ansi \deff0 {\fonttbl {\f0\fnil Arial;} {\f1\fnil Symbol;} } {\plain {The transfer function:}} } .EQN 0 20 9 0 {0:K}NAME({0:u}NAME,{0:\g}NAME):({0:H}NAME({0:u}NAME))/((({0:H}NAME({0:u}NAME)))^(2)+{0:\g}NAME*(({0:Pn}NAME({0:u}NAME))/({0:Pf}NAME({0:u}NAME)))) .EQN 5 -3 11 0 &&(_n_u_l_l_&_n_u_l_l_)&{0:K}NAME({0:u}NAME,0.2),{0:H}NAME({0:u}NAME),{0:Pn}NAME({0:u}NAME),{0:Pf}NAME({0:u}NAME)@&&(_n_u_l_l_&_n_u_l_l_)&{0:u}NAME 0 0 1 0 5 0 1 0 0 1 1 0 0 1 0 1 0 0 NO-TRACE-STRING 0 1 1 0 NO-TRACE-STRING 0 1 2 0 NO-TRACE-STRING 0 1 3 0 NO-TRACE-STRING 0 1 4 0 NO-TRACE-STRING 0 2 5 0 NO-TRACE-STRING 0 3 6 0 NO-TRACE-STRING 0 4 0 0 NO-TRACE-STRING 0 1 1 0 NO-TRACE-STRING 0 2 2 0 NO-TRACE-STRING 0 3 3 0 NO-TRACE-STRING 0 4 4 0 NO-TRACE-STRING 0 1 5 0 NO-TRACE-STRING 0 2 6 0 NO-TRACE-STRING 0 3 0 0 NO-TRACE-STRING 0 4 1 0 NO-TRACE-STRING 0 1 34 18 0 3 .TXT 11 47 25 0 Cg a5.125000,20.500000,14 {\rtf1\ansi \deff0 {\fonttbl {\f0\fnil Arial;} {\f1\fnil Symbol;} } {\plain {\f1 g} = 0.2} } .EQN 17 -64 13 0 5&0&(_n_u_l_l_&4)&{0:H}NAME({0:u}NAME),{0:K}NAME({0:u}NAME,1),{0:K}NAME({0:u}NAME,.5),{0:K}NAME({0:u}NAME,.2),{0:K}NAME({0:u}NAME,.02),{0:K}NAME({0:u}NAME,0)@&&(_n_u_l_l_&_n_u_l_l_)&{0:u}NAME 0 0 1 0 5 0 1 0 0 1 0 5 0 1 0 1 0 0 NO-TRACE-STRING 0 1 2 0 NO-TRACE-STRING 0 1 1 0 NO-TRACE-STRING 0 1 3 0 NO-TRACE-STRING 0 1 4 0 NO-TRACE-STRING 0 1 5 0 NO-TRACE-STRING 0 3 6 0 NO-TRACE-STRING 0 4 0 0 NO-TRACE-STRING 0 1 1 0 NO-TRACE-STRING 0 2 2 0 NO-TRACE-STRING 0 3 3 0 NO-TRACE-STRING 0 4 4 0 NO-TRACE-STRING 0 1 5 0 NO-TRACE-STRING 0 2 6 0 NO-TRACE-STRING 0 3 0 0 NO-TRACE-STRING 0 4 1 0 NO-TRACE-STRING 0 1 60 32 0 3 .TXT 10 73 14 0 Cg a6.500000,7.000000,87 {\rtf1\ansi \deff0 {\fonttbl {\f0\fnil Arial;} {\f1\fnil Symbol;} } {\plain {\f1 g}{ = 0\par }{\f1 g}{ = 0.05\par }{\f1 g}{ = 0.2\par }{\f1 g}{ = 0.5\par }{\f1 g}{ = 1\par H(u)}} } .TXT 33 -48 15 0 Cg a30.625000,63.625000,45 {\rtf1\ansi \deff0 {\fonttbl {\f0\fnil Arial;} {\f1\fnil Symbol;} } {\plain {Parametric Wiener Filter Transfer Function}} } .TXT 5 -27 22 0 Cg b86.500000,86.500000,348 {\rtf1\ansi \deff0 {\fonttbl {\f0\fnil Arial;} {\f1\fnil Symbol;} } {\plain Note that {\f1 g} = 0 is straight deconvolution, {\f1 g} = 1 is Wiener deconvolution, and 0 < {\f1 g} < 1 is everything in between.\par In this case the signal dominates the noise over almost the entire spectrum, but the blurring is rather severe.\par {\f1 g} = 0.02 gives a pretty good looking restoration MTF. Wiener deconvolution is usually too mild.} }