DSP实验报告

DSP第六、七次实验报告1. 实验目的:（1）进一步熟悉Matlab实验环境和语言。

（2）熟悉各种滤波器的结构及Matlab实现语言。

（3）掌握用冲击响应不变法和双线性变换法设计IIR滤波器的方法。

（4）掌握用窗函数法和频率抽样法设计FIR滤波器的方法。

2. 实验内容及总结:1.滤波器结构:（1）IIR滤波器各种结构1、直接型结构例如直接型滤波器系统函数, 则有系数向量a=[1,a1,a2,a3],b=[b0,b1,b2], 利用:Y=filter[b,a,x]求信号x(n)通过此滤波器的输出。

2、由系统函数或差分方程求系统的二阶分式（含一阶分式）的级联结构将例如的系统函数重写为二阶分式节的级联型, 利用:[sos,G]=tf2sos(b,a)3、由二阶分式的级联结构转换成系统函数的直接结构是第二步的逆运算, 调用函数:[b,a] = sos2tf(sos)可以求得系数向量a,b, 从而得到H(z)4、由系统函数求部分分式展开（留数及其极点计算）即求z反变换的部分分式展开法, 利用:[r,p,c]=residuez(b,a)其中极点为p, 留数为r, 直接项系数为c。

5、由r,p,c求系统函数即第4步的逆运算, 利用:[b,a]=residuez(r,p,c)6、由直接型结构转换为并联型结构需开发函数:[C,B,A]=tf2par(b,a)其中, b,a为直接型的系数向量, C,B,A为并联型实系数向量, 基本思想是: 1.反复调用[r,p,c]=residuez(b,a)求出极点及留数；2.利用cplxpair函数把极点、留数对按复共轭极点-留数对, 实极点-留数对的顺序排列；3.开发cplxcomp函数, 保证极点和留数相互对应；4.调用[b,a]=residuez(r,p,c)计算并联二阶节的分子分母。

7、由并联型结构转换成直接型结构开发函数:[b,a]=par2tf(C,B,A)为[C,B,A]=tf2par(b,a)的逆函数。

实验一数字IO应用实验—、实验目的1. 了解DSP开发系统的组成和结构2. 在实验设备上完成I/O硬件连接，编写I/O实验程序并运行验证。

3. 内存观察工具的使用二、实验设备计算机，CCS3.1版本软件，DSP仿真器，教学实验箱三、实验原理2.键值读取程序：该部分有两种方法进行键值的判断。

方法1：利用内存观察工具进行观察方法2：利用LED1-LED8的亮灭对应显示键值。

a)外部中断1的应用参照实验五；b)内存观察键值：程序中定义了三个变量“W”“row”和“col”。

“W”代表是CPLD中键盘的扫描数值，“row”和“col”分别代表键盘的行和列，由行和列可以判定按键的位置。

上述三个变量可以在观察窗口中观察的。

c)利用LED灯显示键值原理，参看实验一。

具体的LED灯显示值以查表的形式读出，请参看“”库文件。

本实验的CPLD地址译码说明：基地址：0x0000，当底板片选CS0为低时，分配有效。

CPU的IO空间：基地址+0x0200 LED灯output 8位外部中断用XINT1：由CPLD分配，中断信号由键盘按键产生。

中断下降沿触发。

KEY_DAT_REG(R)：基地址+0x0004;四、实验步骤和内容1.2407CPU板JUMP1的1和2脚短接，拨码开关S1的第一位置ON,其余置OFF；2.E300板上的开关SW4的第一位置ON，其余OFF；SW3的第四位置ON其余的SW置OFF“DEBUG→Connect”）4.打开系统项目文件 \e300.test\ normal \05_key interface \；“\Debug\”文件“Debug\Go Main”跳到主程序的开始；7.指定位置设置断点；8.View--〉Watch Window打开变量观察窗口；9. 将变量“w”“row”和“col”添加到观察窗口中，改变变量观察窗口的显示方式为HEX显示。

“Debug--〉Animate”全速运行，然后点击E300板上键盘按键，观察窗口中变量变化，同时LED1-LED8灯也相应变化，指示键值。

实验一 CCS编程环境及实验箱使用入门一、实验目的1、熟悉DSP集成开发环境CCS，熟悉CCS中工程管理和程序编辑、汇编和链接的过程。

2、熟悉SEED-DTK系列DSP实验平台，掌握CCS中源文件管理、编辑和调试的方法。

