(* Content-type: application/mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 6.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[      6701,        208]
NotebookOptionsPosition[      5924,        177]
NotebookOutlinePosition[      6268,        192]
CellTagsIndexPosition[      6225,        189]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["\<\
Group testing matrices
\
\>", "Title",
 CellChangeTimes->{{3.440930658059*^9, 3.4409306911470003`*^9}, {
  3.4409308995109997`*^9, 3.440930910046*^9}, {3.441217566064286*^9, 
  3.441217573158286*^9}, {3.4412176894802856`*^9, 3.441217733173286*^9}, {
  3.441220781460286*^9, 3.441220784191286*^9}, {3.441220887138286*^9, 
  3.441220898752286*^9}, {3.441221676395286*^9, 3.441221702452286*^9}, {
  3.4423098479874*^9, 3.4423098535474*^9}, {3.4423113103924*^9, 
  3.4423113210624*^9}}],

Cell["\<\
From Optical Imaging and Spectroscopy (Wiley - OSA 2009)
By David J. Brady
www.opticalimaging.org
Figure 7.21\
\>", "Subtitle",
 CellChangeTimes->{{3.4409307282320004`*^9, 3.440930761152*^9}, {
   3.440930913613*^9, 3.440930913829*^9}, {3.441217577128286*^9, 
   3.441217577638286*^9}, {3.441217735493286*^9, 3.441217744740286*^9}, {
   3.441220815238286*^9, 3.441220815693286*^9}, 3.4412208544512863`*^9, {
   3.441221706495286*^9, 3.4412217089142857`*^9}, {3.4423098590344*^9, 
   3.4423098760734*^9}, {3.4423113274954*^9, 3.4423113292064*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"h", "[", 
    RowBox[{"i_", ",", "j_", ",", "d_"}], "]"}], "=", 
   RowBox[{"Mod", "[", 
    RowBox[{
     RowBox[{"Floor", "[", 
      RowBox[{"j", "/", 
       RowBox[{"(", 
        RowBox[{"d", "^", "i"}], ")"}]}], "]"}], ",", "d"}], "]"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.41086817854*^9, 3.410868263935*^9}, {
  3.4108683434440002`*^9, 3.41086834873*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"a", "=", 
  RowBox[{"ArrayPlot", "[", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"h", "[", 
       RowBox[{"i", ",", "j", ",", "2"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "0", ",", "6"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"j", ",", "0", ",", "127"}], "}"}]}], "]"}], ",", 
    RowBox[{"AspectRatio", "\[Rule]", "1"}], ",", 
    RowBox[{"PlotLabel", "\[Rule]", "\"\<(a)\>\""}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.410868474723*^9, 3.410868564894*^9}, {3.410868702166*^9,
   3.410868704734*^9}, {3.410869008093*^9, 3.410869008804*^9}, {
  3.410870036491*^9, 3.410870046993*^9}, {3.410870107718*^9, 
  3.4108701435810003`*^9}, {3.4109587500699997`*^9, 3.410958750509*^9}, {
  3.410958837566*^9, 3.410958952*^9}}],

Cell[BoxData[
 GraphicsBox[RasterBox[CompressedData["
1:eJxTTMoPSmJiYGBgB+IGIGakEDBQCCg1n9r6B0p+IPkDSw8wHmAIAOP/A/A=

   "], {{0, 0}, {128, 7}}, {0, 1}],
  AspectRatio->1,
  Frame->Automatic,
  FrameLabel->{None, None},
  FrameTicks->{{None, None}, {None, None}},
  PlotLabel->FormBox["\"(a)\"", TraditionalForm]]], "Output",
 CellChangeTimes->{3.4109589527790003`*^9}]
}, Open  ]],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.41095889351*^9, 3.410958932552*^9}}],

