音乐合成_matlab.docVIP

MATLAB音乐合成 一、简单音乐合成 1、根据音乐简谱和十二平均律计算出每个乐音的频率，在 MATLAB中生成幅度为 1，抽样频率为 8kHz 的正弦信号表示这些乐音，用 sound 播放合成的音乐 A=[174.61 196 220 246.94 261.62 293.66 329.63 349.23 392 440 493.88 523.25 587.33 659.25 698.45 739.99 783.99]; % 定义各调频率 fs=8000; %修改频率 pu=[1 1 1 3 5 5 5 5 6 6 6 8 5 4 4 6 6 3 3 3 3 2 2 2 2 5 5 1 1 1 3 5 5 5 5 6 6 6 8 5 4 4 4 6 3 3 3 3 3 2 2 2 3 1 ]; %在pu中只需输入乐谱数字，在下面B=A-什么什么那部分调整音调，省去了写频率的麻烦 p=0.5; last_time=[p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p p p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p ]; %各乐音持续时间为节省空间，用p代替 N=length(pu); f=zeros(1,N); B=A-5; for i=1:N f(i)=B(pu(i)); %f为各个乐音对应的频率，修改其它调 end point=fs*last_time; %各个乐音的抽样点数 total_point=sum(point); %这段音乐的总抽样点数 store=zeros(1,total_point); %用store向量来储存抽样点 m=1; for num=1:N %利用循环产生抽样数据，num表示乐音编号 t=1/fs:1/fs:point(num)/fs; %产生第num个乐音的抽样点 store(m:m+point(num)-1)=sin(2*pi*f(num)*t); %抽样点对应的幅值 m=m+point(num); end sound(store,8000); %播放出store plot(store); %画出store图像 2、除噪音，加包络 在1.1中听到有“啪”的杂声，这是由于相位不连续产生了高频分量，下面通过加包络来消噪 音。比较科学的包络如下图所示，每个乐音都经过冲激、衰减、持续、消失四个阶段。 由上图可以看出这个包络是四段直线段构成的，因此只要确定了每段线段的端点，即可用 端点数据写出直线方程，因为直线方程可以用通式写出（可以用斜截式），因此这段包络可以 用简单的循环来完成。例如认为包络线上的数据如下图所示： (1.2) A=[174.61 196 220 246.94 261.62 293.66 329.63 349.23 392 440 493.88 523.25 587.33 659.25 698.45 739.99 783.99]; % 定义各调频率 fs=8000; %修改频率 pu=[1 1 1 3 5 5 5 5 6 6 6 8 5 4 4 6 6 3 3 3 3 2 2 2 2 5 5 1 1 1 3 5 5 5 5 6 6 6 8 5 4 4 4 6 3 3 3 3 3 2 2 2 3 1 ]; p=0.5; last_time=[p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p p p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p/2 p ]; N=length