二、实验要求1、掌握CCS环境下软件工作及编译方法。

2、了解SEED-DTK实验平台的基本组成，熟悉SEED-DEC5416实验环境。

三、实验内容1、利用SEED-DTK实验箱来完成工程管理，以及源文件管理、编译、汇编、链接和调试。

（1）安装、配置SEED-XDSUSB2.0仿真器。

（2）建立DSP工程（3）建立DSP源文件（4）学习使用CCS的调试工具2、在CCS下开发DSP软件项目需要先建立一个工程文件，扩展名缺省为*.pjt，所有关于项目的信息都会存储在这个工程文件中；若用C语言开发，需要使用Project → Add Files to Project…命令将C语言的标准支持库rts.lib 或rts_ext.lib添加到工程文件中。

此外，还需要将链接器命令文件*.cmd添加到工程文件中。

执行Project → Rebuild All命令，对工程进行编译、汇编和链接，在Output 窗口中将显示相关信息。

目标文件*.out生成成功后，执行菜单命令File →Load Program，选择*.out文件并打开，将生成的可执行程序加载到DSP中，CCS 将自动打开一个“反汇编”窗口，显示加载程序的反汇编指令。

执行菜单命令Debug → Run运行程序。

四、实验步骤1、将DSP仿真器与计算机连接好；2、将DSP仿真器的JTAG插头与SEED-DEC54xx单元的J1连接；3、启动计算机后，打开SEED-DTK实验箱的电源。

观察SEED-DTK-101单元的+5V、+3.3V、+15V和-15V的电源指示灯，以及SEED-DEC54xx的D1以及SEED-DSK2812的D2是否均亮；若有不亮，断开电源，检查电源。

实验一: 闪灯实验熟悉DSP 软硬件测试系统实验目的1.了解SHARC 系列高性能数字信号处理器的程序开发过程和编程语言；2.熟悉集成开发工具VisualDSP++, 学会使用VisualDSP++进行SHARC 系列ADSP 的程序开发、编译与调试；3.掌握SHARC 系列ADSP 的程序加载设计和加载过程。

实验内容利用波形产生信号板, 结合FPGA 编程技术和程序编程器, 编写测试ADSP21065L 和FPGA 之间硬件连接的应用程序, 同时完成应用程序的加载和脱机操作, 在信号指示灯“HL2”上产生可调周期的脉冲信号, “点亮”与“熄灭”指示灯HL2。

实验要求通过DSP 编程, 在其FLAG11引脚上模拟如下波形的周期信号:要求:(1) 500H T ms >,500L T ms >. (2) 并用示波器查看波形, 测量信号周期。

实验步骤1. 熟悉电路图, 清楚波形产生电路板ADSP21065L 与可编程FPGA 器件之间的连接关系；2. 编写FPGA 程序。

在FPGA 内部将ADSP21065L 的标志引脚FLAG11（引脚号26）设置为输出, 作为FPGA 的输入信号, 在FPGA 内部编程将该信号直接输出在发FPGA 的37引脚号上, 设置37引脚为输出信号, 驱动板上的HL2 LED 指示灯；3. 启动VisualDsp++4.5,选择project 工程选项菜单, 创建一个名称为Test.dpj 的工程文件, 选择处理器的型号为ADSP-21065L ；4．弹出一个对话框, 选择是否需要加入VDSP kernel ,选择“NO ”；5. 在工程中加入以下参考源文件:\exp1\test(boot)\ boot1.asm 和boot1.ldf 6．编译, 链接调试, 生成可执行文件。

7．运行程序, 可以看到波形发生电路板上的指示灯“HL2”不断闪动。

8. 利用示波器观测系统时钟,并测量产生信号的波形和周期。

一、综合实验内容和目的1、实验目的(1) 通过实验学习掌握TMS320F28335的浮点处理； (2) 学习并掌握A/D 模块的使用方法；(3) 学习并掌握中断方式和查询方式的相关知识及其相互之间的转换； (4) 学习信号时域分析的方法，了解相关电量参数的计算方法； (5) 了解数字滤波的一些基本方法。

2、实验内容要求1：对给定的波形信号，采用TMS320F28335的浮点功能计算该信号的以下时域参数：信号的周期T ，信号的均方根大小V rms 、平均值V avg 、峰-峰值V pp 。

