Найдено 18 результатов

protsiv
Пт фев 05, 2010 2:20 pm
Форум: Mathematica
Тема: блокирование строк в коде Mathematica
Ответы: 2
Просмотры: 1671

Есть тая возможность, нежно вставить нечитаемое выражение в скобки со звездочками, например, слово как
(*как*) не будет восприниматься программой как обрабатываемое выражение.
protsiv
Чт фев 04, 2010 5:18 pm
Форум: Mathematica
Тема: Как изметить производную неизвестной в процессе вычислений
Ответы: 0
Просмотры: 1271

Как изметить производную неизвестной в процессе вычислений

Подскажите, знатоки, могу ли я изменить в процессе вычисления NDSolve производную одного из неизвестных f2 с помощью
EvaluationMonitor:>If[(f2'[t]<0.,f2'[t]=0.,f2'[t]=f2'[t]]
или это нужно делать по другому?
protsiv
Пт янв 29, 2010 12:12 pm
Форум: Mathematica
Тема: Помогите найти ошибку
Ответы: 11
Просмотры: 3794

Да я ошибся, правильное время 5.26002 с, тормозной путь 3.36571 м, скорость 1.48047 м/с. И именно оно появляется в файле test3-a. При этом счет останавливался при n=5000 (условная защита от зацикливания). В файле test3-1 я уменьшил общее время вічислений до 2 с поєтому он стал останавливаться раньше...
protsiv
Чт янв 28, 2010 11:46 pm
Форум: Mathematica
Тема: Помогите найти ошибку
Ответы: 11
Просмотры: 3794

В последнем файле test4 Вы оставили круглые скобки пустыми после Block[{}, data = Append[data, {t, Mdata[t]}]]. Насколько я разбираюсь, это означает, что локальные переменные процедуры не определены? Если в ранее присланном Вами файле test3 сделать так же, т.е. Block[{}, data = Append[data, {t, M[t]...
protsiv
Чт янв 28, 2010 9:04 pm
Форум: Mathematica
Тема: Помогите найти ошибку
Ответы: 11
Просмотры: 3794

Уважаемый, Gordon. Я сравнил результаты расчетов двух вариантов решения с формированием data и без. Они почему то разные. Без: Время до остановки 5.26002 с, тормозной путь 3.36571 м, скорость 1.48047 м/с. С: Время до остановки 3.88568 с, тормозной путь 1.75993 м, скорость 0.853614 м/с. Исходные данн...
protsiv
Чт янв 28, 2010 7:41 pm
Форум: Mathematica
Тема: Помогите найти ошибку
Ответы: 11
Просмотры: 3794

Еще раз спасибо Вам,Gordon, Вы профи. Вы так красиво подредактировали мою программу, что я теперь могу идти дальше в своих расчетах.
protsiv
Чт янв 28, 2010 6:40 pm
Форум: Mathematica
Тема: Помогите найти ошибку
Ответы: 11
Просмотры: 3794

Спасибо, Gordon. Моя ошибка состояла в переприсвоении, как я понял.
Block[{t = t, M = M[t]}, data = Append[data, {t, M}]]
Теперь попробую построить новый график M[t].
protsiv
Чт янв 28, 2010 11:56 am
Форум: Mathematica
Тема: Помогите найти ошибку
Ответы: 11
Просмотры: 3794

Вот пример. ClearAll["Global`*"] Clear[n,DM,data] nvag=0;mvag=6.3*10^3; mл=12*10^3;mc=nvag*mvag;m3=3.5* 10^2;m4=3.5* 10^2;g=9.81; j2=1.2* 10^3;j3=10^3;j4=10^3; c3=8.2* 10^6;c4=8.2* 10^6; b3=100;b4=100; cf3=4.2* 10^4;cf4=4.2* 10^4; bf3=10;bf4=10; r=0.265; data {}; T=45.;V0=1.0;M0=7516.0;ip=...
protsiv
Чт янв 28, 2010 12:11 am
Форум: Mathematica
Тема: Помогите найти ошибку
Ответы: 11
Просмотры: 3794

Помогите найти ошибку

Помогите найти ошибку. Решаю систему с помощью NDSolve. Пытаюсь накопить и сохранить промежуточные результаты по одному из параметров M[t] в списке data, а затем вывести их на печать. Вставляю в систему опцию StepMonitor :> Block[{t, M}, data = Append[data, {t, M[t]}]], а после системы Print[data]. ...
protsiv
Вт янв 26, 2010 9:00 pm
Форум: Mathematica
Тема: прекращение вычислений NDSolve[]
Ответы: 18
Просмотры: 6927

Во-первых, хочется остановки при выполнения определенного условия.
Во-вторых,"Alt+." (то есть Abort Evalution)не всегда останавливает, почему-то!
protsiv
Вт янв 26, 2010 1:06 am
Форум: Mathematica
Тема: Подскажите, что означает прерывание NDSolve?
Ответы: 0
Просмотры: 1358

Подскажите, что означает прерывание NDSolve?

Подскажите,пожалуйста, что означает прерывание NDSolve::ndsz: At t == 18.283501766806012`, step size is effectively zero; singularity or stiff system suspected. >>
protsiv
Вс янв 24, 2010 12:15 pm
Форум: Mathematica
Тема: прекращение вычислений NDSolve[]
Ответы: 18
Просмотры: 6927

Спасибо за подсказку, надо попробовать.
protsiv
Вс янв 24, 2010 12:03 am
Форум: Mathematica
Тема: прекращение вычислений NDSolve[]
Ответы: 18
Просмотры: 6927

В четвертой строчке последнюю единицу можно менять на 2, 3, 4, 5. Выдаются списки промежуточных результатов?
protsiv
Сб янв 23, 2010 11:27 pm
Форум: Mathematica
Тема: прекращение вычислений NDSolve[]
Ответы: 18
Просмотры: 6927

У меня в 7-й версии такого примера для оператора Part нет. Возможно у Вас более ранняя версия. Но если вернуться к sol[[1, 1, 2, 1, 1, 2]], то очень хотелось бы знать, для всех ли результатов решения NDSolve будут вызываться одни и те же параметры при использовании [[1, 1, 2, 1, 1, 2]]? Если мы не з...
protsiv
Сб янв 23, 2010 9:25 pm
Форум: Mathematica
Тема: прекращение вычислений NDSolve[]
Ответы: 18
Просмотры: 6927

Пример, о123, как раз не на эту тему. Функция Part может использоваться для выбора заданного элемента выражения из списка. В этом случае вместо одного надо указать три числа — номер выражения как элемента списка, уровень выражения и порядковый номер извлекаемого из выражения объекта. А в нашем случа...