标签: phase

之间有什么区别吗

cmath.phase()

模块中的函数cmath，以及

np.angle()

函数来自numpy.

我有在相关点具有等效间隔和相应测量值的数据。例如，这是我拥有的数据的摘录：

y =[2.118, 2.1289, 2.1374, 2.1458, 2.1542, 2.1615, 2.1627, 2.165 2.1687...]

点之间的间隔为 0.1

所以，我需要从数据中获取的是幅度谱（幅度与频率）和相位谱（相位角与频率）。此外，我应该将数据的相位偏移负 90 度 (-pi/2)。

在移动相位并保持幅度不变时，我需要执行逆 fft 并获得新信号。我想在 Python 中做到这一点。

你能给我一个执行此操作的示例吗？

我使用的代码取自另一个 SO 问题，但我做了一些修改

## Perform FFT WITH SCIPY
signalFFT = np.fft.fft(y)

## Get Power Spectral Density
signalPSD = np.abs(signalFFT) ** 2
signalPhase = np.angle(signalFFT)

## Shift the phase py +90 degrees
new_signalPhase =(180/np.pi)*np.angle(signalFFT)+90 


## Get frequencies corresponding to signal 
fftFreq = np.fft.fftfreq(len(signalPSD), 0.1)

## Get positive half of frequencies
i = fftFreq>0

##
plt.figurefigsize=(8,4)
#plt.plot(fftFreq[i], 10*np.log10(signalPSD[i]));

plt.plot(fftFreq[i], …

假设我有一系列类似于此的阶段:

import numpy as np
import matplotlib.pyplot as plt
phase = np.linspace(0., 100., 1000) % np.pi
plt.plot(phase)
plt.show()

(有很多这样的不连续性)

如何从中获得一系列更"连续"的阶段？

当然,我已经尝试过np.unwrap:

plt.plot(np.unwrap(phase))

要么

plt.plot(np.unwrap(phase),discont=0.1)

但它保持完全相似:

我的期望是这样的展开:

下面是比较 fft 相位图与 2 种不同方法的代码：

import numpy as np
import matplotlib.pyplot as plt
import scipy.fftpack

phase = np.pi / 4
f = 1
fs = f*20
dur=10
t = np.linspace(0, dur, num=fs*dur, endpoint=False)
y = np.cos(2 * np.pi * t + phase)
Y = scipy.fftpack.fftshift(scipy.fftpack.fft(y))
f = scipy.fftpack.fftshift(scipy.fftpack.fftfreq(len(t)))

p = np.angle(Y)
p[np.abs(Y) < 1] = 0

fig, ax = plt.subplots(2, 1)
ax[0].plot(t, y)
ax[1].plot(f*fs, p, label='from fft')
ax[1].phase_spectrum(y, fs, window=None, label='from phase_spectrum')
plt.legend()
plt.show()

这是结果：

这是当信号周期数不是整数时的结果：

我有几个问题：

• 为什么使用phase_spectrum 或使用fft 的相位图和角度如此不同？使用 …

我将FFT应用于每个图像.我提取了图像1和2的角度或相位.之后,我想交换两个图像之间的相位响应,但保持大小完整.我到目前为止的代码是:

Image1=double(imread('/home/anelmad/Desktop/ml/signal_images_processing/tp6/barbara.png'))/255;
Image2=double(imread('/home/anelmad/Desktop/ml/signal_images_processing/tp6/lena.bmp'))/255;
Image1=fft2(Image1);
Image2=fft2(Image2);

figure(1);
imshow(angle(Image1));
x=(angle(Image1));
angle(Image2)=x;

figure(2);
imshow(Image2);

我收到此错误:

error: assignment failed, or no method for '<unknown type> = matrix'

我认为问题是我们应该修改傅里叶变换矩阵(fft2)中的复数.那是对的吗？