使用Octave音频处理（二）：基本数学信号处理_写一个octave程序,把4段不同采样率的单通道音频,统一重采样成16khz,并且把这四段

2401_86359338

于 2024-09-08 15:53:52 发布

阅读量385

点赞数 3

文章标签：音视频信号处理

本文链接：https://blog.csdn.net/2401_86359338/article/details/142028141

版权

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

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

外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传

圣母歌曲

赞美诗曲(Hymnus Track)

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

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

外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传

赞美诗曲

圣母颂曲 + 赞美诗曲

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

外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传

圣母歌曲 + 赞美诗曲

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

两个信号的乘积

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

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 矩阵操作产品手册。

乘积生成信号的图像

外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传

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

第一步：

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

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

载波

外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传

第二步：

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

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

调试中

外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传

第三步：

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

>> plot(y1.*y2);

调制后的信号

外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传

一个信号和一个标量相乘

一个函数和一个标量相乘的效果等于更改它的值域，在某些情况下，更改的是相标志。给定一个标量 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);

原始信号的图像