Cell[BoxData[""], "Input",
 CellChangeTimes->{3.410958751583*^9, {3.410958850166*^9, 3.410958871725*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"b", "=", 
  RowBox[{"ArrayPlot", "[", 
   RowBox[{
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"h", "[", 
       RowBox[{"i", ",", "j", ",", "4"}], "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", "0", ",", "3"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"j", ",", "0", ",", "127"}], "}"}]}], "]"}], ",", 
    RowBox[{"AspectRatio", "\[Rule]", "1"}], ",", 
    RowBox[{"PlotLabel", "\[Rule]", "\"\<(b)\>\""}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.410868474723*^9, 3.410868564894*^9}, {3.410868702166*^9,
   3.410868704734*^9}, {3.410869008093*^9, 3.410869008804*^9}, {
  3.410870036491*^9, 3.410870046993*^9}, {3.410870107718*^9, 
  3.4108701869700003`*^9}, {3.410958813196*^9, 3.410958813619*^9}, {
  3.410958963516*^9, 3.410958967421*^9}}],

Cell[BoxData[
 GraphicsBox[RasterBox[CompressedData["
1:eJxTTMoPSmJiYGBgAeIGIGamEDBRCAiZx4gGGNAANfTD1A0In4lxQDEAa+kF
rQ==
   "], {{0, 0}, {128, 4}}, {0, 3}],
  AspectRatio->1,
  Frame->Automatic,
  FrameLabel->{None, None},
  FrameTicks->{{None, None}, {None, None}},
  PlotLabel->FormBox["\"(b)\"", TraditionalForm]]], "Output",
 CellChangeTimes->{3.4109589794779997`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"GraphicsArray", "[", 
   RowBox[{"{", 
    RowBox[{"a", ",", "b"}], "}"}], "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.41095881935*^9, 3.4109588292349997`*^9}, 
   3.410958988493*^9}],

Cell[BoxData[
 GraphicsBox[{{}, {InsetBox[
     GraphicsBox[RasterBox[CompressedData["
1:eJxTTMoPSmJiYGBgB+IGIGakEDBQCCg1n9r6B0p+IPkDSw8wHmAIAOP/A/A=

       "], {{0, 0}, {128, 7}}, {0, 1}],
      AspectRatio->1,
      Frame->Automatic,
      FrameLabel->{None, None},
      FrameTicks->{{None, None}, {None, None}},
      PlotLabel->FormBox["\"(a)\"", TraditionalForm]], {192., -189.}, 
     ImageScaled[{0.5, 0.5}], {360., 360.}], InsetBox[
     GraphicsBox[RasterBox[CompressedData["
1:eJxTTMoPSmJiYGBgAeIGIGamEDBRCAiZx4gGGNAANfTD1A0In4lxQDEAa+kF
rQ==
       "], {{0, 0}, {128, 4}}, {0, 3}],
      AspectRatio->1,
      Frame->Automatic,
      FrameLabel->{None, None},
      FrameTicks->{{None, None}, {None, None}},
      PlotLabel->FormBox["\"(b)\"", TraditionalForm]], {576., -189.}, 
     ImageScaled[{0.5, 0.5}], {360., 360.}]}, {}},
  ContentSelectable->True,
  ImageSize->{641., Automatic},
  PlotRangePadding->{6, 5}]], "Output",
 CellChangeTimes->{3.410958989618*^9}]
}, Open  ]]
}, Open  ]]
},
WindowSize->{856, 688},
WindowMargins->{{49, Automatic}, {0, Automatic}},
FrontEndVersion->"7.0 for Microsoft Windows (32-bit) (November 10, 2008)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[567, 22, 492, 10, 137, "Title"],
Cell[1062, 34, 556, 11, 133, "Subtitle"],
Cell[1621, 47, 418, 13, 31, "Input"],
Cell[CellGroupData[{
Cell[2064, 64, 800, 18, 31, "Input"],
Cell[2867, 84, 369, 10, 390, "Output"]
}, Open  ]],
Cell[3251, 97, 85, 1, 31, "Input"],
Cell[3339, 100, 105, 1, 31, "Input"],
Cell[CellGroupData[{
Cell[3469, 105, 798, 18, 31, "Input"],
Cell[4270, 125, 373, 10, 390, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[4680, 140, 233, 6, 31, "Input"],
Cell[4916, 148, 980, 25, 350, "Output"]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)