unwrapMultipart

创建一个相位角向量，其取值范围为 0 到 4 $\pi$，且该向量被一个 NaN 的值分为两部分。利用相位角来定义螺旋线的 $\mathit{xy}$ 坐标。

t1 = linspace(0,1.75*pi);
t2 = linspace(2.25*pi,4*pi);
t = [t1 NaN t2];

x = t/pi.*cos(t);
y = t/pi.*sin(t);
plot(x,y)

使用 atan2 函数，根据 $\mathit{xy}$ 坐标求出螺旋线的相位角。atan2 函数返回闭区间 $-\pi$ 到 $\pi$ 内的角度值。

P = atan2(y,x);
plot(t,P)

请注意，该图存在不连续点。使用 unwrapMultipart 函数消除不连续点。unwrapMultipart 函数会独立地解包向量中由 NaN 分隔的部分。

Q1 = unwrapMultipart(P);
plot(t,Q1)

将结果与 unwrap 函数进行比较，该函数不会独立地拆分由 NaN 分隔的部分。

Q2 = unwrap(P);
plot(t,Q2)

创建一个经度向量，范围为 0 到 720 度，该向量由 NaN 值分为两部分。

lon = [0:20:340 NaN 380:20:720];

使用 wrapTo180 函数将经度值回绕到从 –180 度到 180 度的闭区间内。

P = wrapTo180(lon);
plot(P)

使用 unwrapMultipart 函数对经度角进行展开。该函数会独立地解包向量中以 NaN 分隔的部分。

Q = unwrapMultipart(P,"degrees");
plot(Q)

将结果与 unwrap 函数进行比较，该函数仅支持弧度制，且不会独立地对 NaN 分割的部分进行解卷。

rad = deg2rad(P);
Q2 = unwrap(rad);
plot(Q2)