================================================================================ Now Frequency(Hz): 0.00146484, Period(s): 682.667 ================================================================================ -------------------------------------------------------------------------------- Calculate response functions for output variable 0 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 83 Squared coherence: 0.88469 Estimated response function: ( 3.7554e-01, 4.2087e-01), ( 3.3244e-01, 4.4139e-01) Amplitude of the estimated response function: 5.6406e-01, 5.5257e-01 Phase(deg.) of the estimated response function: 48.3, 53.0 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.022806 Sum of weights: 82.9476 Squared coherence: 0.885022 Estimated response function: ( 3.7573e-01, 4.2087e-01), ( 3.3240e-01, 4.4119e-01) Amplitude of the estimated response function: 5.6419e-01, 5.5240e-01 Phase(deg.) of the estimated response function: 48.2, 53.0 Weighted residual power: 0.000981678 Iteration number = 1 Scale factor: 0.0227546 Sum of weights: 82.9452 Squared coherence: 0.885037 Estimated response function: ( 3.7574e-01, 4.2087e-01), ( 3.3240e-01, 4.4118e-01) Amplitude of the estimated response function: 5.6419e-01, 5.5239e-01 Phase(deg.) of the estimated response function: 48.2, 53.0 Weighted residual power: 0.000981555 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.0227546 Parameter c: 2.72115 Sum of weights: 52.3421 Squared coherence: 0.936934 Estimated response function: ( 3.6478e-01, 4.3058e-01), ( 3.1974e-01, 4.3805e-01) Amplitude of the estimated response function: 5.6433e-01, 5.4233e-01 Phase(deg.) of the estimated response function: 49.7, 53.9 Weighted residual power: 0.000472609 Iteration number = 1 Scale factor: 0.0227546 Parameter c: 2.72115 Sum of weights: 52.4908 Squared coherence: 0.936921 Estimated response function: ( 3.6093e-01, 4.3574e-01), ( 3.1606e-01, 4.3815e-01) Amplitude of the estimated response function: 5.6581e-01, 5.4025e-01 Phase(deg.) of the estimated response function: 50.4, 54.2 Weighted residual power: 0.000466372 Iteration number = 2 Scale factor: 0.0227546 Parameter c: 2.72115 Sum of weights: 52.4937 Squared coherence: 0.937038 Estimated response function: ( 3.5945e-01, 4.3869e-01), ( 3.1506e-01, 4.3901e-01) Amplitude of the estimated response function: 5.6715e-01, 5.4036e-01 Phase(deg.) of the estimated response function: 50.7, 54.3 Weighted residual power: 0.000464009 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Calculate response functions for output variable 1 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 83 Squared coherence: 0.943059 Estimated response function: ( -6.7784e-01, -8.3556e-01), ( -1.5476e-01, -1.5345e-01) Amplitude of the estimated response function: 1.0759e+00, 2.1794e-01 Phase(deg.) of the estimated response function: -129.1, -135.2 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.0222735 Sum of weights: 82.9436 Squared coherence: 0.943228 Estimated response function: ( -6.7752e-01, -8.3556e-01), ( -1.5541e-01, -1.5301e-01) Amplitude of the estimated response function: 1.0757e+00, 2.1809e-01 Phase(deg.) of the estimated response function: -129.0, -135.4 Weighted residual power: 0.00091751 Iteration number = 1 Scale factor: 0.0223012 Sum of weights: 82.9424 Squared coherence: 0.943231 Estimated response function: ( -6.7752e-01, -8.3556e-01), ( -1.5542e-01, -1.5300e-01) Amplitude of the estimated response function: 1.0757e+00, 2.1809e-01 Phase(deg.) of the estimated response function: -129.0, -135.4 Weighted residual power: 0.000917455 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.0223012 Parameter c: 2.72115 Sum of weights: 52.5864 Squared coherence: 0.96767 Estimated response function: ( -6.7183e-01, -8.2903e-01), ( -1.6088e-01, -1.5195e-01) Amplitude of the estimated response function: 1.0671e+00, 2.2129e-01 Phase(deg.) of the estimated response function: -129.0, -136.6 Weighted residual power: 0.000494981 Iteration number = 1 Scale factor: 0.0223012 Parameter c: 2.72115 Sum of weights: 52.6135 Squared coherence: 0.967771 Estimated response function: ( -6.6963e-01, -8.2683e-01), ( -1.6317e-01, -1.5197e-01) Amplitude of the estimated response function: 1.0640e+00, 2.2298e-01 Phase(deg.) of the estimated response function: -129.0, -137.0 Weighted residual power: 0.000492262 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Estimate errors by fixed-weights bootstrap -------------------------------------------------------------------------------- ================================================================================ Now Frequency(Hz): 0.00195312, Period(s): 512 ================================================================================ -------------------------------------------------------------------------------- Calculate response functions for output variable 0 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 83 Squared coherence: 0.881175 Estimated response function: ( 4.5545e-01, 4.7352e-01), ( 3.8421e-01, 5.4449e-01) Amplitude of the estimated response function: 6.5701e-01, 6.6640e-01 Phase(deg.) of the estimated response function: 46.1, 54.8 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.0218496 Sum of weights: 82.6197 Squared coherence: 0.882984 Estimated response function: ( 4.5500e-01, 4.7346e-01), ( 3.8344e-01, 5.4390e-01) Amplitude of the estimated response function: 6.5665e-01, 6.6547e-01 Phase(deg.) of the estimated response function: 46.1, 54.8 Weighted residual power: 0.000963873 Iteration number = 1 Scale factor: 0.0216392 Sum of weights: 82.5946 Squared coherence: 0.883104 Estimated response function: ( 4.5497e-01, 4.7346e-01), ( 3.8339e-01, 5.4387e-01) Amplitude of the estimated response function: 6.5663e-01, 6.6541e-01 Phase(deg.) of the estimated response function: 46.1, 54.8 Weighted residual power: 0.000962433 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.0216392 Parameter c: 2.72115 Sum of weights: 51.9848 Squared coherence: 0.94671 Estimated response function: ( 4.5796e-01, 4.7390e-01), ( 3.9126e-01, 5.4359e-01) Amplitude of the estimated response function: 6.5902e-01, 6.6975e-01 Phase(deg.) of the estimated response function: 46.0, 54.3 Weighted residual power: 0.000405453 Iteration number = 1 Scale factor: 0.0216392 Parameter c: 2.72115 Sum of weights: 51.9775 Squared coherence: 0.947136 Estimated response function: ( 4.5955e-01, 4.7370e-01), ( 3.9451e-01, 5.4318e-01) Amplitude of the estimated response function: 6.5998e-01, 6.7133e-01 Phase(deg.) of the estimated response function: 45.9, 54.0 Weighted residual power: 0.000403474 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Calculate response functions for output variable 1 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 83 Squared coherence: 0.952068 Estimated response function: ( -8.0225e-01, -1.0040e+00), ( -1.8291e-01, -2.1313e-01) Amplitude of the estimated response function: 1.2852e+00, 2.8086e-01 Phase(deg.) of the estimated response function: -128.6, -130.6 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.0202773 Sum of weights: 82.6902 Squared coherence: 0.953928 Estimated response function: ( -8.0554e-01, -1.0066e+00), ( -1.8058e-01, -2.1369e-01) Amplitude of the estimated response function: 1.2892e+00, 2.7978e-01 Phase(deg.) of the estimated response function: -128.7, -130.2 Weighted residual power: 0.000718665 Iteration number = 1 Scale factor: 0.0200992 Sum of weights: 82.6802 Squared coherence: 0.953988 Estimated response function: ( -8.0565e-01, -1.0067e+00), ( -1.8050e-01, -2.1371e-01) Amplitude of the estimated response function: 1.2894e+00, 2.7974e-01 Phase(deg.) of the estimated response function: -128.7, -130.2 Weighted residual power: 0.000717807 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.0200992 Parameter c: 2.72115 Sum of weights: 54.5999 Squared coherence: 0.977877 Estimated response function: ( -8.1360e-01, -1.0059e+00), ( -1.8360e-01, -1.9836e-01) Amplitude of the estimated response function: 1.2937e+00, 2.7029e-01 Phase(deg.) of the estimated response function: -129.0, -132.8 Weighted residual power: 0.000366205 Iteration number = 1 Scale factor: 0.0200992 Parameter c: 2.72115 Sum of weights: 54.666 Squared coherence: 0.978069 Estimated response function: ( -8.1499e-01, -1.0054e+00), ( -1.8576e-01, -1.9361e-01) Amplitude of the estimated response function: 1.2943e+00, 2.6831e-01 Phase(deg.) of the estimated response function: -129.0, -133.8 Weighted residual power: 0.000362375 Iteration number = 2 Scale factor: 0.0200992 Parameter c: 2.72115 Sum of weights: 54.6546 Squared coherence: 0.978125 Estimated response function: ( -8.1503e-01, -1.0052e+00), ( -1.8687e-01, -1.9200e-01) Amplitude of the estimated response function: 1.2941e+00, 2.6793e-01 Phase(deg.) of the estimated response function: -129.0, -134.2 Weighted residual power: 0.000361139 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Estimate errors by fixed-weights bootstrap -------------------------------------------------------------------------------- ================================================================================ Now Frequency(Hz): 0.00292969, Period(s): 341.333 ================================================================================ -------------------------------------------------------------------------------- Calculate response functions for output variable 0 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 167 Squared coherence: 0.810574 Estimated response function: ( 4.7326e-01, 5.9456e-01), ( 4.3491e-01, 7.6086e-01) Amplitude of the estimated response function: 7.5992e-01, 8.7639e-01 Phase(deg.) of the estimated response function: 51.5, 60.2 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.0288886 Sum of weights: 166.415 Squared coherence: 0.814823 Estimated response function: ( 4.7388e-01, 5.9538e-01), ( 4.3979e-01, 7.6089e-01) Amplitude of the estimated response function: 7.6094e-01, 8.7885e-01 Phase(deg.) of the estimated response function: 51.5, 60.0 Weighted residual power: 0.00175102 Iteration number = 1 Scale factor: 0.0285009 Sum of weights: 166.298 Squared coherence: 0.81525 Estimated response function: ( 4.7421e-01, 5.9551e-01), ( 4.4048e-01, 7.6065e-01) Amplitude of the estimated response function: 7.6125e-01, 8.7898e-01 Phase(deg.) of the estimated response function: 51.5, 59.9 Weighted residual power: 0.00174595 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.0285009 Parameter c: 2.69063 Sum of weights: 100.247 Squared coherence: 0.894387 Estimated response function: ( 4.9245e-01, 5.9734e-01), ( 4.8112e-01, 7.4061e-01) Amplitude of the estimated response function: 7.7416e-01, 8.8317e-01 Phase(deg.) of the estimated response function: 50.5, 57.0 Weighted residual power: 0.00083162 Iteration number = 1 Scale factor: 0.0285009 Parameter c: 2.69063 Sum of weights: 100.909 Squared coherence: 0.897317 Estimated response function: ( 5.0209e-01, 5.9962e-01), ( 5.0016e-01, 7.3271e-01) Amplitude of the estimated response function: 7.8207e-01, 8.8714e-01 Phase(deg.) of the estimated response function: 50.1, 55.7 Weighted residual power: 0.000825873 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Calculate response functions for output variable 1 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 167 Squared coherence: 0.917953 Estimated response function: ( -9.3048e-01, -1.3085e+00), ( -1.8189e-01, -2.2954e-01) Amplitude of the estimated response function: 1.6056e+00, 2.9286e-01 Phase(deg.) of the estimated response function: -125.4, -128.4 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.0333414 Sum of weights: 166.919 Squared coherence: 0.918131 Estimated response function: ( -9.3035e-01, -1.3084e+00), ( -1.8163e-01, -2.2931e-01) Amplitude of the estimated response function: 1.6054e+00, 2.9253e-01 Phase(deg.) of the estimated response function: -125.4, -128.4 Weighted residual power: 0.00170763 Iteration number = 1 Scale factor: 0.0332804 Sum of weights: 166.913 Squared coherence: 0.918144 Estimated response function: ( -9.3035e-01, -1.3084e+00), ( -1.8161e-01, -2.2931e-01) Amplitude of the estimated response function: 1.6054e+00, 2.9252e-01 Phase(deg.) of the estimated response function: -125.4, -128.4 Weighted residual power: 0.00170733 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.0332804 Parameter c: 2.69063 Sum of weights: 114.527 Squared coherence: 0.95728 Estimated response function: ( -9.3352e-01, -1.3030e+00), ( -1.7784e-01, -2.2987e-01) Amplitude of the estimated response function: 1.6029e+00, 2.9064e-01 Phase(deg.) of the estimated response function: -125.6, -127.7 Weighted residual power: 0.000879729 Iteration number = 1 Scale factor: 0.0332804 Parameter c: 2.69063 Sum of weights: 114.55 Squared coherence: 0.957274 Estimated response function: ( -9.3460e-01, -1.3012e+00), ( -1.7667e-01, -2.2967e-01) Amplitude of the estimated response function: 1.6020e+00, 2.8976e-01 Phase(deg.) of the estimated response function: -125.7, -127.6 Weighted residual power: 0.00087964 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Estimate errors by fixed-weights bootstrap -------------------------------------------------------------------------------- ================================================================================ Now Frequency(Hz): 0.00390625, Period(s): 256 ================================================================================ -------------------------------------------------------------------------------- Calculate response functions for output variable 0 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 167 Squared coherence: 0.795293 Estimated response function: ( 6.0472e-01, 7.0462e-01), ( 4.3766e-01, 9.1245e-01) Amplitude of the estimated response function: 9.2854e-01, 1.0120e+00 Phase(deg.) of the estimated response function: 49.4, 64.4 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.0310511 Sum of weights: 166.784 Squared coherence: 0.796375 Estimated response function: ( 6.0490e-01, 7.0408e-01), ( 4.3770e-01, 9.1133e-01) Amplitude of the estimated response function: 9.2824e-01, 1.0110e+00 Phase(deg.) of the estimated response function: 49.3, 64.3 Weighted residual power: 0.00163652 Iteration number = 1 Scale factor: 0.0309407 Sum of weights: 166.773 Squared coherence: 0.796417 Estimated response function: ( 6.0491e-01, 7.0403e-01), ( 4.3769e-01, 9.1124e-01) Amplitude of the estimated response function: 9.2821e-01, 1.0109e+00 Phase(deg.) of the estimated response function: 49.3, 64.3 Weighted residual power: 0.00163594 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.0309407 Parameter c: 2.69063 Sum of weights: 109.585 Squared coherence: 0.890475 Estimated response function: ( 6.2530e-01, 6.8916e-01), ( 4.4520e-01, 9.1482e-01) Amplitude of the estimated response function: 9.3056e-01, 1.0174e+00 Phase(deg.) of the estimated response function: 47.8, 64.0 Weighted residual power: 0.000825154 Iteration number = 1 Scale factor: 0.0309407 Parameter c: 2.69063 Sum of weights: 109.718 Squared coherence: 0.890703 Estimated response function: ( 6.3352e-01, 6.8379e-01), ( 4.4827e-01, 9.1508e-01) Amplitude of the estimated response function: 9.3215e-01, 1.0190e+00 Phase(deg.) of the estimated response function: 47.2, 63.9 Weighted residual power: 0.000824492 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Calculate response functions for output variable 1 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 167 Squared coherence: 0.898233 Estimated response function: ( -1.2107e+00, -1.5316e+00), ( -2.0797e-01, -2.9546e-01) Amplitude of the estimated response function: 1.9524e+00, 3.6131e-01 Phase(deg.) of the estimated response function: -128.3, -125.1 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.0320828 Sum of weights: 166.624 Squared coherence: 0.899563 Estimated response function: ( -1.2101e+00, -1.5314e+00), ( -2.0696e-01, -2.9203e-01) Amplitude of the estimated response function: 1.9518e+00, 3.5793e-01 Phase(deg.) of the estimated response function: -128.3, -125.3 Weighted residual power: 0.00174138 Iteration number = 1 Scale factor: 0.0320228 Sum of weights: 166.616 Squared coherence: 0.899589 Estimated response function: ( -1.2101e+00, -1.5314e+00), ( -2.0694e-01, -2.9196e-01) Amplitude of the estimated response function: 1.9517e+00, 3.5786e-01 Phase(deg.) of the estimated response function: -128.3, -125.3 Weighted residual power: 0.00174085 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.0320228 Parameter c: 2.69063 Sum of weights: 111.539 Squared coherence: 0.956214 Estimated response function: ( -1.2112e+00, -1.5584e+00), ( -1.9635e-01, -3.2332e-01) Amplitude of the estimated response function: 1.9737e+00, 3.7827e-01 Phase(deg.) of the estimated response function: -127.9, -121.3 Weighted residual power: 0.000787953 Iteration number = 1 Scale factor: 0.0320228 Parameter c: 2.69063 Sum of weights: 111.669 Squared coherence: 0.957058 Estimated response function: ( -1.2109e+00, -1.5662e+00), ( -1.9305e-01, -3.3339e-01) Amplitude of the estimated response function: 1.9797e+00, 3.8525e-01 Phase(deg.) of the estimated response function: -127.7, -120.1 Weighted residual power: 0.000778579 Iteration number = 2 Scale factor: 0.0320228 Parameter c: 2.69063 Sum of weights: 111.647 Squared coherence: 0.957266 Estimated response function: ( -1.2106e+00, -1.5685e+00), ( -1.9204e-01, -3.3643e-01) Amplitude of the estimated response function: 1.9813e+00, 3.8739e-01 Phase(deg.) of the estimated response function: -127.7, -119.7 Weighted residual power: 0.000776143 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Estimate errors by fixed-weights bootstrap -------------------------------------------------------------------------------- ================================================================================ Now Frequency(Hz): 0.00585938, Period(s): 170.667 ================================================================================ -------------------------------------------------------------------------------- Calculate response functions for output variable 0 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 336 Squared coherence: 0.69714 Estimated response function: ( 8.2872e-01, 9.2629e-01), ( 6.3990e-01, 1.1135e+00) Amplitude of the estimated response function: 1.2429e+00, 1.2843e+00 Phase(deg.) of the estimated response function: 48.2, 60.1 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.0475882 Sum of weights: 334.961 Squared coherence: 0.70271 Estimated response function: ( 8.2515e-01, 9.2992e-01), ( 6.4011e-01, 1.1171e+00) Amplitude of the estimated response function: 1.2432e+00, 1.2875e+00 Phase(deg.) of the estimated response function: 48.4, 60.2 Weighted residual power: 0.00363032 Iteration number = 1 Scale factor: 0.0477639 Sum of weights: 334.979 Squared coherence: 0.702626 Estimated response function: ( 8.2517e-01, 9.2988e-01), ( 6.4016e-01, 1.1171e+00) Amplitude of the estimated response function: 1.2432e+00, 1.2875e+00 Phase(deg.) of the estimated response function: 48.4, 60.2 Weighted residual power: 0.00363163 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.0477639 Parameter c: 2.67584 Sum of weights: 229.978 Squared coherence: 0.830192 Estimated response function: ( 8.3095e-01, 9.1333e-01), ( 6.5443e-01, 1.1158e+00) Amplitude of the estimated response function: 1.2348e+00, 1.2935e+00 Phase(deg.) of the estimated response function: 47.7, 59.6 Weighted residual power: 0.00176792 Iteration number = 1 Scale factor: 0.0477639 Parameter c: 2.67584 Sum of weights: 229.971 Squared coherence: 0.830423 Estimated response function: ( 8.3299e-01, 9.0717e-01), ( 6.5910e-01, 1.1151e+00) Amplitude of the estimated response function: 1.2316e+00, 1.2954e+00 Phase(deg.) of the estimated response function: 47.4, 59.4 Weighted residual power: 0.00176226 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Calculate response functions for output variable 1 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 336 Squared coherence: 0.858605 Estimated response function: ( -1.5595e+00, -1.8754e+00), ( -3.0034e-01, -3.2833e-01) Amplitude of the estimated response function: 2.4391e+00, 4.4498e-01 Phase(deg.) of the estimated response function: -129.7, -132.5 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.043179 Sum of weights: 334.663 Squared coherence: 0.862036 Estimated response function: ( -1.5605e+00, -1.8813e+00), ( -3.0019e-01, -3.2961e-01) Amplitude of the estimated response function: 2.4443e+00, 4.4583e-01 Phase(deg.) of the estimated response function: -129.7, -132.3 Weighted residual power: 0.00348942 Iteration number = 1 Scale factor: 0.0434217 Sum of weights: 334.699 Squared coherence: 0.861959 Estimated response function: ( -1.5605e+00, -1.8812e+00), ( -3.0017e-01, -3.2963e-01) Amplitude of the estimated response function: 2.4442e+00, 4.4582e-01 Phase(deg.) of the estimated response function: -129.7, -132.3 Weighted residual power: 0.00349146 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.0434217 Parameter c: 2.67584 Sum of weights: 218.773 Squared coherence: 0.934817 Estimated response function: ( -1.5645e+00, -1.9067e+00), ( -2.9279e-01, -3.4999e-01) Amplitude of the estimated response function: 2.4664e+00, 4.5631e-01 Phase(deg.) of the estimated response function: -129.4, -129.9 Weighted residual power: 0.00156197 Iteration number = 1 Scale factor: 0.0434217 Parameter c: 2.67584 Sum of weights: 218.759 Squared coherence: 0.93531 Estimated response function: ( -1.5655e+00, -1.9144e+00), ( -2.9153e-01, -3.5669e-01) Amplitude of the estimated response function: 2.4730e+00, 4.6067e-01 Phase(deg.) of the estimated response function: -129.3, -129.3 Weighted residual power: 0.00155742 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Estimate errors by fixed-weights bootstrap -------------------------------------------------------------------------------- ================================================================================ Now Frequency(Hz): 0.0078125, Period(s): 128 ================================================================================ -------------------------------------------------------------------------------- Calculate response functions for output variable 0 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 336 Squared coherence: 0.662431 Estimated response function: ( 1.0339e+00, 1.0460e+00), ( 8.5128e-01, 1.3658e+00) Amplitude of the estimated response function: 1.4707e+00, 1.6093e+00 Phase(deg.) of the estimated response function: 45.3, 58.1 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.0460144 Sum of weights: 335.168 Squared coherence: 0.666596 Estimated response function: ( 1.0353e+00, 1.0504e+00), ( 8.5002e-01, 1.3582e+00) Amplitude of the estimated response function: 1.4749e+00, 1.6022e+00 Phase(deg.) of the estimated response function: 45.4, 58.0 Weighted residual power: 0.00337387 Iteration number = 1 Scale factor: 0.0459654 Sum of weights: 335.155 Squared coherence: 0.666657 Estimated response function: ( 1.0354e+00, 1.0505e+00), ( 8.4997e-01, 1.3581e+00) Amplitude of the estimated response function: 1.4750e+00, 1.6021e+00 Phase(deg.) of the estimated response function: 45.4, 58.0 Weighted residual power: 0.00337296 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.0459654 Parameter c: 2.67584 Sum of weights: 231.347 Squared coherence: 0.814872 Estimated response function: ( 1.0069e+00, 1.0782e+00), ( 8.7246e-01, 1.3195e+00) Amplitude of the estimated response function: 1.4753e+00, 1.5818e+00 Phase(deg.) of the estimated response function: 47.0, 56.5 Weighted residual power: 0.00158587 Iteration number = 1 Scale factor: 0.0459654 Parameter c: 2.67584 Sum of weights: 231.639 Squared coherence: 0.815162 Estimated response function: ( 9.9871e-01, 1.0857e+00), ( 8.7903e-01, 1.3108e+00) Amplitude of the estimated response function: 1.4752e+00, 1.5783e+00 Phase(deg.) of the estimated response function: 47.4, 56.2 Weighted residual power: 0.00158008 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Calculate response functions for output variable 1 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 336 Squared coherence: 0.81653 Estimated response function: ( -1.7951e+00, -2.2401e+00), ( -4.1262e-01, -3.9622e-01) Amplitude of the estimated response function: 2.8706e+00, 5.7205e-01 Phase(deg.) of the estimated response function: -128.7, -136.2 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.0455373 Sum of weights: 335.227 Squared coherence: 0.818577 Estimated response function: ( -1.7943e+00, -2.2358e+00), ( -4.1218e-01, -3.9641e-01) Amplitude of the estimated response function: 2.8668e+00, 5.7187e-01 Phase(deg.) of the estimated response function: -128.7, -136.1 Weighted residual power: 0.0034975 Iteration number = 1 Scale factor: 0.0453679 Sum of weights: 335.203 Squared coherence: 0.818638 Estimated response function: ( -1.7943e+00, -2.2357e+00), ( -4.1220e-01, -3.9638e-01) Amplitude of the estimated response function: 2.8667e+00, 5.7186e-01 Phase(deg.) of the estimated response function: -128.7, -136.1 Weighted residual power: 0.00349601 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.0453679 Parameter c: 2.67584 Sum of weights: 224.586 Squared coherence: 0.912459 Estimated response function: ( -1.8245e+00, -2.2802e+00), ( -3.8474e-01, -4.1682e-01) Amplitude of the estimated response function: 2.9203e+00, 5.6724e-01 Phase(deg.) of the estimated response function: -128.7, -132.7 Weighted residual power: 0.00160985 Iteration number = 1 Scale factor: 0.0453679 Parameter c: 2.67584 Sum of weights: 224.574 Squared coherence: 0.913729 Estimated response function: ( -1.8321e+00, -2.2944e+00), ( -3.7612e-01, -4.2190e-01) Amplitude of the estimated response function: 2.9361e+00, 5.6521e-01 Phase(deg.) of the estimated response function: -128.6, -131.7 Weighted residual power: 0.00159654 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Estimate errors by fixed-weights bootstrap -------------------------------------------------------------------------------- ================================================================================ Now Frequency(Hz): 0.0117188, Period(s): 85.3333 ================================================================================ -------------------------------------------------------------------------------- Calculate response functions for output variable 0 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 674 Squared coherence: 0.562555 Estimated response function: ( 1.3891e+00, 1.2224e+00), ( 1.0533e+00, 1.8187e+00) Amplitude of the estimated response function: 1.8504e+00, 2.1017e+00 Phase(deg.) of the estimated response function: 41.3, 59.9 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.0742573 Sum of weights: 673.05 Squared coherence: 0.565359 Estimated response function: ( 1.3889e+00, 1.2208e+00), ( 1.0530e+00, 1.8221e+00) Amplitude of the estimated response function: 1.8491e+00, 2.1045e+00 Phase(deg.) of the estimated response function: 41.3, 60.0 Weighted residual power: 0.00743701 Iteration number = 1 Scale factor: 0.0744731 Sum of weights: 673.077 Squared coherence: 0.565286 Estimated response function: ( 1.3889e+00, 1.2207e+00), ( 1.0531e+00, 1.8221e+00) Amplitude of the estimated response function: 1.8491e+00, 2.1045e+00 Phase(deg.) of the estimated response function: 41.3, 60.0 Weighted residual power: 0.00743906 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.0744731 Parameter c: 2.66852 Sum of weights: 489.075 Squared coherence: 0.742511 Estimated response function: ( 1.3675e+00, 1.2292e+00), ( 1.0951e+00, 1.8011e+00) Amplitude of the estimated response function: 1.8387e+00, 2.1079e+00 Phase(deg.) of the estimated response function: 42.0, 58.7 Weighted residual power: 0.00353595 Iteration number = 1 Scale factor: 0.0744731 Parameter c: 2.66852 Sum of weights: 489.254 Squared coherence: 0.74254 Estimated response function: ( 1.3622e+00, 1.2308e+00), ( 1.1061e+00, 1.7975e+00) Amplitude of the estimated response function: 1.8358e+00, 2.1105e+00 Phase(deg.) of the estimated response function: 42.1, 58.4 Weighted residual power: 0.00353299 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Calculate response functions for output variable 1 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 674 Squared coherence: 0.758964 Estimated response function: ( -2.4524e+00, -2.9277e+00), ( -4.3863e-01, -5.8231e-01) Amplitude of the estimated response function: 3.8192e+00, 7.2903e-01 Phase(deg.) of the estimated response function: -130.0, -127.0 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.0736573 Sum of weights: 672.775 Squared coherence: 0.761626 Estimated response function: ( -2.4472e+00, -2.9333e+00), ( -4.4179e-01, -5.8236e-01) Amplitude of the estimated response function: 3.8200e+00, 7.3097e-01 Phase(deg.) of the estimated response function: -129.8, -127.2 Weighted residual power: 0.00768526 Iteration number = 1 Scale factor: 0.0735713 Sum of weights: 672.759 Squared coherence: 0.761659 Estimated response function: ( -2.4471e+00, -2.9334e+00), ( -4.4183e-01, -5.8239e-01) Amplitude of the estimated response function: 3.8201e+00, 7.3102e-01 Phase(deg.) of the estimated response function: -129.8, -127.2 Weighted residual power: 0.00768386 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.0735713 Parameter c: 2.66852 Sum of weights: 477.002 Squared coherence: 0.869332 Estimated response function: ( -2.4661e+00, -2.9406e+00), ( -4.6564e-01, -6.1156e-01) Amplitude of the estimated response function: 3.8378e+00, 7.6865e-01 Phase(deg.) of the estimated response function: -130.0, -127.3 Weighted residual power: 0.003871 Iteration number = 1 Scale factor: 0.0735713 Parameter c: 2.66852 Sum of weights: 476.746 Squared coherence: 0.86986 Estimated response function: ( -2.4713e+00, -2.9427e+00), ( -4.7131e-01, -6.2091e-01) Amplitude of the estimated response function: 3.8428e+00, 7.7953e-01 Phase(deg.) of the estimated response function: -130.0, -127.2 Weighted residual power: 0.00386602 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Estimate errors by fixed-weights bootstrap -------------------------------------------------------------------------------- ================================================================================ Now Frequency(Hz): 0.015625, Period(s): 64 ================================================================================ -------------------------------------------------------------------------------- Calculate response functions for output variable 0 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 674 Squared coherence: 0.488236 Estimated response function: ( 1.5185e+00, 1.3183e+00), ( 1.3435e+00, 2.1202e+00) Amplitude of the estimated response function: 2.0109e+00, 2.5100e+00 Phase(deg.) of the estimated response function: 41.0, 57.6 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.0724904 Sum of weights: 673.116 Squared coherence: 0.491499 Estimated response function: ( 1.5222e+00, 1.3148e+00), ( 1.3426e+00, 2.1181e+00) Amplitude of the estimated response function: 2.0114e+00, 2.5078e+00 Phase(deg.) of the estimated response function: 40.8, 57.6 Weighted residual power: 0.00708811 Iteration number = 1 Scale factor: 0.0723684 Sum of weights: 673.104 Squared coherence: 0.491541 Estimated response function: ( 1.5223e+00, 1.3148e+00), ( 1.3426e+00, 2.1181e+00) Amplitude of the estimated response function: 2.0114e+00, 2.5077e+00 Phase(deg.) of the estimated response function: 40.8, 57.6 Weighted residual power: 0.00708714 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.0723684 Parameter c: 2.66852 Sum of weights: 484.607 Squared coherence: 0.673544 Estimated response function: ( 1.5379e+00, 1.3279e+00), ( 1.3537e+00, 2.0805e+00) Amplitude of the estimated response function: 2.0319e+00, 2.4821e+00 Phase(deg.) of the estimated response function: 40.8, 57.0 Weighted residual power: 0.00356131 Iteration number = 1 Scale factor: 0.0723684 Parameter c: 2.66852 Sum of weights: 484.737 Squared coherence: 0.673597 Estimated response function: ( 1.5427e+00, 1.3317e+00), ( 1.3569e+00, 2.0699e+00) Amplitude of the estimated response function: 2.0379e+00, 2.4750e+00 Phase(deg.) of the estimated response function: 40.8, 56.8 Weighted residual power: 0.00356028 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Calculate response functions for output variable 1 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 674 Squared coherence: 0.671128 Estimated response function: ( -2.9422e+00, -3.2299e+00), ( -5.5561e-01, -5.7310e-01) Amplitude of the estimated response function: 4.3691e+00, 7.9821e-01 Phase(deg.) of the estimated response function: -132.3, -134.1 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.0683443 Sum of weights: 671.868 Squared coherence: 0.676499 Estimated response function: ( -2.9375e+00, -3.2372e+00), ( -5.6057e-01, -5.7689e-01) Amplitude of the estimated response function: 4.3713e+00, 8.0439e-01 Phase(deg.) of the estimated response function: -132.2, -134.2 Weighted residual power: 0.00826111 Iteration number = 1 Scale factor: 0.0686162 Sum of weights: 671.933 Squared coherence: 0.676362 Estimated response function: ( -2.9375e+00, -3.2371e+00), ( -5.6056e-01, -5.7699e-01) Amplitude of the estimated response function: 4.3713e+00, 8.0445e-01 Phase(deg.) of the estimated response function: -132.2, -134.2 Weighted residual power: 0.00826535 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.0686162 Parameter c: 2.66852 Sum of weights: 453.048 Squared coherence: 0.839968 Estimated response function: ( -2.9356e+00, -3.2364e+00), ( -5.7037e-01, -6.1945e-01) Amplitude of the estimated response function: 4.3695e+00, 8.4204e-01 Phase(deg.) of the estimated response function: -132.2, -132.6 Weighted residual power: 0.00349775 Iteration number = 1 Scale factor: 0.0686162 Parameter c: 2.66852 Sum of weights: 453.093 Squared coherence: 0.840113 Estimated response function: ( -2.9349e+00, -3.2365e+00), ( -5.7423e-01, -6.3253e-01) Amplitude of the estimated response function: 4.3690e+00, 8.5430e-01 Phase(deg.) of the estimated response function: -132.2, -132.2 Weighted residual power: 0.00349658 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Estimate errors by fixed-weights bootstrap -------------------------------------------------------------------------------- ================================================================================ Now Frequency(Hz): 0.0234375, Period(s): 42.6667 ================================================================================ -------------------------------------------------------------------------------- Calculate response functions for output variable 0 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 1349 Squared coherence: 0.285299 Estimated response function: ( 1.8094e+00, 1.5365e+00), ( 1.7851e+00, 2.7188e+00) Amplitude of the estimated response function: 2.3738e+00, 3.2525e+00 Phase(deg.) of the estimated response function: 40.3, 56.7 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.115513 Sum of weights: 1347.73 Squared coherence: 0.288686 Estimated response function: ( 1.8099e+00, 1.5348e+00), ( 1.7857e+00, 2.7240e+00) Amplitude of the estimated response function: 2.3730e+00, 3.2572e+00 Phase(deg.) of the estimated response function: 40.3, 56.8 Weighted residual power: 0.0161098 Iteration number = 1 Scale factor: 0.115815 Sum of weights: 1347.76 Squared coherence: 0.288633 Estimated response function: ( 1.8098e+00, 1.5348e+00), ( 1.7857e+00, 2.7240e+00) Amplitude of the estimated response function: 2.3730e+00, 3.2572e+00 Phase(deg.) of the estimated response function: 40.3, 56.8 Weighted residual power: 0.0161126 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.115815 Parameter c: 2.66487 Sum of weights: 1007.23 Squared coherence: 0.522855 Estimated response function: ( 1.9083e+00, 1.5870e+00), ( 1.8507e+00, 2.7172e+00) Amplitude of the estimated response function: 2.4819e+00, 3.2876e+00 Phase(deg.) of the estimated response function: 39.7, 55.7 Weighted residual power: 0.0082682 Iteration number = 1 Scale factor: 0.115815 Parameter c: 2.66487 Sum of weights: 1005.03 Squared coherence: 0.526725 Estimated response function: ( 1.9289e+00, 1.5975e+00), ( 1.8642e+00, 2.7140e+00) Amplitude of the estimated response function: 2.5045e+00, 3.2926e+00 Phase(deg.) of the estimated response function: 39.6, 55.5 Weighted residual power: 0.0082418 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Calculate response functions for output variable 1 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 1349 Squared coherence: 0.541271 Estimated response function: ( -3.7694e+00, -3.9053e+00), ( -8.2834e-01, -9.5551e-01) Amplitude of the estimated response function: 5.4277e+00, 1.2646e+00 Phase(deg.) of the estimated response function: -134.0, -130.9 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.123226 Sum of weights: 1347.51 Squared coherence: 0.545925 Estimated response function: ( -3.7673e+00, -3.9069e+00), ( -8.2413e-01, -9.6010e-01) Amplitude of the estimated response function: 5.4274e+00, 1.2653e+00 Phase(deg.) of the estimated response function: -134.0, -130.6 Weighted residual power: 0.0190019 Iteration number = 1 Scale factor: 0.123251 Sum of weights: 1347.51 Squared coherence: 0.545922 Estimated response function: ( -3.7673e+00, -3.9069e+00), ( -8.2410e-01, -9.6011e-01) Amplitude of the estimated response function: 5.4274e+00, 1.2653e+00 Phase(deg.) of the estimated response function: -134.0, -130.6 Weighted residual power: 0.0190021 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.123251 Parameter c: 2.66487 Sum of weights: 1000.12 Squared coherence: 0.747648 Estimated response function: ( -3.7331e+00, -3.8562e+00), ( -8.1722e-01, -1.0234e+00) Amplitude of the estimated response function: 5.3671e+00, 1.3097e+00 Phase(deg.) of the estimated response function: -134.1, -128.6 Weighted residual power: 0.00904234 Iteration number = 1 Scale factor: 0.123251 Parameter c: 2.66487 Sum of weights: 1001.08 Squared coherence: 0.746773 Estimated response function: ( -3.7264e+00, -3.8462e+00), ( -8.1508e-01, -1.0360e+00) Amplitude of the estimated response function: 5.3553e+00, 1.3182e+00 Phase(deg.) of the estimated response function: -134.1, -128.2 Weighted residual power: 0.00905133 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Estimate errors by fixed-weights bootstrap -------------------------------------------------------------------------------- ================================================================================ Now Frequency(Hz): 0.03125, Period(s): 32 ================================================================================ -------------------------------------------------------------------------------- Calculate response functions for output variable 0 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 1349 Squared coherence: 0.182882 Estimated response function: ( 2.4038e+00, 1.7105e+00), ( 2.2720e+00, 3.4053e+00) Amplitude of the estimated response function: 2.9502e+00, 4.0937e+00 Phase(deg.) of the estimated response function: 35.4, 56.3 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.117692 Sum of weights: 1347.7 Squared coherence: 0.18659 Estimated response function: ( 2.3990e+00, 1.7126e+00), ( 2.2691e+00, 3.4123e+00) Amplitude of the estimated response function: 2.9476e+00, 4.0978e+00 Phase(deg.) of the estimated response function: 35.5, 56.4 Weighted residual power: 0.0178857 Iteration number = 1 Scale factor: 0.117671 Sum of weights: 1347.69 Squared coherence: 0.186601 Estimated response function: ( 2.3990e+00, 1.7126e+00), ( 2.2691e+00, 3.4123e+00) Amplitude of the estimated response function: 2.9476e+00, 4.0978e+00 Phase(deg.) of the estimated response function: 35.5, 56.4 Weighted residual power: 0.0178853 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.117671 Parameter c: 2.66487 Sum of weights: 989.787 Squared coherence: 0.467743 Estimated response function: ( 2.3768e+00, 1.7658e+00), ( 2.3306e+00, 3.4004e+00) Amplitude of the estimated response function: 2.9610e+00, 4.1224e+00 Phase(deg.) of the estimated response function: 36.6, 55.6 Weighted residual power: 0.00848306 Iteration number = 1 Scale factor: 0.117671 Parameter c: 2.66487 Sum of weights: 989.147 Squared coherence: 0.468887 Estimated response function: ( 2.3703e+00, 1.7797e+00), ( 2.3475e+00, 3.3996e+00) Amplitude of the estimated response function: 2.9640e+00, 4.1313e+00 Phase(deg.) of the estimated response function: 36.9, 55.4 Weighted residual power: 0.00847388 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Calculate response functions for output variable 1 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 1349 Squared coherence: 0.416897 Estimated response function: ( -4.5305e+00, -4.3722e+00), ( -1.0212e+00, -1.1243e+00) Amplitude of the estimated response function: 6.2961e+00, 1.5189e+00 Phase(deg.) of the estimated response function: -136.0, -132.2 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.131509 Sum of weights: 1348.11 Squared coherence: 0.419911 Estimated response function: ( -4.5371e+00, -4.3670e+00), ( -1.0142e+00, -1.1239e+00) Amplitude of the estimated response function: 6.2973e+00, 1.5139e+00 Phase(deg.) of the estimated response function: -136.1, -132.1 Weighted residual power: 0.0197379 Iteration number = 1 Scale factor: 0.131457 Sum of weights: 1348.11 Squared coherence: 0.41992 Estimated response function: ( -4.5371e+00, -4.3670e+00), ( -1.0141e+00, -1.1239e+00) Amplitude of the estimated response function: 6.2973e+00, 1.5138e+00 Phase(deg.) of the estimated response function: -136.1, -132.1 Weighted residual power: 0.0197375 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.131457 Parameter c: 2.66487 Sum of weights: 1019.39 Squared coherence: 0.653343 Estimated response function: ( -4.4342e+00, -4.3396e+00), ( -9.1178e-01, -1.0439e+00) Amplitude of the estimated response function: 6.2044e+00, 1.3861e+00 Phase(deg.) of the estimated response function: -135.6, -131.1 Weighted residual power: 0.00995038 Iteration number = 1 Scale factor: 0.131457 Parameter c: 2.66487 Sum of weights: 1022.78 Squared coherence: 0.651157 Estimated response function: ( -4.4130e+00, -4.3338e+00), ( -8.9339e-01, -1.0216e+00) Amplitude of the estimated response function: 6.1852e+00, 1.3571e+00 Phase(deg.) of the estimated response function: -135.5, -131.2 Weighted residual power: 0.00998484 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Estimate errors by fixed-weights bootstrap -------------------------------------------------------------------------------- ================================================================================ Now Frequency(Hz): 0.046875, Period(s): 21.3333 ================================================================================ -------------------------------------------------------------------------------- Calculate response functions for output variable 0 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 2699 Squared coherence: 0 Estimated response function: ( 3.4951e+00, 2.1127e+00), ( 2.7520e+00, 4.2758e+00) Amplitude of the estimated response function: 4.0841e+00, 5.0849e+00 Phase(deg.) of the estimated response function: 31.2, 57.2 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.190613 Sum of weights: 2696.85 Squared coherence: 0 Estimated response function: ( 3.4905e+00, 2.1061e+00), ( 2.7290e+00, 4.2591e+00) Amplitude of the estimated response function: 4.0767e+00, 5.0584e+00 Phase(deg.) of the estimated response function: 31.1, 57.4 Weighted residual power: 0.0415517 Iteration number = 1 Scale factor: 0.190219 Sum of weights: 2696.83 Squared coherence: 0 Estimated response function: ( 3.4905e+00, 2.1061e+00), ( 2.7287e+00, 4.2589e+00) Amplitude of the estimated response function: 4.0767e+00, 5.0581e+00 Phase(deg.) of the estimated response function: 31.1, 57.4 Weighted residual power: 0.0415476 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.190219 Parameter c: 2.66305 Sum of weights: 2039.99 Squared coherence: 0.25142 Estimated response function: ( 3.4017e+00, 2.0860e+00), ( 2.6299e+00, 4.2306e+00) Amplitude of the estimated response function: 3.9903e+00, 4.9814e+00 Phase(deg.) of the estimated response function: 31.5, 58.1 Weighted residual power: 0.0203781 Iteration number = 1 Scale factor: 0.190219 Parameter c: 2.66305 Sum of weights: 2046.23 Squared coherence: 0.246995 Estimated response function: ( 3.3798e+00, 2.0824e+00), ( 2.6063e+00, 4.2245e+00) Amplitude of the estimated response function: 3.9698e+00, 4.9638e+00 Phase(deg.) of the estimated response function: 31.6, 58.3 Weighted residual power: 0.0204523 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Calculate response functions for output variable 1 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 2699 Squared coherence: 0.137491 Estimated response function: ( -5.7019e+00, -4.8532e+00), ( -1.0653e+00, -1.0137e+00) Amplitude of the estimated response function: 7.4877e+00, 1.4705e+00 Phase(deg.) of the estimated response function: -139.6, -136.4 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.208203 Sum of weights: 2697.61 Squared coherence: 0.141308 Estimated response function: ( -5.6961e+00, -4.8610e+00), ( -1.0740e+00, -1.0233e+00) Amplitude of the estimated response function: 7.4883e+00, 1.4834e+00 Phase(deg.) of the estimated response function: -139.5, -136.4 Weighted residual power: 0.0454558 Iteration number = 1 Scale factor: 0.208441 Sum of weights: 2697.62 Squared coherence: 0.141288 Estimated response function: ( -5.6961e+00, -4.8610e+00), ( -1.0740e+00, -1.0233e+00) Amplitude of the estimated response function: 7.4883e+00, 1.4834e+00 Phase(deg.) of the estimated response function: -139.5, -136.4 Weighted residual power: 0.0454571 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.208441 Parameter c: 2.66305 Sum of weights: 2082.6 Squared coherence: 0.444755 Estimated response function: ( -5.6530e+00, -4.9602e+00), ( -8.8944e-01, -1.1365e+00) Amplitude of the estimated response function: 7.5206e+00, 1.4432e+00 Phase(deg.) of the estimated response function: -138.7, -128.0 Weighted residual power: 0.0238622 Iteration number = 1 Scale factor: 0.208441 Parameter c: 2.66305 Sum of weights: 2081.68 Squared coherence: 0.44573 Estimated response function: ( -5.6433e+00, -4.9820e+00), ( -8.4906e-01, -1.1529e+00) Amplitude of the estimated response function: 7.5278e+00, 1.4318e+00 Phase(deg.) of the estimated response function: -138.6, -126.4 Weighted residual power: 0.0238195 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Estimate errors by fixed-weights bootstrap -------------------------------------------------------------------------------- ================================================================================ Now Frequency(Hz): 0.0625, Period(s): 16 ================================================================================ -------------------------------------------------------------------------------- Calculate response functions for output variable 0 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 2699 Squared coherence: 0 Estimated response function: ( 3.8661e+00, 2.4595e+00), ( 3.7929e+00, 5.6612e+00) Amplitude of the estimated response function: 4.5821e+00, 6.8144e+00 Phase(deg.) of the estimated response function: 32.5, 56.2 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.214947 Sum of weights: 2696.72 Squared coherence: 0 Estimated response function: ( 3.8557e+00, 2.4604e+00), ( 3.7434e+00, 5.6688e+00) Amplitude of the estimated response function: 4.5738e+00, 6.7933e+00 Phase(deg.) of the estimated response function: 32.5, 56.6 Weighted residual power: 0.049531 Iteration number = 1 Scale factor: 0.21527 Sum of weights: 2696.76 Squared coherence: 0 Estimated response function: ( 3.8559e+00, 2.4600e+00), ( 3.7437e+00, 5.6683e+00) Amplitude of the estimated response function: 4.5738e+00, 6.7930e+00 Phase(deg.) of the estimated response function: 32.5, 56.6 Weighted residual power: 0.0495366 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.21527 Parameter c: 2.66305 Sum of weights: 2080.02 Squared coherence: 0.0308249 Estimated response function: ( 3.9277e+00, 2.4448e+00), ( 3.2823e+00, 5.5237e+00) Amplitude of the estimated response function: 4.6264e+00, 6.4254e+00 Phase(deg.) of the estimated response function: 31.9, 59.3 Weighted residual power: 0.0247656 Iteration number = 1 Scale factor: 0.21527 Parameter c: 2.66305 Sum of weights: 2092.65 Squared coherence: 0.0208152 Estimated response function: ( 3.9482e+00, 2.4427e+00), ( 3.1831e+00, 5.4907e+00) Amplitude of the estimated response function: 4.6427e+00, 6.3466e+00 Phase(deg.) of the estimated response function: 31.7, 59.9 Weighted residual power: 0.0249937 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Calculate response functions for output variable 1 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 2699 Squared coherence: 0 Estimated response function: ( -5.7741e+00, -5.9493e+00), ( -1.3897e+00, -1.3007e+00) Amplitude of the estimated response function: 8.2906e+00, 1.9034e+00 Phase(deg.) of the estimated response function: -134.1, -136.9 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.219695 Sum of weights: 2698.47 Squared coherence: 0 Estimated response function: ( -5.7752e+00, -5.9441e+00), ( -1.4083e+00, -1.3112e+00) Amplitude of the estimated response function: 8.2876e+00, 1.9242e+00 Phase(deg.) of the estimated response function: -134.2, -137.0 Weighted residual power: 0.0498645 Iteration number = 1 Scale factor: 0.219928 Sum of weights: 2698.48 Squared coherence: 0 Estimated response function: ( -5.7751e+00, -5.9442e+00), ( -1.4081e+00, -1.3112e+00) Amplitude of the estimated response function: 8.2877e+00, 1.9241e+00 Phase(deg.) of the estimated response function: -134.2, -137.0 Weighted residual power: 0.0498661 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.219928 Parameter c: 2.66305 Sum of weights: 2091.2 Squared coherence: 0.248979 Estimated response function: ( -6.0605e+00, -5.8105e+00), ( -1.4275e+00, -1.3844e+00) Amplitude of the estimated response function: 8.3959e+00, 1.9886e+00 Phase(deg.) of the estimated response function: -136.2, -135.9 Weighted residual power: 0.0263899 Iteration number = 1 Scale factor: 0.219928 Parameter c: 2.66305 Sum of weights: 2085.17 Squared coherence: 0.254814 Estimated response function: ( -6.1214e+00, -5.7715e+00), ( -1.4449e+00, -1.3872e+00) Amplitude of the estimated response function: 8.4132e+00, 2.0030e+00 Phase(deg.) of the estimated response function: -136.7, -136.2 Weighted residual power: 0.0262236 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Estimate errors by fixed-weights bootstrap -------------------------------------------------------------------------------- ================================================================================ Now Frequency(Hz): 0.09375, Period(s): 10.6667 ================================================================================ -------------------------------------------------------------------------------- Calculate response functions for output variable 0 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 5399 Squared coherence: 0 Estimated response function: ( 3.9019e+00, 2.2656e+00), ( 4.4292e+00, 7.3770e+00) Amplitude of the estimated response function: 4.5120e+00, 8.6045e+00 Phase(deg.) of the estimated response function: 30.1, 59.0 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.3301 Sum of weights: 5394.64 Squared coherence: 0 Estimated response function: ( 3.8949e+00, 2.2811e+00), ( 4.3922e+00, 7.3913e+00) Amplitude of the estimated response function: 4.5137e+00, 8.5979e+00 Phase(deg.) of the estimated response function: 30.4, 59.3 Weighted residual power: 0.117904 Iteration number = 1 Scale factor: 0.32998 Sum of weights: 5394.63 Squared coherence: 0 Estimated response function: ( 3.8950e+00, 2.2811e+00), ( 4.3922e+00, 7.3915e+00) Amplitude of the estimated response function: 4.5138e+00, 8.5980e+00 Phase(deg.) of the estimated response function: 30.4, 59.3 Weighted residual power: 0.117905 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.32998 Parameter c: 2.66214 Sum of weights: 4164.24 Squared coherence: 0 Estimated response function: ( 3.8730e+00, 2.1104e+00), ( 4.5287e+00, 7.1763e+00) Amplitude of the estimated response function: 4.4107e+00, 8.4858e+00 Phase(deg.) of the estimated response function: 28.6, 57.7 Weighted residual power: 0.0551405 Iteration number = 1 Scale factor: 0.32998 Parameter c: 2.66214 Sum of weights: 4179.34 Squared coherence: 0 Estimated response function: ( 3.8739e+00, 2.0798e+00), ( 4.5622e+00, 7.1418e+00) Amplitude of the estimated response function: 4.3969e+00, 8.4746e+00 Phase(deg.) of the estimated response function: 28.2, 57.4 Weighted residual power: 0.0556177 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Calculate response functions for output variable 1 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 5399 Squared coherence: 0 Estimated response function: ( -6.8066e+00, -7.2849e+00), ( -4.3713e-01, -2.5180e+00) Amplitude of the estimated response function: 9.9699e+00, 2.5557e+00 Phase(deg.) of the estimated response function: -133.1, -99.8 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.304303 Sum of weights: 5391.17 Squared coherence: 0 Estimated response function: ( -6.7845e+00, -7.3230e+00), ( -4.5213e-01, -2.6111e+00) Amplitude of the estimated response function: 9.9827e+00, 2.6500e+00 Phase(deg.) of the estimated response function: -132.8, -99.8 Weighted residual power: 0.121071 Iteration number = 1 Scale factor: 0.305432 Sum of weights: 5391.25 Squared coherence: 0 Estimated response function: ( -6.7849e+00, -7.3232e+00), ( -4.5196e-01, -2.6117e+00) Amplitude of the estimated response function: 9.9832e+00, 2.6505e+00 Phase(deg.) of the estimated response function: -132.8, -99.8 Weighted residual power: 0.121094 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.305432 Parameter c: 2.66214 Sum of weights: 3996.97 Squared coherence: 0.0985905 Estimated response function: ( -6.8657e+00, -6.9912e+00), ( -1.3267e+00, -2.5057e+00) Amplitude of the estimated response function: 9.7987e+00, 2.8352e+00 Phase(deg.) of the estimated response function: -134.5, -117.9 Weighted residual power: 0.0485671 Iteration number = 1 Scale factor: 0.305432 Parameter c: 2.66214 Sum of weights: 4013.31 Squared coherence: 0.0861207 Estimated response function: ( -6.9019e+00, -6.9385e+00), ( -1.5094e+00, -2.4943e+00) Amplitude of the estimated response function: 9.7866e+00, 2.9154e+00 Phase(deg.) of the estimated response function: -134.8, -121.2 Weighted residual power: 0.0492692 Iteration number = 2 Scale factor: 0.305432 Parameter c: 2.66214 Sum of weights: 4012.28 Squared coherence: 0.0853595 Estimated response function: ( -6.9144e+00, -6.9310e+00), ( -1.5498e+00, -2.4941e+00) Amplitude of the estimated response function: 9.7902e+00, 2.9364e+00 Phase(deg.) of the estimated response function: -134.9, -121.9 Weighted residual power: 0.0493184 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Estimate errors by fixed-weights bootstrap -------------------------------------------------------------------------------- ================================================================================ Now Frequency(Hz): 0.125, Period(s): 8 ================================================================================ -------------------------------------------------------------------------------- Calculate response functions for output variable 0 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 5399 Squared coherence: 0 Estimated response function: ( 5.2436e+00, 1.8981e+00), ( 5.8367e+00, 7.6169e+00) Amplitude of the estimated response function: 5.5766e+00, 9.5961e+00 Phase(deg.) of the estimated response function: 19.9, 52.5 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.361732 Sum of weights: 5397.06 Squared coherence: 0 Estimated response function: ( 5.2345e+00, 1.8737e+00), ( 5.8302e+00, 7.6385e+00) Amplitude of the estimated response function: 5.5598e+00, 9.6092e+00 Phase(deg.) of the estimated response function: 19.7, 52.6 Weighted residual power: 0.136418 Iteration number = 1 Scale factor: 0.362392 Sum of weights: 5397.09 Squared coherence: 0 Estimated response function: ( 5.2343e+00, 1.8739e+00), ( 5.8296e+00, 7.6385e+00) Amplitude of the estimated response function: 5.5596e+00, 9.6089e+00 Phase(deg.) of the estimated response function: 19.7, 52.6 Weighted residual power: 0.136419 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.362392 Parameter c: 2.66214 Sum of weights: 4190.59 Squared coherence: 0 Estimated response function: ( 4.7311e+00, 1.9972e+00), ( 5.8180e+00, 7.9060e+00) Amplitude of the estimated response function: 5.1354e+00, 9.8160e+00 Phase(deg.) of the estimated response function: 22.9, 53.7 Weighted residual power: 0.0674299 Iteration number = 1 Scale factor: 0.362392 Parameter c: 2.66214 Sum of weights: 4193.77 Squared coherence: 0 Estimated response function: ( 4.6529e+00, 2.0178e+00), ( 5.7916e+00, 7.9785e+00) Amplitude of the estimated response function: 5.0716e+00, 9.8589e+00 Phase(deg.) of the estimated response function: 23.4, 54.0 Weighted residual power: 0.0674119 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Calculate response functions for output variable 1 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 5399 Squared coherence: 0 Estimated response function: ( -9.2962e+00, -8.5980e+00), ( 2.5354e-02, -1.6449e+00) Amplitude of the estimated response function: 1.2663e+01, 1.6451e+00 Phase(deg.) of the estimated response function: -137.2, -89.1 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.34741 Sum of weights: 5390.19 Squared coherence: 0 Estimated response function: ( -9.2735e+00, -8.6384e+00), ( -3.3656e-02, -1.5436e+00) Amplitude of the estimated response function: 1.2674e+01, 1.5440e+00 Phase(deg.) of the estimated response function: -137.0, -91.2 Weighted residual power: 0.156267 Iteration number = 1 Scale factor: 0.348718 Sum of weights: 5390.42 Squared coherence: 0 Estimated response function: ( -9.2731e+00, -8.6408e+00), ( -3.3343e-02, -1.5414e+00) Amplitude of the estimated response function: 1.2675e+01, 1.5418e+00 Phase(deg.) of the estimated response function: -137.0, -91.2 Weighted residual power: 0.156339 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.348718 Parameter c: 2.66214 Sum of weights: 4025.8 Squared coherence: 0 Estimated response function: ( -8.1755e+00, -8.0965e+00), ( -7.7671e-01, -2.2409e+00) Amplitude of the estimated response function: 1.1506e+01, 2.3717e+00 Phase(deg.) of the estimated response function: -135.3, -109.1 Weighted residual power: 0.0570727 Iteration number = 1 Scale factor: 0.348718 Parameter c: 2.66214 Sum of weights: 4121.07 Squared coherence: 0 Estimated response function: ( -8.0524e+00, -7.9570e+00), ( -9.3388e-01, -2.3522e+00) Amplitude of the estimated response function: 1.1321e+01, 2.5308e+00 Phase(deg.) of the estimated response function: -135.3, -111.7 Weighted residual power: 0.0604893 Iteration number = 2 Scale factor: 0.348718 Parameter c: 2.66214 Sum of weights: 4134.36 Squared coherence: 0 Estimated response function: ( -8.0428e+00, -7.9240e+00), ( -9.6379e-01, -2.3677e+00) Amplitude of the estimated response function: 1.1291e+01, 2.5563e+00 Phase(deg.) of the estimated response function: -135.4, -112.1 Weighted residual power: 0.0610442 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Estimate errors by fixed-weights bootstrap -------------------------------------------------------------------------------- ================================================================================ Now Frequency(Hz): 0.1875, Period(s): 5.33333 ================================================================================ -------------------------------------------------------------------------------- Calculate response functions for output variable 0 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 10799 Squared coherence: 0 Estimated response function: ( 4.7388e+00, 3.6659e+00), ( 8.5635e+00, 9.6804e+00) Amplitude of the estimated response function: 5.9912e+00, 1.2925e+01 Phase(deg.) of the estimated response function: 37.7, 48.5 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.502386 Sum of weights: 10757.9 Squared coherence: 0 Estimated response function: ( 4.6041e+00, 3.4384e+00), ( 8.5797e+00, 9.8077e+00) Amplitude of the estimated response function: 5.7463e+00, 1.3031e+01 Phase(deg.) of the estimated response function: 36.8, 48.8 Weighted residual power: 0.372829 Iteration number = 1 Scale factor: 0.502108 Sum of weights: 10757.4 Squared coherence: 0 Estimated response function: ( 4.6032e+00, 3.4316e+00), ( 8.5797e+00, 9.8081e+00) Amplitude of the estimated response function: 5.7415e+00, 1.3031e+01 Phase(deg.) of the estimated response function: 36.7, 48.8 Weighted residual power: 0.372609 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.502108 Parameter c: 2.66169 Sum of weights: 7881.37 Squared coherence: 0 Estimated response function: ( 4.4361e+00, 2.6423e+00), ( 7.6949e+00, 1.0975e+01) Amplitude of the estimated response function: 5.1634e+00, 1.3404e+01 Phase(deg.) of the estimated response function: 30.8, 55.0 Weighted residual power: 0.126534 Iteration number = 1 Scale factor: 0.502108 Parameter c: 2.66169 Sum of weights: 7871.52 Squared coherence: 0 Estimated response function: ( 4.4089e+00, 2.5843e+00), ( 7.4182e+00, 1.1148e+01) Amplitude of the estimated response function: 5.1105e+00, 1.3391e+01 Phase(deg.) of the estimated response function: 30.4, 56.4 Weighted residual power: 0.124191 Iteration number = 2 Scale factor: 0.502108 Parameter c: 2.66169 Sum of weights: 7877.21 Squared coherence: 0 Estimated response function: ( 4.3979e+00, 2.5802e+00), ( 7.3493e+00, 1.1174e+01) Amplitude of the estimated response function: 5.0989e+00, 1.3374e+01 Phase(deg.) of the estimated response function: 30.4, 56.7 Weighted residual power: 0.124257 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Calculate response functions for output variable 1 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 10799 Squared coherence: 0 Estimated response function: ( -1.1045e+01, -1.1288e+01), ( -2.2411e-01, -5.9487e-01) Amplitude of the estimated response function: 1.5793e+01, 6.3568e-01 Phase(deg.) of the estimated response function: -134.4, -110.6 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.478591 Sum of weights: 10711.5 Squared coherence: 0 Estimated response function: ( -1.0700e+01, -1.0857e+01), ( -2.9335e-01, -3.0862e-01) Amplitude of the estimated response function: 1.5243e+01, 4.2580e-01 Phase(deg.) of the estimated response function: -134.6, -133.5 Weighted residual power: 0.386061 Iteration number = 1 Scale factor: 0.469278 Sum of weights: 10715.1 Squared coherence: 0 Estimated response function: ( -1.0695e+01, -1.0861e+01), ( -2.9250e-01, -3.1651e-01) Amplitude of the estimated response function: 1.5243e+01, 4.3097e-01 Phase(deg.) of the estimated response function: -134.6, -132.7 Weighted residual power: 0.386992 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.469278 Parameter c: 2.66169 Sum of weights: 7624.59 Squared coherence: 0 Estimated response function: ( -9.5926e+00, -1.0178e+01), ( -1.4037e-01, -9.6849e-01) Amplitude of the estimated response function: 1.3986e+01, 9.7861e-01 Phase(deg.) of the estimated response function: -133.3, -98.2 Weighted residual power: 0.0962821 Iteration number = 1 Scale factor: 0.469278 Parameter c: 2.66169 Sum of weights: 7821.6 Squared coherence: 0 Estimated response function: ( -9.4018e+00, -1.0021e+01), ( 3.1420e-02, -8.2711e-01) Amplitude of the estimated response function: 1.3741e+01, 8.2770e-01 Phase(deg.) of the estimated response function: -133.2, -87.8 Weighted residual power: 0.103991 Iteration number = 2 Scale factor: 0.469278 Parameter c: 2.66169 Sum of weights: 7862.78 Squared coherence: 0 Estimated response function: ( -9.3796e+00, -1.0002e+01), ( 8.9541e-02, -7.3514e-01) Amplitude of the estimated response function: 1.3712e+01, 7.4057e-01 Phase(deg.) of the estimated response function: -133.2, -83.1 Weighted residual power: 0.105888 Iteration number = 3 Scale factor: 0.469278 Parameter c: 2.66169 Sum of weights: 7868.34 Squared coherence: 0 Estimated response function: ( -9.3809e+00, -1.0003e+01), ( 1.0283e-01, -7.0506e-01) Amplitude of the estimated response function: 1.3714e+01, 7.1252e-01 Phase(deg.) of the estimated response function: -133.2, -81.7 Weighted residual power: 0.106198 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Estimate errors by fixed-weights bootstrap -------------------------------------------------------------------------------- ================================================================================ Now Frequency(Hz): 0.25, Period(s): 4 ================================================================================ -------------------------------------------------------------------------------- Calculate response functions for output variable 0 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 10799 Squared coherence: 0 Estimated response function: ( 5.8663e+00, 4.0237e+00), ( 1.1449e+01, 8.7057e+00) Amplitude of the estimated response function: 7.1136e+00, 1.4383e+01 Phase(deg.) of the estimated response function: 34.4, 37.2 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.542211 Sum of weights: 10757.9 Squared coherence: 0 Estimated response function: ( 6.1989e+00, 3.9568e+00), ( 1.1643e+01, 8.5442e+00) Amplitude of the estimated response function: 7.3541e+00, 1.4442e+01 Phase(deg.) of the estimated response function: 32.6, 36.3 Weighted residual power: 0.445821 Iteration number = 1 Scale factor: 0.549359 Sum of weights: 10759.4 Squared coherence: 0 Estimated response function: ( 6.2000e+00, 3.9423e+00), ( 1.1628e+01, 8.5284e+00) Amplitude of the estimated response function: 7.3473e+00, 1.4420e+01 Phase(deg.) of the estimated response function: 32.5, 36.3 Weighted residual power: 0.445345 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.549359 Parameter c: 2.66169 Sum of weights: 7867.71 Squared coherence: 0 Estimated response function: ( 3.6134e+00, 5.8625e+00), ( 9.3774e+00, 1.2662e+01) Amplitude of the estimated response function: 6.8866e+00, 1.5756e+01 Phase(deg.) of the estimated response function: 58.4, 53.5 Weighted residual power: 0.169381 Iteration number = 1 Scale factor: 0.549359 Parameter c: 2.66169 Sum of weights: 7723.02 Squared coherence: 0 Estimated response function: ( 3.0765e+00, 6.4204e+00), ( 8.6928e+00, 1.3783e+01) Amplitude of the estimated response function: 7.1195e+00, 1.6295e+01 Phase(deg.) of the estimated response function: 64.4, 57.8 Weighted residual power: 0.162365 Iteration number = 2 Scale factor: 0.549359 Parameter c: 2.66169 Sum of weights: 7630.66 Squared coherence: 0 Estimated response function: ( 2.9543e+00, 6.5618e+00), ( 8.4242e+00, 1.4246e+01) Amplitude of the estimated response function: 7.1962e+00, 1.6550e+01 Phase(deg.) of the estimated response function: 65.8, 59.4 Weighted residual power: 0.158462 Iteration number = 3 Scale factor: 0.549359 Parameter c: 2.66169 Sum of weights: 7590.53 Squared coherence: 0 Estimated response function: ( 2.9383e+00, 6.5768e+00), ( 8.3066e+00, 1.4406e+01) Amplitude of the estimated response function: 7.2033e+00, 1.6629e+01 Phase(deg.) of the estimated response function: 65.9, 60.0 Weighted residual power: 0.156174 Iteration number = 4 Scale factor: 0.549359 Parameter c: 2.66169 Sum of weights: 7579.46 Squared coherence: 0 Estimated response function: ( 2.9426e+00, 6.5689e+00), ( 8.2604e+00, 1.4446e+01) Amplitude of the estimated response function: 7.1979e+00, 1.6641e+01 Phase(deg.) of the estimated response function: 65.9, 60.2 Weighted residual power: 0.15535 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Calculate response functions for output variable 1 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 10799 Squared coherence: 0 Estimated response function: ( -1.5791e+01, -8.7705e+00), ( -2.1661e+00, 6.6073e-01) Amplitude of the estimated response function: 1.8063e+01, 2.2646e+00 Phase(deg.) of the estimated response function: -151.0, 163.0 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.537948 Sum of weights: 10714.4 Squared coherence: 0 Estimated response function: ( -1.4965e+01, -9.4064e+00), ( -1.6618e+00, -2.1359e-01) Amplitude of the estimated response function: 1.7675e+01, 1.6754e+00 Phase(deg.) of the estimated response function: -147.8, -172.7 Weighted residual power: 0.490929 Iteration number = 1 Scale factor: 0.527053 Sum of weights: 10714.2 Squared coherence: 0 Estimated response function: ( -1.4927e+01, -9.3861e+00), ( -1.5579e+00, -1.9915e-01) Amplitude of the estimated response function: 1.7633e+01, 1.5705e+00 Phase(deg.) of the estimated response function: -147.8, -172.7 Weighted residual power: 0.488524 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.527053 Parameter c: 2.66169 Sum of weights: 7617.27 Squared coherence: 0 Estimated response function: ( -1.3752e+01, -1.0902e+01), ( -1.3903e+00, -4.5545e-01) Amplitude of the estimated response function: 1.7549e+01, 1.4630e+00 Phase(deg.) of the estimated response function: -141.6, -161.9 Weighted residual power: 0.134292 Iteration number = 1 Scale factor: 0.527053 Parameter c: 2.66169 Sum of weights: 7639.36 Squared coherence: 0 Estimated response function: ( -1.3453e+01, -1.1256e+01), ( -1.2581e+00, -3.2005e-01) Amplitude of the estimated response function: 1.7541e+01, 1.2982e+00 Phase(deg.) of the estimated response function: -140.1, -165.7 Weighted residual power: 0.134753 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Estimate errors by fixed-weights bootstrap -------------------------------------------------------------------------------- ================================================================================ Now Frequency(Hz): 0.375, Period(s): 2.66667 ================================================================================ -------------------------------------------------------------------------------- Calculate response functions for output variable 0 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 21599 Squared coherence: 0 Estimated response function: ( -3.4466e+00, -5.2926e+00), ( 2.3350e+01, 8.4248e+00) Amplitude of the estimated response function: 6.3159e+00, 2.4823e+01 Phase(deg.) of the estimated response function: -123.1, 19.8 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.901226 Sum of weights: 21306.4 Squared coherence: 0 Estimated response function: ( -3.2957e-01, -5.7979e+00), ( 2.4154e+01, 1.1911e+01) Amplitude of the estimated response function: 5.8072e+00, 2.6931e+01 Phase(deg.) of the estimated response function: -93.3, 26.2 Weighted residual power: 1.98007 Iteration number = 1 Scale factor: 0.942643 Sum of weights: 21285.8 Squared coherence: 0 Estimated response function: ( -3.0956e-02, -5.6683e+00), ( 2.3931e+01, 1.2061e+01) Amplitude of the estimated response function: 5.6684e+00, 2.6799e+01 Phase(deg.) of the estimated response function: -90.3, 26.7 Weighted residual power: 1.94338 Iteration number = 2 Scale factor: 0.938307 Sum of weights: 21286.5 Squared coherence: 0 Estimated response function: ( -1.7188e-02, -5.6557e+00), ( 2.3901e+01, 1.2075e+01) Amplitude of the estimated response function: 5.6557e+00, 2.6777e+01 Phase(deg.) of the estimated response function: -90.2, 26.8 Weighted residual power: 1.94071 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.938307 Parameter c: 2.66146 Sum of weights: 14548.6 Squared coherence: 0 Estimated response function: ( 6.1592e+00, -5.2026e-01), ( 1.3440e+01, 1.8263e+01) Amplitude of the estimated response function: 6.1811e+00, 2.2675e+01 Phase(deg.) of the estimated response function: -4.8, 53.6 Weighted residual power: 0.360629 Iteration number = 1 Scale factor: 0.938307 Parameter c: 2.66146 Sum of weights: 15382.8 Squared coherence: 0 Estimated response function: ( 5.7045e+00, 9.2482e-01), ( 1.3924e+01, 1.6274e+01) Amplitude of the estimated response function: 5.7790e+00, 2.1418e+01 Phase(deg.) of the estimated response function: 9.2, 49.4 Weighted residual power: 0.378927 Iteration number = 2 Scale factor: 0.938307 Parameter c: 2.66146 Sum of weights: 15703.7 Squared coherence: 0 Estimated response function: ( 6.1688e+00, 1.4135e+00), ( 1.4146e+01, 1.6292e+01) Amplitude of the estimated response function: 6.3286e+00, 2.1577e+01 Phase(deg.) of the estimated response function: 12.9, 49.0 Weighted residual power: 0.424679 Iteration number = 3 Scale factor: 0.938307 Parameter c: 2.66146 Sum of weights: 15617.3 Squared coherence: 0 Estimated response function: ( 6.3659e+00, 1.7619e+00), ( 1.4129e+01, 1.6405e+01) Amplitude of the estimated response function: 6.6052e+00, 2.1651e+01 Phase(deg.) of the estimated response function: 15.5, 49.3 Weighted residual power: 0.420139 Iteration number = 4 Scale factor: 0.938307 Parameter c: 2.66146 Sum of weights: 15574.6 Squared coherence: 0 Estimated response function: ( 6.4454e+00, 1.9647e+00), ( 1.4078e+01, 1.6482e+01) Amplitude of the estimated response function: 6.7381e+00, 2.1676e+01 Phase(deg.) of the estimated response function: 17.0, 49.5 Weighted residual power: 0.417482 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Calculate response functions for output variable 1 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 21599 Squared coherence: 0 Estimated response function: ( -1.3650e+01, -3.9352e+00), ( -4.6476e-02, -1.1209e+01) Amplitude of the estimated response function: 1.4206e+01, 1.1209e+01 Phase(deg.) of the estimated response function: -163.9, -90.2 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.697875 Sum of weights: 21353.1 Squared coherence: 0 Estimated response function: ( -1.3453e+01, -4.7179e+00), ( -9.4150e-01, -1.2448e+01) Amplitude of the estimated response function: 1.4256e+01, 1.2484e+01 Phase(deg.) of the estimated response function: -160.7, -94.3 Weighted residual power: 1.07156 Iteration number = 1 Scale factor: 0.729792 Sum of weights: 21357.2 Squared coherence: 0 Estimated response function: ( -1.3539e+01, -4.6822e+00), ( -8.1343e-01, -1.2587e+01) Amplitude of the estimated response function: 1.4326e+01, 1.2613e+01 Phase(deg.) of the estimated response function: -160.9, -93.7 Weighted residual power: 1.08562 Iteration number = 2 Scale factor: 0.731445 Sum of weights: 21353.4 Squared coherence: 0 Estimated response function: ( -1.3548e+01, -4.6895e+00), ( -7.9792e-01, -1.2619e+01) Amplitude of the estimated response function: 1.4337e+01, 1.2644e+01 Phase(deg.) of the estimated response function: -160.9, -93.6 Weighted residual power: 1.08713 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.731445 Parameter c: 2.66146 Sum of weights: 14750.6 Squared coherence: 0 Estimated response function: ( -1.4603e+01, -1.2206e+01), ( 2.1207e+00, -5.3262e+00) Amplitude of the estimated response function: 1.9032e+01, 5.7328e+00 Phase(deg.) of the estimated response function: -140.1, -68.3 Weighted residual power: 0.329703 Iteration number = 1 Scale factor: 0.731445 Parameter c: 2.66146 Sum of weights: 14740.9 Squared coherence: 0 Estimated response function: ( -1.1163e+01, -1.3313e+01), ( -6.8935e-01, -2.9053e+00) Amplitude of the estimated response function: 1.7373e+01, 2.9860e+00 Phase(deg.) of the estimated response function: -130.0, -103.3 Weighted residual power: 0.220446 Iteration number = 2 Scale factor: 0.731445 Parameter c: 2.66146 Sum of weights: 15448.6 Squared coherence: 0 Estimated response function: ( -1.0483e+01, -1.2945e+01), ( -2.2372e+00, -3.3493e+00) Amplitude of the estimated response function: 1.6657e+01, 4.0278e+00 Phase(deg.) of the estimated response function: -129.0, -123.7 Weighted residual power: 0.248637 Iteration number = 3 Scale factor: 0.731445 Parameter c: 2.66146 Sum of weights: 15555.1 Squared coherence: 0 Estimated response function: ( -1.0766e+01, -1.2790e+01), ( -2.3568e+00, -3.8572e+00) Amplitude of the estimated response function: 1.6718e+01, 4.5202e+00 Phase(deg.) of the estimated response function: -130.1, -121.4 Weighted residual power: 0.26058 Iteration number = 4 Scale factor: 0.731445 Parameter c: 2.66146 Sum of weights: 15486.1 Squared coherence: 0 Estimated response function: ( -1.0888e+01, -1.2874e+01), ( -2.1337e+00, -3.8963e+00) Amplitude of the estimated response function: 1.6860e+01, 4.4422e+00 Phase(deg.) of the estimated response function: -130.2, -118.7 Weighted residual power: 0.257918 Iteration number = 5 Scale factor: 0.731445 Parameter c: 2.66146 Sum of weights: 15455.4 Squared coherence: 0 Estimated response function: ( -1.0861e+01, -1.2923e+01), ( -2.0779e+00, -3.8331e+00) Amplitude of the estimated response function: 1.6881e+01, 4.3601e+00 Phase(deg.) of the estimated response function: -130.0, -118.5 Weighted residual power: 0.25523 Iteration number = 6 Scale factor: 0.731445 Parameter c: 2.66146 Sum of weights: 15458.6 Squared coherence: 0 Estimated response function: ( -1.0833e+01, -1.2927e+01), ( -2.0986e+00, -3.8127e+00) Amplitude of the estimated response function: 1.6866e+01, 4.3521e+00 Phase(deg.) of the estimated response function: -130.0, -118.8 Weighted residual power: 0.255003 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Estimate errors by fixed-weights bootstrap -------------------------------------------------------------------------------- ================================================================================ Now Frequency(Hz): 0.5, Period(s): 2 ================================================================================ -------------------------------------------------------------------------------- Calculate response functions for output variable 0 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 21599 Squared coherence: 0 Estimated response function: ( 1.2138e+00, -1.2061e+00), ( 2.1695e+01, 6.6378e+00) Amplitude of the estimated response function: 1.7112e+00, 2.2687e+01 Phase(deg.) of the estimated response function: -44.8, 17.0 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.777726 Sum of weights: 21266.4 Squared coherence: 0 Estimated response function: ( 6.1250e+00, -2.1352e+00), ( 2.5598e+01, 9.8701e+00) Amplitude of the estimated response function: 6.4865e+00, 2.7435e+01 Phase(deg.) of the estimated response function: -19.2, 21.1 Weighted residual power: 1.99314 Iteration number = 1 Scale factor: 0.938241 Sum of weights: 21267.9 Squared coherence: 0 Estimated response function: ( 6.4251e+00, -2.3332e+00), ( 2.6330e+01, 1.0466e+01) Amplitude of the estimated response function: 6.8356e+00, 2.8333e+01 Phase(deg.) of the estimated response function: -20.0, 21.7 Weighted residual power: 2.11245 Iteration number = 2 Scale factor: 0.9598 Sum of weights: 21258.9 Squared coherence: 0 Estimated response function: ( 6.5193e+00, -2.3577e+00), ( 2.6417e+01, 1.0673e+01) Amplitude of the estimated response function: 6.9325e+00, 2.8492e+01 Phase(deg.) of the estimated response function: -19.9, 22.0 Weighted residual power: 2.13004 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.9598 Parameter c: 2.66146 Sum of weights: 14409.3 Squared coherence: 0 Estimated response function: ( 7.9981e+00, 1.2492e+00), ( 1.7627e+01, 1.8127e+01) Amplitude of the estimated response function: 8.0951e+00, 2.5284e+01 Phase(deg.) of the estimated response function: 8.9, 45.8 Weighted residual power: 0.377248 Iteration number = 1 Scale factor: 0.9598 Parameter c: 2.66146 Sum of weights: 14938.6 Squared coherence: 0 Estimated response function: ( 8.0010e+00, 2.2643e+00), ( 1.7061e+01, 1.7211e+01) Amplitude of the estimated response function: 8.3152e+00, 2.4234e+01 Phase(deg.) of the estimated response function: 15.8, 45.2 Weighted residual power: 0.406347 Iteration number = 2 Scale factor: 0.9598 Parameter c: 2.66146 Sum of weights: 15126.6 Squared coherence: 0 Estimated response function: ( 8.3339e+00, 2.6289e+00), ( 1.7571e+01, 1.6561e+01) Amplitude of the estimated response function: 8.7387e+00, 2.4146e+01 Phase(deg.) of the estimated response function: 17.5, 43.3 Weighted residual power: 0.432913 Iteration number = 3 Scale factor: 0.9598 Parameter c: 2.66146 Sum of weights: 15106.5 Squared coherence: 0 Estimated response function: ( 8.5018e+00, 2.6752e+00), ( 1.8038e+01, 1.6063e+01) Amplitude of the estimated response function: 8.9128e+00, 2.4154e+01 Phase(deg.) of the estimated response function: 17.5, 41.7 Weighted residual power: 0.432157 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Calculate response functions for output variable 1 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 21599 Squared coherence: 0 Estimated response function: ( -5.3388e+00, -5.2437e+00), ( -6.0756e+00, 1.4369e+00) Amplitude of the estimated response function: 7.4833e+00, 6.2432e+00 Phase(deg.) of the estimated response function: -135.5, 166.7 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.499754 Sum of weights: 21434 Squared coherence: 0 Estimated response function: ( -5.1478e+00, -4.8192e+00), ( -7.8560e+00, 1.9903e+00) Amplitude of the estimated response function: 7.0516e+00, 8.1042e+00 Phase(deg.) of the estimated response function: -136.9, 165.8 Weighted residual power: 0.481852 Iteration number = 1 Scale factor: 0.522253 Sum of weights: 21421.5 Squared coherence: 0 Estimated response function: ( -5.4162e+00, -4.4032e+00), ( -8.1627e+00, 1.7402e+00) Amplitude of the estimated response function: 6.9802e+00, 8.3461e+00 Phase(deg.) of the estimated response function: -140.9, 168.0 Weighted residual power: 0.487055 Iteration number = 2 Scale factor: 0.523792 Sum of weights: 21417 Squared coherence: 0 Estimated response function: ( -5.4480e+00, -4.3095e+00), ( -8.3763e+00, 1.7677e+00) Amplitude of the estimated response function: 6.9464e+00, 8.5608e+00 Phase(deg.) of the estimated response function: -141.7, 168.1 Weighted residual power: 0.494324 Iteration number = 3 Scale factor: 0.527559 Sum of weights: 21415.8 Squared coherence: 0 Estimated response function: ( -5.4748e+00, -4.2584e+00), ( -8.4422e+00, 1.7271e+00) Amplitude of the estimated response function: 6.9360e+00, 8.6170e+00 Phase(deg.) of the estimated response function: -142.1, 168.4 Weighted residual power: 0.496214 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.527559 Parameter c: 2.66146 Sum of weights: 15218.2 Squared coherence: 0 Estimated response function: ( -8.8000e+00, -6.7967e+00), ( -1.1073e+01, -1.0097e+00) Amplitude of the estimated response function: 1.1119e+01, 1.1119e+01 Phase(deg.) of the estimated response function: -142.3, -174.8 Weighted residual power: 0.235204 Iteration number = 1 Scale factor: 0.527559 Parameter c: 2.66146 Sum of weights: 13765.9 Squared coherence: 0 Estimated response function: ( -9.7411e+00, -1.0796e+01), ( -8.8000e+00, -2.6436e+00) Amplitude of the estimated response function: 1.4541e+01, 9.1886e+00 Phase(deg.) of the estimated response function: -132.1, -163.3 Weighted residual power: 0.168926 Iteration number = 2 Scale factor: 0.527559 Parameter c: 2.66146 Sum of weights: 13408.6 Squared coherence: 0 Estimated response function: ( -1.0260e+01, -1.2521e+01), ( -7.6798e+00, -5.3105e+00) Amplitude of the estimated response function: 1.6187e+01, 9.3371e+00 Phase(deg.) of the estimated response function: -129.3, -145.3 Weighted residual power: 0.156163 Iteration number = 3 Scale factor: 0.527559 Parameter c: 2.66146 Sum of weights: 13042 Squared coherence: 0 Estimated response function: ( -1.0032e+01, -1.3651e+01), ( -5.4167e+00, -5.9975e+00) Amplitude of the estimated response function: 1.6941e+01, 8.0815e+00 Phase(deg.) of the estimated response function: -126.3, -132.1 Weighted residual power: 0.133816 Iteration number = 4 Scale factor: 0.527559 Parameter c: 2.66146 Sum of weights: 13027.9 Squared coherence: 0 Estimated response function: ( -1.0665e+01, -1.3885e+01), ( -3.8490e+00, -6.4626e+00) Amplitude of the estimated response function: 1.7508e+01, 7.5219e+00 Phase(deg.) of the estimated response function: -127.5, -120.8 Weighted residual power: 0.134361 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Estimate errors by fixed-weights bootstrap -------------------------------------------------------------------------------- ================================================================================ Now Frequency(Hz): 0.75, Period(s): 1.33333 ================================================================================ -------------------------------------------------------------------------------- Calculate response functions for output variable 0 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 43199 Squared coherence: 0 Estimated response function: ( -2.1007e+01, 2.5791e+00), ( 1.6554e+01, 1.0975e+01) Amplitude of the estimated response function: 2.1164e+01, 1.9862e+01 Phase(deg.) of the estimated response function: 173.0, 33.5 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 1.26013 Sum of weights: 42352.1 Squared coherence: 0 Estimated response function: ( -2.1390e+01, 6.2375e+00), ( 2.1976e+01, 1.2097e+01) Amplitude of the estimated response function: 2.2281e+01, 2.5085e+01 Phase(deg.) of the estimated response function: 163.7, 28.8 Weighted residual power: 5.25165 Iteration number = 1 Scale factor: 1.406 Sum of weights: 42283 Squared coherence: 0 Estimated response function: ( -2.0743e+01, 6.1785e+00), ( 2.1674e+01, 1.1687e+01) Amplitude of the estimated response function: 2.1644e+01, 2.4624e+01 Phase(deg.) of the estimated response function: 163.4, 28.3 Weighted residual power: 4.97399 Iteration number = 2 Scale factor: 1.38322 Sum of weights: 42298.4 Squared coherence: 0 Estimated response function: ( -2.0761e+01, 6.1442e+00), ( 2.1642e+01, 1.1682e+01) Amplitude of the estimated response function: 2.1651e+01, 2.4594e+01 Phase(deg.) of the estimated response function: 163.5, 28.4 Weighted residual power: 4.97881 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 1.38322 Parameter c: 2.66135 Sum of weights: 28051.8 Squared coherence: 0 Estimated response function: ( -1.8390e+01, -6.3793e+00), ( 2.1348e+01, 2.1292e+01) Amplitude of the estimated response function: 1.9465e+01, 3.0151e+01 Phase(deg.) of the estimated response function: -160.9, 44.9 Weighted residual power: 1.11169 Iteration number = 1 Scale factor: 1.38322 Parameter c: 2.66135 Sum of weights: 27139.4 Squared coherence: 0 Estimated response function: ( -1.3550e+01, -1.1726e+01), ( 1.5351e+01, 2.3291e+01) Amplitude of the estimated response function: 1.7919e+01, 2.7895e+01 Phase(deg.) of the estimated response function: -139.1, 56.6 Weighted residual power: 0.733976 Iteration number = 2 Scale factor: 1.38322 Parameter c: 2.66135 Sum of weights: 27965.4 Squared coherence: 0 Estimated response function: ( -1.3161e+01, -1.1265e+01), ( 1.3629e+01, 2.2692e+01) Amplitude of the estimated response function: 1.7324e+01, 2.6470e+01 Phase(deg.) of the estimated response function: -139.4, 59.0 Weighted residual power: 0.780166 Iteration number = 3 Scale factor: 1.38322 Parameter c: 2.66135 Sum of weights: 28464.5 Squared coherence: 0 Estimated response function: ( -1.3633e+01, -1.0610e+01), ( 1.3431e+01, 2.2749e+01) Amplitude of the estimated response function: 1.7275e+01, 2.6418e+01 Phase(deg.) of the estimated response function: -142.1, 59.4 Weighted residual power: 0.854842 Iteration number = 4 Scale factor: 1.38322 Parameter c: 2.66135 Sum of weights: 28487 Squared coherence: 0 Estimated response function: ( -1.3833e+01, -1.0441e+01), ( 1.3673e+01, 2.2901e+01) Amplitude of the estimated response function: 1.7332e+01, 2.6673e+01 Phase(deg.) of the estimated response function: -143.0, 59.2 Weighted residual power: 0.869199 Iteration number = 5 Scale factor: 1.38322 Parameter c: 2.66135 Sum of weights: 28404.2 Squared coherence: 0 Estimated response function: ( -1.3804e+01, -1.0485e+01), ( 1.3776e+01, 2.2929e+01) Amplitude of the estimated response function: 1.7334e+01, 2.6749e+01 Phase(deg.) of the estimated response function: -142.8, 59.0 Weighted residual power: 0.858698 Iteration number = 6 Scale factor: 1.38322 Parameter c: 2.66135 Sum of weights: 28382.5 Squared coherence: 0 Estimated response function: ( -1.3762e+01, -1.0515e+01), ( 1.3762e+01, 2.2909e+01) Amplitude of the estimated response function: 1.7319e+01, 2.6725e+01 Phase(deg.) of the estimated response function: -142.6, 59.0 Weighted residual power: 0.853772 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Calculate response functions for output variable 1 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 43199 Squared coherence: 0 Estimated response function: ( 4.6419e+00, -4.1036e+00), ( -1.1865e+01, -1.0362e+01) Amplitude of the estimated response function: 6.1958e+00, 1.5753e+01 Phase(deg.) of the estimated response function: -41.5, -138.9 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 0.887056 Sum of weights: 42643.1 Squared coherence: 0 Estimated response function: ( 6.0624e+00, -5.0822e+00), ( -1.5007e+01, -1.1896e+01) Amplitude of the estimated response function: 7.9108e+00, 1.9150e+01 Phase(deg.) of the estimated response function: -40.0, -141.6 Weighted residual power: 2.30498 Iteration number = 1 Scale factor: 1.00498 Sum of weights: 42561.4 Squared coherence: 0 Estimated response function: ( 5.2540e+00, -5.4387e+00), ( -1.4979e+01, -1.0909e+01) Amplitude of the estimated response function: 7.5620e+00, 1.8531e+01 Phase(deg.) of the estimated response function: -46.0, -143.9 Weighted residual power: 2.14692 Iteration number = 2 Scale factor: 0.983764 Sum of weights: 42578.3 Squared coherence: 0 Estimated response function: ( 5.3811e+00, -5.4998e+00), ( -1.4975e+01, -1.0914e+01) Amplitude of the estimated response function: 7.6944e+00, 1.8530e+01 Phase(deg.) of the estimated response function: -45.6, -143.9 Weighted residual power: 2.16148 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 0.983764 Parameter c: 2.66135 Sum of weights: 29005.5 Squared coherence: 0 Estimated response function: ( -5.7355e+00, -8.8745e-01), ( -9.5320e+00, -1.2645e+01) Amplitude of the estimated response function: 5.8037e+00, 1.5835e+01 Phase(deg.) of the estimated response function: -171.2, -127.0 Weighted residual power: 0.554733 Iteration number = 1 Scale factor: 0.983764 Parameter c: 2.66135 Sum of weights: 30815.8 Squared coherence: 0 Estimated response function: ( -6.1502e+00, -6.0600e+00), ( -5.1246e+00, -1.1874e+01) Amplitude of the estimated response function: 8.6342e+00, 1.2932e+01 Phase(deg.) of the estimated response function: -135.4, -113.3 Weighted residual power: 0.450484 Iteration number = 2 Scale factor: 0.983764 Parameter c: 2.66135 Sum of weights: 31610.2 Squared coherence: 0 Estimated response function: ( -4.8719e+00, -8.2765e+00), ( -3.7022e+00, -9.2094e+00) Amplitude of the estimated response function: 9.6040e+00, 9.9257e+00 Phase(deg.) of the estimated response function: -120.5, -111.9 Weighted residual power: 0.399293 Iteration number = 3 Scale factor: 0.983764 Parameter c: 2.66135 Sum of weights: 32704 Squared coherence: 0 Estimated response function: ( -4.2901e+00, -8.9273e+00), ( -5.4396e+00, -8.1713e+00) Amplitude of the estimated response function: 9.9046e+00, 9.8163e+00 Phase(deg.) of the estimated response function: -115.7, -123.7 Weighted residual power: 0.455879 Iteration number = 4 Scale factor: 0.983764 Parameter c: 2.66135 Sum of weights: 32621.4 Squared coherence: 0 Estimated response function: ( -4.4083e+00, -8.8343e+00), ( -6.6191e+00, -7.3481e+00) Amplitude of the estimated response function: 9.8731e+00, 9.8898e+00 Phase(deg.) of the estimated response function: -116.5, -132.0 Weighted residual power: 0.447719 Iteration number = 5 Scale factor: 0.983764 Parameter c: 2.66135 Sum of weights: 32603.7 Squared coherence: 0 Estimated response function: ( -4.2917e+00, -8.5306e+00), ( -7.4222e+00, -6.9884e+00) Amplitude of the estimated response function: 9.5493e+00, 1.0194e+01 Phase(deg.) of the estimated response function: -116.7, -136.7 Weighted residual power: 0.445879 Iteration using Tukey's biweights weight converged -------------------------------------------------------------------------------- Estimate errors by fixed-weights bootstrap -------------------------------------------------------------------------------- ================================================================================ Now Frequency(Hz): 1, Period(s): 1 ================================================================================ -------------------------------------------------------------------------------- Calculate response functions for output variable 0 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 43199 Squared coherence: 0 Estimated response function: ( 1.1221e+01, 1.0270e+01), ( 2.1550e+01, -2.5333e+01) Amplitude of the estimated response function: 1.5211e+01, 3.3259e+01 Phase(deg.) of the estimated response function: 42.5, -49.6 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 1.49203 Sum of weights: 42288.2 Squared coherence: 0 Estimated response function: ( 1.4983e+01, 7.4984e+00), ( 1.2376e+01, -2.6529e+01) Amplitude of the estimated response function: 1.6755e+01, 2.9274e+01 Phase(deg.) of the estimated response function: 26.6, -65.0 Weighted residual power: 5.15193 Iteration number = 1 Scale factor: 1.41619 Sum of weights: 42337.4 Squared coherence: 0 Estimated response function: ( 1.4879e+01, 6.6271e+00), ( 1.0962e+01, -2.6287e+01) Amplitude of the estimated response function: 1.6288e+01, 2.8481e+01 Phase(deg.) of the estimated response function: 24.0, -67.4 Weighted residual power: 4.92398 Iteration number = 2 Scale factor: 1.38763 Sum of weights: 42351.6 Squared coherence: 0 Estimated response function: ( 1.4863e+01, 6.6775e+00), ( 1.1056e+01, -2.6207e+01) Amplitude of the estimated response function: 1.6294e+01, 2.8444e+01 Phase(deg.) of the estimated response function: 24.2, -67.1 Weighted residual power: 4.92551 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 1.38763 Parameter c: 2.66135 Sum of weights: 28107.2 Squared coherence: 0 Estimated response function: ( 6.7594e+01, 3.6278e+00), ( -3.2167e+00, -1.2341e+02) Amplitude of the estimated response function: 6.7692e+01, 1.2345e+02 Phase(deg.) of the estimated response function: 3.1, -91.5 Weighted residual power: 15.6297 Iteration number = 1 Scale factor: 1.38763 Parameter c: 2.66135 Sum of weights: 16463.4 Squared coherence: 0 Estimated response function: ( -1.7109e+01, -1.4646e+01), ( 5.9694e+00, 9.2518e+01) Amplitude of the estimated response function: 2.2522e+01, 9.2711e+01 Phase(deg.) of the estimated response function: -139.4, 86.3 Weighted residual power: 1.04341 Iteration number = 2 Scale factor: 1.38763 Parameter c: 2.66135 Sum of weights: 19832.2 Squared coherence: 0 Estimated response function: ( 8.8412e+00, -1.1025e+00), ( -7.4878e+00, 4.2728e+01) Amplitude of the estimated response function: 8.9097e+00, 4.3379e+01 Phase(deg.) of the estimated response function: -7.1, 99.9 Weighted residual power: 0.367189 Iteration number = 3 Scale factor: 1.38763 Parameter c: 2.66135 Sum of weights: 25368.9 Squared coherence: 0 Estimated response function: ( 2.5536e+01, -2.6414e+01), ( -3.2464e+01, 1.3315e+02) Amplitude of the estimated response function: 3.6740e+01, 1.3705e+02 Phase(deg.) of the estimated response function: -46.0, 103.7 Weighted residual power: 8.65379 Iteration number = 4 Scale factor: 1.38763 Parameter c: 2.66135 Sum of weights: 16790.4 Squared coherence: 0 Estimated response function: ( 8.9621e+00, -9.8033e+00), ( 1.2525e+01, 6.3150e+01) Amplitude of the estimated response function: 1.3283e+01, 6.4380e+01 Phase(deg.) of the estimated response function: -47.6, 78.8 Weighted residual power: 0.34843 Iteration number = 5 Scale factor: 1.38763 Parameter c: 2.66135 Sum of weights: 22355.6 Squared coherence: 0 Estimated response function: ( 1.4553e+01, -1.2539e+01), ( 3.8589e+00, 6.4558e+01) Amplitude of the estimated response function: 1.9210e+01, 6.4674e+01 Phase(deg.) of the estimated response function: -40.7, 86.6 Weighted residual power: 0.779038 Iteration number = 6 Scale factor: 1.38763 Parameter c: 2.66135 Sum of weights: 22168.1 Squared coherence: 0 Estimated response function: ( 2.5982e+01, -7.3970e+00), ( -2.3741e+01, 9.8036e+01) Amplitude of the estimated response function: 2.7014e+01, 1.0087e+02 Phase(deg.) of the estimated response function: -15.9, 103.6 Weighted residual power: 1.72513 Iteration number = 7 Scale factor: 1.38763 Parameter c: 2.66135 Sum of weights: 19243.8 Squared coherence: 0 Estimated response function: ( 1.0013e+01, -7.7338e+00), ( 1.0924e+01, 6.1710e+01) Amplitude of the estimated response function: 1.2652e+01, 6.2669e+01 Phase(deg.) of the estimated response function: -37.7, 80.0 Weighted residual power: 0.398859 Iteration number = 8 Scale factor: 1.38763 Parameter c: 2.66135 Sum of weights: 22539.9 Squared coherence: 0 Estimated response function: ( 1.3997e+01, -1.1760e+01), ( 1.1347e+00, 6.2641e+01) Amplitude of the estimated response function: 1.8281e+01, 6.2651e+01 Phase(deg.) of the estimated response function: -40.0, 89.0 Weighted residual power: 0.775522 Iteration number = 9 Scale factor: 1.38763 Parameter c: 2.66135 Sum of weights: 22387.6 Squared coherence: 0 Estimated response function: ( 2.6026e+01, -9.5294e+00), ( -2.3714e+01, 9.6517e+01) Amplitude of the estimated response function: 2.7715e+01, 9.9387e+01 Phase(deg.) of the estimated response function: -20.1, 103.8 Weighted residual power: 1.77157 Iteration using Tukey's biweights weight does not converge -------------------------------------------------------------------------------- Calculate response functions for output variable 1 -------------------------------------------------------------------------------- Calculate response functions by the ordinary least square method Sum of weights: 43199 Squared coherence: 0 Estimated response function: ( -1.1946e+01, -5.6520e+00), ( -2.0656e+01, -1.9537e+01) Amplitude of the estimated response function: 1.3216e+01, 2.8432e+01 Phase(deg.) of the estimated response function: -154.7, -136.6 Calculate response functions by iteratively reweighted remote reference using Huber weight Iteration number = 0 Scale factor: 1.30468 Sum of weights: 42325.6 Squared coherence: 0 Estimated response function: ( -1.3853e+01, -5.1899e+00), ( -1.2184e+01, -6.4621e+00) Amplitude of the estimated response function: 1.4793e+01, 1.3792e+01 Phase(deg.) of the estimated response function: -159.5, -152.1 Weighted residual power: 2.13191 Iteration number = 1 Scale factor: 0.976564 Sum of weights: 42535.1 Squared coherence: 0 Estimated response function: ( -1.5634e+01, -8.1530e+00), ( -1.3950e+01, -6.8908e+00) Amplitude of the estimated response function: 1.7632e+01, 1.5559e+01 Phase(deg.) of the estimated response function: -152.5, -153.7 Weighted residual power: 2.76019 Iteration number = 2 Scale factor: 1.07101 Sum of weights: 42453.6 Squared coherence: 0 Estimated response function: ( -1.5975e+01, -7.7011e+00), ( -1.3058e+01, -6.8235e+00) Amplitude of the estimated response function: 1.7734e+01, 1.4733e+01 Phase(deg.) of the estimated response function: -154.3, -152.4 Weighted residual power: 2.6312 Iteration number = 3 Scale factor: 1.05938 Sum of weights: 42472.3 Squared coherence: 0 Estimated response function: ( -1.6052e+01, -7.9389e+00), ( -1.3311e+01, -6.8264e+00) Amplitude of the estimated response function: 1.7908e+01, 1.4959e+01 Phase(deg.) of the estimated response function: -153.7, -152.8 Weighted residual power: 2.69327 Iteration number = 4 Scale factor: 1.06686 Sum of weights: 42463.8 Squared coherence: 0 Estimated response function: ( -1.6083e+01, -7.8747e+00), ( -1.3233e+01, -6.8082e+00) Amplitude of the estimated response function: 1.7907e+01, 1.4882e+01 Phase(deg.) of the estimated response function: -153.9, -152.8 Weighted residual power: 2.67971 Iteration using Huber weight converged Calculate response functions by iteratively reweighted remote reference using Tukey's biweights weight Iteration number = 0 Scale factor: 1.06686 Parameter c: 2.66135 Sum of weights: 28716.2 Squared coherence: 0 Estimated response function: ( -1.2142e+01, -7.5734e+00), ( -3.4025e+00, -6.5192e-01) Amplitude of the estimated response function: 1.4310e+01, 3.4644e+00 Phase(deg.) of the estimated response function: -148.0, -169.2 Weighted residual power: 0.349296 Iteration number = 1 Scale factor: 1.06686 Parameter c: 2.66135 Sum of weights: 33530.6 Squared coherence: 0 Estimated response function: ( -1.4746e+01, -6.4551e+00), ( -8.4020e-01, 2.0779e+00) Amplitude of the estimated response function: 1.6097e+01, 2.2413e+00 Phase(deg.) of the estimated response function: -156.4, 112.0 Weighted residual power: 0.610527 Iteration number = 2 Scale factor: 1.06686 Parameter c: 2.66135 Sum of weights: 32681.9 Squared coherence: 0 Estimated response function: ( -1.6438e+01, -5.8912e+00), ( 9.0023e-01, 6.4834e+00) Amplitude of the estimated response function: 1.7462e+01, 6.5456e+00 Phase(deg.) of the estimated response function: -160.3, 82.1 Weighted residual power: 0.710637 Iteration number = 3 Scale factor: 1.06686 Parameter c: 2.66135 Sum of weights: 31158.1 Squared coherence: 0 Estimated response function: ( -1.8141e+01, -4.8866e+00), ( 7.8586e+00, 1.1977e+01) Amplitude of the estimated response function: 1.8788e+01, 1.4325e+01 Phase(deg.) of the estimated response function: -164.9, 56.7 Weighted residual power: 0.901353 Iteration number = 4 Scale factor: 1.06686 Parameter c: 2.66135 Sum of weights: 28638 Squared coherence: 0 Estimated response function: ( -1.9969e+01, -1.7484e+00), ( 1.9546e+01, 4.5856e+00) Amplitude of the estimated response function: 2.0045e+01, 2.0077e+01 Phase(deg.) of the estimated response function: -175.0, 13.2 Weighted residual power: 0.851249 Iteration number = 5 Scale factor: 1.06686 Parameter c: 2.66135 Sum of weights: 26826.6 Squared coherence: 0 Estimated response function: ( -1.3861e+01, 1.0239e+00), ( 1.7565e+01, -1.6266e+01) Amplitude of the estimated response function: 1.3899e+01, 2.3940e+01 Phase(deg.) of the estimated response function: 175.8, -42.8 Weighted residual power: 0.745024 Iteration number = 6 Scale factor: 1.06686 Parameter c: 2.66135 Sum of weights: 27003.1 Squared coherence: 0 Estimated response function: ( -9.4145e+00, -7.9354e+00), ( -1.9777e+01, -2.9175e+00) Amplitude of the estimated response function: 1.2313e+01, 1.9991e+01 Phase(deg.) of the estimated response function: -139.9, -171.6 Weighted residual power: 1.13915 Iteration number = 7 Scale factor: 1.06686 Parameter c: 2.66135 Sum of weights: 28646.6 Squared coherence: 0 Estimated response function: ( -1.5411e+01, -4.9800e+00), ( 7.4761e+00, 1.6842e+01) Amplitude of the estimated response function: 1.6195e+01, 1.8426e+01 Phase(deg.) of the estimated response function: -162.1, 66.1 Weighted residual power: 1.20402 Iteration number = 8 Scale factor: 1.06686 Parameter c: 2.66135 Sum of weights: 28200.4 Squared coherence: 0 Estimated response function: ( -1.7984e+01, 5.4339e+00), ( 2.3078e+01, -7.9002e+00) Amplitude of the estimated response function: 1.8787e+01, 2.4393e+01 Phase(deg.) of the estimated response function: 163.2, -18.9 Weighted residual power: 1.13709 Iteration number = 9 Scale factor: 1.06686 Parameter c: 2.66135 Sum of weights: 25998.5 Squared coherence: 0 Estimated response function: ( -4.2807e+00, -7.6892e+00), ( -1.0323e+01, -3.3470e+01) Amplitude of the estimated response function: 8.8005e+00, 3.5026e+01 Phase(deg.) of the estimated response function: -119.1, -107.1 Weighted residual power: 1.86695 Iteration using Tukey's biweights weight does not converge -------------------------------------------------------------------------------- Estimate errors by fixed-weights bootstrap --------------------------------------------------------------------------------