cepstral.m
View Options

			files = {"audio1.wav", "audio4.wav"}
			exclude = 5

			figure
			hax1 = axes
			grid on
			title('Cepstr')
			xlabel('Quefrency')
			ylabel('Power Cepstrum')
			hold(hax1)
			lgd = legend


			for i=1:length(files)
			[sf, t, k, S, FS] = load_audio(files{i});
			S = S(:,1);
			cepstr = mypowercepstrum(S);

			figure
			hax2 = axes
			title(hax2, sprintf('Fourier transform of %s', files{i}))
			xlabel(hax2, 'Frequency')
			ylabel(hax2, 'Signal sample')
			semilogy(hax2, k, abs(FS));

			% normalize cepstr
			t = t(exclude:end-exclude);
			cepstr = cepstr(exclude:end-exclude);

			norm = discr_integral(t, cepstr);
			cepstr = cepstr / norm;

			plot(hax1, t, cepstr, 'DisplayName', files{i});

			foutname = sprintf("output/%s_ceps.out", files{i});
			fout = fopen(foutname, 'w');
			fprintf(fout, "t C\n");
			for i=1:length(t)
			fprintf(fout, "%g %g\n", t(i), cepstr(i));
			end
			fclose(fout);
			end

compphys/exo2/cepstral.m4b851d2c5111master