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

实验 双音多频拨号系统matlab实现xxxxxxxxxx班xx号大学霸实验内容： 把八位电话号码通过DTMF编码生成调制信号，并在接收端检测。 技术指标： 根据ITU Q.23建议，DTMF信号的技术指标是：传送/接收率为每秒10个号码，或每个号码100ms。 每个号码传送过程中，信号存在时间至少45ms，且不多于55ms，100ms的其余时间是静音。 在每个频率点上允许有不超过1.5%的频率误差。任何超过给定频率3.5%的信号，均被认为是无效的，拒绝接收。（其中关键是不同频率的正弦波的产生。可以使用查表方式模拟产生两个不同频率的正弦波。正弦表的制定要保证合成信号的频率误差在1.5%以内，同时是取样点数尽量少）dtmf信号产生:双音多频调制DTMF实质是将两个不同频率的正弦波予以混合相加即可得到，所以只要能够产生两个不同点的正弦波的函数值，再相加就可以了Z-1Z-1Z-1Z-1-b11-b21-b21-b22Z-1a20(n)Z-1a10DAC低通滤波器按键音产生低频音的IIR DF产生高频音的IIR DFDTMF生成器其中 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算法）VKyK(n)z-1z-1X(n)2cos2kN-exp(2kjN)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; %低频门限检测 end end d_out(l)=tm(r,s-4); %还原信号enddisp(双音多频信号已生成并发出)disp(接收端检测到的号码: )disp(setstr(d_out)程序调试截图信号频谱图：观察上面的频谱图，说明接收信号正好是78934620，输入输出一致，程序运行无误实验总结这次matlab实验明显比之前的有难度，在掌握了第四章