对函数调用
(艰难)= tfstability (afz,该组织)

输入:
afz -时间/频率振幅轮廓模型，通过Timbre Modeler函数计算我的文件交换
ffl—基频归一化频率集(如Timbre Modeler中使用的)

输出:
Slog -输出稳定向量

来源:
[1]https://qr.ae/TaISSq

示例:使用3个级别分析在900hz下演奏的长笛样本:一个平滑滤波器、一个一阶瞬态滤波器和一个二阶瞬态滤波器。
1.绘制输入样例并播放声音。暂停。
2.对于每个过滤器，绘制模型并播放模型的声音。
3.对于每个模型，计算并绘制音乐稳定性。暂停。
4.最后，回放原始声音进行比较。暂停。然后对每个模型进行加权组合，进行更精确的再创造。

% %用户输入
%%得到采样音乐音频
易= audioread (flutea5.mp3) ';Fs = 44100;
指定hermite滤波器结构，基频，
频率设置和带宽
ff = 900;该组织= (1:25);fd = 1.5 * ff;
Hc = 3;yobulk = 0 (Hc、长度(yi));%选择3级分析Hc
w = [1,1,1].* 1 (1,Hc);%选择输出为等权重
% %非用户输入
声音(咦,Fs);播放原声并暂停
图(1);次要情节(Hc + 2 1 1);t =(0:长度(yi) 1) / Fs;
情节(t,易);xlim ([0, t(结束)));ylabel (AMP);标题(“样本仪器”)
次要情节(Hc + 2 1 Hc + 2);cla重置%清除稳定保持
legtxt =细胞(Hc);启动传奇细胞以获得稳定性
暂停
for it=1:Hc %for所有过滤级别
跳频= [0 (1,1),1];%指定埃尔米特的结构
生成标准高斯小波模型
[哟,afz1, t, G1, tg) = waveletmodel(咦,Fs,跳频,ff, fd,该组织);%叫音色modeler
%%情节模式在指定的子情节和播放声音
声音(哟,Fs);
次要情节(Hc + 2, 1 + 1)
achtxt = ['Model H'，num2str(it-1)];
区(t, afz1”);xlim ([0, t(结束)));标题(achtxt)
yobulk(,) =哟;
计算稳定性和标绘
Sk = tfstability (afz1,该组织);%的电话音乐稳定
次要情节(Hc Hc + 2, 1 + 2)
情节(t, Sk);xlim ([0, t(结束)));抓住;
标题(“稳定”);包含(时间(s));ylabel (' r ')
legtxt{它}= [' H ', num2str (1)];
暂停%暂停
结束
传奇(legtxt)
yplay = 0(1、长度(yi));
播放原声
声音(咦,Fs);暂停重放原声
播放加权过滤输出
为它= 1:Hc
yplay = yplay + w(it)*yobulk(it，:);%
结束
声音(max (abs (yi)) * yplay / max (abs (yplay)), Fs)