Решение систем алгебраических и трансцендентных уравнений, страница 5

·  При формировании матрицы М достаточно большого размера следует иметь в виду, что через панельку Matrix в Mathcad 13можно сформировать матрицу размером только 24х24.

Для большего размера следует использовать таблицу данных, шаблон которой вставляется в РДМ командами меню Insert-Data-Table, вставляющими таблицу практически неограниченного размера (проверялась на матрице 700х700). Этой таблице нужно присвоить имя (например, М) и заполнить ячейки коэффициентами системы уравнений вручную или через буфер при наличии уордовского или эксель-файла с таблицей или векторами коэффициентов.

Если наборы коэффициентов и свободных членов большой системы уравнений можно аппроксимировать статистическим распределением, то решатель большой системы линейных уравнений можно формировать, как показано на рис. 6.1 или в усовершенствованном виде – на рис.6. 3. Этапы ввода пронумерованы; в формирователях n-компонентного векторов с нормальным распределением величин компонент rnorm(n,1,2) и rnorm(n,0.5,2) можно


                                                        

Рис. 6.3. Решатель системы линейных уравнений с n неизвестными (для определенности n = 70)                                    с нормальным распределением  коэффициентов уравнений системы (матрица Z)                                                                    и свободных членов (вектор v1); m – вектор координат найденного решения                                                                          (5 – начало вектора, 6 – конец)

изменять параметры нормального распределения (среднее значение 1 и 0.5, среднеквадратическое отклонение 2), кроме того, можно вместо нормальных векторов использовать векторы с другим распределением величин компонент, выбирая их из 17 подразделов распределений (от бета-распределения до распределения Вейбулла) в разделе справочника (хелпа) Mathcad «Density and Distributions Functions» [45].

Если компоненты М и v можно вычислять по определенным правилам ai,j = fM(i,j) и bi = fv(i), то NxN-компонентную матрицу М и N-компонентный вектор v можно создать, используя функции формирования матриц в Mathcad 13: matrix(N,N, fM(i,j)) и matrix(N,1, fv(i)).