Как сделать полярную интерполяцию в matlab

Question

У меня есть скалярная величина Q(r, theta), которую я хочу интерполировать в полярных координатах. Проблема состоит в том, что theta работает от 0 до 2π. Таким образом, существует разрыв на интерфейсах 0 и 2π, что создает проблему.

Мои исходные данные находятся на неструктурированной сетке, поэтому я использую griddata в MATLAB.

Answer 1

Я бы триангулировал домен, по существу преобразовывая в декартовы координаты. Конечно, триангуляция будет завершена, поскольку она больше не заботится о разрыве. Теперь интерполяция внутри триангуляции будет тривиальной, и она будет правильно использовать информацию, которая лежит через эту периодическую границу. (TriScatteredInterp хватит для этой проблемы, как только вы перейдете в декартовы координаты.)