基于Python怎么把实数矩阵合并成复数矩阵,有几种方法
Admin 2022-06-29 群英技术资讯 687 次浏览
有时需要把两个实数矩阵,一个作为实部,一个作为虚部,合并为一个复数矩阵,该如何操作?
假如是在第二个维度上进行合并(real: Data[:, 0, :, :] imag: Data[:, 1, :, :]),有两种方法
result = Data[:, 0, :, :] + 1j*Data[:, 1, :, :]
result = 1j*Data[:, 1, :, :] result += Data[:, 0, :, :]
第二种方法更节省内存~
补充:python numpy 分离与合并复数矩阵实部虚部的方法
在进行数字信号处理的过程中,我们往往有对短时傅里叶变换频谱(spectrogram)进行分析的需求。
常见的分析手段对应欧拉公式分为两种,要么使用模与相位的形式,要么使用实部虚部。
本文分享一个简单的将复数光谱图分解为实部与虚部以及将两个部分重新合并为一个复数矩阵的过程,以下为python代码。
import numpy as np import librosa # load the original wav test_wave, _ = librosa.load("../RecFile_1_20200617_153719_Sound_Capture_DShow_5_monoOutput1.wav", sr=44100) # calculate the complex spectrogram stft spectrogram_test_wav = librosa.stft(test_wave, n_fft=735*2, win_length=735*2, hop_length=735) # calculate the real part of the spectrogram real_spectrogram = spectrogram_test_wav.real # calculate the imaginary part of the spectrogram imaginary_spectrogram = spectrogram_test_wav.imag # combine these two parts reconstruction_spectrogram = real_spectrogram + 1j * imaginary_spectrogram print(np.array_equal(spectrogram_test_wav, reconstruction_spectrogram))
其中librosa库为常用的音频处理库。
上述代码实现了对wavfile进行短时傅里叶变换,分离出实部虚部并重新合并的过程。
最终的输出为True, 证明了经过这些步骤过后,重构的复数矩阵与初始的光谱图是一致的。
免责声明:本站发布的内容(图片、视频和文字)以原创、转载和分享为主,文章观点不代表本网站立场,如果涉及侵权请联系站长邮箱:mmqy2019@163.com进行举报,并提供相关证据,查实之后,将立刻删除涉嫌侵权内容。
猜你喜欢
本文主要介绍了python实现自动抢课脚本的示例代码,文中通过示例代码介绍的非常详细,具有一定的参考价值,感兴趣的小伙伴们可以参考一下
内容介绍1.关于箱线图及plt.boxplot()方法2.绘制一幅简单的箱线图3.绘制一幅更精致的图像4.异常值的标准5.异常值的输出1.关于箱线图及plt.boxplot()方法
每天上班最痛苦的事情就是早起早起早起!这是大部分上班族的痛苦,但是不上班又是不可能的啦,因为都是为了搞钱。本文用Python制作了三款有趣的闹钟,感兴趣的可以学习一下
这篇文章主要介绍了python 实现二维数组的索引、删除、拼接操作,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧
这篇文章主要为大家介绍了python机器学习tensorflow构建长短时记忆网络LSTM,有需要的朋友可以借鉴参考下,希望能够有所帮助,祝大家多多进步,早日升职加薪
成为群英会员,开启智能安全云计算之旅
立即注册Copyright © QY Network Company Ltd. All Rights Reserved. 2003-2020 群英 版权所有
增值电信经营许可证 : B1.B2-20140078 粤ICP备09006778号 域名注册商资质 粤 D3.1-20240008