第一实验室doc-成都信息工程学院电子工程学院

第一实验室：基础实验篇第Ⅰ部分基本训练题目第Ⅱ部分简介各题目的原理、程序、效果第Ⅲ部分基研训练程序软件压缩文件第Ⅰ部分基本训练题目1-1-1序列的图示方法DSP11011-1-2连续信号及采样信号的图示方法DSP11021-1-3单位冲激序列函数impseq单位冲激序列图示DSP11031-1-4单位阶跃序列函数stepseq单位阶跃序列图示DSP11041-1-5矩形序列)(nRN及图示DSP11051-1-6实指数序列)(nan及图示DSP11061-1-7正弦序列)sin(n及图示DSP11071-1-8复指数序列njme)(及图示DSP11081-1-9周期序列)()(Nnxnx及图示DSP11091-1-10常用5种连续信号及图示DSP11101-1-11离散序列的运算DSP11111-1-12输入序列)(nx与系统冲激响应)(nh的卷积),(hxconvDSP11121-1-13非零起点时两信号的卷积),(hxconvmDSP11131-2-1指数序列)(8.0)(nnxn的离散时间傅立叶变换DSP12011-2-2矩形序列)(nRN的离散时间傅立叶变换DSP12021-2-3离散时间傅立叶变换的性质DSP12031-2-4正弦序列输入，输出为正弦序列，幅度相位因)(jeH变化DSP12041-2-5模拟信号taetx1000)(付氏变换与采样信号的离散时间傅立叶变换DSP12051-3-1N点离散傅立叶变换dft(xn,N)1-3-2N点离散傅立叶反变换idft(xn,N)1-3-3DFT与)(nx的Z变换关系DSP13031-3-4DFT与)(nx的离散时间傅立叶变换的关系DSP13041-3-5有限长序列添零填充，得高密度DFT，离散时间付氏频谱不变DSP13051-3-6采样点增多的高分辨率DFT，采样点数少仅添零的高密度DFTDSP13061-3-7DFT的圆周移位函数cirshftt1-3-8DFT圆周移位实例DSP13081-3-9圆周卷积DSP13091-3-10复共轭序列的DFTDSP13101-3-11DFT的共轭对称性DSP13111-3-12补零填充实现线性卷积DSP13121-3-13重迭保留法实现线性卷积DSP13131-3-14重迭保留实现函数ovrlpsav1-3-15DET对连续信号作近似谱分析：滤高频，避免混迭频谱；截高时；变有限长序列，避免泄漏频谱DSP13151-3-16采样点为100，进行200点DFT，对)(tet进行谱分析DSP13161-3-17实序列的奇偶分解及DFT的虚实分量DSP13171-3-18实序列的奇偶分解函数DSP13181-3-19用FFT分析信号频率成分DSP13191-3-20用FFT分析语言信号的频谱DSP13201-3-21DCT变换DSP13211-3-22用DCT变换进行语言压缩DSP13221-3-23线性调频Z变换DSP13231-3-24利用CZT计算滤波器100—150HZ频率特性的细节DSP13242-1-1直接型实现系统函数H（Z）的IIR数字滤波器DSP21012-1-2级联型实现系统H（Z）的IIR数字滤波器DSP21022-1-3级联型实现H（Z）的IIR数字滤波器DSP21032-1-4直接型实现H（Z）的IIR数字滤波器DSP21042-1-5并联型实现H（Z）的IIR数字滤波器DSP21052-1-6并联型DSP21062-1-7直接型DSP21072-1-8最终的级联，并联DSP21082-1-9直接型级联型dir2cas(b,a)2-1-10级联型直接型cas2par(b0,B,A)2-1-11直接型并联型dir2par2-1-12并联型直接型par2dire2-1-13并联型级联型casfilter2-1-14级联型并联型parfilter2-2-1FIR直接型滤波器DSP22012-2-2FIR级联型滤波器DSP22022-2-3FIR的频率取样形式结构DSP22032-2-4（原例11）由频率样本2,1,0315,......5,415.00)(kkkkH求频率采样形式，及单位冲激响应)(nhDSP22042-2-5窄带滤波器中的频率采取滤波器是由直接型转换为频率采样型dir2fs(n)3-1-1偶对称奇序列的⒈型FIR滤波器的振幅响应hr_type13-1-2偶对称奇序列的)(H及零极点分布DSP31023-1-3偶对称偶序列的Ⅱ型FIR滤波器的振幅响应hr_type23-1-4偶对称偶序列的)(H及零极点分布DSP31043-1-5奇对称奇序列的Ⅲ型FIR滤波器的振幅响应hr_type33-1-6奇对称奇序列的)(H及零极点分布DSP31063-1-7奇对称偶序列的Ⅳ型FIR滤波器的振幅响应hr_type43-1-8奇对称偶序列的)(H及零极点分布DSP31083-1-9线性相位FIR滤波器的零点位置有4种可能DSP31093-1-10常用加窗函数DSP32103-1-11对信号用加窗函数的DFT分析频谱DSP32113-2-1计算理想低通滤波器的)(nhaDSP32013-2-2计算FIR滤波器的绝对和相对的幅度响应DSP32023-2-3提取大于50dB衰减的汉明窗FIR低通滤波器DSP32033-2-4理想高通，偶对称因果序列，N为奇的窗函数，滤波器的单位冲激响应)(nhdDSP32043-2-5汉宁窗，44dB最小阻带衰减，过度带N2.6DSP32053-2-6理想高通，奇对称因果序列，N为偶的窗函数，滤波器的单位冲激响应)(nhdDSP32063-2-7汉宁窗，44dB最小阻带衰减，过度带N2.6DSP32073-2-8理想高通，偶对称因果序列，N为奇的窗函数，滤波器的单位冲激响应)(nhdDSP32083-2-9设计一个数字FIR带通滤波器DSP32093-2-10理想带通数字滤波器的频率响应)(eHjwdDSP32103-2-11设计一个具有2相移的数字FIR带通滤波器DSP32113-2-12理想带阻，偶对称因果序列，N为奇的窗函数，滤波器的单位冲激响应)(nhdideal-be()3-2-13设计一个数字FIR带阻滤波器DSP32133-3-1采样点=0处的频率采样法DSP33013-3-2在过渡带上加两个T1和T2DSP33023-3-3设计2型FIR低通滤波器DSP33033-3-4设计1型FIR高通滤波器DSP33043-3-5设计4型FIR高通滤波器DSP33053-3-6设计2型FIR带通滤波器DSP33063-3-7设计1型FIR带阻滤波器DSP33073-3-8设计1型FIR低通滤波器DSP33083-3-9设计1型FIR高通滤波器DSP33093-3-10设计4型FIR高通滤波器DSP33103-3-11设计3型FIR带通滤波器DSP33113-4-1用频率响应采样法1设计具有线性相位DSP34013-4-2用窗函数法设计具有线性相位DSP34023-4-3用频率采样法1设计低通滤波器对其进行除噪DSP34034-1-1在MATLAB中用afd_butt(Omegap,Omegar,Ap,Ar)函数来设计巴特沃斯模拟低通滤波器DSP41014-1-2若设计非归一化（c≠1）巴特沃斯模拟低通滤波器原型DSP41024-1-3freqs_m(b,a,Omega_max)函数DSP41034-1-4sdir2cas函数DSP41044-1-5设计一个巴特沃斯模拟滤波器DSP41054-2-1用来实现N阶、通带波动为δ的归一化切比学夫1型模拟低通滤波器DSP42014-2-2根据技术指标设计切比学夫1型模拟滤波器DSP42024-2-3设计一个低通切比学夫1型滤波器DSP42034-2-4设计归一化切比学夫2型模拟滤波器DSP42044-2-5根据给定指标设计切比学夫2型模拟滤波器DSP42054-2-6设计一个切比学夫2型低通滤波器DSP42064-3-1用imp_invr函数实现脉冲响应不变法DSP43014-3-2设计一个巴特沃斯模拟滤波器DSP43024-3-3设计低通数字滤波器DSP43034-3-4设计低通数字滤波器DSP43044-4-1双线性变换法设计低通数字滤波器DSP44014-4-2切比雪夫滤波器原型用双线性变换法设计低通数字滤波器DSP44024-5-14-5-24-5-34-5-4分别设计一个巴特沃斯滤波器和切比雪夫高通滤波器DSP45044-5-5分别设计一个巴特沃斯滤波器和切比雪夫高通滤波器DSP45054-5-6设计一个巴特沃斯带通滤波器DSP45064-5-7设计一个切比雪夫带通滤波器DSP45074-5-8设计一个滤波器DSP45084-5-9设计一个滤波器DSP45094-5-10设计一个滤波器DSP45104-6-1zamppingDSP46014-6-2用zmapping函数实现例11中的高通滤波器DSP46024-6-3切比雪夫1型高通数字滤波器，上述过程由chebhpf函数实现DSP46034-6-4用数字频域变换法，设计一个切比雪夫1型高通数字滤波器DSP46044-6-5用双线性变换法设计低通滤波器DSP46054-6-6用脉冲响应不变法设计的低通滤波器对其除噪DSP46064-6-7模拟信号DSP46075-1-1下采样DSP51015-1-2例题DSP51025-1-3上采样DSP51035-1-4程序DSP51045-1-5采样率的非整数倍转换DSP51055-1-6程序DSP51065-1-7例题DSP51075-1-8用傅立叶变换对信号进行消噪声处理DSP51085-1-9信号特定频率的提取DSP51095-1-10例题DSP51105-1-11信号特定频率区间的抑制DSP5111第Ⅱ部分简介各题目的原理、程序、效果1-1-1序列的图示方法DSP1101原理：数字信号处理中，所有信号都是离散时间信号——序列。x(n)={…，x(-1),x(0),x(1),…}如果x(n)={0,5,7,9,6,3,2,1}，-1=n=6。程序：n=-1:6;%序列的序号依兴次取1至6的各整数x=[0,5,7,9,6,3,2,1];%对应序号的序列各值stem(n,x);%调绘离散序列函数图形如下：1-1-2连续信号及采样信号的图示方法DSP1102原理：例)2sin(2)2sin()(21tftfty,当f1=50HZ,f2=120HZ,fs=1000HZ时的信号为)1000240sin(2)1000100sin()(nnny。程序：f1=50;f2=120;fs=1000;%fs表采样频率t=0:1/fs:1;n=t*fs;%时刻t从0至1，步长为1/fsy=sin(2*pi*f1*t)+2*sin(1*pi*f2*t);subplot(211);plot(t(1:50),y(1:50));title('y(t)');subplot(212);stem(n(1:50),y(1:50));title('y(n)');图形如下：1-1-3单位冲激序列函数impseq单位冲激序列图示DSP1103原理：产生)(0nn,可以用函数impseqfunction[x,n]=impseq(n0,n1,n2)n=[n1:n2];%n取从n1至n2的各整数x=[(n-n0)==0];%n仅当n=n0时x值为1，其它x值为0单位冲激序列图示DSP1103n1=-4;%指定参数n2=6;n0=2;impseq(n0,n1,n2);%调用函数stem(n,x)1-1-4单位阶跃序列stepseq单位阶跃序列图示DSP1104function[x,n]=stepseq(n0,n1,n2)n=[n1:n2];%n取从n1至n2的各整数x=[(n-n0)=0];%n当n=n0时x值为1，其它x值为0单位阶跃序列图示DSP114n1=-4;n2=1