Документ взят из кэша поисковой машины. Адрес оригинального документа : http://www.mmonline.ru/forum/read/7/42839/
Дата изменения: Mon Apr 11 15:07:40 2016
Дата индексирования: Mon Apr 11 15:07:40 2016
Кодировка: Windows-1251
MMOnline | Форумы | Разное | Чтение матрицы большого размера в Mathematica

Чтение матрицы большого размера в Mathematica

Автор темы Наташа 
15.11.2004 23:01
Наташа
Чтение матрицы большого размера в Mathematica
Поделитесь, пожалуйста, кто знает, как прочитать из файла матрицу в Mathematica? Размер матрицы 1000*1000 и больше. какой формат должен иметь файл? И еще, осилит ли Mathematica расчет собственных векторов и значений для такой матрицы? Если не осилит, может, кто-нибуь знает, какая справится??? Спасибо!
15.11.2004 23:16
Fractal
Maple попробуйте
Создайте тестовый наподобие следующего:
> with(linalg):
> A:=matrix(
[[a_{11}, a_{12}, ..., a_{1N}], [...], ..., [a_{N1}, a_{N2}, ..., a_{NN}]]):
> eigenvectors(A);
Откройте его в Maple и выполните...
Считать будет, вероятно, несколько часов, если памяти вообще хватит. Возможно, проще написать свою программу для приближенного расчета с учетом особенностей матриц (разреженные, блочные и тд...).
16.11.2004 01:08
или Matlab
Если речь идет о числовых (не символьных) матрицах, то удобно с помощью Matlab. Загрузка из файла командой load:
A=load('mymatrix.txt');
Файл может быть обычный текстовый (в каждой строчке элементы разделяются пробелом) либо специальных форматов (я не пробовал).

Процедура eig (с двумя выходными параметрами) находит собственные числа и векторы:
[eigvec,eigval]=eig(A);
Для матрицы порядка 1000 думает менее минуты.

16.11.2004 01:31
Да по-простому:
посмотреть, как у нее в каком-нибудь файле *.nb сделано, да самому написать так же.
Я вот сейчас глянул -
RowBox[{"(", GridBox[{
{"0", \(-1\), "0", "0", "0", \(-1\), \(-1\), "0", "0", "0",
"0", "0", "0", "0"},
и так далее, вроде понятно все.
16.11.2004 22:08
Наташа
спасисбо
Спасибо всем большое! Своя прога есть, но результаты странные получаются если матрица большого размера. Вот я и решила проверить :)
Извините, только зарегистрированные пользователи могут публиковать сообщения в этом форуме.

Кликните здесь, чтобы войти