其中，均方根V rms 的计算公式如下：()21N rms iV u i N =∑ 式中N 为采样点数，()u i 为采样序列中的第i 个采样点。

要求2：所设计软件需要计算采样的波形周期个数，并控制采样点数大于1个波形周期，且小于3个波形周期大小。

要求3：对采集的数据需要加一定的数字滤波。

二、硬件电路相关硬件：TMS320F28335DSP 实验箱，仿真器。

硬件结构图三、程序流程图1、主程序流程图程序的主流程图2、子程序流程图参数计算的流程图四、实验结果和分析1、实验过程分析(1) 使用的函数原型声明对ADC模件相关参数进行定义：ADC时钟预定标，使外设时钟HSPCLK 为25MHz，ADC模块时钟为12.5MHz，采样保持周期为16个ADC时钟。

(2) 定义全局变量根据程序需要，定义相关变量。

主要有：ConversionCount、Voltage[1024]、Voltage1[1024]、Voltage2[1024]、filter_buf[N]、filter_i、Max、Min、T、temp、temp1、temp2、temp3、Num、V、Vav、Vpp、Vrm、fre。

这些变量的声明请见报告后所附的源程序。

(3) 编写主函数完成系统寄存器及GPIO初始化；清除所有中断，初始化PIE向量表，将程序中使用的ADC中断重新映射到ISR；对外部设备ADC模块进行初始化；使能PIE中ADCINT中断；配置ADC模块；等待ADC中断。

