s1_8.m

来自「能够实现音乐不同乐器的转换,在MatLab上实现」· M 代码 · 共 34 行

clear all;
load('music.mat');

ff=ones(2,1);
ff=kron(ff,f);   %将信号重复2倍

T=1910/8000;    %2430是重复10倍后的长度，8000采样率，T是rw持续的时间
N=1910;
t=linspace(0,T-T/N,N)';

OMG=10000;
K=2000;
omg = linspace(-5000,5000-OMG/K,K)';
FT = T/N*exp(-j*kron(omg,t.'));

FM=zeros(K,30);  %存30个频谱
FBA=zeros(30,1); %存30个频谱中最大的幅度
FB=zeros(30,1);  %存30个最大频率

k=[-5:5];   %确定系数的起止下标
FF=zeros(11,30); %存30套级数系数

for n=1:30
    FM(:,n)=FT*ff(:,n);   %做傅里叶变换
    FBA(n)=max(abs(FM(:,n)));   %找频谱中的最大幅度
    for nf=1:2000
        if abs(FM(nf,n))==FBA(n)
            FB(n)=nf*5-5005;    %找对应最大幅度的频率
        end
    end
    FF(:,n)=1/N*exp(-j*kron(FB(n)*k.',t.'))*ff(:,n); %级数系数(复数)
end

FBm=floor(12*log2(FB/2/pi/220)); %各音对应的音符