matlab 画短时平均幅度谱

和上一节 matlab 画短时能量图的区别在于，上一节是求平方，现在是求绝对值。但是得到的图形和上一节一样，短时平均幅度谱相比短时能量谱的好处在于节省计算量，而且不会出现计算溢出错误。

fid=fopen('zqq.txt', 'rt');
x=fscanf(fid, '%f');
fclose(fid);%计算窗长为50，帧移=0时的语音能量
s=fra(50, 50, x);%帧长为50个数字，帧移为50个数字, 得到s为1319行，50列的矩阵。也就是每一行有50个数字，这50个数字为一帧
%disp(length(s2));
%disp(size(s2, 1));
%disp(size(s2, 2));
s2=abs(s);
%b=sum(a,dim); a表示矩阵；dim等于1或者2，1表示每一列进行求和，2表示每一行进行求和；表示每列求和还是每行求和；b表示求得的行向量
energy=sum(s2, 2);%对矩阵s2每一行求和, energy为1319行，1列，也就是每个数字代表一帧能量
%disp(size(energy, 1));
%disp(size(energy, 2));
subplot(2, 2, 1);
plot(energy);
xlabel('帧数');
ylabel('短时能量E');
legend('窗长N=50');
axis([0, 1500, 0, 10*10^5]);%计算窗长为100，帧移=0时的语音能量
s=fra(100, 100, x);
s2=abs(s);
energy=sum(s2, 2);
subplot(2, 2, 2);
plot(energy);
xlabel('帧数');
ylabel('短时能量E');
legend('窗长N=100');
axis([0, 750, 0, 2*10^6]);%计算窗长为400，帧移=0时的语音能量
s=fra(400, 400, x);
s2=abs(s);
energy=sum(s2, 2);
subplot(2, 2, 3);
plot(energy);
xlabel('帧数');
ylabel('短时能量E');
legend('窗长N=400');
axis([0, 190, 0, 7*10^6]);%计算窗长为800，帧移=0时的语音能量
s=fra(800, 800, x);
s2=abs(s);
energy=sum(s2, 2);
subplot(2, 2, 4);
plot(energy);
xlabel('帧数');
ylabel('短时能量E');
legend('窗长N=800');
axis([0, 95, 0, 14*10^6]);function f=fra(len, inc, x)%len为帧长，inc为帧长一帧移%size(X,1),返回矩阵X的行数；%size(X,2),返回矩阵X的列数；%disp(size(x, 1));fh=fix(((size(x, 1) - len)/inc) + 1);%计算帧数，fh=1319%disp(fh);f=zeros(fh, len);%fh行，len列的0矩阵，这里是1319行，50列的0矩阵i=1; n=1;while i<=fhj=1;%下面的循环是每次循环一次填满一行50个数字，也就是一帧，填满一行以后，将行数加1后再次填while j<=lenf(i, j)=x(n);j=j+1;n=n+1;endn=n-len+inc;% n-len代表将恢复到上帧的末尾，+inc的意思是在上一帧的末尾再移动inc长度 i=i+1;end
end

matlab 画短时平均幅度谱相关推荐

【MATLAB图像处理】傅里叶变换--幅度谱、相位谱、逆变换
fft2() 傅里叶正变换 fftshift() 频谱搬移-直流量(f=0)搬移至频谱中心 I=imread('exp2.tif'); %读入原图像 I1=I(:,:,1:3); %四通道转为三通 ...
用MATLAB画出离散幅度和相位谱
Eg:画出x(t)=cos(πt/8),|t|<=4 信号的x(t)的离散幅度和相位谱 main.m echo on fnct='Fcos'; a=-4; b=4; n=24; tol=0.1; ...
matlab 画短时能量图
使用的 zqq.txt 是一个语音文件,只有1列,数字大概有60000多个 fid=fopen('zqq.txt', 'rt'); x=fscanf(fid, '%f'); fclose(fid);% ...
语音信号短时域分析之短时平均能量（四）
由于语音信号的能量随时间而变化,清音和浊音之间的能量差别相当显著,因此对短时能量和短时平均幅度进行分析,可以描述语音的这种特征变换情况.定义n时刻某语音信号的短时平均能量E为: 式中,N为窗长,可见短 ...
【数字语音处理】Part3 语音信号的短时时域分析：短时平均、短时自相关、语音端点检测、基音周期估值
Part3 语音信号的短时时域分析一.帧和加窗的概念二.短时平均能量三.短时平均幅度函数四.短时平均过零率五.短时自相关分析六.基于能量和过零率的语音端点检测七.基音周期估值八.总结 ...
matlab如何将函数做成幅度谱,好用的画包络谱和幅值谱matlab函数 | 学步园
画包络函数: function hua_baol (y,fs,style,varargin) %画包络函数是hua_baoluo的代码简化版 %输入参数要么是3个,要么是5个 %当三个的时候 y_hh ...
matlab画傅立叶变换后相位谱,对一幅图像进行傅里叶变换后，包含频谱（也叫幅度谱）和相位谱两部分，请问那一部分更重要？...
对一幅图像进行傅里叶变换后,包含频谱(也叫幅度谱)和相位谱两部分,请问那一部分更重要? 更多相关问题 [判断题] 抗盐产品切换到中分时,生产线不需要酸洗. [判断题] 机械密封中的载荷系数不是越大越好 ...
音频信号过零率 matlab,MATLAB计算信号短时平均过零率
MATLAB计算信号短时平均过零率 MATLAB计算信号短时平均过零率一.原理: 短时平均过零率表示一帧语音中语音信号波形穿过横轴的次数二.计算过零率代码 clear all; clc; clos ...
MATLAB计算信号短时平均过零率
一.原理: 短时平均过零率表示一帧语音中语音信号波形穿过横轴的次数二.计算过零率代码 clear all; clc; close all; filedir=[]; % 设置路径 filename=' ...
使用matlab读取图像并通过matlab自带的fft2，fftshift，ifftshift，ifft2等函数获取该图像的傅里叶频谱图，幅度图，相位谱图，幅度图，相位谱图，幅度谱重建图以及相位谱重建
%% 读取图像.预处理 Picture1 = imread('Cameraman.tif');tryPicture1=rgb2gray(Picture1); % 如果是RGB图像,转成灰度图 catc ...

matlab 画短时平均幅度谱

matlab 画短时平均幅度谱相关推荐

最新文章

热门文章