实验报告实验一、1、A=[1.2,3,5.1,-1.5,2.4],n=[-2,-1,0,1,2];B=0.5m,m=0~10;请编程完成下面操作并画图表示结果A×B，A-B，A.*B，A（n-2）程序A=[1.2,3,5.1,-1.5,2.4];ns1=[-2,-1,0,1,2];m=0:10;B=(0.5).^m;nf1=-2+length(A)-1;nf2=0+length(B)-1;n1=-2:nf1;n2=0:nf2;n=min(ns1,0):max(nf1,nf2);y1=zeros(1,length(n));y2=y1;y1(find((n>=-2)&(n<=nf1)==1))=A; %给y1赋值Ay2(find((n>=0)&(n<=nf2)==1))=B;y3=y2';ya=y1-y2;ym=y1.*y2;ys=y1*y3;ny=ns1+2;y=A;subplot(321);stem(n1,A,'.');ylabel('A(n)');grid;subplot(322);stem(n2,B,'.');xlabel('n');ylabel('B(n)');grid; subplot(323);stem(n,ya,'.');ylabel('A-B');grid;subplot(324);stem(n,ym,'.');xlabel('n');ylabel('A.*B');grid;subplot(325);stem(ny,y,'.');xlabel('ny');ylabel('A(n-2)');grid; subplot(326);stem(1,ys,'.');xlabel('n');ylabel('A*B');grid;2、p138 2.28n=0:50;x1=0.9.^(abs(n));x2=3*cos(0.5*pi*n+0.25*pi)+2*sin(0.3*pi*n);a=1;N=1000;[X,w]=freqz(x1,a,N);magX=abs(X);angX=angle(X);%subplot(251);stem(n,x1);grid;%title('原序列');ylabel('h(n)');subplot(241);plot(w/pi,magX);grid;title('幅频特性');ylabel('模值');subplot(242);plot(w/pi,angX);grid;title('相频特性');ylabel('弧度');%输入xn,求输出ynny1=n(1)+n(1);ny2=n(end)+n(end);ny=ny1:ny2;y1=conv(x1,x2);subplot(243);stem(n,x2);grid;title('输入序列x(n)');ylabel('y(n)');subplot(244);stem(ny,y1);grid;title('h2(n)*x(n)');ylabel('y(n)');n=0:50;x3=0.3.^n+0.5.^n;x2=3*cos(0.5*pi*n+0.25*pi)+2*sin(0.3*pi*n);a=1;N=1000;[X,w]=freqz(x3,a,N);magX=abs(X);angX=angle(X);%subplot(256);stem(n,x3);grid;%title('原序列');ylabel('h(n)'); subplot(245);plot(w/pi,magX);grid; title('幅频特性');ylabel('模值'); subplot(246);plot(w/pi,angX);grid; title('相频特性');ylabel('弧度');%输入xn,求输出ynny1=n(1)+n(1);ny2=n(end)+n(end);ny=ny1:ny2;y2=conv(x3,x2);subplot(247);stem(n,x2);grid;title('输入序列x(n)');ylabel('y(n)'); subplot(248);stem(ny,y2);grid;title('h2(n)*x(n)');ylabel('y(n)');2.34.clear all;n=0:0.5:39;x=2*cos(pi*n/5);subplot(331);stem(n,x);B1=[1,1];A1=[1,-0.6];y1=filter(B1,A1,x);subplot(332);stem(n,y1);subplot(333);zplane(B1,A1);subplot(334);stem(n,x);B2=[1,1,1];A2=[1,0.5,-0.25]; y2=filter(B2,A2,x);subplot(335);stem(n,y2);subplot(336);zplane(B2,A2);subplot(337);stem(n,x);B3=[1,0,1];A3=[1,-6,9];y3=filter(B3,A3,x);subplot(338);stem(n,y3);subplot(339);zplane(B3,A3);实验二、p211 3.6clc; clear all;n=0:4;x=[2,1,4,2,3];[H, w]=freqz(x,1,'whole'); figure(1);plot(abs(H));hold on;[Y]=fft(x);stem(n*512/5,abs(Y)); grid;n=0:7;x=[2,1,4,2,3,0,0,0];[H, w]=freqz(x,1,'whole'); figure(2);plot(abs(H));hold on;[Y]=fft(x);stem(n*512/8,abs(Y)); grid;Xk =Columns 1 through 77.3223 + 0.0000i -2.3347 - 0.4963i 0.4792 + 0.2134i -0.2926 - 0.2126i 0.1167 + 0.1296i -0.0806 - 0.1397i 0.0357 + 0.1100iColumns 8 through 14-0.0120 - 0.1139i -0.0120 + 0.1139i 0.0357 - 0.1100i -0.0806 + 0.1397i 0.1167 - 0.1296i -0.2926 + 0.2126i 0.4792 - 0.2134iColumn 15-2.3347 + 0.4963i程序：N=8;n=0:1:N-1;xn=3*cos(0.25*pi*n);[Xk]=fft(xn,N);magXk=abs(Xk);angXk=angle(Xk);subplot(221);stem(n,magXk,'.');grid;xlabel('n');title('幅频特性');ylabel('模值');subplot(222);stem(n,angXk,'.');grid;xlabel('n');title('相频特性');ylabel('模值');N=15;n=-7:1:7;for i=1:1:Nxn(i)=0.8.^(abs(i-8));[Xk]=fft(xn,N)magXk=abs(Xk);angXk=angle(Xk);subplot(223);stem(n,magXk,'.');grid;xlabel('n');title('幅频特性');ylabel('模值');subplot(224);stem(n,angXk,'.');grid;xlabel('n');title('相频特性');ylabel('模值');3.32N=11;n=0:1:N-1;k=0:1:N-1;xn=cos(0.48*pi*n)+cos(0.52*pi*n);WN=exp(-j*2*pi/N);nk=n'*k;WNnk=WN.^nk;Xk=xn*WNnk;subplot(221);plot(k,abs(Xk));gridxlabel('k');title('频谱曲线');ylabel('x(k)');N=101;n=0:1:N-1;k=0:1:N-1;n1=0:1:10;xn=cos(0.48*pi*n1)+cos(0.52*pi*n1);xn=[xn,zeros(1,N-length(n1))];%将序列补90个0所需长度WN=exp(-j*2*pi/N);nk=n'*k;WNnk=WN.^nk;Xk=xn*WNnk;subplot(222);plot(k,abs(Xk));gridxlabel('k');title('频谱曲线');ylabel('x(k)');N=101;n=0:1:N-1;k=0:1:N-1;xn=cos(0.48*pi*n)+cos(0.52*pi*n);WN=exp(-j*2*pi/N);nk=n'*k;WNnk=WN.^nk;Xk=xn*WNnk;subplot(223);plot(k,abs(Xk));grid;xlabel('k');title('幅频特性曲线');ylabel('|x(k)|');实验三、5.18 （1）（2）（3）（4）b =1.0000 1.6667 1.3333 1.3333 0.3333 -0.3333a =1.0000 -1.2500 1.2400 -0.8465 0.3675 -0.1147C =8.0587 - 2.9703i8.0587 + 2.9703i31.6073 + 0.0000i-29.5838 -10.1490i-29.5838 +10.1490iB =-0.0000 + 0.7000i-0.0000 - 0.7000i0.6500 + 0.0000i0.3000 + 0.5196i0.3000 - 0.5196iA =17.4429程序：sos=[1,0,1,1,-0.6,0.36;1,-1/3,0,1,-0.65,0;1,2,1,1,0,0.49]; [b,a]=sos2tf(sos)%级联型到直接型的转换[C,B,A]=residuez(6*b,a)%直接型到并联型的转换直接I型直接II型级联型：实验四、P464 7.17椭圆函数数字带通滤波器程序clear,clcFs=25000;wp1=2*5000/Fs;wp2=2*7000/Fs;ws1=2*3500/Fs;ws2=2*8500/Fs; wp=[wp1,wp2];ws=[ws1,ws2];Rp=0.5;As=45;[N,wc]=ellipord(wp,ws,Rp,As);[b,a]=ellip(N,Rp,As,wc);hn=dimpulse(b,a);w0=[ws1*pi,wp1*pi,wp2*pi,ws2*pi];Hx=freqz(b,a,w0);[H,w]=freqz(b,a);dbHx=-20*log10(abs(H)/max(abs(H)));[db,mag,pha,grd,w]=freqz_m(b,a);subplot(2,1,1),plot(12.5*w/pi,db,'linewidth',2);title('椭圆函数数字带通幅度响应（dB)');xlabel('f(KHz)');ylabel('dB');axis([0,12.5,-50,3]);set(gca,'xtickmode','manual','xtick',[1.25,3.5,5,7,8.5,10,12.5]); set(gca,'ytickmode','manual','ytick',[-50,-As,-10,0]);grid; subplot(2,1,2),plot(w/pi,pha,'linewidth',2);title('椭圆函数数字带通相位响应');xlabel('\omega/\pi');ylabel('rad');grid;axis([0,1,-4,4]);figure(2),stem(hn,'.','linewidth',2);title('系统冲激响应');xlabel('n');axis([0,121,-0.15,0.15]);N=4;切比雪夫I型带通滤波器clear,clcFs=25000;wp1=2*5000/Fs;wp2=2*7000/Fs;ws1=2*3500/Fs;ws2=2*8500/Fs; wp=[wp1,wp2];ws=[ws1,ws2];Rp=0.5;As=45;[N,wc]=cheb1ord(wp,ws,Rp,As);[b,a]=cheby1(N,Rp,wc)hn=dimpulse(b,a)w0=[ws1*pi,wp1*pi,wp2*pi,ws2*pi];Hx=freqz(b,a,w0);[H,w]=freqz(b,a);dbHx=-20*log10(abs(H)/max(abs(H)));[db,mag,pha,grd,w]=freqz_m(b,a);subplot(2,1,1),plot(12.5*w/pi,db,'linewidth',2);title('切比雪夫I型数字带通幅度响应（dB)');xlabel('f(KHz)');ylabel('dB');axis([0,12.5,-50,3]);set(gca,'xtickmode','manual','xtick',[1.25,3.5,5,7,8.5,10,12.5]); set(gca,'ytickmode','manual','ytick',[-50,-As,-10,0]);grid; subplot(2,1,2),plot(w/pi,pha,'linewidth',2);title('切比雪夫I型数字带通相位响应');xlabel('\omega/\pi');ylabel('rad');grid;axis([0,1,-4,4]);figure(2),stem(hn,'.','linewidth',2);title('系统冲激响应');xlabel('n');axis([0,151,-0.16,0.16]);N=5;切比雪夫II型带通滤波器clear,clcFs=25000;wp1=2*5000/Fs;wp2=2*7000/Fs;ws1=2*3500/Fs;ws2=2*8500/Fs; wp=[wp1,wp2];ws=[ws1,ws2];Rp=0.5;As=45;[N,wc]=cheb2ord(wp,ws,Rp,As);[b,a]=cheby2(N,As,wc);hn=dimpulse(b,a);w0=[ws1*pi,wp1*pi,wp2*pi,ws2*pi];Hx=freqz(b,a,w0);[H,w]=freqz(b,a);dbHx=-20*log10(abs(H)/max(abs(H)));[db,mag,pha,grd,w]=freqz_m(b,a);subplot(2,1,1),plot(12.5*w/pi,db,'linewidth',2);title('切比雪夫II型数字带通幅度响应（dB)');xlabel('f(KHz)');ylabel('dB');axis([0,12.5,-50,3]);set(gca,'xtickmode','manual','xtick',[1.25,3.5,5,7,8.5,10,12.5]); set(gca,'ytickmode','manual','ytick',[-50,-As,-10,0]);grid; subplot(2,1,2),plot(w/pi,pha,'linewidth',2);title('切比雪夫II型数字带通相位响应');xlabel('\omega/\pi');ylabel('rad');grid;axis([0,1,-4,4]);figure(2),stem(hn,'.','linewidth',2);title('系统冲激响应');xlabel('n');axis([0,60,-0.22,0.22]);N=5;巴特沃斯型带通滤波器clear,clcFs=25000;wp1=2*5000/Fs;wp2=2*7000/Fs;ws1=2*3500/Fs;ws2=2*8500/Fs; wp=[wp1,wp2];ws=[ws1,ws2];Rp=0.5;As=45;[N,wc]=buttord(wp,ws,Rp,As);[b,a]=butter(N,wc);hn=dimpulse(b,a);w0=[ws1*pi,wp1*pi,wp2*pi,ws2*pi];Hx=freqz(b,a,w0);[H,w]=freqz(b,a);dbHx=-20*log10(abs(H)/max(abs(H)));[db,mag,pha,grd,w]=freqz_m(b,a);subplot(2,1,1),plot(12.5*w/pi,db,'linewidth',2);title('巴特沃斯型数字带通幅度响应（dB)');xlabel('f(KHz)');ylabel('dB');axis([0,12.5,-50,3]);set(gca,'xtickmode','manual','xtick',[1.25,3.5,5,7,8.5,10,12.5]); set(gca,'ytickmode','manual','ytick',[-50,-As,-10,0]);grid; subplot(2,1,2),plot(w/pi,pha,'linewidth',2);title('巴特沃斯型数字带通相位响应');xlabel('\omega/\pi');ylabel('rad');grid;axis([0,1,-4,4]);figure(2),stem(hn,'.','linewidth',2);title('系统冲激响应');xlabel('n');axis([0,60,-0.2,0.2]);N=7;实验五、P561 8.11，设计线性相位FIR低通数字滤波器dev =0.0575 0.0002h =Columns 1 through 13-0.0028 -0.0087 -0.0086 0.0044 0.0153 0.0066 -0.0058 0.0032 0.0148 0.0005 -0.0126 0.0073 0.0194Columns 14 through 26-0.0096 -0.0201 0.0199 0.0253 -0.0309 -0.0268 0.0537 0.0302 -0.0991 -0.0307 0.3165 0.5320 0.3165Columns 27 through 39-0.0307 -0.0991 0.0302 0.0537 -0.0268 -0.0309 0.0253 0.0199 -0.0201 -0.0096 0.0194 0.0073 -0.0126Columns 40 through 490.0005 0.0148 0.0032 -0.0058 0.0066 0.0153 0.0044 -0.0086 -0.0087 -0.0028程序：clc,clear all;Rp=1;As=75;Fs=100*10^3;f=[25*10^3,30*10^3];a=[1,0];dev=[(10^(Rp/20)-1)/(10^(Rp/20)+1),10^(-As/20)][M,f0,a0,weights]=firpmord(f,a,dev,Fs);h=firpm(M,f0,a0,weights)[H,f]=freqz(h,1,1024,Fs);subplot(211);plot(f/Fs,20*log10(abs(H)),'linewidth',2);title('幅度响应(dB)');xlabel('f/Fs');ylabel('20log|H(e^j^\omega)|(dB)');axis([0,0.5,-140,20]);set(gca,'xtickmode','manual','xtick',[0,0.1,0.2,0.25,0.3,0.4,0.5]); set(gca,'ytickmode','manual','ytick',[-140,-110,-70,-30,0,20]);grid; subplot(212);plot(f/Fs,angle(H),'linewidth',2);grid;title('相位响应');xlabel('f/Fs');ylabel('arg[H(e^j^\omega)]');8.18用firpm,firpmord设计一个FIR线性相位双带通滤波器M =36f0 =0.10000.20000.40000.50000.60000.70000.80000.90001.0000 a0 =11110 weights =17.19471.000017.19471.000017.1947h =Columns 1 through 130.0159 -0.0002 -0.0071 -0.0027 -0.0281 0.0135 -0.0019 -0.0273 -0.0010 0.0271 0.0458 -0.0003 0.0441Columns 14 through 26-0.0461 0.0363 0.0835 -0.2605 -0.0810 -0.0732 0.0333 0.4606 0.0333 -0.0732 -0.0810 -0.2605 0.0835Columns 27 through 390.0363 -0.0461 0.0441 -0.0003 0.0458 0.0271 -0.0010 -0.0273 -0.0019 0.0135 -0.0281 -0.0027 -0.0071Columns 40 through 41-0.0002 0.0159程序：clc,clear all;Rp=1;As=50;f=[0.1,0.2,0.4,0.5,0.6,0.7,0.8,0.9,];a=[0,1,0,1,0];delta1=(10^(Rp/20)-1)/(10^(Rp/20)+1);delta2=(1+delta1)*(10^(-As/20));dev=[delta2,delta1,delta2,delta1,delta2];[M,f0,a0,weights]=firpmord(f,a,dev)h=firpm(M+4,f0,a0,weights)[H,f]=freqz(h,1,1024);subplot(211);plot(f/pi,20*log10(abs(H)),'linewidth',2);grid;title('幅度响应（dB)');xlabel('\omega/\pi');ylabel('20log|H(e^j^\omega)|(dB)');axis([0,1,-100,10]);set(gca,'xtickmode','manual','xtick',[0:0.1:1]);grid on;set(gca,'ytickmode','manual','ytick',[-100,-50,-1,10]);grid onsubplot(212);plot(f/pi,angle(H),'linewidth',2);gridtitle('相位响应');xlabel('\omega/\pi');ylabel('arg[H(e^j^\omega)]'); figure(2),stem(h,'.','linewidth',2);title('系统冲激响应');xlabel('n');axis([0,40,-0.6,0.6]);grid;N=5;实验六频域抽样定理和音频信号的处理一、频域抽样定理的验证给定信号1, 013 ()27, 14260,n nx n n n+≤≤⎧⎪=-≤≤⎨⎪⎩其它。

