文档章节

用 Octave 对音频文件进行基本数学的信号处理

linuxprobe
 linuxprobe
发布于 2016/10/14 17:56
字数 1180
阅读 26
收藏 0

信号叠加

两个信号 S1(t)和 S2(t)相加形成一个新的信号 R(t),这个信号在任何瞬间的值等于构成它的两个信号在那个时刻的值之和。就像下面这样:

R(t) = S1(t) + S2(t)

我们将用 Octave 重新产生两个信号的和并通过图表看达到的效果。首先,我们生成两个不同频率的信号,看一看它们的叠加信号是什么样的。

第一步:产生两个不同频率得信号(oog文件)

>> sig1='cos440.ogg';                  %creating the audio file @440 Hz
>> sig2='cos880.ogg';                  %creating the audio file @880 Hz
>> fs=44100;                           %generating the parameters values (Period, sampling frequency and angular frequency)
>> t=0:1/fs:0.02;
>> w1=2*pi*440*t;
>> w2=2*pi*880*t;
>> audiowrite(sig1,cos(w1),fs);        %writing the function cos(w) on the files created
>> audiowrite(sig2,cos(w2),fs);

然后我们绘制出两个信号的图像。

信号 1 的图像(440 赫兹)

>> [y1, fs] = audioread(sig1);
>> plot(y1)

用 Octave 对音频文件进行基本数学的信号处理用 Octave 对音频文件进行基本数学的信号处理
<
信号 2 的图像(880 赫兹)

>> [y2, fs] = audioread(sig2);
>> plot(y2)

用 Octave 对音频文件进行基本数学的信号处理用 Octave 对音频文件进行基本数学的信号处理

第二步:把两个信号叠加

现在我们展示一下前面步骤中产生的两个信号的和。

>> sumres=y1+y2;
>> plot(sumres)

叠加信号的图像

用 Octave 对音频文件进行基本数学的信号处理用 Octave 对音频文件进行基本数学的信号处理

在 Octaver 中,这个效果产生的声音是独特的,因为它可以仿真音乐家弹奏的低八度或者高八度音符(取决于内部程序设计),仿真音符和原始音符成对,也就是两个音符发出相同的声音。

第三步:把两个真实的信号相加(比如两首音乐歌曲)

为了实现这个目的,我们使用格列高利圣咏中的两首歌曲(声音采样)。
圣母颂曲

首先,我们看一下圣母颂曲并绘出它的图像:

>> [y1,fs]=audioread('avemaria_.ogg');
>> plot(y1)

用 Octave 对音频文件进行基本数学的信号处理用 Octave 对音频文件进行基本数学的信号处理
赞美诗曲

现在我们看一下赞美诗曲并绘出它的图像。

>> [y2,fs]=audioread('hymnus.ogg');
>> plot(y2)

用 Octave 对音频文件进行基本数学的信号处理用 Octave 对音频文件进行基本数学的信号处理

圣母颂曲 + 赞美诗曲

>> y='avehymnus.ogg';
>> audiowrite(y, y1+y2, fs);
>> [y, fs]=audioread('avehymnus.ogg');
>> plot(y)

用 Octave 对音频文件进行基本数学的信号处理用 Octave 对音频文件进行基本数学的信号处理

结果,从音频的角度来看,两个声音信号混合在了一起。

两个信号的乘积

对于求两个信号的乘积,我们可以使用类似求和的方法。我们使用之前生成的相同文件。

R(t) = S1(t) * S2(t)
>> sig1='cos440.ogg';       %creating the audio file @440 Hz
>> sig2='cos880.ogg';                  %creating the audio file @880 Hz
>> product='prod.ogg';                 %creating the audio file for product
>> fs=44100;                           %generating the parameters values (Period, sampling frequency and angular frequency)
>> t=0:1/fs:0.02;
>> w1=2*pi*440*t;
>> w2=2*pi*880*t;
>> audiowrite(sig1, cos(w1), fs);      %writing the function cos(w) on the files created
>> audiowrite(sig2, cos(w2), fs);>> [y1,fs]=audioread(sig1);>> [y2,fs]=audioread(sig2);
>> audiowrite(product, y1.*y2, fs);    %performing the product
>> [yprod,fs]=audioread(product);
>> plot(yprod);                        %plotting the product

注意:我们必须使用操作符 ‘.*’,因为在参数文件中,这个乘积是值与值相乘。更多信息,请参考 Octave 矩阵操作产品手册。

乘积生成信号的图像

用 Octave 对音频文件进行基本数学的信号处理用 Octave 对音频文件进行基本数学的信号处理

两个基本频率相差很大的信号相乘后的图表效果(调制原理)。

第一步:

生成两个频率为 220 赫兹的声音信号。

>> fs=44100;
>> t=0:1/fs:0.03;
>> w=2*pi*220*t;
>> y1=cos(w);
>> plot(y1);

用 Octave 对音频文件进行基本数学的信号处理用 Octave 对音频文件进行基本数学的信号处理

第二步:

生成一个 22000 赫兹的高频调制信号。

>> y2=cos(100*w);
>> plot(y2);

用 Octave 对音频文件进行基本数学的信号处理用 Octave 对音频文件进行基本数学的信号处理

第三步:

把两个信号相乘并绘出图像。

>> plot(y1.*y2);

用 Octave 对音频文件进行基本数学的信号处理用 Octave 对音频文件进行基本数学的信号处理

一个信号和一个标量相乘

