annotate toolboxes/MIRtoolbox1.3.2/MIRToolbox/@mirscalar/mirsave.m @ 0:cc4b1211e677 tip

initial commit to HG from Changeset: 646 (e263d8a21543) added further path and more save "camirversion.m"
author Daniel Wolff
date Fri, 19 Aug 2016 13:07:06 +0200
parents
children
rev   line source
Daniel@0 1 function mirsave(a)
Daniel@0 2
Daniel@0 3 d = get(a,'Data');
Daniel@0 4 n = get(a,'Name');
Daniel@0 5 t = get(a,'Title');
Daniel@0 6 c = get(a,'Channels');
Daniel@0 7 fp = get(a,'FramePos');
Daniel@0 8 out = cell(1,length(d));
Daniel@0 9
Daniel@0 10 for k = 1:length(d)
Daniel@0 11 dk = d{k};
Daniel@0 12 if not(iscell(dk))
Daniel@0 13 dk = {dk};
Daniel@0 14 end
Daniel@0 15 out = [];
Daniel@0 16 for l = 1:size(dk{1},3)
Daniel@0 17 for i = 1:length(dk)
Daniel@0 18 di = dk{i};
Daniel@0 19 synth = zeros(1,ceil((fp{k}{i}(end)-fp{k}{i}(1))*44100)+1);
Daniel@0 20 for j = 1:size(di,2)
Daniel@0 21 if iscell(di)
Daniel@0 22 dj = di{j};
Daniel@0 23 else
Daniel@0 24 dj = di(j);
Daniel@0 25 end
Daniel@0 26 if not(isempty(dj))
Daniel@0 27 k1 = floor((fp{k}{i}(1,j)-fp{k}{i}(1))*44100)+1;
Daniel@0 28 k2 = floor((fp{k}{i}(2,j)-fp{k}{i}(1))*44100)+1;
Daniel@0 29 synth(k1:k2) = synth(k1:k2) ...
Daniel@0 30 + sum(sin(2*pi*dj*(0:k2-k1)/44100),1).*hann(k2-k1+1)';
Daniel@0 31 end
Daniel@0 32 end
Daniel@0 33 out = [out synth];
Daniel@0 34 if size(dk{1},3)>1
Daniel@0 35 out = [out rand(1,10)];
Daniel@0 36 end
Daniel@0 37 end
Daniel@0 38 end
Daniel@0 39 fout = miraudio(out,44100);
Daniel@0 40 mirsave(fout,[n{k},'.',t]);
Daniel@0 41 end
Daniel@0 42
Daniel@0 43
Daniel@0 44 function oldmirsave % not used anymore
Daniel@0 45 d = get(a,'Data');
Daniel@0 46 nf = length(d);
Daniel@0 47 fp = get(a,'FramePos');
Daniel@0 48 nm = get(a,'Name');
Daniel@0 49 t = get(a,'Title');
Daniel@0 50 for i = 1:nf
Daniel@0 51 nmi = nm{i};
Daniel@0 52 di = d{i}{1};
Daniel@0 53 fpi = fp{i}{1};
Daniel@0 54
Daniel@0 55 %Let's remove the extension from the original files
Daniel@0 56 if length(nmi)>3 && strcmpi(nmi(end-3:end),'.wav')
Daniel@0 57 nmi(end-3:end) = [];
Daniel@0 58 elseif length(nmi)>2 && strcmpi(nmi(end-2:end),'.au')
Daniel@0 59 nmi(end-2:end) = [];
Daniel@0 60 end
Daniel@0 61 n = [nmi,'.',lower(t),'.txt'];
Daniel@0 62
Daniel@0 63 fid = fopen(n, 'wt');
Daniel@0 64 fprintf(fid,'Frame_start Frame_end Data \n');
Daniel@0 65
Daniel@0 66 for j = 1:length(di)
Daniel@0 67 fprintf(fid,'%g %g %g \n',fpi(1,j),fpi(2,j),di(j));
Daniel@0 68 end
Daniel@0 69 fclose(fid);
Daniel@0 70 disp([n,' saved.']);
Daniel@0 71 end