dsp实验报告实验一:CCS入门实验实验目的:1. 熟悉CCS集成开发环境，掌握工程的生成方法;熟悉SEED-DEC643实验环境; 掌握CCS集成开发环境的调试方法。

2.学习用标准C 语言编写程序;了解TI CCS开发平台下的C 语言程序设计方法和步骤; 熟悉使用软件仿真方式调试程序。

3. 学习用汇编语言编写程序; 了解汇编语言与 C 语言程序的区别和在设置上的不同;了解TMS320C6000 汇编语言程序结果和一些简单的汇编语句用法学习在CCS 环境中调试汇编代码。

4. 在了解纯C 语言程序工程和汇编语言程序工程结构的基础上，学习在C 工程中加入汇编编程的混合编程方法; 了解混合编程的注意事项;理解混合编程的必要性和在什么情况下要采用混合编程5. 熟悉CCS集成开发环境，掌握工程的生成方法; 熟悉SEED-DEC643实验环境;掌握CCS集成开发环境的调试方法。

实验原理:CCS 提供了配置、建立、调试、跟踪和分析程序的工具，它便于实时、嵌入式信号处理程序的编制和测试，它能够加速开发进程，提高工作效率。

CCS 提供了基本的代码生成工具，它们具有一系列的调试、分析能力序。