一个函数和一个标量相乘的效果等于更改它的值域,在某些情况下,更改的是相标志。给定一个标量 K ,一个函数 F(t) 和这个标量相乘定义为:

R(t) = K*F(t)
>> [y,fs]=audioread('cos440.ogg');        %creating the work files
>> res1='coslow.ogg';                
>> res2='coshigh.ogg';>> res3='cosinverted.ogg';
>> K1=0.2;                                %values of the scalars
>> K2=0.5;>> K3=-1;
>> audiowrite(res1, K1*y, fs);            %product function-scalar
>> audiowrite(res2, K2*y, fs);
>> audiowrite(res3, K3*y, fs);

原始信号的图像

>> plot(y)

用 Octave 对音频文件进行基本数学的信号处理用 Octave 对音频文件进行基本数学的信号处理

信号振幅减为原始信号振幅的 0.2 倍后的图像

>> plot(res1)

用 Octave 对音频文件进行基本数学的信号处理用 Octave 对音频文件进行基本数学的信号处理

信号振幅减为原始振幅的 0.5 倍后的图像

>> plot(res2)

用 Octave 对音频文件进行基本数学的信号处理用 Octave 对音频文件进行基本数学的信号处理

倒相后的信号图像

>> plot(res3)

用 Octave 对音频文件进行基本数学的信号处理用 Octave 对音频文件进行基本数学的信号处理

结论

基本数学运算比如代数和、乘,以及函数与常量相乘是更多高级运算比如谱分析、振幅调制,角调制等的支柱和基础。在下一个教程中,我们来看一看如何进行这样的运算以及它们对声音文件产生的效果。

免费提供最新Linux技术教程书籍,为开源技术爱好者努力做得更多更好:http://www.linuxprobe.com/

本文转载自:http://www.linuxprobe.com/audio-octave.html

共有 人打赏支持
linuxprobe
粉丝 21
博文 257
码字总数 45072
作品 0
朝阳
私信 提问
使用Octave音频处理(一):对文件进行读写操作

环境准备 首先是安装 octave,在 Ubuntu 终端运行下面的命令添加 Octave PPA,然后安装 Octave 。 sudo apt-add-repository ppa:octave/stablesudo apt-get updatesudo apt-get install oct......

linux小陶
2016/11/13
42
0
如何使用Octave 在Ubuntu上科学处理音频

本教程中关注的不是安装和学习使用安装好的音频处理软件,而是从设计和音频工程的角度理解它是如何工作的。 环境准备 首先是安装 octave,在 Ubuntu 终端运行下面的命令添加 Octave PPA,然后...

linuxprobe
2016/09/16
46
0
使用Octave音频处理(二):基本数学信号处理

信号叠加 两个信号 S1(t)和 S2(t)相加形成一个新的信号 R(t),这个信号在任何瞬间的值等于构成它的两个信号在那个时刻的值之和。就像下面这样: R(t) = S1(t) + S2(t) 我们将用 Octave 重新产...

linux小陶
2016/11/13
44
0
数学、信号处理、通讯类C++库--IT++

IT++是一个数学,信号处理和通信类功能的C++库。它的主要用途是模拟通信系统,并进行研究相关领域的通信。它综合了Matlab的功能和C的速度,适合于信号处理等领域。 功能和特点 模板阵列和堆叠...

匿名
2008/12/05
11.5K
1
Octave 3.4.3 发布,数学软件包

Octave 3.4.3 发布了,该版本主要是小bug的修复。 Octave 是一个类似matlab和Scilab的数学软件包,可以进行各种运算,编程。它还有丰富的C++接口可以让用户编程时调用。它绘图使用gnuplot。 ...

红薯
2011/10/24
1K
1

没有更多内容

加载失败,请刷新页面

加载更多

总结:volatile关键字

实现内存可见性原理: 对volatile变量执行写操作时,会在写操作之后加入一条store指令,将CPU缓存数据强制刷新到主内存中 对volatile变量执行读操作的时候,会在读操作前加入一条load指令,重...

浮躁的码农
30分钟前
0
0
OSChina 周六乱弹 —— 看见这花臂了么?赶紧叫大佬!

Osc乱弹歌单(2019)请戳(这里) 【今日歌曲】 @宇辰OSC :分享周华健的单曲《有没有一首歌会让你想起我》 《有没有一首歌会让你想起我》- 周华健 手机党少年们想听歌,请使劲儿戳(这里) ...

小小编辑
今天
119
4
Confluence 6 升级中的一些常见问题

升级的时候遇到了问题了吗? 如果你想尝试重新进行升级的话,你需要首先重新恢复老的备份。不要尝试再次对 Confluence 进行升级或者在升级失败后重新启动老的 Confluence。 在升级过程中的一...

honeymoose
今天
2
0
C++随笔(四)Nuget打包

首先把自己编译好的包全部准备到一个文件夹 像这样 接下来新建一个文本文档,后缀名叫.nuspec 填写内容 <?xml version="1.0"?><package xmlns="http://schemas.microsoft.com/packaging/201......

Pulsar-V
今天
3
0
再谈使用开源软件搭建数据分析平台

三年前,我写了这篇博客使用开源软件快速搭建数据分析平台, 当时收到了许多的反馈,有50个点赞和300+的收藏。到现在我还能收到一些关于dataplay2的问题。在过去的三年,开源社区和新技术的发...

naughty
今天
23
0

没有更多内容

加载失败,请刷新页面

加载更多

返回顶部
顶部