压缩感知处理语音信号代码（python）

先安装convex优化包

conda install -c conda-forge cvxpy

语音信号的生成

import numpy as np
import matplotlib as mpl
import matplotlib.pyplot as plt
import scipy.optimize as spopt
import scipy.fftpack as spfft
import scipy.ndimage as spimg
import cvxpy as cvx

T= 0.000025
N = 5000
fn = np.linspace(0,1-1/N,int(N)) 

t = np.linspace(0, T*N, N)
y = np.sin(1394 * np.pi * t) + np.sin(3266 * np.pi * t)
yt = spfft.dct(y, norm='ortho')


fig, (axs1,axs2) = plt.subplots(2,1)
axs1.plot(t, y)
axs1.set_xlabel('t')
axs1.set_ylabel('y')
axs1.grid(True)

axs2.plot(fn, yt)
axs2.set_xlabel('k')
axs2.set_ylabel('F(y)')
axs2.grid(True)

局部放大