Pull to refresh
19
Цыбулин Иван@Uranix

User

1
Subscribers
Send message
А как обстоит дело с векторизованными операциями, когда одновременно можно проводить несколько умножений и сложений?
SS — это производная x по a.
f=diff(x(t),'a');
S1=simplify(f);
SS=eval(S1);

Можете объяснить, почему вы приравниваете производную от x1 по a к значению x1 из модели? У них даже размерность разная. Почему вы игнорируете вторую компоненту решения, хотя она явно фигурирует в МНК?
У вас нет шума, метод Рунге-Кутты, который вы используете, дает 6-7 верных знаков при стандартных настройках. В моем понимании, оценка для а должна с самого начала совпадать с точным значением (с теми же 6-7 верными знаками). Возможно (пальцем в небо), проблема в том, что вы сравниваете численное решение с аналитическим в разных точках. Метод сравнения мне кажется очень странным, вы сравниваете два решения всего лишь в одной точке (текущей) и по этой точке оцениваете параметр. На мой взгляд, с шумом это будет работать крайне плохо. Более логично сравнивать, например, интеграл квадрата отклонения за весь промежуток времени. Хотя я бы сравнивал локальные куски решения (например, в окрестности текущего момента времени), либо вообще подставлял траекторию в диффур и искал МНК оценку для а.
Совершенно не понимаю фразу
примерно на 3,5 секунде процесс стабилизируется
. У вас параметр a оценивается непрерывно во времени? Или когда уже известны обе траектории — модельная и фактическая? Почему такая значительная ошибка в определении а, если шума у вас нет? Самую интересную часть спрятали в код.
А бесплатная специализация только для студентов или аспирантам и сотрудникам тоже доступна?
А с чем в Матлабе сравнивали? Я помню там есть PDEToolbox, но он конечно-элементный, а не конечно-разностный. У вас же на картинках из Матлаба явно результаты расчета какого-то конечно-объемного метода.
А для проверки критичности кросс-диффузии попробуйте сетку сильно растянуть по одной координате и запустить просто с точечным начальным условием. В моем понимании, вместо кругового пятна у вас получится эллиптическое.
Ваше утверждение про точность в четвертом знаке неверно. Ошибка аппроксимации данной схемы — image, константа устойчивости — единица. Очевидно, что ошибка не фиксирована, она зависит от гладкости решения. Похоже, что о зануднсти вычислительной математики вы судите исключительно из своего опыта.
Сечение плоскостью, параллельной образующим — квадрат. Площадь элементарно выражается через растояние от плоскости до точки пересечения. Итого image
А как обстоит дело с более сложными областями, в которых расщепление переменных не работает? Например, решить уравнение Лапласа в области Polygon[{{0, 0}, {0, 1}, {3, 3}, {1, 0}}]. А в качестве общего решения уравнения Лапласа Математика по-прежнему выдает C1[x + I y] + C2[x - I y]?
С удовольствием бы почитал про обращение преобразования Радона с точечным источником и о реализации этого всего с использованием CUDA
Я лишь говорю, что ваш алгоритм соответствует по сложности алгоритму работы с денормализованными числами в IEEE. В любом случае, без сравнения конкретных алгоритмов или их железной реализации, это лишь впечатление, которое может быть неверным.
Поставили бы они лучше зеркальный модулятор в фокальную плоскость вместо рассеивающей пластинки — было бы больше похоже на оптический сопроцессор, а это непонятно что. Много ли задач полагается на умножение на случайную матрицу?
А чего стоят ссылки картинкой
Вот, а у меня после прочтения сложилось впечатление, что вы этот алгоритм предлагаете использовать всегда.
я про реализацию на уровне железа
По-моему логично, что работа с денормализованными числами в IEEE не сложнее работы с предлагаемыми ненормализованными числами
Даже с уравнениями с непостоянными коэффициентами
В одномерной теплопроводности у вас не осцилляции, а обычная неустойчивость от большого шага по времени. Проверьте условие Куранта, скорее всего, оно у вас нарушено. Осцилляциями обычно называют немонотонность численного решения при сохранении устойчивости.
Можно про серидины хорд поподробнее? Похоже, что этот способ только для окружностей работает

Information

Rating
Does not participate
Location
Долгопрудный, Москва и Московская обл., Россия
Date of birth
Registered
Activity