--- Saving session to: ECE481_16-Apr-2002.txt --- Processed startup.m --- doc ; ; ; ; ; ; ; fs=10e3 fs = 10000 dur=2 dur = 2 ss=randn(1,dur*fs); soundsc(ss,fs) soundsc(ss,fs) calspec(ss,[],fs) n=5 n = 5 fc= 1000 / (fs/2) fc = 0.2000 [b,a] = butter(n,fc) b = Columns 1 through 4 0.0013 0.0064 0.0128 0.0128 Columns 5 through 6 0.0064 0.0013 a = Columns 1 through 4 1.0000 -2.9754 3.8060 -2.5453 Columns 5 through 6 0.8811 -0.1254 y = filter (b,a,ss); soundsc(ss,fs) soundsc(y,fs) calspec(y,[],fs) [H,F] = freqz (b,a,fs/2,fs); plot(F,abs(H)) b b = Columns 1 through 4 0.0013 0.0064 0.0128 0.0128 Columns 5 through 6 0.0064 0.0013 a a = Columns 1 through 4 1.0000 -2.9754 3.8060 -2.5453 Columns 5 through 6 0.8811 -0.1254 clear edit lpfsweep lpfsweep Original sound ... (press key to continue) Corner frequency trajectory... (press key to continue) Filtered sound. fschange('c:\personal\class\2001-02\spring\ece481\matlab\lpfsweep.m'); clear lpfsweep lpfsweep fc = Columns 1 through 4 0.1000 0.1127 0.1271 0.1433 Columns 5 through 8 0.1616 0.1822 0.2054 0.2316 Columns 9 through 12 0.2611 0.2944 0.3319 0.3742 Columns 13 through 16 0.4220 0.4757 0.5364 0.6048 Columns 17 through 20 0.6818 0.7688 0.8668 0.9772 Original sound ... (press key to continue) Corner frequency trajectory... (press key to continue) Filtered sound. exit