如何使用matlab模拟钢琴音色？

有网友说想用MATLAB模拟钢琴音，于是我尝试写了一段，供大家参考。复制代码可直接运行。 

function Dian_Zi_Qin()
clc;clear;
P=[1  2  3  4 5 6 7 8 9 10 11 12 13 14]; % 对应的数字音调
M=[0 -5 -6 -7 1 2 3 4 5  6  7 +5 +6 +7]; % 对应的曲谱音调

f=360+80*(P-1); % 440为最低音

f(1)=1; % 让f(1)不发声
fs=3200; % 太低听不见

K=[ 1 2,3 4 5  6  7  8 9 10,11,12 13 14 15 16,17 18 19 20 21, 22 23, 24]; % 计数器
Y=[-5 1 0, 1 1 1 -5 -6 -7 1 0 1, 0 0 3 0 1  2  3  5 0 5, 0 0 3  3  1  3  5  3  2  2 0 0 ... % 义勇军进行曲
    6 5 2 3 5 3 0 0 5 3 2 3 1 3 0 -5 -6 1 1 3 3 5 5 2 2 2 -6 2 0 -5 1 0 1 3 0 3 5 0 0 ... 
    1 3 5 5 6 5 0 0 3 1 5 5 5 3 1 0 -5 1 0 3 1 5 5 5 3 1 0 -5 1 0 -5 1 0 -5 1 1 0];
Q=[ 1 0 0 4 0 0 0 0 ... % 将义勇军进行曲转换为数字，0表示不发音
    4 0 4 0 4 0 1 2 3 0 4 0 0 4 0 0 0 ...
    6 0 4 5 6 0 8 0 0 8 0 0 0 ...
    6 0 6 0 4 0 6 0 8 0 0 6 5 0 0 0 5 0 0 0 0 ... 
    9 0 0 8 0 0 5 0 0 6 0 0 8 0 6 0 0 0 8 0 6 0 5 6 4 0 0 6 0 0 0 ...
    1 0 2 4 0 4 0 6 0 6 0 8 0 8 0 5 5 5 0 2 0 0 5 0 0 0 ...
    1 0 4 0 0 0 4 0 6 0 0 0 6 0 8 0 0 0 0 ...
    4 0 6 0 8 0 8 0 9 9 0 8 0 0 0 ...
    6 0 4 0 8 8 8 0 6 0 0 4 0 0 0 ...
    1 0 0 4 0 0 0 ...
    6 0 4 0 8 8 8 0 6 0 0 4 0 0 0 ...
    1 0 0 4 0 0 0 1 0 0 4 0 0 0 1 0 0 4 0 0 4 0 0 0];

L=length(Q);
for i=1:L
    t=1/fs:1/fs:0.5;
    x=sin(2*pi*f(Q(i))*t); % 音频转换
    sound(x,fs);    % 让电脑发声
    pause(0.16)   
end

如有疑问，欢迎交流讨论。