Latest revision as of 16:22, 21 April 2018

Manipulate piano sounnd with Matlab

Put your content here . . . As we all know, the timbre is determined by the waveform of the sound. To manipulate the sound of the piano, first need to find the Fourier Analysis of piano.

From this link Fourier_analysis_in_Music

• 

  Piano

From this page, can get

a0=1 a1=0.1 a2=0.33 a3=0.06 a4=0.05 a5=0.045 a6=0 a7=0.02 a8=0.005 a9=0.005 a10=0 a11=0.005 a12=0.01

Then find the frequency of sound in music from this page

https://pages.mtu.edu/~suits/notefreqs.html

Use the frequency to make every sound of the piano. The basic signal is a cos signal.

To make the sound perfect, add the exp function before every sound.

Finally find a music score and group these sounds together and use the sound command to play it.

The music score is like this page https://www.bethsnotesplus.com/2015/07/scarborough-fair.html

The final matlab code is like this.

======================================================

clc;

clear;

Fs = 22050;

delta = 1/Fs;

FC5 = 523.25;

FD5 = 587.33;

FE5 = 659.25;

FF5 = 698.46;

FG5 = 783.99;

FA5 = 880;

FB5 = 987.77;

d1 = 2;

d05 = 1;

d15 = 3;

t1 = delta:delta:d1;

t05 = delta:delta:d05;

t15 = delta:delta:d15;

PC4_1 = zeros(1,Fs*d1);

PD5_1 = zeros(1,Fs*d1);

PD5_05 = zeros(1,Fs*d05);

PA5_1 = zeros(1,Fs*d1);

PA5_05 = zeros(1,Fs*d05);

PE5_05 = zeros(1,Fs*d05);

PF5_05 = zeros(1,Fs*d05);

PD5_15 = zeros(1,Fs*d15);

p = [1 0.1 0.33 0.06 0.05 0.045 0 0.02 0.005 0.005 0 0.005 0.01];

for n=1:length(p);

 %PC4_1 = PC4_1 + p(n)*exp((-1)*t1*4).*cos(2*pi*n*FC5*t1); % sythesize waveform

 PD5_1 = PD5_1 + p(n)*exp((-1)*t1*3).*cos(2*pi*n*FD5*t1);

 PD5_05 = PD5_05 + p(n)*exp((-1)*t05*3).*cos(2*pi*n*FD5*t05);

 PD5_15 = PD5_15 + p(n)*exp((-1)*t15*3).*cos(2*pi*n*FD5*t15);

 PA5_1 = PA5_1 + p(n)*exp((-1)*t1*3).*cos(2*pi*n*FA5*t1);

 PA5_05 = PA5_05 + p(n)*exp((-1)*t05*3).*cos(2*pi*n*FA5*t05); 

 PE5_05 = PE5_05 + p(n)*exp((-1)*t05*3).*cos(2*pi*n*FE5*t05); 

 PF5_05 = PF5_05 + p(n)*exp((-1)*t05*3).*cos(2*pi*n*FF5*t05);

end

y = [PD5_1 PD5_05 PA5_1 PA5_05 PE5_05 PF5_05 PE5_05 PD5_15];

sound(y, Fs);

======================================================

Copy this code into matlab can play a piece of music in Matlab.

Back to 2018 Spring ECE 301 Boutin