MATLAB演奏音乐

演奏原理：sound(y,Fs) 以采样率 Fs 向扬声器发送音频信号 y。
Fs 指单位时间对音频信号采样的次数。采样频率的选择应该遵循奈奎斯特（Harry Nyquist）采样理论：如果对某一模拟信号进行采样，则采样后可还原的最高信号频率只有采样频率的一半，或者说只要采样频率高于输入信号最高频率的两倍，就能从采样信号系列重构原始信号。正常人听觉的频率范围大约在20Hz~20kHz之间，根据奈奎斯特采样理论，为了保证声音不失真，采样频率应该在40kHz左右。
y 一般指正弦波信号，对应的频率可以在网上找到数据。

由以上可写出自己喜欢的音乐的曲子了，以生日快乐歌为例，找到对应谱子
main.m

clear;clc;
Fs = 40000;
y = [yy(1,5,0.2),yy(1,5,0.2),yy(1,6,0.3),yy(1,5,0.3),yy(2,1,0.3),yy(1,7,0.3),zeros(1,2478),...
    yy(1,5,0.2),yy(1,5,0.2),yy(1,6,0.3),yy(1,5,0.3),yy(2,2,0.3),yy(2,1,0.3),zeros(1,2478),...
    yy(1,5,0.2),yy(1,5,0.2),yy(2,5,0.3),yy(2,3,0.3),yy(2,1,0.3),yy(1,7,0.2),yy(1,6,0.2),zeros(1,2478),...
    yy(2,4,0.2),yy(2,4,0.2),yy(2,3,0.3),yy(2,1,0.3),yy(2,2,0.3),yy(2,1,0.3)];
sound(y,Fs);

注：y构成一个行向量，sound按顺序演奏

yy.m

function y = yy(type,tone,rythm)
%type为1表示低音，为2表示中音，为3表示高音
%tone为音节序号，1为哆
%rythm为持续时间长短
%y为所输出声波信号所对应离散值
    Fs = 40000;
    freqs = zeros(3,7);
    freqs(1,:) = [220,246.94,277.19,293.67,329.63,370,415.31];
    freqs(2,:) = [440,493.88,554.37,587.33,659.33,739.99,830.61];
    freqs(3,:) = [880,987.76,1108.73,1174.66,1318.52,1479.98,1661.22];
    x = linspace(0, 2 * pi * rythm, floor(Fs * rythm));
    y = sin(freqs(type,tone) * x);
end

注：两个函数里的Fs应保持一致。