wdn_md 2019-06-20 20:04
浏览 747

matlab 音频信号频谱分析中的周期问题

1.问题
用MATLAB做了一个声音信号的时域和频域分析系统 当导入音频文件出来的波形是非周期信号 但是时域分析却能计算出周期频率 频域分析就会显示周期无限大 该怎么解决或者避免这个情况(时域分析周期计算用的是过零检测法)

2.下面是程序截图

图片说明

图片说明

这两幅图是分析一段音频得出的结果,本来是非周期的,可是计算出了周期

3.代码

Fs=str2double(get(findobj('Tag','sampfreq'),'String'));
N=str2double(get(findobj('Tag','sampnum'),'String'));
if handles.inputtype==0
        msgbox('No wave exist! Please choose a input type!');
        return;
end
n=1;
ymax=max([handles.y(1) handles.y(2)]);
ymin=min([handles.y(1) handles.y(2)]);
from=str2double(get(handles.pointfrom,'String'));
to=str2double(get(handles.pointto,'String'));
global ti
global amp
global T
if from<1 || to-from<5
    msgbox('Error range!');
    return;
end
for i=(from+2):(to-1)
    if handles.y(i-1)<0 && handles.y(i-2)<0 && handles.y(i)>=0 && handles.y(i+1)>0
        if handles.y(i)==0
            ti(n)=i;
        else
            ti(n)=i-handles.y(i)/(handles.y(i)-handles.y(i-1));
        end
        amp(n)=(ymax-ymin)/2;
        ymax=0;
        ymin=0;
        n=n+1;
    else
        if ymax<handles.y(i)
            ymax=handles.y(i);
        end
        if ymin>handles.y(i)
            ymin=handles.y(i);
        end
    end
end
n=n-1;
for i=1:n-1
    T=ti(i+1)-ti(i);
end

请大神帮忙看一下,本人萌新,若有不清楚的或错误的地方,谢谢大佬指导!

  • 写回答

0条回答 默认 最新

    报告相同问题?

    悬赏问题

    • ¥15 乘性高斯噪声在深度学习网络中的应用
    • ¥15 运筹学排序问题中的在线排序
    • ¥15 关于docker部署flink集成hadoop的yarn,请教个问题 flink启动yarn-session.sh连不上hadoop,这个整了好几天一直不行,求帮忙看一下怎么解决
    • ¥30 求一段fortran代码用IVF编译运行的结果
    • ¥15 深度学习根据CNN网络模型,搭建BP模型并训练MNIST数据集
    • ¥15 C++ 头文件/宏冲突问题解决
    • ¥15 用comsol模拟大气湍流通过底部加热(温度不同)的腔体
    • ¥50 安卓adb backup备份子用户应用数据失败
    • ¥20 有人能用聚类分析帮我分析一下文本内容嘛
    • ¥30 python代码,帮调试,帮帮忙吧