Как объединить таблицы

@Jennan · Регистрация: 02.10.2017

Author24 — интернет-сервис помощи студентам

Добрый день, господа.

Столкнулся с небольшой проблемой: попытка объединения таблиц привела к "некрасивому" результату - пришлось вручную обходными путями создавать новый пустой массив, который поэлементно затем заполнялся. Нельзя ли как-то эффективнее сделать данный процесс?

u[1..5][t] - численные решения имеющегося ОДУ 1 порядка, при этом значения функции f1[t] (других, соответственно, тоже) выводятся Математикой в фигурных скобках, а не просто числом. В результате весь массив B тоже состоит из кучи скобок, а ведь другой программе (тому же языку C++/фортран) надо будет прочитать этот файл. Можно как-то упростить жизнь?

C++

s = NDSolve[{eqs, u0}, fs, xs]
f1[t_] = u[1][t] /. s;
f2[t_] = u[2][t] /. s;
f3[t_] = u[3][t] /. s;
f4[t_] = u[4][t] /. s;
f5[t_] = u[5][t] /. s;
 
file = File[CreateFile[]]
B = RandomInteger[{0, 0}, {n, n}];
For[i = 1, i < n + 1, i++,
 B[[i, 1]] = f1[i];
 B[[i, 2]] = f2[i];
 B[[i, 3]] = f3[i];
 B[[i, 4]] = f4[i];
 B[[i, 5]] = f5[i];]
B // MatrixForm
 
Write[file, B];
Close[file];

@Mysterious Light · 17.05.2018, 21:09

Сообщение от Jennan

при этом значения функции f1[t] (других, соответственно, тоже) выводятся Математикой в фигурных скобках, а не просто числом

Это из-за того, что NDSolve, как и другие подобные ей функции-решатели, возвращают список решений, каждое из которых для многомерных задач являются списком компонент типа Rule.
Поэтому expr /. NDSolve[...] будет возвращать список.
Если заранее известно, что решение единственно (т-ма Коши), имеет смысл сразу писать NDSolve[...][[1]], тогда "лишних" скобок не будет.

Сообщение от Jennan

В результате весь массив B тоже состоит из кучи скобок [...] Можно как-то упростить жизнь?

Если вдруг всё же мы получили список A списков чего-то (напр., ещё списков), можно "снять" один уровень скобок несколькими способами, которые могут давать различный эффект при нарушении вышеупомянутого предположения о единственности.
1. B = A[[1, All]]. Убираем скобочки с верхнего уровня, оставляя только первый элемент.
2. B = A[[All, 1]]. Убираем скобочки со второго уровня, оставляя только первый элемент.
3. B = Join @@ A. Убираем скобки с верхнего уровня, объединяя (конкатенируя).
4. B = Join @@@ A. Убираем скобки со второго уровня, объединяя (конкатенируя) на втором уровне.

Короче, поиграйтесь.

Lisp

A = RandomInteger[{0, 10}, {2, 3, 1}];
A // MatrixForm
A[[1]] // MatrixForm
A[[All, 1]] // MatrixForm
A[[All, All, 1]] // MatrixForm
Join @@ A // MatrixForm
Join @@@ A // MatrixForm

@Jennan · 17.05.2018, 22:12 **[ТС]**

Да, решение моего уравнения единственно. Просто количество u[1,....] зависит от количества точек разбиения промежутка.
Оба варианта действительно успешно справляются с моими кривыми руками, спасибо!
P.S. Случайно нет волшебной функции, выводящей матрицу в файл в первозданном виде, а не на языке Математики со скобками? Как вариант можно и построчно, но в таком случае придется помучиться со считыванием
Было несколько вариантов:

C++

1
2
3

For[i = 1, i < n + 1, i++, 
 For[j = 1, j < n + 1, j++
   Write[file, B[[1,j]] ] ] ]

C++

1	For[i = 1, i < n + 1, i++, Write[file, B[[i,j]] ] ]

@Mysterious Light · 17.05.2018, 23:17

Сообщение от Jennan

Случайно нет волшебной функции, выводящей матрицу в файл в первозданном виде, а не на языке Математики со скобками?

Для WM первозданный (каноничный) вид со скобками. Остальное уже альтернативные форматы.
Обычно для перевода в другой формат используется Export/Import.
См. ref/Export, там вторым-третьим аргументом нужно будет передать "Table" (в кавычках). Также поищи "Table" (с кавычками) в документации, как формат.

@Jennan · 18.05.2018, 01:30 **[ТС]**

Большое спасибо за пояснения!

@Mysterious Light $Эксперт по математике/физике$ 4300 / 2091 / 431 Регистрация: 19.07.2009 Сообщений: 3,163 Записей в блоге: 24
	17.05.2018, 23:17	4
	Сообщение от Jennan Случайно нет волшебной функции, выводящей матрицу в файл в первозданном виде, а не на языке Математики со скобками? Для WM первозданный (каноничный) вид со скобками. Остальное уже альтернативные форматы. Обычно для перевода в другой формат используется Export/Import. См. ref/Export, там вторым-третьим аргументом нужно будет передать "Table" (в кавычках). Также поищи "Table" (с кавычками) в документации, как формат. 1

@Jennan 0 / 0 / 0 Регистрация: 02.10.2017 Сообщений: 7
	18.05.2018, 01:30 [ТС]	5
	Большое спасибо за пояснения! 0

Как объединить таблицы

Решение