zoukankan      html  css  js  c++  java
  • 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;
    
    


  • 相关阅读:
    Centos7安装部署openstack--Keystone认证服务
    Centos7安装部署openstack----基础服务安装
    Centos7 k8s dns
    集中式存储3apr
    Centos7 k8s部署dahsboard
    Centos7 k8s tomcat-app项目持久化
    Centos7 k8s 容器的网络访问service
    Centos 7 k8s Deployment新副本控制器
    模型层中模型的基本了解
    程序员必知必会Git的小知识
  • 原文地址:https://www.cnblogs.com/llguanli/p/8723494.html
Copyright © 2011-2022 走看看