--- Saving session to:
    ECE481_02-Apr-2002.txt
--- Processed startup.m ---
; ; ; ; ; ; ; 
edit fm6
fm6
Warning: Data clipped during write to file:\temp\junk.wav
> In C:\MatLab6p1\toolbox\matlab\audio\wavwrite.m (PCM_Quantize) at line 239
  In C:\MatLab6p1\toolbox\matlab\audio\wavwrite.m (write_wavedat) at line 257
  In C:\MatLab6p1\toolbox\matlab\audio\wavwrite.m at line 104
  In c:\personal\class\2001-02\spring\ece481\matlab\fm6.m at line 33
type fm6

%function fm6(fc,h,mi)
% FM synthesis demo #6: Double-carrier FM (Dodge & Jerse p. 127)

fs=20e3;
dur=0.3;
tt=0:1/fs:dur;

% parameters
fc1 = 300;		%carrier frequenciese
fc2 = 1500;
fm = fc1;		%modulation frequency
I1 = 2.66;		%modulation indices
I2 = 1.8;
a1 = 1;		% amplitudes
a2 = 0.2;

% Morrill double-carrier instrument
if 1
   a1=gen4([0 0 -4; .03 1 0; .85 1 -4; 1 0 0],dur,fs);
   a2=0.2*gen4([0 0 -4; .03 1 0; .7 1 -4; 1 0 0],dur,fs);
   I1=2.66*gen4([0 0 -4; .03 1 0; .99 1 -4; 1 0 0],dur,fs);
   I2=(1.8/2.66)*I1;
end


% output waveform
y = a1.*sin(2*pi*fc1*tt + I1.*sin(2*pi*fm*tt)) +...
    a2.*sin(2*pi*fc2*tt + I2.*sin(2*pi*fm*tt));
y=y/max(abs(y));
calspec(y,[],fs)

soundsc(y,fs) 
wavwrite(y,fs,16,'\temp\junk.wav')


edit fm5
fm5
??? Error using ==> wavwrite
Sorry. No help in figuring out the problem . . .

Error in ==> c:\personal\class\2001-02\spring\ece481\matlab\fm5.m
On line 41  ==> wavwrite(y,fs,16,'\temp\junk.wav')

wavwrite(y,fs,16,'\temp\junk.wav')
type fm5

%function fm4(fc,h,mi)
% FM synthesis demo #5: Cascade FM (Moore p. 332)

fs=40e3;
dur=3;
tt=0:1/fs:dur;

% parameters
fc = 60;	%carrier frequency
fm1 = sqrt(2)*fc;	%modulating frequency
fm2 = pi*fm1;
mi1 = 2;		%modulation index
mi2 = 4;

% amplitude envelope
a=gen4([0 0 -10; 0.01 1 -1; 1 0 0],dur,fs);

% I1 envelope
p6=sqrt(2);
i1=gen4([0 1*p6 1; 0.045 12*p6 -2; 1 0.5*p6 0],dur,fs);

% I2 envelope
p7=pi;
i2=gen4([0 1*p7 1; 0.03 8*p7 -1; 1 1*p7 0],dur,fs);


% output waveform
y = a.*sin(2*pi*fc*tt + i1.*sin(2*pi*fm1*tt + i2.*sin(2*pi*fm2*tt)));

calspec(y,512,fs)
if 0
[b,f,t]=specgram(y,[],fs);
imagesc(t,f/1000,abs(b))
axis('xy')
xlabel('time (s)')
ylabel('frequency (kHz)')
title(sprintf('FM: fc=%dHz, H=%d, I=%d',fc,h,mi))
end

sound(y,fs) 
wavwrite(y,fs,16,'\temp\junk.wav')


exit