一种可用作晚期混响脉冲响应的衰减曲线

clc;
close all;
clear all;
t = 1;
fs = 16000;
n = (0:1/fs:t)';
N = length(n);
b = normrnd(0,0.1,N,1); % 生成Nx1维服从正态分布的随机数矩阵 
rt = 0.4;
delta = 3*log10(10)/rt;
h = b.*exp(-delta*n);
plot;
plot(h);