%天空の城

fs = 44100; % sample rate

dt = 1/fs;

T16 = 0.125;

t16 = 0:dt:T16;

[temp k] = size(t16);

t1 = linspace(0,16*T16,16*k);

t34 = linspace(0, 12*T16, 12*k);

t2 = linspace(0,8*T16,8*k);

t38 = linspace(0,6*T16,6*k);

t4 = linspace(0,4*T16,4*k);

t8 = linspace(0,2*T16,2*k);

[temp i] = size(t4);

[temp j] = size(t8);

% Modification functions

mod1 = sin(pi*t1/t1(end));

mod34 = sin(pi*t34/t34(end));

mod38 = sin(pi*t38/t38(end));

mod2 = sin(pi*t2/t2(end));

mod4 = sin(pi*t4/t4(end));

mod8 = sin(pi*t8/t8(end));

f0 = 2*146.8; % reference frequency

ScaleTable = [1 9/8 5/4 4/3 3/2 5/3 15/8 ...

2 9/4 5/2 8/3 3 10/3 15/4];

%1 notes

la04 = mod1.*cos(2*pi*ScaleTable(6)*f0*t1);

mi14 = mod1.*cos(2*pi*ScaleTable(10)*f0*t1);

%1/2 notes

do0d = mod2.*cos(2*pi*ScaleTable(1)*f0*t2);

re0d = mod2.*cos(2*pi*ScaleTable(2)*f0*t2);

mi0d = mod2.*cos(2*pi*ScaleTable(3)*f0*t2);

fa0d = mod2.*cos(2*pi*ScaleTable(4)*f0*t2);

so0d = mod2.*cos(2*pi*ScaleTable(5)*f0*t2);

la0d = mod2.*cos(2*pi*ScaleTable(6)*f0*t2);

si0d = mod2.*cos(2*pi*ScaleTable(7)*f0*t2);

do1d = mod2.*cos(2*pi*ScaleTable(8)*f0*t2);

re1d = mod2.*cos(2*pi*ScaleTable(9)*f0*t2);

mi1d = mod2.*cos(2*pi*ScaleTable(10)*f0*t2);

fa1d = mod2.*cos(2*pi*ScaleTable(11)*f0*t2);

so1d = mod2.*cos(2*pi*ScaleTable(12)*f0*t2);

la1d = mod2.*cos(2*pi*ScaleTable(13)*f0*t2);

si1d = mod2.*cos(2*pi*ScaleTable(14)*f0*t2);

%3/8 notes

do0t = mod38.*cos(2*pi*ScaleTable(1)*f0*t38);

re0t = mod38.*cos(2*pi*ScaleTable(2)*f0*t38);

mi0t = mod38.*cos(2*pi*ScaleTable(3)*f0*t38);

fa0t = mod38.*cos(2*pi*ScaleTable(4)*f0*t38);

so0t = mod38.*cos(2*pi*ScaleTable(5)*f0*t38);

la0t = mod38.*cos(2*pi*ScaleTable(6)*f0*t38);

si0t = mod38.*cos(2*pi*ScaleTable(7)*f0*t38);

do1t = mod38.*cos(2*pi*ScaleTable(8)*f0*t38);

re1t = mod38.*cos(2*pi*ScaleTable(9)*f0*t38);

mi1t = mod38.*cos(2*pi*ScaleTable(10)*f0*t38);

fa1t = mod38.*cos(2*pi*ScaleTable(11)*f0*t38);

so1t = mod38.*cos(2*pi*ScaleTable(12)*f0*t38);

la1t = mod38.*cos(2*pi*ScaleTable(13)*f0*t38);

si1t = mod38.*cos(2*pi*ScaleTable(14)*f0*t38);

%3/4 notes

so0q = mod34.*cos(2*pi*ScaleTable(5)*f0*t34);

si0q = mod34.*cos(2*pi*ScaleTable(7)*f0*t34);

%1/4 notes

do0f = mod4.*cos(2*pi*ScaleTable(1)*f0*t4);

re0f = mod4.*cos(2*pi*ScaleTable(2)*f0*t4);

mi0f = mod4.*cos(2*pi*ScaleTable(3)*f0*t4);

fa0f = mod4.*cos(2*pi*ScaleTable(4)*f0*t4);

so0f = mod4.*cos(2*pi*ScaleTable(5)*f0*t4);

la0f = mod4.*cos(2*pi*ScaleTable(6)*f0*t4);

si0f = mod4.*cos(2*pi*ScaleTable(7)*f0*t4);

do1f = mod4.*cos(2*pi*ScaleTable(8)*f0*t4);

re1f = mod4.*cos(2*pi*ScaleTable(9)*f0*t4);

mi1f = mod4.*cos(2*pi*ScaleTable(10)*f0*t4);

fa1f = mod4.*cos(2*pi*ScaleTable(11)*f0*t4);

so1f = mod4.*cos(2*pi*ScaleTable(12)*f0*t4);

la1f = mod4.*cos(2*pi*ScaleTable(13)*f0*t4);

si1f = mod4.*cos(2*pi*ScaleTable(14)*f0*t4);

blkf = zeros(1,i);

%1/8 notes

do0e = mod8.*cos(2*pi*ScaleTable(1)*f0*t8);

re0e = mod8.*cos(2*

總結

以上是生活随笔為你收集整理的matlab 天空之城,matlab演奏天空之城的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內容還不錯，歡迎將生活随笔推薦給好友。