m-files compile clean
This commit is contained in:
@@ -72,10 +72,10 @@ end;
|
|||||||
close all
|
close all
|
||||||
plot(1:nsin, blepm(fix(nharm/2), :)); grid;
|
plot(1:nsin, blepm(fix(nharm/2), :)); grid;
|
||||||
|
|
||||||
wavwrite(0.5*vtri, fs, 'tri.wav');
|
wavw16(0.5*vtri, fs, 'tri.wav');
|
||||||
wavwrite(0.5*vsqr, fs, 'sqr.wav');
|
wavw16(0.5*vsqr, fs, 'sqr.wav');
|
||||||
wavwrite(0.5*vsaw, fs, 'saw.wav');
|
wavw16(0.5*vsaw, fs, 'saw.wav');
|
||||||
wavwrite(0.5*vblep, fs, 'blep.wav');
|
wavw16(0.5*vblep, fs, 'blep.wav');
|
||||||
|
|
||||||
function h = lagrange(N, delay)
|
function h = lagrange(N, delay)
|
||||||
%LAGRANGE h=lagrange(N,delay) returns order N FIR
|
%LAGRANGE h=lagrange(N,delay) returns order N FIR
|
||||||
|
|||||||
@@ -82,8 +82,8 @@ close all
|
|||||||
figure;
|
figure;
|
||||||
plot(1:L, vsaw, 1:L, vblit); grid;
|
plot(1:L, vsaw, 1:L, vblit); grid;
|
||||||
|
|
||||||
wavwrite(0.5*vtri, fs, 'tri.wav');
|
wavw16(0.5*vtri, fs, 'tri.wav');
|
||||||
wavwrite(0.5*vsqr, fs, 'sqr.wav');
|
wavw16(0.5*vsqr, fs, 'sqr.wav');
|
||||||
wavwrite(0.5*vsaw, fs, 'saw.wav');
|
wavw16(0.5*vsaw, fs, 'saw.wav');
|
||||||
wavwrite(0.5*vblit, fs, 'blit.wav');
|
wavw16(0.5*vblit, fs, 'blit.wav');
|
||||||
|
|
||||||
|
|||||||
@@ -1,4 +1,4 @@
|
|||||||
function eval_blit()
|
function eval_blit2()
|
||||||
L = 64000;
|
L = 64000;
|
||||||
Nb = 65;
|
Nb = 65;
|
||||||
Nos = 8;
|
Nos = 8;
|
||||||
@@ -9,7 +9,7 @@ df = (fend-fstart)/L
|
|||||||
fcut = 16000;
|
fcut = 16000;
|
||||||
|
|
||||||
x = (0:Nb*Nos-1) - (Nb*Nos-1)/2;
|
x = (0:Nb*Nos-1) - (Nb*Nos-1)/2;
|
||||||
H = sinc(2*fcut/fs.*x/Nos).*Kaiser(Nos*Nb, 8)';
|
H = sinc(2*fcut/fs.*x/Nos).*kaiser(Nos*Nb, 8)';
|
||||||
for ii=0:Nos-1,
|
for ii=0:Nos-1,
|
||||||
for jj=0:Nb-1
|
for jj=0:Nb-1
|
||||||
Hs(ii+1,jj+1) = H(Nos-ii+Nos*jj);
|
Hs(ii+1,jj+1) = H(Nos-ii+Nos*jj);
|
||||||
@@ -58,7 +58,7 @@ for i = 1:L,
|
|||||||
f = f + df;
|
f = f + df;
|
||||||
|
|
||||||
end;
|
end;
|
||||||
wavwrite(0.7*blit, fs, 'blit.wav');
|
wavw16(0.7*blit, fs, 'blit.wav');
|
||||||
|
|
||||||
for ii = 1:Nos
|
for ii = 1:Nos
|
||||||
step(ii,:) = k(ii).*filter(Hs(ii,:), 1, [ones(1,Nb)]);
|
step(ii,:) = k(ii).*filter(Hs(ii,:), 1, [ones(1,Nb)]);
|
||||||
|
|||||||
@@ -1,4 +1,4 @@
|
|||||||
function eval_blit()
|
function eval_blit3()
|
||||||
L = 64000;
|
L = 64000;
|
||||||
Nb = 137;
|
Nb = 137;
|
||||||
Nphases = 256;
|
Nphases = 256;
|
||||||
@@ -10,7 +10,7 @@ fcut = 16000;
|
|||||||
|
|
||||||
x = (0:Nb*Nphases-1) - (Nb*Nphases-1)/2;
|
x = (0:Nb*Nphases-1) - (Nb*Nphases-1)/2;
|
||||||
H = sinc(2*fcut/fs.*x/Nphases);%.*Kaiser(Nphases*Nb, 8)';
|
H = sinc(2*fcut/fs.*x/Nphases);%.*Kaiser(Nphases*Nb, 8)';
|
||||||
Hw = Kaiser(Nb, 18)';
|
Hw = kaiser(Nb, 18)';
|
||||||
|
|
||||||
for ii=0:Nphases-1,
|
for ii=0:Nphases-1,
|
||||||
for jj=0:Nb-1
|
for jj=0:Nb-1
|
||||||
@@ -62,7 +62,7 @@ plot(saw); grid;
|
|||||||
figure;
|
figure;
|
||||||
plot(saw_nodc); grid;
|
plot(saw_nodc); grid;
|
||||||
|
|
||||||
wavwrite(0.8*blit, fs, 'blit.wav');
|
wavw16(0.8*blit, fs, 'blit.wav');
|
||||||
wavwrite(0.7*saw_nodc, fs, 'saw.wav');
|
wavw16(0.7*saw_nodc, fs, 'saw.wav');
|
||||||
wavwrite(0.5*sqr_nodc, fs, 'sqr.wav');
|
wavw16(0.5*sqr_nodc, fs, 'sqr.wav');
|
||||||
|
|
||||||
|
|||||||
@@ -80,9 +80,9 @@ plot(vc4); grid;
|
|||||||
figure;
|
figure;
|
||||||
plot(vc3); grid;
|
plot(vc3); grid;
|
||||||
|
|
||||||
wavwrite(0.5*vc4, fs, 'c4.wav');
|
wavw16(0.5*vc4, fs, 'c4.wav');
|
||||||
wavwrite(0.5*vc3, fs, 'c3.wav');
|
wavw16(0.5*vc3, fs, 'c3.wav');
|
||||||
wavwrite(0.5*vsqr, fs, 'sqr.wav');
|
wavw16(0.5*vsqr, fs, 'sqr.wav');
|
||||||
wavwrite(0.5*vsaw, fs, 'saw.wav');
|
wavw16(0.5*vsaw, fs, 'saw.wav');
|
||||||
wavwrite(0.5*vblit, fs, 'blit.wav');
|
wavw16(0.5*vblit, fs, 'blit.wav');
|
||||||
|
|
||||||
|
|||||||
Reference in New Issue
Block a user