Форум программистов, компьютерный форум, киберфорум
Wolfram Mathematica
Войти
Регистрация
Восстановить пароль
Блоги Сообщество Поиск Заказать работу  
 
Рейтинг 4.75/4: Рейтинг темы: голосов - 4, средняя оценка - 4.75
0 / 0 / 1
Регистрация: 17.03.2012
Сообщений: 188

Чему равны внутренние парамеры процедуры

13.08.2013, 02:48. Показов 875. Ответов 3
Метки нет (Все метки)

Студворк — интернет-сервис помощи студентам
Доброе время суток форумчане (по Wolfram mathematica),
Есть рабочая процедура:
Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
Interpol2[Lp_, Tzfp_, Hzfp_, Ltp_, Tdp_, tz_, noz_] := 
  Module[{Nmp, tn, tt, i, kt, kta, 
    ktb, \[Rho]0, \[Rho]1, \[Rho]2, \[Rho]3, \[Rho]4, \[Alpha]1, \
\[Alpha]0, \[Alpha]\[Alpha], \[Beta]1, \[Beta]0, \[Beta]\[Beta], c0, 
    c1, c2, k0, xt, Tw, Hw, k},
   Nmp = 2*Ltp + 1; tn = tz - Ltp*Tdp;
   Tw = {}; Hw = {};
   For[i = 1, i <= Nmp, i++, tt = tn + Tdp*(i - 1);
    kt = 1; While[(Tzfp[[kt]] < tt), kt++]; kta = kt - noz; 
    ktb = kt + noz - 1;
    If[kta < 1, kta = 1; ktb = 2*noz]; 
    If[ktb > Lp, kta = Lp - 2*noz + 1; ktb = Lp];
    \[Rho]4 = Sum[(Tzfp[[k]] - tt)^4, {k, kta, ktb}]; \[Rho]3 = 
     Sum[(Tzfp[[k]] - tt)^3, {k, kta, ktb}];
    \[Rho]2 = Sum[(Tzfp[[k]] - tt)^2, {k, kta, ktb}]; \[Rho]1 = 
     Sum[Tzfp[[k]] - tt, {k, kta, ktb}]; \[Rho]0 = 2*noz;
    c2 = Sum[Hzfp[[k]]*(Tzfp[[k]] - tt)^2, {k, kta, ktb}]; 
    c1 = Sum[Hzfp[[k]]*(Tzfp[[k]] - tt), {k, kta, ktb}]; 
    c0 = Sum[Hzfp[[k]], {k, kta, ktb}];
    \[Alpha]1 = (\[Rho]3^2 - \[Rho]2*\[Rho]4); \[Alpha]0 = (\[Rho]2*\
\[Rho]3 - \[Rho]1*\[Rho]4); \[Alpha]\[Alpha] = c2*\[Rho]3 - c1*\[Rho]4;
    \[Beta]1 = (\[Rho]2^2 - \[Rho]1*\[Rho]3); \[Beta]0 = (\[Rho]1*\
\[Rho]2 - \[Rho]0*\[Rho]3); \[Beta]\[Beta] = c1*\[Rho]2 - c0*\[Rho]3;
    xt = (\[Alpha]\[Alpha]*\[Beta]1 - \[Beta]\[Beta]*\[Alpha]1)/(\
\[Alpha]0*\[Beta]1 - \[Beta]0*\[Alpha]1);
    AppendTo[Tw, tt]; AppendTo[Hw, xt]
    ]; {Tw, Hw}];
она формирует 2 массива (Tw, Hw).
Возник такой вопрос:
как посмотреть, чему будут равнятся внутренние параметры процедуры, например: "kta", "ktb" и т.п.?

Заранее спасибо

P.S.
Я ее перевел на Delphi, и мне нужно ее протестировать
0
cpp_developer
Эксперт
20123 / 5690 / 1417
Регистрация: 09.04.2010
Сообщений: 22,546
Блог
13.08.2013, 02:48
Ответы с готовыми решениями:

Чему равны сутки?
Сутки на планете Земля в настоящее время равны 24 часа. Однако солнечные и лунные приливы замедляют вращение Земли и есть...

Чему равны пункт и милл?
В долях дюйма, пожалуйста.

чему равны элементы массива
В программе описан одномерный целочисленный массив A с индексами от 0 до 10. Ниже представлен фрагмент этой программы, в котором значения...

3
Эксперт функциональных языков программированияЭксперт по математике/физике
4313 / 2105 / 431
Регистрация: 19.07.2009
Сообщений: 3,204
Записей в блоге: 24
13.08.2013, 12:20
Страшная у Вас процедура, не понятная без поллитры =)

Самый простой из известных мне способов вывести внутренние переменные с целью узнать, как меняются их значения, заключается в исполнении команды Print[kta] в конце тела цикла; также можно создать локальный список (ktas = {}) и в теле цикла делать ktas = Append[ktas, kta], а перед последней строкой написать Print[ktas], тогда все значения выведутся единочасно и внутри одного блока Out.
Есть ещё один вариант, но Вам, как я думаю, он не пригодится: использовать Dynamic[], чтоб видеть изменение параметров в реальном времени.
Я ее перевел на Delphi, и мне нужно ее протестировать
была необходимость?
0
0 / 0 / 1
Регистрация: 17.03.2012
Сообщений: 188
13.08.2013, 16:32  [ТС]
Цитата Сообщение от Mysterious Light Посмотреть сообщение

была необходимость?
Ага, вообще, мне давно страшно, хоть я уже перевел 90% одной мат.модели, написаной одним к.м.н. на mathematica

