北邮双音多频MATLAB实验报告.docx

实验 双音多频拨号系统matlab实现xxxxxxxxxx班xx号大学霸实验内容：把八位电话号码通过DTMF编码生成调制信号，并在接收端检测。技术指标：根据ITU Q.23建议，DTMF信号的技术指标是：传送/接收率为每秒10个号码，或每个号码100ms。每个号码传送过程中，信号存在时间至少45ms，且不多于55ms，100ms的其余时间是静音。在每个频率点上允许有不超过1.5%的频率误差。任何超过给定频率3.5%的信号，均被认为是无效的，拒绝接收。（其中关键是不同频率的正弦波的产生。可以使用查表方式模拟产生两个不同频率的正弦波。正弦表的制定要保证合成信号的频率误差在1.5%以内，同时是取样点数尽量少）dtmf信号产生:双音多频调制DTMF实质是将两个不同频率的正弦波予以混合相加即可得到，所以只要能够产生两个不同点的正弦波的函数值，再相加就可以了DTMF生成器其中 b21=b22=1Matlab代码m=zeros(1,205);y=[m;m;m;m;m;m;m;m;m;m];a10=[0.6736 0.5204 0.5204 0.5204 0.5686 0.5686 0.5686 0.6203 0.6203 0.6203];a20=[0.8671 0.8131 0.8671 0.9168 0.8131 0.8671 0.9168 0.8131 0.8671 0.9168];b11=[-1.4782 -1.7077 -1.7077 -1.7077 -1.6453 -1.6453 -1.6453 -1.5687 -1.5687 -1.5687];b21=[-0.9964 -1.1641 -0.9964 -0.7986 -1.1641 -0.9964 -0.7986 -1.1641 -0.9964 -0.7986];h=impseq(0,0,204); %长度为205的单位冲激函数n=0:204;for i=1:10 a1=[0,a10(i)];a2=[0,a20(i)];b1=[1,b11(i),1];b2=[1,b21(i),1]; y(i,:)=filter(a1,b1,h)+filter(a2,b2,h); %两个冲激响应相叠加endtone=y(d+1,:); %tone为相应的按键音（正弦函数）figure(1)for t=1:10 subplot(5,2,t) plot(n/205,y(t,:));title(strcat(数字键,num2str(t-1)));endDTMF信号检测（GOERTZEL算法）z-1MATLAB代码f=[697 770 852 941 1209 1336 1477 1633];fs=8000;fn=fs/N;k=f/fn;x=[ysin,zeros(1,N)];tic;%timer startVk(1)=0;Vk(2)=0;for i1=1:8 w=2*cos(2*pi*k(i1)/N); for j1=3:N Vk(j1)=w.*Vk(j1-1)-Vk(j1-2)+x(j1); end Xk(i1)=Vk(N).^2+Vk(N-1).^2-w*Vk(N)*Vk(N-1);endX=Xk/100;主程序N=205;tm=[49,50,51;52,53,54;55,56,57;0,48,0];f=[697 770 852 941 1209 1336 1477 1633];d_in=input(请输入电话号码: ,s); %输入信号d_infor l=1:8 d=d_in(l)-48; d_tone=dsin(d); %产生拨号音 sound(d_tone,8000); %播放拨号音 pause(0.2) X=GOERT(d_tone,N); %GOERT算法求出接收端信号的频率分量 value=abs(X); figure(2) subplot(4,2,l) stem(f,value,.);grid; set(gca,XTick,f); title(FFT x(n));xlabel(频率k/Hz);ylabel(幅度|X(k)|); limit=50; %设置门限值 for s=5:8 if value(s)limit break; %高频门限检测 end end for r=1:4 if value(r)limit break;