• fft分析前后频谱数据


    正弦信号输入

    input

    输入的原始信号

    short [128]


    fir

    滤波后的输出信号

    SHORT [128]


    fft

    傅里叶变换后的freq数据  float [128]


    rmroise

    去除底噪后的fft数据  float[128]


    snr

    wav信号输入


    input

    这里有8位和32位。你研究下16位呗
    For eight bit data use
    MOD(HEX2DEC(A7)+2^7,2^8)-2^7;
     
    For 32 bit data
    MOD(HEX2DEC(A7)+2^31,2^32)-2^31
    
    能够
    MOD(HEX2DEC(A7)+2^15,2^16)-2^15
    -72 
    56 
    56 
    56 
    -72 
    56 
    56 
    56 
    -72 
    -72 
    56 
    -72 
    -200 
    56 
    -72 
    -72 
    -72 
    -72 
    56 
    56 
    -72 
    56 
    56 
    56 
    56 
    -72 
    -72 
    56 
    -72 
    -72 
    56 
    56 
    56 
    -72 
    -72 
    56 
    56 
    -72 
    -72 
    -72 
    56 
    56 
    -72 
    56 
    -72 
    56 
    -72 
    -72 
    56 
    56 
    -72 
    56 
    56 
    56 
    -72 
    56 
    56 
    -72 
    -72 
    56 
    56 
    -72 
    -72 
    -72 
    56 
    56 
    56 
    56 
    -72 
    -72 
    56 
    56 
    -72 
    -72 
    56 
    -72 
    56 
    -72 
    56 
    -72 
    -72 
    56 
    56 
    -72 
    -72 
    56 
    56 
    56 
    -72 
    56 
    56 
    -72 
    56 
    56 
    -72 
    -72 
    56 
    56 
    -72 
    56 
    56 
    -72 
    -72 
    -72 
    -72 
    -72 
    56 
    56 
    56 
    -72 
    -72 
    56 
    56 
    56 
    56 
    56 
    -72 
    56 
    -72 
    -200 
    -72 
    56 
    56 
    -72 
    56 
    56 
    56 
    -72 
    




    fir


    fft


    bn


    snr


    matlab仿真


    仿真的代码

    clf;
    fs=8000;N=128;   %採样频率和数据点数
    n=0:N-1;%t=(n+128)/fs;   %时间序列
    %x=20*sin(2*pi*500*t)+35*sin(2*pi*3000*t); %信号
    x=[-72 ,56 ,56 ,56 ,-72 ,56 ,56 ,56 ,-72 ,-72 ,56 ,-72 ,-200 ,56 ,-72,-72,... 
    -72 ,...
    -72 ,...
    56 ,...
    56 ,...
    -72 ,...
    56 ,...
    56 ,...
    56 ,...
    56 ,...
    -72 ,...
    -72 ,...
    56 ,...
    -72 ,...
    -72 ,...
    56 ,...
    56 ,...
    56 ,...
    -72 ,...
    -72 ,...
    56 ,...
    56 ,...
    -72 ,...
    -72 ,...
    -72 ,...
    56 ,...
    56 ,...
    -72 ,...
    56 ,...
    -72 ,...
    56 ,...
    -72 ,...
    -72 ,...
    56 ,...
    56 ,...
    -72 ,...
    56 ,...
    56 ,...
    56 ,...
    -72 ,...
    56 ,...
    56 ,...
    -72 ,...
    -72 ,...
    56 ,...
    56 ,...
    -72 ,...
    -72 ,...
    -72 ,...
    56 ,...
    56 ,...
    56 ,...
    56 ,...
    -72 ,...
    -72 ,...
    56 ,...
    56 ,...
    -72 ,...
    -72 ,...
    56 ,...
    -72 ,...
    56 ,...
    -72 ,...
    56 ,...
    -72 ,...
    -72 ,...
    56 ,...
    56 ,...
    -72 ,...
    -72 ,...
    56 ,...
    56 ,...
    56 ,...
    -72 ,...
    56 ,...
    56 ,...
    -72 ,...
    56 ,...
    56 ,...
    -72 ,...
    -72 ,...
    56 ,...
    56 ,...
    -72 ,...
    56 ,...
    56 ,...
    -72 ,...
    -72 ,...
    -72 ,...
    -72 ,...
    -72 ,...
    56 ,...
    56 ,...
    56 ,...
    -72 ,...
    -72 ,...
    56 ,...
    56 ,...
    56 ,...
    56 ,...
    56 ,...
    -72 ,...
    56 ,...
    -72 ,...
    -200 ,...
    -72 ,...
    56 ,...
    56 ,...
    -72 ,...
    56 ,...
    56 ,...
    56 ,...
    -72 ,...
    ];
    y=fft(x,N);    %对信号进行高速Fourier变换
    mag=abs(y);     %求得Fourier变换后的振幅
    f=n*fs/N;    %频率序列
    
     figure(1);plot(f,mag);
    figure(2);
    subplot(2,2,1),plot(f,mag);   %绘出随频率变化的振幅
    subplot(2,2,2),plot(n,x);
    xlabel('频率/Hz');
    ylabel('振幅');title('N=128');grid on;
    
    


  • 相关阅读:
    ubuntu 的 软件管理工具包管理(离线安装dpkg、在线安装apt、源码安装(适用于github下载的源码程序))
    前端本地 Nginx 反向代理
    CSS图片高斯模糊方式
    前端开发最喜欢的30个工具
    前端开发网站
    Sql server 根据不同的查询条件拼接SQL
    npm源管理
    工作流引擎
    git常用命令大全
    leaderlinevue移动svg(指示线元素)到另外一个容器中(定位基准)
  • 原文地址:https://www.cnblogs.com/llguanli/p/8723494.html
Copyright © 2020-2023  润新知