WAVファイルから得​られた減衰波形に対応​した包絡線をプロット​する方法をご教授いた​だけないでしょうか。​またもっといいWAV​ファイルの読み込み方​はないでしょうか。

filename1="1th.cut230622-T00"+a+"_Tr1.WAV";

filename2="1th.cut230622-T00"+a+"_Tr2.WAV";

filename3="1th.cut230622-T00"+a+"_Tr3.WAV";

filename4="1th.cut230622-T00"+a+"_Tr4.WAV";

[y1,fs1]=audioread(filename1);

[y2,fs2]=audioread(filename2);

[y3,fs3]=audioread(filename3);

[y4,fs4]=audioread(filename4);

[Jiro Doke]

envelope 関数はいかがでしょう？