спасибо

Добавлено через 18 минут
Извеняюсь за назойливость, это примерно выглядит так:

Matlab M
1
2
3
4
5
Interpol2[...] :=  Module[{...., ktas = {}},
...
For[...., 
      ; AppendTo[Hw, xt]; ktas = Append[ktas, kta] 
      ]; {Tw, Hw} Print[ktas]]
может, я что-то упустил, но списка не получается

Добавлено через 25 минут
А..а, вот как:

Matlab M
1
2
3
4
5
Interpol2[...] := Module[{..., k, ktas},
....Hw = {}; ktas = {};
For[i = 1.... 
...AppendTo[Hw, xt]; ktas = Append[ktas, kta]
]; {Tw, Hw, ktas} ];

вот так выводит (думаю правильно), хотя без "Print" и в массиве вместе с Tw, Hw, доп.массивом, ну и ладушки, но надеюсь правильно?
0
Эксперт функциональных языков программированияЭксперт по математике/физике
4313 / 2105 / 431
Регистрация: 19.07.2009
Сообщений: 3,204
Записей в блоге: 24
13.08.2013, 16:39
Мой вариант:
Matlab M
1
2
3
4
5
6
7
8
Interpol2[...] := Module[{..., k, ktas},
....Hw = {}; ktas = {};
For[i = 1.... 
...AppendTo[Hw, xt]; ktas = Append[ktas, kta]
];
Print[ktas];
{Tw, Hw}
];
«Хорошо то, что нас качает» (c) Рома Михайлов
Правильно так, как нам удобнее )

Не по теме:

По-моему, Do[..., {i,1,Nmp}] выглядит чуть красивее For[i=1, i<=Nmp, ++i, ...]

0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
raxper
Эксперт
30234 / 6612 / 1498
Регистрация: 28.12.2010
Сообщений: 21,154
Блог
13.08.2013, 16:39
Помогаю со студенческими работами здесь

Импульс. Чему равны скорости лодок?
Две груженых лодки,общей массой М=200кг каждая,идут по озеру встречными курсами с одинаковой скоростью V0=2м/с. Когда они поравнялись, из...

чему равны a.b.c не могу понять тк мы их не вводим
procedure f(a,b,c:Real; var x:Real); begin x:=((2*a)-b-sin(c))/(5+Abs(c)) end; procedure TForm1.btn1Click(Sender: TObject); ...

Чему равны нормальное и тангенциальное ускорение
От самолёта летящего горизонтально со скоростью 500 м/с, оторвался предмет. Чему равны нормальное и тангенциальное ускорение предмета через...

Чему будут равны элементы массива
В программе описан одномерный целочисленный массив с индексами от 0 до 10 и целочисленная переменная i. В приведённом ниже фрагменте...

Чему будут равны Фортран выражения?
Чему будут равны фортран выражения? 1/1, 1/2, 3/2 3./2, 3/2. cos(1/3)/3 Добавлено через 40 минут cos(1/3)/3


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
4
Ответ Создать тему
Новые блоги и статьи
SDL3 для Desktop (MinGW): Создаём пустое окно с нуля для 2D-графики на SDL3, Си и C++
8Observer8 10.03.2026
Содержание блога Финальные проекты на Си и на C++: hello-sdl3-c. zip hello-sdl3-cpp. zip Результат:
Установка CMake и MinGW 13.1 для сборки С и C++ приложений из консоли и из Qt Creator в EXE
8Observer8 10.03.2026
Содержание блога MinGW - это коллекция инструментов для сборки приложений в EXE. CMake - это система сборки приложений. Здесь описаны базовые шаги для старта программирования с помощью CMake и. . .
Как дизайн сайта влияет на конверсию: 7 решений, которые реально повышают заявки
Neotwalker 08.03.2026
Многие до сих пор воспринимают дизайн сайта как “красивую оболочку”. На практике всё иначе: дизайн напрямую влияет на то, оставит человек заявку или уйдёт через несколько секунд. Даже если у вас. . .
Модульная разработка через nuget packages
DevAlt 07.03.2026
Сложившийся в . Net-среде способ разработки чаще всего предполагает монорепозиторий в котором находятся все исходники. При создании нового решения, мы просто добавляем нужные проекты и имеем. . .
Модульный подход на примере F#
DevAlt 06.03.2026
В блоге дяди Боба наткнулся на такое определение: В этой книге («Подход, основанный на вариантах использования») Ивар утверждает, что архитектура программного обеспечения — это структуры,. . .
Управление камерой с помощью скрипта OrbitControls.js на Three.js: Вращение, зум и панорамирование
8Observer8 05.03.2026
Содержание блога Финальная демка в браузере работает на Desktop и мобильных браузерах. Итоговый код: orbit-controls-threejs-js. zip. Сканируйте QR-код на мобильном. Вращайте камеру одним пальцем,. . .
SDL3 для Web (WebAssembly): Синхронизация спрайтов SDL3 и тел Box2D
8Observer8 04.03.2026
Содержание блога Финальная демка в браузере. Итоговый код: finish-sync-physics-sprites-sdl3-c. zip На первой гифке отладочные линии отключены, а на второй включены:. . .
SDL3 для Web (WebAssembly): Идентификация объектов на Box2D v3 - использование userData и событий коллизий
8Observer8 02.03.2026
Содержание блога Финальная демка в браузере. Итоговый код: finish-collision-events-sdl3-c. zip Сканируйте QR-код на мобильном и вы увидите, что появится джойстик для управления главным героем. . . .
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin
Copyright ©2000 - 2026, CyberForum.ru