万方数据
多采样率数字信号处理及其MATLAB仿真
作者：黄硕， 魏亚楠， 安永丽
作者单位：唐山钢铁股份有限公司,唐山,063016
刊名：
科技资讯
英文刊名：SCIENCE & TECHNOLOGY INFORMATION
年，卷(期)：2009，(23)
引用次数：0次
参考文献(3条)
1.杨小牛.楼才义.徐建良软件无线电原理与应用 2005
2.李忠琦.凌翔.胡剑浩软件无线电架构研究[期刊论文]-电信科学 2007(7)
3.尹健华试论软件无线电技术及其应用[期刊论文]-企业技术开发(学术版) 2007(8)
相似文献(10条)
1.学位论文赵启敏中频采样技术的分析与研究2004
该课题结合数字软件化雷达的研制,研究了数字软件化雷达中频采样技术的实现以及对雷达主要技术指标的影响.该论文针对传统模拟相参正交采样技术存在的不足,论述了基于A/D变换和数字下变频的中频采样方法,并在此基础上设计了中频采样数据采集卡,并对该硬件进行了调试和试验,试验结果证明,中频采样技术比传统模拟相参正交采样技术更具优势,较好的解决了传统模拟相参正交采样中存在的幅相误差问题,以及该采集卡具有小的孔径抖动,可以满足中频采样的要求.该论文在中频采样技术中首先研究了数据采集技术对雷达性能的影响,接着根据目前数字下变频器件自身的限制不能适应高速数据流的问题,详细研究了利用欠采样技术的镜频加数字下变频实现解调的方法,以及一种利用多速率信号处理技术将抽取和滤波提前的数字下变频的高效结构,通过仿真证明此两种方法都能较好的解决硬件本身限制与高速数据流不匹配的问题,并通过分析得出此数字下变频的高效结构的运算量大大低于传统数字下变频的运算量.此外该论文还着重讨论了孔径抖动对雷达各项性能的影响.
2.期刊论文张明珊.孟利民.ZHANG Ming-shan.MENG Li-min基于频域采样技术的软件无线电接收机-浙江工业大学学报2005,33(1)
目前软件无线电面临的一个难题是如何对高工作频带内的射频信号进行直接模/数转换.利用频域采样技术提出了一种接收信号进行处理的方法,并用数学理论证明了它的可行性,最后还给出了软件无线电接收机模型.其关键思想是提取接收信号的频域成份,然后在频域中对信号进行处理.这种方法大大降低了A/D转换器的要求,从而使得实现软件无线电接收机成为可能,对当前微电子工艺下的软件无线电系统设计带来很大的理论意义和实用价值,而且克服了传统Rake接收机的一些缺点,特别适合于多径丰富的无线环境.
3.学位论文杨清海软件无线电的功能实现2001
1992年，JeoMitola提出了软件无线电的概念，很快引起了国际通信界的关注。软件无线电结构的关键是在尽可能靠近天线的地方使用宽带A／D和
D／A变换器，将尽可能多的无线电功能用软件来定义，从而实现电台在各种网络中的通用性及电台功能升级换代的连续性，软件无线电已成为无线通信的一个主要发展方向。特别是近年来，软件无线电已经不再仅仅局限于军事方面，在GSMMOU会议中，软件无线电被描述成GSM继续发展进步的基础，甚至被称为第三代（3G）全球移动通信实现的技术基础。本文主要探讨软件无线电思想在接收机设计中的应用，论证了系统硬件实现方案和软件实现方案。重点讨论了用到的信号采样技术和数字信号处理技术，包括多速率信号处理、FIR滤波器的多相结构、低通滤波、免混频正交解调和信号的带通采样技术。最后优化了解调算法，利用我们的试验平台实现了AM、FM、SSB和ASK、FSK、PSK信号解调。
4.学位论文洪亮高速并行交替采样ADC系统的研究与实现2009
模数转换器(ADC)是数字信号处理系统的关键组成部分，广泛应用于通信、雷达、测试仪器等领域。随着超宽带雷达技术研究的深入和软件无线电技术的发展，对ADC的速度和精度的要求越来越高，ADC已经成为现代信号处理的瓶颈。在给定的工艺下，ADC工作的最大采样速率受限于它的分辨率，单片ADC芯片很难同时满足高速高精度的要求，而并行交替采样ADC(TIADC)结构是突破这一瓶颈的有效方法之一。
这种方法在前端利用M片采样率为fs/M的ADC并行交替采样，在后端进行拼接使得整个系统的采样率达到fs。然而受到制造工艺的局限，通道失配误差如偏置误差、增益误差、时间偏差和带宽失配误差的存在，将严重降低系统的信纳比(SINAD)和无杂散动态范围(SFDR)。
本论文主要包括三方面的工作。首先，深入研究了并行交替采样技术，对TIADC结构的通道失配误差进行了全面的分析，特别是对带宽失配误差进行了建模分析，给出了四种通道失配误差联合作用于信号的信号频谱，以及系统设计时误差的容忍范围。
其次，通过合理的近似，提出了通道失配误差的测量算法和联合校正算法，其中关键的是时间偏差和带宽失配误差的联合估算与校正，它是在周期非均匀采样信号完美重构基础上提出来的，并通过仿真验证了算法的有效性。
最后，设计了一个基于并行交替采样技术的12bit420MSPS的高速数据采集系统，该系统由两片12bit210MSPS的AD9430组成。其中，结合系统设计进行的信号完整性分析对高速电路的设计具有一定的指导意义。
5.期刊论文王宏.刘丽.宋晓峰.WANG Hong.LIU Li.SONG Xiaofeng基于频域采样技术的软件无线电接收机-现代电子技术2006,29(23)
目前软件无线电面临的一个难题是如何对高工作频带内的射频信号进行直接模/数转换.利用频域采样技术提出了一种接收信号进行处理的方法,并用数学理论证明了他的可行性,最后还给出了软件无线电接收机模型.其关键思想是提取接收信号的频域成份,然后在频域中对信号进行处理.这种方法大大降低了A/D转换器的要求,对当前微电子工艺下的软件无线电系统设计有很大的理论意义和实用价值,而且克服了传统Rake接收机的一些缺点,特别适合于多径丰富的无线环境.
6.学位论文李裕多信道软件无线电接收机实现技术研究2003
软件无线电的基本思想是将宽带A/D及D/A尽可能靠近天线,将无线电台的各种功能在一个开放性、模块化的通用硬件平台上尽可能多的用软件来实现.软件无线电已成为移动通信中的关键技术之一.本文主要研究了软件无线电接收机中的相关理论及实现方案并进行了相应的系统仿真.本文首先深入讨论了软件无线电接收机的基本理论:采样技术、多速率信号处理和调制解调算法.在此基础上研究了下变频技术和带通采样技术在并行多信道接收机中的应用,提出了利用CIC,HBF和FIR级联设计下变频器的方案,并完成了系统仿真.然后深入研究了多相滤波技术在信道化接收机中的应用,推导和建立了实信号接收机的数学模型,给出了真实信道中心频率和带宽的计算公式,简要分析了算法复杂度,最后完成了基于此模型的4信道软件无线电接收机的系统仿真.本文所建立的两个系统作为后续研究的基础平台,可以利用其移植各种通信系统,并分析系统性能,具有一定的应用价值.
7.学位论文徐建良软件无线电理论与实践1999
该文通过对软件无线电的深入研究,提出了理想的软件无线电结构以及随器件技术发展而能逐步实现的三种折衷的软件无线电结构,提出并论证了实现软件无线电的具体方法和数学模型,从而构筑起整个软件无线电理论与实践的框架,为软件无线电的工程实现和进一步研究奠定了坚实的基础.其中,射频宽开采样技术和内插上变频技术是当前实现软件无线电的关键,也是这一研究领域首次提出的新思想,同时,对多速率数字信号处理与正交检测技术应用于软件无线电进行了详细的论述,另外,还对正交调制方法进行了介绍,这些技术将对软件无线电的发展产生积极的影响.该文最后是根据上核定是论所研制的对射频直接进行采样的软件无线电试验平台以验证理论的正确性和可行性.
8.学位论文曹丽君软件无线电在直放站中的应用研究2008
传统无线电系统对硬件依赖性强，而且信号的适应能力很差，软件无线电的出现改变了这种局面。其基本思想就是让A/D、D/A尽可能靠近天线，把硬件作为无线通信的基本平台，让尽可能多的无线通信功能用软件来实现，这样无线通信系统具有很好的通用性和灵活性，使系统的更新和升级变的更加方便。受目前器件水平的限制，现阶段尚无法对射频直接采样，而只能在中频进行采样，再利用数字信号处理技术来完成中频模块的设计。
论文讨论了软件无线电技术在直放站中的应用。首先从阐述软件无线电的理论基础出发，对数字上/下变频技术、多速率采样技术、几种常用的数字滤波器进行了分析和探讨。然后建立系统的实现方案，选择合适芯片设计系统的硬件电路，进而介绍了用FPGA实现中频信号处理的流程。最后利用M盯LAB中的simulink进行系统级仿真，验证了系统方案设计的可行性，并介绍了整个系统调试的过程和结果。
9.期刊论文张沫阳.卢继华.Zhang Moyang.Lu Jihua射频带通采样技术在软件无线电电台中的应用-军民两用技
术与产品2006(8)
理想的软件无线电电台要求对天线接收的模拟信号经过放大后直接采样,但由于技术所限很难实现,而多采用中频采样的方法.由于中频采样增加了较多的模拟电路环节,违背了软件无线电思想的初衷,进而又提出了射频带通采样的方法.本文对这种采样方法进行了分析,并对应用这一方法的实际系统进行了介绍.
10.学位论文李利软件接收机研究2002
软件无线电具有极大的灵活性和适应性,是通信领域的一个研究热点.为了对软件无线电进行探索性研究,该论文主要研究软件接收机中的中频数据采集技术、数字下变频算法、几种模拟调制信号的解调算法.带通采样技术的应用使得软件无线电接收机在中频采样成为可能,该文详细讨论了软件接收机中的带通均匀采样技术及软件无线电中ADC器件的关键技术特性.数字下变频(DDC)完成下变频、抽取滤波功能,将高速宽带数据流变为低速基带数据流,以便DSP实时处理.该文提出了数字下变频的实现方案,给出了相应的算法和算法分析,特别提出了CIC滤波器与内插的二阶多项式滤波器(ISOP-InterpolatedSecond-OrderPolymomial)组合进行有效抽取滤波的设计.在AM信号的解调算法中给出了内部自动增益控制AGC的算法.对于DSB-SC信号解调特别地提出了抑制载波双边带调幅信号(DSB-SC)载波同步的算法,根据下变频后的基带信号估计出收发载波的频差,并经卡尔曼滤波后结合数字锁相环以达到载波捕获范围宽、跟踪速度快、环路噪声小的特点,仿真表明,该算法有效可行.该论文中还给出了适于DSP器件的FM信号解调算法.论文中的算法不仅用MATLAB5.3进行了仿真并给出全部仿真结果,而且在TMS320C32EVM板中进行了仿真验证.
本文链接：https://www.doczj.com/doc/8712470841.html,/Periodical_kjzx200923009.aspx
下载时间：2010年5月17日
实验6 数字滤波器的网络结构 一、实验目的： 1、加深对数字滤波器分类与结构的了解。 2、明确数字滤波器的基本结构及其相互间的转换方法。 3、掌握用MA TLAB 语言进行数字滤波器结构间相互转换的子函数及程序编写方法。 二、实验原理： 1、数字滤波器的分类 离散LSI 系统对信号的响应过程实际上就是对信号进行滤波的过程。因此，离散LSI 系统又称为数字滤波器。 数字滤波器从滤波功能上可以分为低通、高通、带通、带阻以及全通滤波器；根据单位脉冲响应的特性，又可以分为有限长单位脉冲响应滤波器（FIR ）和无限长单位脉冲响应滤波器（IIR ）。 一个离散LSI 系统可以用系统函数来表示： M -m -1-2-m m m=0 012m N -1-2-k -k 12k k k=1 b z b +b z +b z ++b z Y(z)b(z)H(z)=== =X(z)a(z) 1+a z +a z ++a z 1+a z ∑∑ 也可以用差分方程来表示： N M k m k=1 m=0 y(n)+a y(n-k)=b x(n-m)∑∑ 以上两个公式中，当a k 至少有一个不为0时，则在有限Z 平面上存在极点，表达的是以一个IIR 数字滤波器；当a k 全都为0时，系统不存在极点，表达的是一个FIR 数字滤波器。FIR 数字滤波器可以看成是IIR 数字滤波器的a k 全都为0时的一个特例。 IIR 数字滤波器的基本结构分为直接Ⅰ型、直接Ⅱ型、直接Ⅲ型、级联型和并联型。 FIR 数字滤波器的基本结构分为横截型（又称直接型或卷积型）、级联型、线性相位型及频率采样型等。本实验对线性相位型及频率采样型不做讨论，见实验10、12。 另外，滤波器的一种新型结构——格型结构也逐步投入应用，有全零点FIR 系统格型结构、全极点IIR 系统格型结构以及全零极点IIR 系统格型结构。 2、IIR 数字滤波器的基本结构与实现 （1）直接型与级联型、并联型的转换 例6-1 已知一个系统的传递函数为 -1-2-3 -1-2-3 8-4z +11z -2z H(z)=1-1.25z +0.75z -0.125z 将其从直接型（其信号流图如图6-1所示）转换为级联型和并联型。
现代数字信号处理仿真作业 1.仿真题3.17 仿真结果及图形： 图 1 基于FFT的自相关函数计算
图 3 周期图法和BT 法估计信号的功率谱 图 2 基于式3.1.2的自相关函数的计算
图 4 利用LD迭代对16阶AR模型的功率谱估计16阶AR模型的系数为： a1=-0.402637623107952-0.919787323662670i; a2=-0.013530139693503+0.024214641171318i; a3=-0.074241889634714-0.088834852915013i; a4=0.027881022353997-0.040734794506749i; a5=0.042128517350786+0.068932699075038i; a6=-0.0042799971761507 + 0.028686095385146i; a7=-0.048427890183189 - 0.019713457742372i; a8=0.0028768633718672 - 0.047990801912420i a9=0.023971346213842+ 0.046436389191530i; a10=0.026025963987732 + 0.046882756497113i; a11= -0.033929397784767 - 0.0053437929619510i; a12=0.0082735406293574 - 0.016133618316269i; a13=0.031893903622978 - 0.013709547028453i ; a14=0.0099274520678052 + 0.022233240051564i; a15=-0.0064643069578642 + 0.014130696335881i; a16=-0.061704614407581- 0.077423818476583i. 仿真程序（3_17）: clear all clc %% 产生噪声序列 N=32; %基于FFT的样本长度
数字信号处理作业
DFT 习题 1. 如果)(~n x 是一个周期为N 的周期序列，那么它也是周期为N 2的周期序列。把)(~ n x 看作周期为N 的周期序列，令)(~1k X 表示)(~n x 的离散傅里叶级数之系数，再把)(~ n x 看作周期为N 2的周期序列，再令)(~2k X 表示)(~n x 的离散傅里叶级数之系数。当然，)(~1k X 是周期性的，周期为N ，而)(~2k X 也是周期性的，周期为N 2。试利用)(~1k X 确定)(~2k X 。（76-4）
2. 研究两个周期序列)(~n x 和)(~n y 。)(~n x 具有周期N ,而)(~ n y 具有周期M 。序列)(~n w 定义为)()()(~ ~~n y n x n w +=。 a. 证明)(~n w 是周期性的，周期为MN 。 b. 由于)(~n x 的周期为N ，其离散傅里叶级数之系数)(~k X 的周期也是N 。类似地， 由于)(~n y 的周期为M ，其离散傅里叶级数之系数)(~k Y 的周期也是M 。)(~n w 的离散傅里叶级数之系数)(~k W 的周期为MN 。试利用)(~k X 和)(~k Y 求)(~k W 。（76-5）
3. 计算下列各有限长度序列DFT （假设长度为N ）: a. )()(n n x δ= b ．N n n n n x <<-=000) ()(δ c ．10)(-≤≤=N n a n x n （78-7） 4. 欲作频谱分析的模拟数据以10千赫速率被取样，且计算了1024个取样的离散傅里叶变换。试求频谱取样之间的频率间隔，并证明你的回答。（79 -10）
实验5 抽样定理 一、实验目的： 1、了解用MA TLAB 语言进行时域、频域抽样及信号重建的方法。 2、进一步加深对时域、频域抽样定理的基本原理的理解。 3、观察信号抽样与恢复的图形，掌握采样频率的确定方法和插公式的编程方法。 二、实验原理： 1、时域抽样与信号的重建 （1）对连续信号进行采样 例5-1 已知一个连续时间信号sin sin(),1Hz 3 ππ=0001f(t)=(2f t)+6f t f ，取最高有限带宽频率f m =5f 0，分别显示原连续时间信号波形和F s >2f m 、F s =2f m 、F s <2f m 三情况下抽样信号的波形。 程序清单如下： %分别取Fs=fm ，Fs=2fm ，Fs=3fm 来研究问题 dt=0.1; f0=1; T0=1/f0; m=5*f0; Tm=1/fm; t=-2:dt:2; f=sin(2*pi*f0*t)+1/3*sin(6*pi*f0*t); subplot(4,1,1); plot(t,f); axis([min(t),max(t),1.1*min(f),1.1*max(f)]); title('原连续信号和抽样信号'); for i=1:3; fs=i*fm;Ts=1/fs; n=-2:Ts:2; f=sin(2*pi*f0*n)+1/3*sin(6*pi*f0*n); subplot(4,1,i+1);stem(n,f,'filled'); axis([min(n),max(n),1.1*min(f),1.1*max(f)]); end 程序运行结果如图5-1所示：
原连续信号和抽样信号 图5-1 （2）连续信号和抽样信号的频谱 由理论分析可知，信号的频谱图可以很直观地反映出抽样信号能否恢复原模拟信号。因此，我们对上述三种情况下的时域信号求幅度谱，来进一步分析和验证时域抽样定理。 例5-2编程求解例5-1中连续信号及其三种抽样频率（F s>2f m、F s=2f m、F s<2f m）下的抽样信号的幅度谱。 程序清单如下： dt=0.1;f0=1;T0=1/f0;fm=5*f0;Tm=1/fm; t=-2:dt:2;N=length(t); f=sin(2*pi*f0*t)+1/3*sin(6*pi*f0*t); wm=2*pi*fm;k=0:N-1;w1=k*wm/N; F1=f*exp(-j*t'*w1)*dt;subplot(4,1,1);plot(w1/(2*pi),abs(F1)); axis([0,max(4*fm),1.1*min(abs(F1)),1.1*max(abs(F1))]); for i=1:3; if i<=2 c=0;else c=1;end fs=(i+c)*fm;Ts=1/fs; n=-2:Ts:2;N=length(n); f=sin(2*pi*f0*n)+1/3*sin(6*pi*f0*n); wm=2*pi*fs;k=0:N-1; w=k*wm/N;F=f*exp(-j*n'*w)*Ts; subplot(4,1,i+1);plot(w/(2*pi),abs(F)); axis([0,max(4*fm),1.1*min(abs(F)),1.1*max(abs(F))]); end 程序运行结果如图5-2所示。 由图可见，当满足F s≥2f m条件时，抽样信号的频谱没有混叠现象；当不满足F s≥2f m 条件时，抽样信号的频谱发生了混叠，即图5-2的第二行F s＜2f m的频谱图，，在f m=5f0的围，频谱出现了镜像对称的部分。
实验名称：多采样率数字信号处理 一．实验目的：1. 掌握信号抽取和插值的基本原理和实现； 2．掌握信号的有理数倍率转换。 二．实验原理： 多采样率数字信号处理共分为3方面的问题：信号的整数倍抽取、信号的整数倍插值和信号的有理数倍速率转换。 Matlab 信号处理工具箱提供了抽取函数decimate 用于信号整数倍抽取，其调用格式为： y=decimate(x,M) y=decimate(x,M,n) y=decimate(x,M,’fir’) y=decimate(x,M,n,’fir’) 其中y=decimate(x,M)将信号x 的采样率降低为原来的 M 1，抽取前缺省地采用8阶Chebyshev Ⅰ型低通滤波器压缩频带。 y=decimate(x,M,n)指定所采用Chebyshev Ⅰ型低通滤波器的阶数，通常13 n 。 y=decimate(x,M,’fir’)指定用FIR 滤波器来压缩频带。 y=decimate(x,M,n,’fir’) 指定所用FIR 滤波器的阶数。 Matlab 信号处理工具箱提供了插值函数interp 用于信号整数倍插值，其调用格式为： y=interp(x,L) y=interp(x,L,n,alpha) [y,b]=interp(x,L,n,alpha) 其中y=interp(x,L)将信号的采样率提高到原来的L 倍。 y=interp(x,L,n,alpha)指定反混叠滤波器的长度n 和截止频率alpha ，缺省值为4和0.5。 [y,b]=interp(x,L,n,alpha)在插值的同时，返回反混叠滤波器的系数向量。 信号的有理数倍速率转换是使信号的采样率经由一个有理因子M L 来改变，可以通过插值和抽取的级联来实现。Matlab 信号处理工具箱提供了重采样函数resample 用于有理倍数速率转换，其调用格式为： y=resample(x,L,M);
数字信号处理上机作业 学院：电子工程学院 班级：021215 组员：
实验一：信号、系统及系统响应 1、实验目的 (1) 熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定理的理解。 (2) 熟悉时域离散系统的时域特性。 (3) 利用卷积方法观察分析系统的时域特性。 (4) 掌握序列傅里叶变换的计算机实现方法，利用序列的傅里叶变换对连续信号、离散信号及系统响应进行频域分析。 2、实验原理与方法 (1) 时域采样。 (2) LTI系统的输入输出关系。 3、实验内容及步骤 (1) 认真复习采样理论、离散信号与系统、线性卷积、序列的傅里叶变换及性质等有关内容，阅读本实验原理与方法。 (2) 编制实验用主程序及相应子程序。 ①信号产生子程序，用于产生实验中要用到的下列信号序列： a. xa(t)=A*e^-at *sin(Ω0t)u(t) b. 单位脉冲序列：xb(n)=δ(n) c. 矩形序列： xc(n)=RN(n), N=10 ②系统单位脉冲响应序列产生子程序。本实验要用到两种FIR系统。 a. ha(n)=R10(n); b. hb(n)=δ(n)+2.5δ(n-1)+2.5δ(n-2)+δ(n-3) ③有限长序列线性卷积子程序 用于完成两个给定长度的序列的卷积。可以直接调用MATLAB语言中的卷积函数conv。 conv 用于两个有限长度序列的卷积，它假定两个序列都从n=0 开始。调用格式如下： y=conv (x, h) 4、实验结果分析 ①分析采样序列的特性。 a. 取采样频率fs=1 kHz,，即T=1 ms。 b. 改变采样频率，fs=300 Hz，观察|X(e^jω)|的变化，并做记录(打印曲线)；进一步降低采样频率，fs=200 Hz，观察频谱混叠是否明显存在，说明原因，并记录(打印)这时的|X(e^j ω)|曲线。 程序代码如下： close all;clear all;clc; A=50; a=50*sqrt(2)*pi; m=50*sqrt(2)*pi; fs1=1000; fs2=300; fs3=200; T1=1/fs1; T2=1/fs2; T3=1/fs3; N=100;
Chapter 2 Solutions 2.1 最小采样频率为两倍的信号最大频率，即44.1kHz 。 2.2 (a)、由ω = 2πf = 20 rad/sec ，信号的频率为f = 3.18 Hz 。信号的奈奎斯特采样频率为6.37 Hz 。 (b)、3 5000π=ω，所以f = 833.3 Hz ，奈奎斯特采样频率为1666.7 Hz 。 (c)、7 3000π=ω，所以f = 214.3 Hz ，奈奎斯特采样频率为428.6 Hz 。 2.3 (a) 1258000 1f 1T S S ===μs (b)、最大还原频率为采样频率的一半，即4000kHz 。 2.4 ω = 4000 rad/sec ，所以f = 4000/(2π) = 2000/π Hz ，周期T = π/2000 sec 。因此，5个周期为5π/2000 = π/400 sec 。对于这个信号，奈奎斯特采样频率为2(2000/π) = 4000/π Hz 。所以采样频率为f S = 4(4000/π) = 16000/π Hz 。因此5个周期收集的采样点为(16000/π samples/sec )(π/400 sec) = 40。 2.5 ω = 2500π rad/sec ，所以f = 2500π/(2π) = 1250 Hz ，T = 1/1250 sec 。因此，5个周期为5/1250 sec 。对于这个信号，奈奎斯特采样频率为2(1250) = 2500 Hz ，所以采样频率为f S = 7/8(2500) = 2187.5 Hz 。采样点数为(2187.5 点/sec)(5/1250 sec) = 8.75。这意味着在模拟信号的五个周期内只有8个点被采样。事实上，对于这个信号来说，在整数的模拟周期中，是不可能采到整数个点的。 2.6 2.7 信号搬移发生在kf S ± f 处，换句话说，频谱搬移发生在每个采样频率的整数倍 0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 频率/kHz
数字信号处理作业 哈尔滨工业大学 2006.10
DFT 习题 1. 如果)(~n x 是一个周期为N 的周期序列，那么它也是周期为N 2的周期序列。把)(~ n x 看作周期为N 的周期序列，令)(~ 1k X 表示)(~n x 的离散傅里叶级数之系数，再把)(~ n x 看作周期为N 2的周期序列，再令)(~ 2k X 表示)(~n x 的离散傅里叶级数之系数。当然，)(~ 1k X 是周期性的，周期为N ，而)(~ 2k X 也是周期性的，周期为N 2。试利用)(~ 1k X 确定)(~ 2k X 。（76-4）
2. 研究两个周期序列)(~ n x 和)(~ n y 。)(~ n x 具有周期N ,而)(~ n y 具有周期M 。序列 )(~n w 定义为)()()(~ ~~n y n x n w +=。 a. 证明)(~ n w 是周期性的，周期为MN 。 b. 由于)(~n x 的周期为N ，其离散傅里叶级数之系数)(~ k X 的周期也是N 。类似地， 由于)(~n y 的周期为M ，其离散傅里叶级数之系数)(~k Y 的周期也是M 。)(~ n w 的离散傅里叶级数之系数)(~ k W 的周期为MN 。试利用)(~ k X 和)(~ k Y 求)(~ k W 。（76-5）
3. 计算下列各有限长度序列DFT （假设长度为N ）: a. )()(n n x δ= b ．N n n n n x <<-=000)()(δ c ．10)(-≤≤=N n a n x n （78-7） 4. 欲作频谱分析的模拟数据以10千赫速率被取样，且计算了1024个取样的离散傅里叶变换。试求频谱取样之间的频率间隔，并证明你的回答。（79 -10）
仿真作业 姓名：李亮 学号：S130101083
4.17程序 clc; clear; for i=1:500 sigma_v1=0.27; b(1)=-0.8458; b(2)=0.9458; a(1)=-(b(1)+b(2)); a(2)=b(1)*b(2); datlen=500; rand('state',sum(100*clock)); s=sqrt(sigma_v1)*randn(datlen,1); x=filter(1,[1,a],s); %% sigma_v2=0.1; u=x+sqrt(sigma_v2)*randn(datlen,1); d=filter(1,[1,-b(1)],s); %% w0=[1;0]; w=w0; M=length(w0); N=length(u); mu=0.005; for n=M:N ui=u(n:-1:n-M+1); y(n)=w'*ui; e(n)=d(n)-y(n); w=w+mu.*conj(e(n)).*ui; w1(n)=w(1); w2(n)=w(2); ee(:,i)=mean(e.^2,2); end end ep=mean(ee'); plot（ep）; xlabel('迭代次数');ylabel('MSE');title('学习曲线'); plot(w1); hold; plot(w2); 仿真结果：
步长0.015仿真结果 0.10.20.30.4 0.50.60.7迭代次数 M S E 学习曲线
步长0.025仿真结果
步长0.005仿真结果 4.18 程序 data_len = 512; %样本序列的长度 trials = 100; %随机试验的次数 A=zeros(data_len,2);EA=zeros(data_len,1); B=zeros(data_len,2);EB=zeros(data_len,1); for m = 1: trials a1 = -0.975; a2 = 0.95; sigma_v_2 =0.0731; v = sqrt(sigma_v_2) * randn(data_len, 1, trials);%产生v(n) u0 = [0 0]; num = 1; den = [1 a1 a2]; Zi = filtic(num, den, u0); %滤波器的初始条件 u = filter(num, den, v, Zi); %产生样本序列u(n) %（2）用LMS滤波器来估计w1和w2 mu1 = 0.05; mu2 = 0.005; w1 = zeros(2, data_len);
数字信号处理作业-答案
数字信号处理作业
DFT 习题 1. 如果)(~ n x 是一个周期为N 的周期序列，那么它也是周期为N 2的周期序列。把)(~ n x 看作周期为N 的周期序列，令)(~ 1 k X 表示)(~ n x 的离散傅里叶级数之系数，再把)(~ n x 看作周期为N 2的周期序列，再令)(~2 k X 表示)(~ n x 的离散傅里叶级数之系数。当然，)(~ 1 k X 是周期性的，周期为N ，而)(~ 2 k X 也是周期性的，周期为N 2。试利用)(~ 1k X 确定)(~ 2 k X 。（76-4）
2. 研究两个周期序列)(~ n x 和)(~ n y 。)(~ n x 具有周期N ,而)(~ n y 具有周期M 。序列)(~ n w 定义为)()()(~~ ~ n y n x n w +=。 a. 证明)(~ n w 是周期性的，周期为MN 。 b. 由于)(~ n x 的周期为N ，其离散傅里叶级数之系数)(~k X 的周期也是N 。类似地，由于)(~ n y 的周期为M ，其离散傅里叶级数之系数)(~ k Y 的周期也是M 。)(~n w 的离散傅里叶级数之系数)(~ k W 的周期为MN 。试利用)(~k X 和)(~k Y 求)(~ k W 。（76-5）
3. 计算下列各有限长度序列DFT （假设长度为N ）: a. )()(n n x δ= b ．N n n n n x <<-=0 0)()(δ c ．10)(-≤≤=N n a n x n （78-7） 4. 欲作频谱分析的模拟数据以10千赫速率被取样，且计算了1024个取样的离散傅里叶变换。试求频谱取样之间的频率间隔，并证明你的回答。（79 -10）
《数字信号处理Ⅰ》作业 姓名： 学号： 学院： 2012 年春季学期
第一章 时域离散信号和时域离散系统 月 日 一 、判断： 1、数字信号处理和模拟信号处理在方法上是一样的。（ ） 2、如果信号的取值和自变量都离散，则称其为模拟信号。（ ） 3、如果信号的取值和自变量都离散，则称其为数字信号。（ ） 4、时域离散信号就是数字信号。（ ） 5、正弦序列都是周期的。（ ） 6、序列)n (h )n (x 和的长度分别为N 和M 时，则)n (h )n (x *的长度为N+M 。（ ） 7、如果离散系统的单位取样响应绝对可和，则该系统稳定。（ ） 8、若满足采样定理，则理想采样信号的频谱是原模拟信号频谱以s Ω（采样频率）为周期进行周期延拓的结果。（ ） 9、序列)n (h )n (x 和的元素个数分别为21n n 和，则)n (h )n (x *有（1n n 21-+）个元素。（ ） 二、选择 1、R N (n)和u(n)的关系为（ ）： A. R N (n)=u(n)-u(n-N) B. R N (n)=u(n)+u(n-N) C. R N (n)=u(n)-u(n-N-1) D. R N (n)=u(n)-u(n-N+1) 2、若f(n)和h(n)的长度为别为N 、M ，则f(n)*h(n)的长度为 （ ）： A.N+M B.N+M-1 C.N-M D.N-M+1 3、若模拟信号的频率范围为[0，1kHz],对其采样，则奈奎斯特速率为（ ）： A.4kHz B. 3kHz C.2kHz D.1kHz 4、LTIS 的零状态响应等于激励信号和单位序列响应的（ ）： A.相乘 B. 相加 C.相减 D.卷积 5、线性系统需满足的条件是（ ）： A.因果性 B.稳定性 C.齐次性和叠加性 D.时不变性 6、系统y(n)=f(n)+2f(n-1)(初始状态为0)是（ ）： A. 线性时不变系统 B. 非线性时不变系统 C. 线性时变系统 D. 非线性时变系统
数字信号处理大作业 班级：021231 学号： 姓名： 指导老师：吕雁
一写出奈奎斯特采样率和和信号稀疏采样的学习报告和体会 1、采样定理 在进行A/D信号的转换过程中，当采样频率fs.max大于信号中最高频 率fmax的2倍时(fs.max>2fmax)，采样之后的数字信号完整地保留了原始信号中的信息，一般实际应用中保证采样频率为信号最高频率的5～10倍；采样定 理又称奈奎斯特定理。 (1)在时域 频带为F的连续信号 f(t)可用一系列离散的采样值f(t1),f(t1±Δt)，f(t1±2Δt)，...来表示,只要这些采样点的时间间隔Δt≤1/2F，便可根据各 采样值完全恢复原始信号。 (2)在频域 当时间信号函数f(t)的最高频率分量为fmax时,f(t)的值可由一系列 采样间隔小于或等于1/2fo的采样值来确定,即采样点的重复频率fs ≥2fmax。 2、奈奎斯特采样频率 (1)概述 奈奎斯特采样定理：要使连续信号采样后能够不失真还原，采样频率必须 大于信号最高频率的两倍（即奈奎斯特频率）。 奈奎斯特频率（Nyquist frequency）是离散信号系统采样频率的一半，因哈里·奈奎斯特（Harry Nyquist）或奈奎斯特－香农采样定理得名。采样定理指出，只要离散系统的奈奎斯特频率高于被采样信号的最高频率或带宽，就可 以真实的还原被测信号。反之，会因为频谱混叠而不能真实还原被测信号。 采样定理指出，只要离散系统的奈奎斯特频率高于采样信号的最高频率或 带宽，就可以避免混叠现象。从理论上说，即使奈奎斯特频率恰好大于信号带宽，也足以通过信号的采样重建原信号。但是，重建信号的过程需要以一个低 通滤波器或者带通滤波器将在奈奎斯特频率之上的高频分量全部滤除，同时还 要保证原信号中频率在奈奎斯特频率以下的分量不发生畸变，而这是不可能实 现的。在实际应用中，为了保证抗混叠滤波器的性能，接近奈奎斯特频率的分 量在采样和信号重建的过程中可能会发生畸变。因此信号带宽通常会略小于奈 奎斯特频率，具体的情况要看所使用的滤波器的性能。需要注意的是，奈奎斯 特频率必须严格大于信号包含的最高频率。如果信号中包含的最高频率恰好为
数字信号处理第三章作业 1．(第三章习题3)在图P3-2中表示了两个周期都为6的周期性序列，确定这个两个序列的周期卷积的结果3()x n ，并画出草图。 2．(第三章习题5)如果()x n 是一个具有周期为N 的周期性序列，它也是具有周期为2N 的周期性序列。令~1()X k 表示当()x n 看做是具有周期为N 的周期性序列的DFS 系数。而~2()X k 表示当()x n 看作是具有周期为2N 的周期性序列的DFS 系数。当然~1()X k 是具有周期为N 的周期性序列，而~2()X k 是具有周期为2N 的周期性序列，试根据~1()X k 确定~2()X k 。 3．(第三章习题6) （a ）试证明下面列出的周期性序列离散傅里叶级数的对称特性。在证明中，可以利用离散傅里叶级数的定义及任何前面的性质，例如在证明性质③时可以利用性质①和②。 序列 离散傅里叶级数 ① *()x n ~*()X k - ②*()x n - ~*()X k ③Re ()x n ???? ~ e ()X k ④Im ()j x n ???? ~()o X k
（b ）根据已在（a ）部分证明的性质，证明对于实数周期序列()x n ，离散傅里叶级数的下列对称性质成立。 ①~~Re ()Re ()X k X k ????=-???????? ②~~Im ()Im ()X k X k ????=--???????? ③~~()()X k X k =- ④~~arg ()arg ()X k X k ????=--???????? 4．(第三章习题7)求下列序列的DFT (a) {}11 1-，，，-1 (b) {}1 j 1j -，，，- (c) ()cn 0n 1x n N =≤≤-， (d) 2n ()sin 0n 1x n N N π??=≤≤- ??? ， 5．(第三章习题8)计算下列各有限长序列的离散傅立叶变换（假设长度为N ） 1 0)()(0) ()()() ()()(00-≤≤=<<-==N n a n x c N n n n n x b n n x a n δδ 6．(第三章习题9)在图P3-4中表示了一有限长序列)(n x ，画出序列)(1n x 和)(2n x 的草图。（注意：)(1n x 是)(n x 圆周移位两个点） )())(()() ())2(()(442441n R n x n x n R n x n x -=-=
前言 《数字信号处理》是信息电子，通信工程等本科专业及其他相近专业的一门专业必修课。通过本课程的学习，学生应掌握以下基本概念、理论和方法：采样定理、离散序列的变换、离散信号的频谱分析；离散系统的传递函数、频率响应、离散系统的基本分析方法；数字滤波器的设计理论、滤波器的软件实现；离散傅立叶变换理论、快速傅立叶变换方法；有限字长效应。 为了使学生更好地理解和深刻地把握这些知识，并在此基础上，训练和培养学生掌握离散系统的基本概念和分析方法，数字滤波器的设计和实现，以及如何利用快速傅立叶变换等DSP技术对数字信号进行分析、滤波等处理，设置了以下三个实验： （1）离散时间序列卷积和MATLAB实现； 内容：使用任意的编程语言编制一个程序，实现两个任意有限序列的卷积和。 目的：理解线性非移变系统I/O关系和实现 要求：掌握使用计算机实现数字系统的方法 （2）FFT算法的MATLAB实现； 内容：使用MATLAB编程语言编制一个程序，实现任意有限序列的FFT。 目的：理解FFT算法的意义和实现 要求：掌握使用计算机实现FFT算法的方法 （3）数字滤波器的设计； 内容：使用MATLAB编程语言编制一个程序，实现FIR或IIR滤波器的设计目的：理解数字滤波器的设计技术 要求：掌握使用计算机进行数字滤波器设计的方法 （4）窗函数设计FIR滤波器； 内容：使用MATLAB编程语言编制一个程序，实现FIR或IIR滤波器的设计目的：理解数字滤波器的设计技术 要求：掌握使用计算机进行数字滤波器设计的方法 该实验指导书是参照该课程的教学大纲而编制的，适合于信息电子工程、通信工程等本科专业及其他相近专业。
1、什么是DSP？简述DSPs的特点？简述DSPs与MCU、FPGA、ARM的区别？学习DSP开发需要哪些知识？学习DSP开发需要构建什么开发环境？(15分) 答：（1）DSP是Digital Signal Processing(数字信号处理的理论和方法）的缩写，同时也是Digital Signal Processor（数字信号处理的可编程微处理器）的缩写。通常流过器件的电压、电流信号都是时间上连续的模拟信号，可以通过A/D器件对连续的模拟信号进行采样，转换成时间上离散的脉冲信号，然后对这些脉冲信号量化、编码，转化成由0和1构成的二进制编码，也就是常说的数字信号。DSP能够对这些数字信号进行变换、滤波等处理，还可以进行各种各样复杂的运算，来实现预期的目标。 （2）DSP既然是特别适合于数学信号处理运算的微处理器，那么根据数字信号处理的要求，DSP芯片一般具有下面所述的主要特点：1）程序空间和数据空间分开，CPU可以同时访问指令和数据； 2）在一个指令周期内可以完成一次乘法和一次加法运算； 3）片内具有快速RAM，通常可以通过独立的数据总线在程序空间和数据空间同时访问； 4）具有低开销和无开销循环及跳转的硬件支持； 5）具有快速的中断处理和硬件I/O支持； 6）可以并行执行多个操作； 7）支持流水线操作，使得取址、译码和执行等操作可以重复执行。（3）DSP采用的是哈佛结构，数据空间和存储空间是分开的，通过
独立的数据总线在数据空间和程序空间同时访问。而MCU采用的是冯·诺依曼结构，数据空间和存储空间共用一个存储器空间，通过一组总线（地址总线和数据总线）连接到CPU）。很显然，在运算处理能力上，MCU不如DSP；但是MCU价格便宜，在对性能要求不是很高的情况下，还是很具有优势的。 ARM是Advanced RISC（精简指令集）Machines的缩写是面向低运算市场的RISC微处理器。ARM具有比较强的事务管理功能，适合用来跑跑界面、操作系统等，其优势主要体现在控制方面，像手持设备90%左右的市场份额均被其占有。而DSP的优势是其强大的数据处理能力和较高的运算速度，例如加密/解密、调制/解调等。 FPGA是Field Programmable Gate Array（现场可编程门阵列）的缩写，它是在PAL、GAL、PLD等可编程器件的基础上进一步发展的产物，是专用集成电路中集成度最高的一种。FPGA采用了逻辑单元阵列LCA（Logical Cell Array）的概念，内部包括了可配置逻辑模块CLB、输入/输出模块IOB、内部连线三个部分。用户可以对FPGA内部的逻辑模块和I/O模块进行重置配置，已实现用户自己的逻辑。它还具有静态可重复编程和动态在系统重构的特性，使得硬件的功能可以像软件一样通过编程来修改。使用FPGA来开发数字电路，可以大大缩短设计时间，减少PCB面积，提高系统的可靠性；同时FPGA可以用VHDL或Verilog HDL来编程，灵活性强。由于FPGA能够进行编程、除错、再编程和重复操作，因此可以充分地进行设计开发和验证。当电路有少量改动时，更能显示出FPGA的优势，其现场编程能力可
现代数字信号处理仿真作业
第三章仿真作业3.17 （1）代码 clear; N=32; m=[-N+1:N-1]; noise=(randn(1,N)+j*randn(1,N))/sqrt(2); f1=0.15; f2=0.17; f3=0.26; SNR1=30; SNR2=30; SNR3=27; A1=10^(SNR1/20); A2=10^(SNR2/20); A3=10^(SNR3/20); signal1=A1*exp(j*2*pi*f1*(0:N-1)); signal2=A2*exp(j*2*pi*f2*(0:N-1)); signal3=A3*exp(j*2*pi*f3*(0:N-1)); un=signal1+signal2+signal3+noise; uk=fft(un,2*N); sk=(1/N) *abs(uk).^2; r0=ifft(sk); r1=[r0(N+2:2*N),r0(1:N)]; r=xcorr(un,N-1,'biased'); figure subplot(2,2,1) stem(m,real(r1)); xlabel('m'); ylabel('FFT估计r1实部'); subplot(2,2,2) stem(m,imag(r1)); xlabel('m'); ylabel('FFT估计r1虚部'); subplot(2,2,3) stem(m,real(r)); xlabel('m'); ylabel('平均估计r实部'); subplot(2,2,4) stem(m,imag(r)); xlabel('m'); ylabel('平均估计r虚部'); 仿真结果
万方数据
多采样率数字信号处理及其MATLAB仿真 作者：黄硕， 魏亚楠， 安永丽 作者单位：唐山钢铁股份有限公司,唐山,063016 刊名： 科技资讯 英文刊名：SCIENCE & TECHNOLOGY INFORMATION 年，卷(期)：2009，(23) 引用次数：0次 参考文献(3条) 1.杨小牛.楼才义.徐建良软件无线电原理与应用 2005 2.李忠琦.凌翔.胡剑浩软件无线电架构研究[期刊论文]-电信科学 2007(7) 3.尹健华试论软件无线电技术及其应用[期刊论文]-企业技术开发(学术版) 2007(8) 相似文献(10条) 1.学位论文赵启敏中频采样技术的分析与研究2004 该课题结合数字软件化雷达的研制,研究了数字软件化雷达中频采样技术的实现以及对雷达主要技术指标的影响.该论文针对传统模拟相参正交采样技术存在的不足,论述了基于A/D变换和数字下变频的中频采样方法,并在此基础上设计了中频采样数据采集卡,并对该硬件进行了调试和试验,试验结果证明,中频采样技术比传统模拟相参正交采样技术更具优势,较好的解决了传统模拟相参正交采样中存在的幅相误差问题,以及该采集卡具有小的孔径抖动,可以满足中频采样的要求.该论文在中频采样技术中首先研究了数据采集技术对雷达性能的影响,接着根据目前数字下变频器件自身的限制不能适应高速数据流的问题,详细研究了利用欠采样技术的镜频加数字下变频实现解调的方法,以及一种利用多速率信号处理技术将抽取和滤波提前的数字下变频的高效结构,通过仿真证明此两种方法都能较好的解决硬件本身限制与高速数据流不匹配的问题,并通过分析得出此数字下变频的高效结构的运算量大大低于传统数字下变频的运算量.此外该论文还着重讨论了孔径抖动对雷达各项性能的影响. 2.期刊论文张明珊.孟利民.ZHANG Ming-shan.MENG Li-min基于频域采样技术的软件无线电接收机-浙江工业大学学报2005,33(1) 目前软件无线电面临的一个难题是如何对高工作频带内的射频信号进行直接模/数转换.利用频域采样技术提出了一种接收信号进行处理的方法,并用数学理论证明了它的可行性,最后还给出了软件无线电接收机模型.其关键思想是提取接收信号的频域成份,然后在频域中对信号进行处理.这种方法大大降低了A/D转换器的要求,从而使得实现软件无线电接收机成为可能,对当前微电子工艺下的软件无线电系统设计带来很大的理论意义和实用价值,而且克服了传统Rake接收机的一些缺点,特别适合于多径丰富的无线环境. 3.学位论文杨清海软件无线电的功能实现2001 1992年，JeoMitola提出了软件无线电的概念，很快引起了国际通信界的关注。软件无线电结构的关键是在尽可能靠近天线的地方使用宽带A／D和 D／A变换器，将尽可能多的无线电功能用软件来定义，从而实现电台在各种网络中的通用性及电台功能升级换代的连续性，软件无线电已成为无线通信的一个主要发展方向。特别是近年来，软件无线电已经不再仅仅局限于军事方面，在GSMMOU会议中，软件无线电被描述成GSM继续发展进步的基础，甚至被称为第三代（3G）全球移动通信实现的技术基础。本文主要探讨软件无线电思想在接收机设计中的应用，论证了系统硬件实现方案和软件实现方案。重点讨论了用到的信号采样技术和数字信号处理技术，包括多速率信号处理、FIR滤波器的多相结构、低通滤波、免混频正交解调和信号的带通采样技术。最后优化了解调算法，利用我们的试验平台实现了AM、FM、SSB和ASK、FSK、PSK信号解调。 4.学位论文洪亮高速并行交替采样ADC系统的研究与实现2009 模数转换器(ADC)是数字信号处理系统的关键组成部分，广泛应用于通信、雷达、测试仪器等领域。随着超宽带雷达技术研究的深入和软件无线电技术的发展，对ADC的速度和精度的要求越来越高，ADC已经成为现代信号处理的瓶颈。在给定的工艺下，ADC工作的最大采样速率受限于它的分辨率，单片ADC芯片很难同时满足高速高精度的要求，而并行交替采样ADC(TIADC)结构是突破这一瓶颈的有效方法之一。 这种方法在前端利用M片采样率为fs/M的ADC并行交替采样，在后端进行拼接使得整个系统的采样率达到fs。然而受到制造工艺的局限，通道失配误差如偏置误差、增益误差、时间偏差和带宽失配误差的存在，将严重降低系统的信纳比(SINAD)和无杂散动态范围(SFDR)。 本论文主要包括三方面的工作。首先，深入研究了并行交替采样技术，对TIADC结构的通道失配误差进行了全面的分析，特别是对带宽失配误差进行了建模分析，给出了四种通道失配误差联合作用于信号的信号频谱，以及系统设计时误差的容忍范围。 其次，通过合理的近似，提出了通道失配误差的测量算法和联合校正算法，其中关键的是时间偏差和带宽失配误差的联合估算与校正，它是在周期非均匀采样信号完美重构基础上提出来的，并通过仿真验证了算法的有效性。 最后，设计了一个基于并行交替采样技术的12bit420MSPS的高速数据采集系统，该系统由两片12bit210MSPS的AD9430组成。其中，结合系统设计进行的信号完整性分析对高速电路的设计具有一定的指导意义。 5.期刊论文王宏.刘丽.宋晓峰.WANG Hong.LIU Li.SONG Xiaofeng基于频域采样技术的软件无线电接收机-现代电子技术2006,29(23) 目前软件无线电面临的一个难题是如何对高工作频带内的射频信号进行直接模/数转换.利用频域采样技术提出了一种接收信号进行处理的方法,并用数学理论证明了他的可行性,最后还给出了软件无线电接收机模型.其关键思想是提取接收信号的频域成份,然后在频域中对信号进行处理.这种方法大大降低了A/D转换器的要求,对当前微电子工艺下的软件无线电系统设计有很大的理论意义和实用价值,而且克服了传统Rake接收机的一些缺点,特别适合于多径丰富的无线环境. 6.学位论文李裕多信道软件无线电接收机实现技术研究2003 软件无线电的基本思想是将宽带A/D及D/A尽可能靠近天线,将无线电台的各种功能在一个开放性、模块化的通用硬件平台上尽可能多的用软件来实现.软件无线电已成为移动通信中的关键技术之一.本文主要研究了软件无线电接收机中的相关理论及实现方案并进行了相应的系统仿真.本文首先深入讨论了软件无线电接收机的基本理论:采样技术、多速率信号处理和调制解调算法.在此基础上研究了下变频技术和带通采样技术在并行多信道接收机中的应用,提出了利用CIC,HBF和FIR级联设计下变频器的方案,并完成了系统仿真.然后深入研究了多相滤波技术在信道化接收机中的应用,推导和建立了实信号接收机的数学模型,给出了真实信道中心频率和带宽的计算公式,简要分析了算法复杂度,最后完成了基于此模型的4信道软件无线电接收机的系统仿真.本文所建立的两个系统作为后续研究的基础平台,可以利用其移植各种通信系统,并分析系统性能,具有一定的应用价值.
数字信号处理大作业
题目一： 利用matlab设计模拟带通巴特沃斯滤波器，要求通带下、上截止频率分别为fhpl=4kHz,fhpu=6kHz，阻带下、上截止频率分别为fhsl=2kHz，fhsu=9kHz，通带最大衰减rp=1dB，阻带最小衰减rs=20dB。 问题分析： （1）该题目中模拟带通滤波器指标：通带下截止频率wpl=2*pi*fhpl，通带上截止频率wpu=2*pi*fhpu通带最大衰减rp=1dB；阻带下截止频率wsl=2*pi*fhsl，阻带上截止频率wsu=2*pi*fhsu；阻带最小衰减rs=20dB。 （2）计算模拟滤波器的阶N和3dB截止频率Wc。 [N,Wc]=buttord(wp,ws,rp,rs,'s'); 其中，wp=[wpl,wpu];ws=[wsl,wsu] （3）计算模拟滤波器系统函数 [Bs,As]=butter(N,Wc,'s'); （4）画图检验所设计出模拟滤波器是否满足指标要求。 原始MATLAB程序： fhpl=4000;fhpu=6000;fhsl=2000;fhsu=9000; wpl=2*pi*fhpl;wpu=2*pi*fhpu;wsl=2*pi*fhsl;wsu=2*pi*fhsu; wp=[wpl,wpu];ws=[wsl,wsu];rp=1;rs=20; [N,Wc]=buttord(wp,ws,rp,rs,'s'); [B,A]=butter(N,Wc,'s'); f=1:20:15000;w=2*pi.*f; [H,WH]=freqs(B,A,w); H=20*log10(abs(H)); plot(f/1000,H),grid on xlabel('频率(KHz)');ylabel('幅度(dB)'); 结果图分析：