我有一个标量 Q(r, theta) ,我想在极坐标中插值 . 问题是θ从0到2π运行 . 因此,在0和2π接口处存在不连续性,这产生了问题 .
Q(r, theta)
我的原始数据是在非结构化网格上,所以我在MATLAB中使用griddata .
我会对域进行三角测量,基本上转换为笛卡尔坐标 . 三角测量当然会包装,因为它不再关心不连续性 . 现在,三角测量中的插值将是微不足道的,并且它将使用位于该周期边界上的正确信息 . (一旦转换为笛卡尔坐标, TriScatteredInterp 就足以解决这个问题了 . )
TriScatteredInterp
1 回答
我会对域进行三角测量,基本上转换为笛卡尔坐标 . 三角测量当然会包装,因为它不再关心不连续性 . 现在,三角测量中的插值将是微不足道的,并且它将使用位于该周期边界上的正确信息 . (一旦转换为笛卡尔坐标,
TriScatteredInterp
就足以解决这个问题了 . )