matlab仿真.rar

matlab仿真.rar

% Most shining national wind//最炫民族风 on Matlab 
% The Modification is from “canon“ not by me 

fs = 44100; % sample rate 
dt = 1/fs; 

T16 = 0.125; 

t16 = [0:dt:T16]; 
[temp k] = size（t16）; 

t4 = linspace（04*T164*k）; 
t8 = linspace（02*T162*k）; 

[temp i] = size（t4）; 
[temp j] = size（t8）; 

% Modification functions 
mod4=（t4.^4）.*exp（-30*（t4.^0.5））; 
mod4=mod4*（1/max（mod4））; 
mod8=（t8.^4）.*exp（-50*（t8.^0.5））; 
mod8=mod8*（1/max（mod8））; 
mod16=（t16.^4）.*exp（-90*（t16.^0.5））; 
mod16=mod16*（1/max（mod16））; 

f0 = 2*146.8; % reference frequency 

ScaleTable = [2/3 3/4 5/6 15/16 ... 
1 9/8 5/4 4/3 3/2 5/3 9/5 15/8 ... 
2 9/4 5/2 8/3 3 10/3 15/4 4 ... 
1/2 9/16 5/8]; 

% 1/4 notes 
do0f = mod4.*cos（2*pi*ScaleTable（21）*f0*t4）; 
re0f = mod4.*cos（2*pi*ScaleTable（22）*f0*t4）; 
mi0f = mod4.*cos（2*pi*ScaleTable（23）*f0*t4）; 

fa0f = mod4.*cos（2*pi*ScaleTable（1）*f0*t4）; 
so0f = mod4.*cos（2*pi*ScaleTable（2）*f0*t4）; 
la0f = mod4.*cos（2*pi*ScaleTable（3）*f0*t4）; 
ti0f = mod4.*cos（2*pi*ScaleTable（4）*f0*t4）; 
do1f = mod4.*cos（2*pi*ScaleTable（5）*f0*t4）; 
re1f = mod4.*cos（2*pi*ScaleTable（6）*f0*t4）; 
mi1f = mod4.*cos（2*pi*ScaleTable（7）*f0*t4）; 
fa1f = mod4.*cos（2*pi*ScaleTable（8）*f0*t4）; 
so1f = mod4.*cos（2*pi*ScaleTable（9）*f0*t4）; 
la1f = mod4.*cos（2*pi*ScaleTable（10）*f0*t4）; 
tb1f = mod4.*cos（2*pi*ScaleTable（11）*f0*t4）; 
ti1f = mod4.*cos（2*pi*ScaleTable（12）*f0*t4）; 
do2f = mod4.*cos（2*pi*ScaleTable（13）*f0*t4）; 
re2f = mod4.*cos（2*pi*ScaleTable（14）*f0*t4）; 
mi2f = mod4.*cos（2*pi*ScaleTable（15）*f0*t4）; 
fa2f = mod4.*cos（2*pi*ScaleTable（16）*f0*t4）; 
so2f = mod4.*cos（2*pi*ScaleTable（17）*f0*t4）; 
la2f = mod4.*cos（2*pi*ScaleTable（18）*f0*t4）; 
ti2f = mod4.*cos（2*pi*ScaleTable（19）*f0*t4）; 
do3f = mod4.*cos（2*pi*ScaleTable（20）*f0*t4）; 
blkf = zeros（1i）; 

% 1/8 notes 
do0e = mod8.*cos（2*pi*ScaleTable（21）*f0*t8）; 
re0e = mod8.*cos（2*pi*ScaleTable（22）*f0*t8）; 
mi0e = mod8.*cos（2*pi*ScaleTable（23）*f0*t8）; 

fa0e = mod8.*cos（2*pi*ScaleTable（1）*f0*t8）; 
so0e = mod8.*cos（2*pi*ScaleTable（2）*f0*t8）; 
la0e = mod8.*cos（2*pi*ScaleTable（3）*f0*t8）; 
ti0e = mod8.*cos（2*pi*ScaleTable（4）*f0*t8）; 
do1e = mod8.*cos（2*pi*ScaleTable（5）*f0*t8）; 
re1e = mod8.*cos（2*pi*ScaleTable（6）*f0*t8）; 
mi1e = mod8.*cos（2*pi*ScaleTable（7）*f0*t8）; 
fa1e = mod8.*cos（2*pi*ScaleTable（8）*f0*t8）; 
so1e = mod8.*cos（2*pi*ScaleTable（9）*f0*t8）; 
la1e = mod8.*cos（2*pi*ScaleTable（10）*f0*t8）; 
tb1e = mod8.*cos（2*pi*ScaleTable（11）*f0*t8）; 
ti1e = mod8.*cos（2*pi*ScaleTable（12）*f0*t8）; 
do2e = mod8.*cos（2*pi*ScaleTable（13）*f0*t8）; 
re2e = mod8.*cos（2*pi*ScaleTable（14）*f0*t8）; 
mi2e = mod8.*cos（2*pi*ScaleTable（15）*f0*t8）; 
fa2e = mod8.*cos（2*pi*ScaleTable（16）*f0*t8）; 
so2e = mod8.*cos（2*pi*ScaleTable（17）*f0*t8）; 
la2e = mod8.*cos（2*pi*ScaleTable（18）*f0*t8）; 
ti2e = mod8.*cos（2*pi*ScaleTable（19）*f0*t8）; 
do3e = mod8.*cos（2*pi*ScaleTable（20）*f0*t8）; 
blke = zeros（1j）; 

% 1/16 notes 
do0s = mod16.*cos（2*pi*ScaleTable（21）*f0*t16）; 
re0s = m

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      39445  2014-06-26 20:30  matlab仿真\boost.mdl

     文件      28760  2014-06-25 20:35  matlab仿真\D.mdl

     文件     231276  2014-09-12 21:32  matlab仿真\MATLAB仿真.docx

     文件      74006  2014-06-25 20:49  matlab仿真\pvv.mdl

     文件      75460  2014-06-26 10:36  matlab仿真\pvv1.mdl

     文件      20992  2014-06-25 20:41  matlab仿真\pvv_acc.mexw32

     文件      45118  2014-06-30 20:14  matlab仿真\pwm.mdl

     文件      30876  2014-06-30 16:59  matlab仿真\raodong.mdl

     文件       6999  2014-06-25 19:59  matlab仿真\Untitled2.m

     文件      34438  2014-06-25 20:19  matlab仿真\X1.mdl

     文件      26883  2014-06-25 20:17  matlab仿真\X2.mdl

     目录          0  2014-09-12 21:36  matlab仿真

----------- ---------  ---------- -----  ----

               614253                    12