есть способ перевести выраж. из Mathematica в читаемый вид?

[salang]

т.е. в виде обычной формулы. Вот например это:

h0approx = 
  Interpolation[
   Flatten[Table[{{sigma = sigmaint (j - 1), 
       tor = (istart + (i - 1) 0.1) ns}, echoes0[[j, 2]][tor]}, {j, 1,
       nsigma + 1}, {i, 1, npoints}], 1]];
zerosig = {{0., g12}}
echoes = Table[{sigma = (2/c) sigmaint j, 
    samf0 = Table[g12[istart ns + 0.1 i ns], {i, 0, npoints}];
    sampan1 = 
     Table[If[i < 2048, rough0[0.1 (i - 1) ns, sigma], 
       rough0[0.1 (i - 4197) ns, sigma]], {i, 1, 4196}];
    zeroes = Table[0., {i, npoints + 2, 4196}];
    samf0extend = Join[samf0, zeroes];
    trsamf0extend = 
     Sqrt[4196] 10^(-10) Fourier[samf0extend, 
       FourierParameters -> {0, 1}];
    trsampan = 
     Sqrt[4196] 10^(-10) Fourier[sampan1, FourierParameters -> {0, 1}];
    prod = Table[trsampan[[i]] trsamf0extend[[i]], {i, 1, 4196}];
    conv = (1/Sqrt[4196]) (1/10^(-10)) Fourier[prod, 
       FourierParameters -> {0, -1}];
    pwr = Table[{istart ns + 0.1 i ns, conv[[i]]}, {i, 0, npoints}];
    g1 = Interpolation[pwr]}, {j, 1, nsigma}];

echoes2 = Join[zerosig, echoes];
Plot[echoes2[[1, 2]][t], {t, -50 ns, 120 ns}]

[Vince Diesel]

Код нерабочий.

[salang]

это неважно. Как перевести-то?

[ewert]

salang в сообщении #880123 писал(а):

это неважно. Как перевести-то?

Действительно неважно. Никак в принципе -- независимо от рабочести или нерабочести алгоритм невозможно представить в виде формулы.

[Vince Diesel]

Что "это"? Последней командой стоит построение графика. Требуется найти формулу для функции, задающей график? Или "это" весь приведенный текст? О чем речь то вообще?

[salang]

ewert в сообщении #880180 писал(а):

Никак в принципе

благодарю, это и хотел уточнить

Vince Diesel в сообщении #880197 писал(а):

Что "это"? Последней командой стоит построение графика. Требуется найти формулу для функции, задающей график? Или "это" весь приведенный текст? О чем речь то вообще?

это только небольшая часть численного расчета, значительно отличающегося от аналитики. Да, требуется найти математическую модель, использованную в расчете.
По каким критериям работает $Flatten$?
Перед интерполированием сглаживание зачем производится? Вроде при интерполяции оно автоматически выполняется?

[Vince Diesel]

По всем командам математики есть подробный хелп. В том числе и по Flatten и Interpolation.

[arseniiv]

salang в сообщении #880247 писал(а):

По каким критериям работает $Flatten$?

Здесь это, скорее, техническая деталь реализации, т. е. для математической модели, если это не модель реализации, это не нужно.

-- Чт июн 26, 2014 18:03:40 --

Кстати, для маленьких кусков кода есть тег tt: Flatten. В виде формулы их можно записывать, но надо тогда использовать \verbsкодs, где s — символ, не встречающийся в коде; типа \verb|Flatten|: $\verb|Flatten|$, но это лишние изощрения при наличии tt. Так, на заметку.

[salang]

А возможно ли не использовать вложенный один в другой массив в одной строке, т.е. сделать это в разных строках

echoes = Table[{sigma = (2/c) sigmaint j, 
   samf0 = Table[g01[istart ns + 0.1 i ns], {i, 0, npoints}];

[Aritaborian]

У вас тут вообще какой-то бред написан. Вы точно ничего не потеряли при копировании?

[salang]

нет, это одна строка из файла.

[Aritaborian]

Там тупо не хватает скобок. Перепроверьте.

[salang]

cкопировал целиком:

echoes = Table[{sigma = (2/c) sigmaint j, 
    samf0 = Table[g01[istart ns + 0.1 i ns], {i, 0, npoints}];
    sampan1 = 
     Table[If[i < 2048, rough0[0.1 (i - 1) ns, sigma], 
       rough0[0.1 (i - 4197) ns, sigma]], {i, 1, 4196}];
    zeroes = Table[0., {i, npoints + 2, 4196}];
    samf0extend = Join[samf0, zeroes];
    trsamf0extend = 
     Sqrt[4196] 10^(-10) Fourier[samf0extend, 
       FourierParameters -> {0, 1}];
    trsampan = 
     Sqrt[4196] 10^(-10) Fourier[sampan1, FourierParameters -> {0, 1}];
    prod = Table[trsampan[[i]] trsamf0extend[[i]], {i, 1, 4196}];
    conv = (1/Sqrt[4196]) (1/10^(-10)) Fourier[prod, 
       FourierParameters -> {0, -1}];
    pwr = Table[{istart ns + 0.1 i ns, conv[[i]]}, {i, 0, npoints}];
    g1 = Interpolation[pwr]}, {j, 1, nsigma}];

[Aritaborian]

Это ужас какой-то, а не код. Чего вы вообще хотите добиться? Даже не знаю, с какого конца к этому подбираться. Вот вы понимаете, например, что переменная nsigma у вас нигде не определена, а вы пытаетесь использовать её в качестве итератора. Ай, что это я о мелочах. Да так вообще никто не пишет. Жуть и мрак.

[salang]

3-мя страницами ранее в этом файле все константы и функции определены. Если требуется, я готов выложить. Мне бы упростить немного для облегчения восприятия и все. Результатом расчета является численное интегрирование и построение графика этой функции.