clc
close all
clear all
t0 = 56000; %只取信号前56000个采样点
[simul1,fs] = audioread('fapb0_sa1.wav'); %提取干声信号和其采样率
t2 = 1/fs; %延时时间,先测试刚好移动一点的结果
signal1 = zeros(t0,1);
signal1(1:t0) = simul1(1:t0) ; %提取信号
signal2 = fft(signal1); %作fft
%反变换
frequency= (linspace(0,fs,t0))';%每点对应的频率,且要转置矩阵,否则为行向量,有可能死机
w = 2*pi*frequency; %模拟角频率表达式
transfer = exp(-i*w.*t2); %时域延时对应频域的系数矩阵
Out_w1 = signal2 .* transfer; %fft结果点乘系数(或者说传递函数)
Out1 = ifft(Out_w1); %反变换输出
figure; %作图,把两个信号放在一张图上,放大对比看延时情况
plot (signal1);
hold on;
plot(real(Out1));
对信号进行极短时间延时
需积分: 9 193 浏览量
2018-07-28
11:27:18
上传
评论
收藏 79KB RAR 举报
FIVEε
- 粉丝: 0
- 资源: 2
最新资源
- javaweb毕业设计-java作业管理系统设计(源代码+lw)(可做课程设计).rar
- 我的自己111111111
- 判断一个数是否为质数以及通过C语言实现.docx
- IPv4-HEADER-03 测试数据
- 高级网络人才培训专家-X00070005 第32章 ISDN和DCC
- 高级网络人才培训专家-X00070004 第31章 配置帧中继
- IPv4-HEADER-02测试数据
- tensorflow-2.6.2-cp38-cp38-manylinux2010-x86-64.whl
- tensorflow-2.6.2-cp37-cp37m-win-amd64.whl
- RV180-Firmware-1.0.5.4.img
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