重庆大学通信学院数字信号处理实验电话号码

2.（1）

clear all;clc;

tm=[1,2,3,65;4,5,6,66;7,8,9,67;42,0,35,68]; n=205;

f1=[697,770,852,941]; f2=[1209,1336,1477,1633]; tn=input('输入6位电话号码='); for m=1:8

d=fix(tn/10^(8-m)); tn=tn-d*10^(8-m); for p=1:4 for q=1:4 if

tm(p,q)==abs(d);break;end end if

tm(p,q)==abs(d);break;end end n=0:1023;

x=sin(2*pi*n*f1(p)/8000)+sin(2*pi*n*f2(q)/8000); xk=fft(x); subplot(4,2,m);

plot(8000*n/1024,abs(xk)); axis([600 1700 0 1000]); sound(x,8000); pause(0.5) end

(2)

clear all;clc;

tm=[1,2,3,65;4,5,6,66;7,8,9,67;42,0,35,68]; n=205;

f1=[697,770,852,941]; f2=[1209,1336,1477,1633]; tn=input('输入6位电话号码='); for m=1:8

d=fix(tn/10^(8-m)); tn=tn-d*10^(8-m); for p=1:4 for q=1:4 if

tm(p,q)==abs(d);break;end end if

tm(p,q)==abs(d);break;end end n=0:204 ;

x=sin(2*pi*n*f1(p)/8000)+sin(2*pi*n*f2(q)/8000);

k=[18,20,22,24,31,34,38,42]; xgk=goertzel(x,k+1); subplot(4,2,m);

stem(8000*k/205,abs(xgk)); axis([600 1700 0 100]); sound(x,8000); pause(0.5) end

第1页