使用此命令后，要重新装载.out 文件后，再执行程序。

使用 CCS常遇见文件简介1. program.c: C 程序源文件;2. program.asm: 汇编程序源文件;3. filename.h: C 程序的头文件，包含DSP/BIOS API模块的头文件;4. filename.lib: 库文件;5. project.cmd: 连接命令文件;6. program.obj: 由源文件编译或汇编而得的目标文件;7. program.out: 经完整的编译、汇编以及连接后生成可执行文件; 8. program.map: 经完整的编译、汇编以及连接后生成空间分配文件; 9.project.wks: 存储环境设置信息的工作区文件。

P.S(CMD文件中常用的程序段名与含义1. .cinit 存放C程序中的变量初值和常量;2. .const 存放C程序中的字符常量、浮点常量和用const声明的常量;3. .text 存放C程序的代码;4. .bss 为C 程序中的全局和静态变量保留存储空间;5. .far 为C 程序中用far声明的全局和静态变量保留空间;6. .stack 为 C 程序系统堆栈保留存储空间，用于保存返回地址、函数间的参数传递、存储局部变量和保存中间结果;7. .sysmem 用于 C 程序中malloc、calloc 和 realloc 函数动态分配存储空间。

dsp实验报告DSP实验报告一、引言数字信号处理（Digital Signal Processing，DSP）是一种对数字信号进行处理和分析的技术。

它在许多领域中被广泛应用，如通信、音频处理、图像处理等。

本实验旨在通过实际操作，探索和理解DSP的基本原理和应用。

二、实验目的1. 理解数字信号处理的基本概念和原理；2. 掌握DSP实验平台的使用方法；3. 进行一系列DSP实验，加深对DSP技术的理解。

三、实验器材和软件1. DSP开发板；2. 电脑；3. DSP开发软件。

四、实验内容1. 实验一：信号采集与重构在此实验中，我们将通过DSP开发板采集模拟信号，并将其转换为数字信号进行处理。

首先，我们需要连接信号源和开发板，然后设置采样频率和采样时间。

接下来，我们将对采集到的信号进行重构，还原出原始模拟信号，并进行观察和分析。

2. 实验二：滤波器设计与实现滤波器是DSP中常用的模块，用于去除或增强信号中的特定频率成分。

在此实验中，我们将学习滤波器的设计和实现方法。

首先，我们将选择合适的滤波器类型和参数，然后使用DSP开发软件进行滤波器设计。

最后，我们将将设计好的滤波器加载到DSP开发板上，并进行实时滤波处理。

3. 实验三：频谱分析与频域处理频谱分析是DSP中常用的方法，用于分析信号的频率成分和能量分布。

在此实验中，我们将学习频谱分析的基本原理和方法，并进行实际操作。

我们将采集一个包含多个频率成分的信号，并使用FFT算法进行频谱分析。

然后，我们将对频谱进行处理，如频率选择、频率域滤波等，并观察处理后的效果。

4. 实验四：音频处理与效果实现音频处理是DSP中的重要应用之一。

在此实验中，我们将学习音频信号的处理方法，并实现一些常见的音频效果。

例如，均衡器、混响、合唱等。

我们将使用DSP开发软件进行算法设计，并将设计好的算法加载到DSP开发板上进行实时处理。

五、实验结果与分析通过以上实验，我们成功完成了信号采集与重构、滤波器设计与实现、频谱分析与频域处理以及音频处理与效果实现等一系列实验。

dsp信号处理实验报告DSP信号处理实验报告一、引言数字信号处理（DSP）是一种将连续信号转换为离散信号，并对其进行处理和分析的技术。

在现代通信、音频处理、图像处理等领域中，DSP技术被广泛应用。

本实验旨在通过对DSP信号处理的实践，加深对该技术的理解与应用。

二、实验目的本实验旨在通过对DSP信号处理的实践，掌握以下内容：1. 学习使用DSP芯片进行信号采集和处理；2. 理解离散信号的采样和重构过程；3. 掌握常见的DSP信号处理算法和方法。

三、实验原理1. 信号采集与重构在DSP信号处理中，首先需要对模拟信号进行采样，将连续信号转换为离散信号。

采样过程中需要注意采样频率的选择，以避免混叠现象的发生。

采样完成后，需要对离散信号进行重构，恢复为连续信号。

2. DSP信号处理算法DSP信号处理涉及到多种算法和方法，如滤波、频谱分析、时域分析等。

其中，滤波是一种常见的信号处理方法，可以通过滤波器对信号进行去噪、增强等处理。

频谱分析可以将信号在频域上进行分析，了解信号的频率成分和能量分布。

时域分析则关注信号的时序特征，如幅值、相位等。

四、实验步骤1. 信号采集与重构在实验中，我们使用DSP芯片进行信号采集与重构。

将模拟信号输入DSP芯片的模拟输入端口，通过ADC（模数转换器）将模拟信号转换为数字信号。

然后，通过DAC（数模转换器）将数字信号转换为模拟信号输出。

2. 滤波处理为了演示滤波处理的效果，我们选择了一个含有噪声的信号进行处理。

首先，使用FIR滤波器对信号进行低通滤波，去除高频噪声。

然后，使用IIR滤波器对信号进行高通滤波，增强低频成分。

3. 频谱分析为了对信号的频率成分和能量分布进行分析，我们使用FFT（快速傅里叶变换）算法对信号进行频谱分析。

通过观察频谱图，可以了解信号的频率特性。

4. 时域分析为了对信号的时序特征进行分析，我们使用时域分析方法对信号进行处理。

通过计算信号的均值、方差、峰值等指标，可以了解信号的幅值、相位等特性。