Power Comparison

All Tables

Initialize

<< d:/math/2004/uta.m

<< d:/math/2004/WriteTable.m

?MergeXYZ

?PowerIA

PowerIA[50, 25, {1}, {0.5}, {0.5}, {}, 
      SigmaSquared -> 0.5]

?SEGain

PowerIA[50, 25, {2}, {}, {0.5}, {}, 
      SigmaSquared -> 0.5]

Compute Tables

ds = {0.25, 0.5, 0.75}; ps = {0.5, 0.75, 1.0};
\(\(tablecols = 
        Prepend[{\ "\omega_0=0.5", \ "\omega_0=0.75", "\omega_0=1.0"}, "\delta_0"];\)\)
rownames = PaddedForm[#, {3, 2}] & /@ 
        Range[0.25, 0.75, 0.25]; 
MySy = 
      StyleForm[#, FontFamily -> Helvetica, 
          FontColor -> RGBColor[0, 0, 1]] &; 
\(\(outDirectory = "d:/math/2004/PowerIA/Simulation/SimpleDynamicIA";\)\)
\(\(outF = {"/outm0p75.tsv", "/outm0p5.tsv", "/outphi0.tsv", "/out0p25.tsv", "/out2.tsv", "/out0p75.tsv", "/out0p9.tsv"};\)\)
phis = {-0.75, -0.5, 0, 0.25, 0.5, 0.75, 0.9}; 
tbl = tdy = tsia = SimOut = 
              Array[0 &, Length[phis]];
\(Do[
phi1 = 
          phis[[
            iPhi]]; 
ω = 
          0.5; \((t1 = 
            Table[{δ, 
                GainTest /. 
                  PowerIA[50, 25, {ω}, {δ}, {phi1}, {}, 
                    SigmaSquared -> 1], 
WaldTest /. 
                  PowerIA[50, 
                    25, {ω/\((1 - δ)\)}, {}, {0.5}, {}, 
                    SigmaSquared -> 1]}, {δ, 0.25, 0.75, 
                0.25}])\); ω = 
          0.75; \((t2 = 
            Table[{δ, 
                GainTest /. 
                  PowerIA[50, 25, {ω}, {δ}, {phi1}, {}, 
                    SigmaSquared -> 1], 
WaldTest /. 
                  PowerIA[50, 
                    25, {ω/\((1 - δ)\)}, {}, {0.5}, {}, 
                    SigmaSquared -> 1]}, {δ, 0.25, 0.75, 
                0.25}])\); ω = 
          1.0; \((t3 = 
            Table[{δ, 
                GainTest /. 
                  PowerIA[50, 25, {ω}, {δ}, {phi1}, {}, 
                    SigmaSquared -> 1], 
WaldTest /. 
                  PowerIA[50, 
                    25, {ω/\((1 - δ)\)}, {}, {0.5}, {}, 
                    SigmaSquared -> 1]}, {δ, 0.25, 0.75, 
                0.25}])\); 
        tdy[[iPhi]] = {\(Transpose[
                t1]\)[[2]], \(Transpose[
                t2]\)[[2]], \(Transpose[
                t3]\)[[2]]}; 
        tsia[[iPhi]] = {\(Transpose[
                t1]\)[[3]], \(Transpose[
                t2]\)[[3]], \(Transpose[
                t3]\)[[3]]}; 
outFil = 
          outDirectory <> 
            outF[[iPhi]]; 
SimOut[[iPhi]] = 
          ReadList[outFil, {Real, Real, Real}]; 
xyz = 
          MergeXYZ[tdy[[iPhi]], 
            tsia[[iPhi]], 
            SimOut[[iPhi]]]; 
tbl[[iPhi]] = 
          TableForm[xyz, 
            TableHeadings -> {Table[
                  StringForm["δ = ``", 
                      NumberForm[
                        ds[[i]], 2]] // 
                    MySy, {i, Length[ds]}], 
                Table[StringForm[\*""(Subscript[ω, 0]) = ``"", 
                      NumberForm[
                        ps[[i]], 3]] // 
                    MySy, {i, Length[ps]}]}]; 
Print[
          StyleForm[
            StringForm[\*""\n\n(Subscript[ϕ, 1]) = ``"", 
              NumberForm[phi1, 2]], FontFamily -> "Helvetica", 
            FontSize -> 14, FontWeight -> "Bold", 
            FontColor -> 
              RGBColor[0, 0.45, 1]]]; 
Print[
          tbl[[iPhi]]], 
{iPhi, Length[phis]}] // Timing\)