s1_7.m

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

clear all;

load('guitar.mat');
s1=ones(10,1);
rw=kron(s1,realwave);   %将信号重复10倍

T=2430/8000;    %2430是重复10倍后的长度，8000采样率，T是rw持续的时间
N=2430;
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.'));

RW=FT*rw;   %做傅里叶变换
plot(omg,abs(RW));

%得到基频为2069rad/s==329.2916Hz