function signal = harmonic(pitch, no)

% function to create harmonic pitch of specified fundamental frequency and number of harmonics

Fs=44100;
t=0:1/Fs:1;
phase=0;
for n=1:no % or n=2:no for missing fundamental
    phase=2*pi*rand;
    z(n,:)=sin(2*pi*n*pitch*t+phase);
end
pitch_comp=sum(z);
signal=0.2*pitch_comp/std(pitch_comp);
wavname=[num2str(pitch) '_' num2str(no)];
wavwrite(signal, 44100, wavname);

Published with MATLAB® 7.1