@@ -4,16 +4,14 @@ discard block |
||
| 4 | 4 | return $q != 1 ? ($b1 * (pow($q, $n) - 1)/($q - 1)) : ($n * $b1); |
| 5 | 5 | } |
| 6 | 6 | |
| 7 | -function sn_floor($value) |
|
| 8 | -{ |
|
| 7 | +function sn_floor($value) { |
|
| 9 | 8 | return $value >= 0 ? floor($value) : ceil($value); |
| 10 | 9 | } |
| 11 | 10 | |
| 12 | 11 | // Эта функция выдает нормально распределенное случайное число с матожиднием $mu и стандартным отклонением $sigma |
| 13 | 12 | // $strict - количество $sigma, по которым идет округление функции. Т.е. $strict = 3 означает, что диапазон значений обрезается по +-3 * $sigma |
| 14 | 13 | // Используется http://ru.wikipedia.org/wiki/Преобразование_Бокса_—_Мюллера |
| 15 | -function sn_rand_gauss($mu = 0, $sigma = 1, $strict = false) |
|
| 16 | -{ |
|
| 14 | +function sn_rand_gauss($mu = 0, $sigma = 1, $strict = false) { |
|
| 17 | 15 | // http://ru.wikipedia.org/wiki/Среднеквадратическое_отклонение |
| 18 | 16 | // При $mu = 0 (график симметричный, цифры только для половины графика) |
| 19 | 17 | // От 0 до $sigma ~ 34.1% |
@@ -43,7 +41,7 @@ discard block |
||
| 43 | 41 | * |
| 44 | 42 | * @return float|int |
| 45 | 43 | */ |
| 46 | -function sn_rand_gauss_range($range_start, $range_end, $round = true, $strict = 4, $cut_extreme = false) { |
|
| 44 | +function sn_rand_gauss_range($range_start, $range_end, $round = true, $strict = 4, $cut_extreme = false) { |
|
| 47 | 45 | if($cut_extreme) { |
| 48 | 46 | $range_start--; |
| 49 | 47 | $range_end++; |
@@ -98,12 +96,10 @@ discard block |
||
| 98 | 96 | * |
| 99 | 97 | * @return float |
| 100 | 98 | */ |
| 101 | -function avg($array) |
|
| 102 | -{ |
|
| 99 | +function avg($array) { |
|
| 103 | 100 | return is_array($array) && count($array) ? array_sum($array) / count($array) : 0; |
| 104 | 101 | } |
| 105 | -function linear_calc(&$linear, $from = 0, $debug = false) |
|
| 106 | -{ |
|
| 102 | +function linear_calc(&$linear, $from = 0, $debug = false) { |
|
| 107 | 103 | for($i = $from; $i < count($linear); $i++) |
| 108 | 104 | { |
| 109 | 105 | $eq = &$linear[$i]; |
@@ -112,7 +108,9 @@ discard block |
||
| 112 | 108 | $eq[$j] /= $eq[$from]; |
| 113 | 109 | } |
| 114 | 110 | } |
| 115 | - if($debug) pdump($linear, 'Нормализовано по х' . $from); |
|
| 111 | + if($debug) { |
|
| 112 | + pdump($linear, 'Нормализовано по х' . $from); |
|
| 113 | + } |
|
| 116 | 114 | |
| 117 | 115 | for($i = $from + 1; $i < count($linear); $i++) |
| 118 | 116 | { |
@@ -122,7 +120,9 @@ discard block |
||
| 122 | 120 | $eq[$j] -= $linear[$from][$j]; |
| 123 | 121 | } |
| 124 | 122 | } |
| 125 | - if($debug) pdump($linear, 'Подставили х' . $from); |
|
| 123 | + if($debug) { |
|
| 124 | + pdump($linear, 'Подставили х' . $from); |
|
| 125 | + } |
|
| 126 | 126 | |
| 127 | 127 | if($from < count($linear) - 1) |
| 128 | 128 | { |
@@ -139,11 +139,14 @@ discard block |
||
| 139 | 139 | $eq[$j] = $eq[$j] - $eq[$from] * $linear[$from][$j]; |
| 140 | 140 | } |
| 141 | 141 | } |
| 142 | - if($debug) pdump($linear, 'Подставили обратно х' . $from); |
|
| 143 | - } |
|
| 144 | - else |
|
| 142 | + if($debug) { |
|
| 143 | + pdump($linear, 'Подставили обратно х' . $from); |
|
| 144 | + } |
|
| 145 | + } else |
|
| 145 | 146 | { |
| 146 | - if($debug) pdump($linear, 'Результат' . $from); |
|
| 147 | + if($debug) { |
|
| 148 | + pdump($linear, 'Результат' . $from); |
|
| 149 | + } |
|
| 147 | 150 | foreach($linear as $index => &$eq) |
| 148 | 151 | { |
| 149 | 152 | pdump($eq[count($linear)], 'x' . $index); |