NDFT的信号频谱的程序

NDFT（Non-Uniform Discrete Fourier Transform）是一种非均匀离散傅里叶变换，它可以用于处理非周期信号的频域分析。下面是一个计算NDFT信号频谱的Python程序示例：

import numpy as np

def ndft(x, omega):
    N = len(x)
    X = 0
    for n in range(N):
        X += x[n] * np.exp(-1j * omega[n] * np.arange(N))
    return X

# 生成非周期信号
t = np.linspace(0, 1, 1000)
x = np.sin(2*np.pi*10*t) + np.sin(2*np.pi*30*t) + np.sin(2*np.pi*50*t)

# 生成非均匀采样点
omega = np.sort(np.random.uniform(low=0, high=2*np.pi, size=100))

# 计算NDFT信号频谱
X = ndft(x, omega)

# 绘制频谱图
import matplotlib.pyplot as plt
plt.plot(omega, np.abs(X))
plt.xlabel('Frequency (rad/sample)')
plt.ylabel('Magnitude')
plt.show()

在这个程序中，我们首先生成一个包含三个正弦波的非周期信号。然后，我们生成100个非均匀采样点，并使用`ndft`函数计算该信号在这些采样点处的NDFT。最后，我们绘制NDFT信号频谱图。