杭电数字信号处理实验10
信号、系统与信号处理实验Ⅱ
实验报告
姓名：王健
学号：14072119
班级：14083413
上课时间：周五-六七八
实验名称：IIR和FIR数字滤波器过滤信号的实现及比较
一、实验目的
1.掌握数字滤波器的计算机仿真方法。
2.通过观察对实际心电图信号的滤波作用，获得数字滤波器的感性知识。
二、实验原理与要求
实验原理
若h(n)是因果FIR滤波器的单位冲激响应，其长度为N，当输入为x(n)时，输出序列y(n)为:
取h(n)为低通滤波器，N=21，截止频率为0.2*pi，采用海明窗设计出h(n)。
实验要求
1.编写FIR数字滤波器实现程序，计算其对心电图信号取样序列x(n)的响应序列。
2.编写IIR数字滤波器实现程序，计算其对心电图信号取样序列x(n)的响应序列。
3.在通用计算机上运行上述两个程序，并画图进行比较
三、实验程序与结果
clear;clc;close all
xn=[ 4 -2 0 -4 -6 -4 -2 -4 -6 -6
-4 -4 -6 -6 -2 6 12 8 0 -16
-38 -60 -84 -90 -66 -32 -4 -2 -4 8
12 12 10 6 6 6 4 0 0 0
0 0 -2 -4 0 0 0 -2 -2 0
0 -2 -2 -2 -2 0];
y1=filter(0.2318.*[1 2 1],[1 0.4164 -0.6006],xn);
y2=filter(0.2318.*[1 2 1],[1 0.3093 -0.1888],y1);
y3=filter(0.2318.*[1 2 1],[1 0.2692 -0.0349],y2);
b=fir1(20,0.2);
yf=filter(b,1,xn)
subplot(3,1,1)
plot(0:length(xn)-1,xn);
title('原信号')
axis([0 60 -100 20])
subplot(3,1,2)
plot(0:length(y3)-1,y3);
title('IIR 滤波器输出')
axis([0 60 -100 20])
subplot(3,1,3)
plot(0:length(yf)-1,yf);
title('FIR 滤波器输出')
axis([0 60 -100 20])
结果：
0102030
405060
-100-50
原信号
0102030
405060-100-50
IIR 滤波器输出
0102030405060-100-50
FIR 滤波器输出
四、 仿真结果分析 两种滤波器都能有效地进行滤波，消除了高频分量使原波形更加的平滑。而且两者都存在延迟。
IIR 滤波器延迟大概3个单位
FIR 滤波器延迟20个单位，和滤波器阶数相同。
五、 实验问题解答与体会
一次数字信号处理实验，虽然题目看起来简单，但是编程的时候却会有陷阱，加上自己的粗心用了好久才完成。以后做实验一定不能大意，一定要预习，特别是例程，例程能很好地表达函数用法，使逻辑更加清楚。
另外，在以后实验的时候一定要带上数字信号处理的教材，因为实验能很好地实践验证教材所教的东西，加深自己的理解，纠正自己的错误观念，结合教材去验证加深知识，而不是一味为完成实验而做实验
信号与系统实验总结1
实验总结 班级：10电子班学号：1039035 姓名：田金龙这学期的实验都有：信号的时域分析、线性时不变系统的时域分析、连续时间信号系统的频域分析、连续时间在连续时间信号的频域LTI系统的复频域分析、连续时间LTI系统的频域分析。在这学期的学习中学习了解到很多关于信号方面的处理方法加上硬件动手的实践能力，让我对课堂上所学到的知识有了更深层次的理解也加深了所学知识的印象。下面则是对每次实验的分析和总结： 实验一：信号的时域分析 在第一次试验中进行信号的时域分析还有的就是学会使用MATLAB软件来利用它实现一些相关的运算并且绘制出相关的信号图。在时域分析中掌握连续时间信号和离散时间信号的描述方法，并能够实现各种信号的时域变化和运算。了解单位阶跃信号和单位冲激信号的拓展函数，以便于熟悉这两种函数在之后的程序中的应用。在能够对简单信号的描述的前提下，通过一些简单的程序，实现信号的分析，时域反相，时域尺度变换和周期信号的描述。 clear, close all dt=0.01; t=-2:dt:2; x=u(t); plot(t,x) title('u signal u(t)') grid on 连续时间信号的时域分析后，则是离散时间信号的仿真。通过对连续时间信号的描述和对离散时间信号的描述，发现它们的不同之处在于对时间的定义和对函数的图形描述。在离散时间信号的图形窗口描述时，使用的是stem（n，x）函数。 在硬件实验中，使用一些信号运算单元，加法器，减法器，倍乘器，反相器，积分器和微分器。输入相应的简单信号，观察通过不同运算单元输出的信号。 实验二：线性时不变系统的时域分析 在线性时不变系统的时域分析中主要研究的就是信号的卷积运算，学会进行信号的卷积
杭电数字信号处理实验7
信号、系统与信号处理实验Ⅱ 实验报告 姓名：王健 学号：14072119 班级：14083413 上课时间：周五-六七八
实验名称：用双线性变换法设计IIR数字滤波器 一、实验目的 熟悉模拟巴特沃兹滤波器设计和用双线性变换法设计IIR数字滤波器的方法 二、实验原理与要求 实验原理 利用双线性变换法设计IIR数字滤波器，首先要设计出满足指标要求的模拟滤波器的传递函数Ha(s),然后由Ha(s)通过双线性变换可得要设计的IIR数字滤波器的系统函数H(z),如果给定的指标为数字滤波器的指标，直接利用模拟滤波器的低通原理，通过式子 到式子 的频率变换关系，可一步完成数字滤波器的设计。式中是低通模拟滤波器的截止频率 实验要求 （1）编写用双线性变换法设计的巴特沃兹低通IIR滤波器的程序，要求通带内频率低于,容许幅度误差在1dB之内，频率在到之间的阻带衰减大于10dB。 （2）用法设计的巴特沃兹低通IIR滤波器，要求使用buttord,butter和biliner函数，滤波器技术指标：取样频率为1Hz;通带内衰减小于1Db; 阻带临界频率0.3Hz,阻带内衰减大于25dB。 （3）以pi/64为取样间隔，在屏幕上打印出数字滤波器的频率区间[0 pi]上的幅频响应特性曲线。 （4）在屏幕上打印出H(z)的分子，分母多项式系数。 三、实验程序与结果 1. 用双线性变换法设计的巴特沃兹低通IIR滤波器的程序，要求通带内频率低于,容许幅度误差在1dB之内，频率在到之间的阻带衰减大于10dB。 clear;clc;close all; Rp=1; Rs=10; Fs=1; Ts=1/Fs
《数字信号处理》课程研究性学习报告解读
《数字信号处理》课程研究性学习报告 指导教师薛健 时间2014.6
【目的】 (1) 掌握IIR 和FIR 数字滤波器的设计和应用； (2) 掌握多速率信号处理中的基本概念和方法 ； (3) 学会用Matlab 计算小波分解和重建。 (4)了解小波压缩和去噪的基本原理和方法。 【研讨题目】 一、 (1)播放音频信号 yourn.wav ，确定信号的抽样频率，计算信号的频谱，确定噪声信号的频率范围； (2)设计IIR 数字滤波器，滤除音频信号中的噪声。通过实验研究s P ,ΩΩ，s P ,A A 的选择对滤波效果及滤波器阶数的影响，给出滤波器指标选择的基本原则，确定你认为最合适的滤波器指标。 （3）设计FIR 数字滤波器，滤除音频信号中的噪声。与（2）中的IIR 数字滤波器，从滤波效果、幅度响应、相位响应、滤波器阶数等方面进行比较。 【设计步骤】 【仿真结果】
【结果分析】 由频谱知噪声频率大于3800Hz。FIR和IIR都可以实现滤波，但从听觉上讲，人对于听觉不如对图像（视觉）明感，没必要要求线性相位，因此，综合来看选IIR滤波器好一点，因为在同等要求下，IIR滤波器阶数可以做的很低而FIR滤波器阶数太高，自身线性相位的良好特性在此处用处不大。【自主学习内容】 MATLAB滤波器设计 【阅读文献】 老师课件，教材 【发现问题】(专题研讨或相关知识点学习中发现的问题)： 过渡带的宽度会影响滤波器阶数N 【问题探究】 通过实验，但过渡带越宽时，N越小，滤波器阶数越低，过渡带越窄反之。这与理论相符合。 【仿真程序】 信号初步处理部分： [x1,Fs,bits] = wavread('yourn.wav'); sound(x1,Fs); y1=fft(x1,1024); f=Fs*(0:511)/1024; figure(1) plot(x1) title('原始语音信号时域图谱'); xlabel('time n'); ylabel('magnitude n'); figure(2) freqz(x1) title('频率响应图') figure(3) subplot(2,1,1); plot(abs(y1(1:512))) title('原始语音信号FFT频谱') subplot(2,1,2); plot(f,abs(y1(1:512))); title(‘原始语音信号频谱') xlabel('Hz'); ylabel('magnitude'); IIR： fp=2500;fs=3500; wp = 2*pi*fp/FS; ws = 2*pi*fs/FS; Rp=1; Rs=15;
杭电《过程控制系统》实验报告
实验时间：5月25号 序号： 杭州电子科技大学 自动化学院实验报告 课程名称：自动化仪表与过程控制 实验名称：一阶单容上水箱对象特性测试实验 实验名称：上水箱液位PID整定实验 实验名称：上水箱下水箱液位串级控制实验 指导教师：尚群立 学生姓名：俞超栋 学生学号：09061821
实验一、一阶单容上水箱对象特性测试实验一．实验目的 （1）熟悉单容水箱的数学模型及其阶跃响应曲线。 （2）根据由实际测得的单容水箱液位的阶跃响应曲线，用相关的方法分别确定它们的参数。二．实验设备 AE2000型过程控制实验装置，PC机，DCS控制系统与监控软件。 三、系统结构框图 单容水箱如图1-1所示： Q2 图1-1、单容水箱系统结构图 四、实验原理 阶跃响应测试法是系统在开环运行条件下，待系统稳定后，通过调节器或其他操作器，手动改变对象的输入信号（阶跃信号），同时记录对象的输出数据或阶跃响应曲线。然后根据已给定对象模型的结构形式，对实验数据进行处理，确定模型中各参数。 图解法是确定模型参数的一种实用方法。不同的模型结构，有不同的图解方法。单容水箱对象模型用一阶加时滞环节来近似描述时，常可用两点法直接求取对象参数。 如图1-1所示，设水箱的进水量为Q1，出水量为Q2，水箱的液面高度为h，出水阀
h1( t ) h1(∞ ) 0.63h1(∞) 0 T V 2固定于某一开度值。根据物料动态平衡的关系，求得： 在零初始条件下，对上式求拉氏变换，得： 式中，T 为水箱的时间常数（注意：阀V 2的开度大小会影响到水箱的时间常数），T=R 2*C ，K=R 2为单容对象的放大倍数，R 1、R 2分别为V 1、V 2阀的液阻，C 为水箱的容量系数。令输入流量Q 1 的阶跃变化量为R 0，其拉氏变换式为Q 1（S ）=R O /S ，R O 为常量，则输出液位高度的拉氏变换式为： 当t=T 时，则有： h(T)=KR 0(1-e -1)=0.632KR 0=0.632h(∞) 即 h(t)=KR 0(1-e -t/T ) 当t —>∞时，h （∞）=KR 0，因而有 K=h （∞）/R0=输出稳态值/阶跃输入 式（1-2）表示一阶惯性环节的响应曲线是一单调上升的指数函数，如图1-2所示。当由实验求得图1-2所示的 阶跃响应曲线后，该曲线上升到稳态值的63%所对应时间，就是水箱的时间常数T ，该时间常数T 也可以通过坐标原点对响应曲线 图 1-2、 阶跃响应曲线
杭电_数字信号处理课程设计_实验5
实验5 IIR和FIR滤波器过滤信号的实现及比较：以心电信号为例 一、实验目的 1、探究心电信号的初步分析。心电信号(频率-般在0.05Hz ~100Hz范围)是一种基本的人体生理信号，体表检测人体心电信号中常带有工频干扰(50HZ)、基线漂移(频率低于0.5Hz)和肌电干扰等各种噪声。 2、为了得到不失真的原始心电信号，需要滤波预处理。设计数字低通滤波器、高通滤波器、带阻滤波器，用MATLAB软件对含噪心电信号分别进行高通、带阻和低通滤波等处理，将心电信号中的低频基线漂移、50Hz 工频高频和高频杂波进行滤除。 3、通过观察对含噪心电图信号的滤波作用，获得数字滤波的感性知识。 二、实验要求及内容 实验题目： 给定一组干净心电信号数据，数据文件存于C盘Ecg.txt。采样频率Fs = 500Hz。 1、编写程序读出心电信号，并在屏幕上打印出其波形。 2、产生模拟高斯白噪声信号，与干净心电混合,设计一个IIR低通滤波器和一个FIR 低通滤波器分别滤除心电信号中的白噪声干扰，调整白噪声信噪比大小，对滤波前后的心电信号的频谱进行分析比较。其中数字低通滤波器指标要求，通带截止频率Wp=0.1π，阻带截止频率 Ws=0.16π，阻带衰减不小于15 dB，通带衰减不大于1 dB。 要求:编写一个IIR低通滤波器和一个FIR低通滤波器仿真程序，在屏幕上打印出数字滤波器的频率区间[0, π]上的幅频响应特性由线(H(e^jw)) ;计算其对含噪心电信号的低通滤波响应序列，并在屏幕上打印出干净心电信号波形，含工频干扰的心电信号波形以及IIR低通滤波和FIR低通后的信号波形，并进行比较;同时对滤波前后的心电信号的频谱进行分析比较，并在屏幕上打印出滤波前后的心电信号的频谱，观察其变化。 3、产生模拟工频信号，与干净心电混合，设计一个带阻滤波器(50Hz 陷波器)滤除心电信号中的电源线干扰，调整工频幅度大小，对滤波前后的心电信号的频谱进行分析比较。其中带阻滤波器指标要求，通带下限频率Wp1=0.18π，阻带下截止频率Ws1=0.192 π，阻带上截止频率Ws2=0.208π，通带上限频率Wp2=0.22π，阻带衰减不小于15 dB, 通带衰减不大于1 dB。 要求:编写IIR带阻滤波器仿真程序,在屏幕上打印出数字滤波器的频率区间[0, π]上的幅频响应特性由线(H(e^jw ));计算其对含工频干扰的心电信号的带阻滤波响应序列，并在屏幕上打印出干净心电信号波形，含工频干扰的心电信号波形以及滤波后的信号波形，并进行比较;同时对滤波前后的心电信号的频谱进行分析比较，并在屏幕上打印出滤波前后的心电信号的频谱，观察其变化。 4、产生模拟基线漂移信号，与干净心电信号混合，设计一个高通滤波器滤除心电信号中的基线低频干扰，调整基线的幅度大小,对滤波前后的心电信号的频谱进行分析比较。其中，高通滤波器指标要求，通带截止频率Wp=0.0028π，阻带截止频率Ws=0.0012π，阻带衰减不小于15 dB，通带衰减不大于1 dB。 要求:编写IIR高通滤波器(或FIR高通滤波器)仿真程序，在屏幕上打印出数字滤波器的频率区间[0,π]上的幅频响应特性由线(H(e^jw);计算其对含基线低频干扰的心电信号的高通滤波响应序
数字信号处理课程实验报告4
数字信号处理课程实验报告 实验名称FIR数字滤 班级姓名 波器设计 教师姓名实验地点实验日期 一、实验内容 1、设计一个最小阶次的低通FIR数字滤波器，性能指标为：通带0Hz~1500Hz，阻带截 止频率2000Hz，通带波动不大于1%，阻带波动不大于1%，采样频率为8000Hz； 2、用一个仿真信号来验证滤波器的正确性（注意：要满足幅度要求和线性相位特性）。 二、实验目的 1、利用学习到的数字信号处理知识解决实际问题； 2、了解线性相位滤波器的特殊结构； 3、熟悉FIR数字滤波器的设计方法。 三、涉及实验的相关情况介绍（包含使用软件或实验设备等情况） 计算机一台(安装MATLAB6.5版本或以上版本) 四、实验记录（以下1~5项必须完成，第6项为选择性试做） 1.原理基础 令希望设计的滤波器的传输函数是H（ejw，hd（n）是与其对应的单位脉冲响应。一般情况下，由Hd（ejw）求出hd（n），然后由Z变换求出滤波器的系统函数。但是通常Hd（ejw）在边界频率处有不连续点，这使得hd（n）是无限长的非因果序列，所以实际是不能实现的。为了构造一个长度为N的线性相位滤波器，可以将hd（n）截取一段来近似，并且根据线性相位的特点，需要保证截取后的序列关于（N-1）/2对称。设截取的一段为h(n)，则 Wr（n）称为矩形窗函数。 当hd（n的对称中心点取值为（N-1)/2时，就可以保证所设计的滤波器具有线性相位。 2 实验流程
1.信号的谱分析 2.信号的采样 3.信号的恢复 3源程序代码 clc; clear all; close all; fs=700;%采样频率 f=[30 40];%截止频率 a=[1 0]; dev=[0.01 0.1]; % dev纹波 [n,fo,ao,w]=remezord(f,a,dev,fs);%n滤波器阶数fo过渡带起止频率ao频带内幅度————firpmord b=remez(n,fo,ao,w);%firpm b=b.*blackman(length(b))'; b=b; a=1; figure(1) % [H,W]=freqz(b,1,1024,Fs); % plot(W,20*log10(abs(H))); freqz(b,1,1024,fs);grid title('滤波器') grid %%%%%%%%%%%%%%%% fc=28; fcl1=50; fcl2=100; fcl3=150; N=1024; n=1:N; % x=2*cos(2*pi*fc/fs*n)+j*2*sin(2*pi*fc/fs*n)+cos(2*pi*fcl/fs*n)+j*sin(2*pi*fcl/fs*n)+1*r and(1,N); xc=2*cos(2*pi*fc/fs*n); x=2*cos(2*pi*fc/fs*n)+2*cos(2*pi*fcl1/fs*n)+2*cos(2*pi*fcl2/fs*n)+0.1*rand(1,N); % x=2*cos(2*pi*fc/fs*n); xfft=abs(fft(x,N));
信号与系统实验实验报告
信号与系统实验实验报 告 文档编制序号：[KKIDT-LLE0828-LLETD298-POI08]
实验五连续系统分析一、实验目的 深刻理解连续时间系统的系统函数在分析连续系统的时域特性、频域特性及稳定性中的重要作用及意义，掌握根据系统函数的零极点设计简单的滤波器的方法。掌握利用MATLAB分析连续系统的时域响应、频响特性和零极点的基本方法。 二、实验原理 MATLAB提供了许多可用于分析线性时不变连续系统的函数，主要包含有系统函数、系统时域响应、系统频域响应等分析函数。 三、实验内容 1.已知描述连续系统的微分方程为，输入，初始状态 ，计算该系统的响应，并与理论结果比较，列出系统响应分析的步骤。 实验代码： a=[1 10]; b=[2]; [A B C D]=tf2ss(b,a); sys=ss(A,B,C,D); t=0: :5; xt=t>0; sta=[1]; y=lsim(sys,xt,t,sta); subplot(3,1,1); plot(t,y); xlabel('t'); title('系统完全响应 y(t)'); subplot(3,1,2); plot(t,y,'-b'); hold on yt=4/5*exp(-10*t)+1/5; plot(t,yt,' : r'); legend('数值计算','理论计算'); hold off xlabel('t'); subplot(3, 1 ,3); k=y'-yt; plot(t,k); k(1) title('误差');
实验结果： 结果分析： 理论值 y(t)=0. 8*exp(-10t)*u(t)+ 程序运行出的结果与理论预期结果相差较大误差随时间增大而变小，初始值相差最大，而后两曲线基本吻合，表明该算法的系统响应在终值附近有很高的契合度，而在初值附近有较大的误差。 2.已知连续时间系统的系统函数为，求输入分别为，， 时，系统地输出，并与理论结果比较。 a=[1,3,2,0]; b=[4,1]; sys=tf(b,a); t=0: :5; x1=t>0; x2=(sin(t)).*(t>0); x3=(exp(-t)).*(t>0); y1=lsim(sys,x1,t); y2=lsim(sys,x2,t); y3=lsim(sys,x3,t); subplot(3,1,1); plot(t,y1); xlabel('t'); title('X(t)=u(t)'); subplot(3,1,2); plot(t,y2); xlabel('t'); title('X(t)=sint*u(t)'); subplot(3, 1 ,3); plot(t,y3); xlabel('t'); title('X(t)=exp(-t)u(t)'); 实验结果： 结果分析： a=[1,3,2,0]; b=[4,1]; sys=tf(b,a); t=0: :5; x1=t>0; x2=(sin(t)).*(t>0); x3=(exp(-t)).*(t>0); y1=lsim(sys,x1,t); y2=lsim(sys,x2,t); y3=lsim(sys,x3,t); subplot(3,1,1); plot(t,y1,'-b');
数字信号处理实验报告(实验1_4)
实验一 MATLAB 仿真软件的基本操作命令和使用方法 实验容 1、帮助命令 使用 help 命令，查找 sqrt （开方）函数的使用方法； 2、MATLAB 命令窗口 （1）在MATLAB 命令窗口直接输入命令行计算3 1)5.0sin(21+=πy 的值； （2）求多项式 p(x) = x3 + 2x+ 4的根； 3、矩阵运算 （1）矩阵的乘法 已知 A=[1 2;3 4]， B=[5 5;7 8]，求 A^2*B
（2）矩阵的行列式 已知A=[1 2 3;4 5 6;7 8 9]，求A （3）矩阵的转置及共轭转置 已知A=[1 2 3;4 5 6;7 8 9]，求A' 已知B=[5+i,2-i,1;6*i,4,9-i], 求B.' , B' （4）特征值、特征向量、特征多项式 已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ，求矩阵A的特征值、特征向量、特征多项式；
（5）使用冒号选出指定元素 已知：A=[1 2 3;4 5 6;7 8 9]；求A 中第3 列前2 个元素；A 中所有列第2，3 行的元素； 4、Matlab 基本编程方法 （1）编写命令文件：计算1+2+…+n<2000 时的最大n 值；
（2）编写函数文件：分别用for 和while 循环结构编写程序，求 2 的0 到15 次幂的和。
5、MATLAB基本绘图命令 （1）绘制余弦曲线 y=cos(t)，t∈[0，2π]
（2）在同一坐标系中绘制余弦曲线 y=cos(t-0.25)和正弦曲线 y=sin(t-0.5)， t∈[0，2π] （3）绘制[0，4π]区间上的 x1=10sint 曲线，并要求： （a）线形为点划线、颜色为红色、数据点标记为加号； （b）坐标轴控制：显示围、刻度线、比例、网络线 （c）标注控制：坐标轴名称、标题、相应文本； >> clear;
杭电通信系统课程设计报告实验报告
通信系统课程设计实验报告 XX：田昕煜 学号：13081405 班级：通信四班 班级号：13083414 基于FSK调制的PC机通信电路设计
一、目的、容与要求 目的: 掌握用FSK调制和解调实现数据通信的方法,掌握FSK调制和解调电路中相关模块的设计方法。初步体验从事通信产品研发的过程. 课程设计任务：设计并制作能实现全双工FSK调制解调器电路，掌握用Orcad Pspice、Protel99se进行系统设计及电路仿真。 要求：合理设计各个电路，尽量使仿真时的频率响应和其他参数达到设计要求。尽量选择符合标称值的元器件构成电路，正确完成电路调试。 二、总体方案设计 信号调制过程如下： 调制数据由信号发生器产生（电平为TTL，波特率不超过9600Baud），送入电平/幅度调整电路完成电平的变换，再经过锁相环（CD4046）,产生两个频率信号分别为30kHz和40kHz（发“1”时产生30kHz方波，发“0”时产生40kHz方波），再经过低通滤波器2，变成平滑的正弦波，最后通过线圈实现单端到差分信号的转换。
信号的解调过程如下： 首先经过带通滤波器1，滤除带外噪声，实现信号的提取。在本设计中FSK 信号的解调方式是过零检测法。所以还要经过比较器使正弦信号变成方波，再经过微分、整流电路和低通滤波器1实现信号的解调，最后经过比较器使解调信号成为TTL电平。在示波器上会看到接收数据和发送数据是一致的。 各主要电路模块作用： 电平/幅度调整电路：完成TTL电平到VCO控制电压的调整； VCO电路：在控制电压作用下，产生30KHz和40KHz方波； 低通2：把30KHz、40KHz方波滤成正弦波； 线圈：完成单端信号和差分信号的相互转换； 带通1：对带外信号抑制，完成带信号的提取； 限放电路：正弦波整形成方波，同时保留了过零点的信息； 微分、整流、脉冲形成电路：完成信号过零点的提取； 低通1：提取基带信号，实现初步解调； 比较器：把初步解调后的信号转换成TTL电平 三、单元电路设计原理与仿真分析 (1)带通1(4阶带通)-- 接收滤波器（对带外信号抑制，完成带信号的提取） 要求通带：26KHz—46KHz，通带波动3dB； 阻带截止频率：fc＝75KHz时，要求衰减大于10dB。经分析，二级四阶巴特沃斯带通滤波器来提取信号。 具体数值和电路见图1仿真结果见图2。
数字信号处理实验及参考程序
数字信号处理实验实验一离散时间信号与系统及MA TLAB实现 1.单位冲激信号： n = -5:5; x = (n==0); subplot(122); stem(n, x); 2.单位阶跃信号： x=zeros(1,11); n0=0; n1=-5; n2=5; n = n1:n2; x(:,n+6) = ((n-n0)>=0); stem(n,x); 3.正弦序列： n = 0:1/3200:1/100; x=3*sin(200*pi*n+1.2); stem(n,x); 4.指数序列 n = 0:1/2:10; x1= 3*(0.7.^n); x2=3*exp((0.7+j*314)*n); subplot(221); stem(n,x1); subplot(222); stem(n,x2); 5.信号延迟 n=0:20; Y1=sin(100*n); Y2=sin(100*(n-3)); subplot(221); stem(n,Y1); subplot(222); stem(n,Y2);
6.信号相加 X1=[2 0.5 0.9 1 0 0 0 0]; X2=[0 0.1 0.2 0.3 0.4 0.5 0.6 0.7]; X=X1+X2; stem(X); 7.信号翻转 X1=[2 0.5 0.9 1]; n=1:4; X2=X1(5-n); subplot(221); stem(n,X1); subplot(222); stem(n,X2); 8.用MATLAB计算序列{-2 0 1 –1 3}和序列{1 2 0 -1}的离散卷积。a=[-2 0 1 -1 3]; b=[1 2 0 -1]; c=conv(a,b); M=length(c)-1; n=0:1:M; stem(n,c); xlabel('n'); ylabel('幅度'); 9．用MA TLAB计算差分方程 当输入序列为时的输出结果。 N=41; a=[0.8 -0.44 0.36 0.22]; b=[1 0.7 -0.45 -0.6]; x=[1 zeros(1,N-1)]; k=0:1:N-1; y=filter(a,b,x); stem(k,y) xlabel('n'); ylabel('幅度') 10.冲激响应impz N=64; a=[0.8 -0.44 0.36 0.22];
07杭电信号与系统期末试题
2007信号卷 一．填空题（每小题3分，10小题，共30分） 1．信号)3 π cos()4πsin()(t t t f +=的基本周期是 。 2．信号)()(n u n x =的功率是 。 3．=+?∞ -ττδd )1(t 。 4．信号)()(t u t f =的傅里叶变换为 。 5．信号)()(n u n x =的算子表示为 。 6．{}{}=-*--2012 112 。 7．已知LTI 系统方程)()()(2)(d d t u t t r t r t +=+δ且1)0(=-r ，则=+)0(r 。 8．无失真传输系统)1(2)(-=t e t r ，其冲激响应为=)(t h 。 9．信号)()1()(t u t t f +=的拉氏变换为 。 10．已知)21(232 3)(22<<+-+=z z z z z X ，则=)(n x 。 解答：1．24；2．0.5 ；3．)1(+t u ；4．ωωδj 1)(π+；5．)(1 n E E δ-； 6．{}21304--；7．2；8．)1(2-t δ；9． )0(12>+σs s ； 10．())1(27)(5)(----n u n u n n δ 二．画图题（每小题5分，4小题，共20分） 1．信号)(t f 的波形如题图2－1，画出)42(+t f 的波形。 题图2－1 解： 2．已知周期函数)(t f 半个周期的波形如题图2-2，根据下列条件画出)(t f 在一个周期()10T t <≤ 的波 形。（1）)(t f 是偶函数； （2）)(t f 是奇函数。
题图2－2 解：（1）()t f 是偶函数，则()()t f t f =-，波形对称于纵轴。 题图2-12 ① 对褶得()t f 1 ②将()t f 1向右平移1T 得()t f 2 ③取10T -的波形得到()t f 在一 个周期()10T t <≤ 的波形。如图（1）所示。 图（1） （2）)(t f 是奇函数，波形对称于原点。过程与（1）相似，如图（2）。 图（2） 3．已知系统的传输算子233 )(2+++=p p p p H ，画出并联结构的信号流图。 解：p p p p p p p p p p H 2 11 1122112233)(2+- ++=+-++=+++= 4．系统方程为)1()2(3)1(2)(-=-+-+n x n y n y n y ，画出信号流图。 解：23211)(E E E E H ++=
杭电通信系统课程设计实验报告
通信系统课程设计实验报告 姓名：田昕煜 学号： 13081405 班级：通信四班 班级号： 13083414
基于FSK调制的PC机通信电路设计 一、目的、内容与要求 目的: 掌握用FSK调制和解调实现数据通信的方法,掌握FSK调制和解调电路中相关模块的设计方法。初步体验从事通信产品研发的过程. 课程设计任务：设计并制作能实现全双工FSK调制解调器电路，掌握用Orcad Pspice、Protel99se进行系统设计及电路仿真。 要求：合理设计各个电路，尽量使仿真时的频率响应和其他参数达到设计要求。尽量选择符合标称值的元器件构成电路，正确完成电路调试。 二、总体方案设计 信号调制过程如下： 调制数据由信号发生器产生（电平为TTL，波特率不超过9600Baud），送入电平/幅度调整电路完成电平的变换，再经过锁相环（CD4046）,产生两个频率信号分别为30kHz和40kHz（发“1”时产生30kHz方波，发“0”时产生40kHz方波），再经过低通滤波器2，变成平滑的正弦波，最后通过线圈实现单端到差分信号的转换。 信号的解调过程如下： 首先经过带通滤波器1，滤除带外噪声，实现信号的提取。在本设计中FSK 信号的解调方式是过零检测法。所以还要经过比较器使正弦信号变成方波，再经过微分、整流电路和低通滤波器1实现信号的解调，最后经过比较器使解调信号成为TTL电平。在示波器上会看到接收数据和发送数据是一致的。 各主要电路模块作用： 电平/幅度调整电路：完成TTL电平到VCO控制电压的调整； VCO电路：在控制电压作用下，产生30KHz和40KHz方波； 低通2：把30KHz、40KHz方波滤成正弦波； 线圈：完成单端信号和差分信号的相互转换；
数字信号处理实验
实验六: 用FFT对信号作频谱分析 一、实验目的 1．了解双音多频信号的产生、检测、包括对双音多频信号进行DFT时的参数选择等。 2．初步了解数字信号处理在是集中的使用方法和重要性。 3．掌握matlab的开发环境。 二、实验原理与方法 1、引言 双音多频（Dual Tone Multi Frequency, DTMF）信号是音频电话中的拨号信号，由美国AT&T贝尔公司实验室研制，并用于电话网络中。这种信号制式具有很高的拨号速度，且容易自动监测识别，很快就代替了原有的用脉冲计数方式的拨号制式。这种双音多频信号制式不仅用在电话网络中，还可以用于传输十进制数据的其它通信系统中，用于电子邮件和银行系统中。这些系统中用户可以用电话发送DTMF信号选择语音菜单进行操作。DTMF信号系统是一个典型的小型信号处理系统，它要用数字方法产生模拟信号并进行传输，其中还用到了D/A变换器；在接收端用A/D变换器将其转换成数字信号，并进行数字信号处理与识别。为了系统的检测速度并降低成本，还开发一种特殊的DFT算法，称为戈泽尔(Goertzel)算法，这种算法既可以用硬件（专用芯片）实现，也可以用软件实现。下面首先介绍双音多频信号的产生方法和检测方法，包括戈泽尔算法，最后进行模拟实验。下面先介绍电话中的DTMF信号的组成。在电话中，数字0~9的中每一个都用两个不同的单音频传输，所用的8个频率分成高频带和低频带两组，低频带有四个频率：679Hz,770Hz,852Hz和941Hz；高频带也有四个频率：1209Hz,1336Hz,1477Hz和1633Hz.。每一个数字均由高、低频带中各一个频率构成，例如1用697Hz和1209Hz两个频率，信号用表示。这样8个频率形成16种不同的双频信号。具体号码以及符号对应的频率如表10.6.1所示。表中最后一列在电话中暂时未用。DTMF信号在电话中有两种作用，一个是用拨号信号去控制交换机接通被叫的用户电话机，另一个作用是控制电话机的各种动作，如播
郑州大学数字信号处理课程设计报告
实验一：基于DFT的数字谱分析以及可能出现的问题 一、实验目的： 1.进一步加深对DFT的基本性质的理解。 2.掌握在MATLAB环境下采用FFT函数编程实现DFT的语句用法。 3.学习用DFT进行谱分析的方法，了解DFT谱分析中出现的频谱泄露和栅栏效应现 象，以便在实际中正确应用DFT。 二、实验步骤: 1.复习DFT的定义、物理含义以及主要性质。 2.复习采用DFT进行谱分析可能出现的三个主要问题以及改善方案。 3.按实验内容要求，上机实验，编写程序。 4.通过观察分析实验结果，回答思考题，加深对DFT相关知识的理解。 三、上机实验内容: 1.编写程序产生下列信号供谱分析用： 离散信号： x1=R10(n) x2={1,2,3,4,4,3,2,1}，n=0,1,2,3,4,5,6,7 x3={4,3,2,1, 1,2,3,4}，n=0,1,2,3,4,5,6,7 连续信号： x4=sin(2πf1t)+sin(2πf2t) f1=100Hz, f2=120Hz，采样率fs=800Hz 2.对10点矩形信号x1分别进行10点、16点、64点和256点谱分析，要求256点 频谱画出连续幅度谱，10点、16点和64点频谱画出离散幅度谱，观察栅栏效应。 3.产生信号x2和x3分别进行8点、16点谱分析，画出离散幅度谱，观察两个信 号的时域关系和幅度谱的关系。 4.对双正弦信号x4以采样率fs=800Hz抽样，生成离散双正弦信号并画出连续波形； 对离散双正弦信号进行时域截断，截取样本数分别为1000、250、50。对不同样本的双正弦信号分别进行1024点谱分析，画出连续幅度谱，观察频谱泄露现象。
杭电信号与系统实验离散时间系统的时域分析
《信号、系统与信号处理实验I》 实验报告 实验名称：离散时间系统的时域分析 姓名： 学号： 专业：通信工程 实验时间 杭州电子科技大学 通信工程学院
一、实验目的 1.通过matlab 仿真一些简单的离散时间系统，并研究它们的时域特性。 2.掌握利用matlab 工具箱求解LTI 系统的单位冲激响应。 二、实验内容 1、离散时间系统的时域分析 1.1 线性与非线性系统 假定系统为y[n]-0.4y[n-1]=2.24x[n]+2.49x[n-1](2.9) 输入三个不同的输入序列x1[n]、x2[n]和，计算并求出相应的三个输出，并判断是否线性。x[n]=a x1[n]+b x2[n] clear all; n=0:40; a=2;b=-3; x1=cos(2*pi*0.1*n); x2=sin(2*pi*0.4*n); x=a*x1+b*x2; num=[2.24 2.49]; den=[1 -0.4]; y1=filter(num.den,x1); y2=filter(num.den,x2); y=filter(num.den,x); yt=a*y1+b*y2; d=y-yt;%计算差值输出d[n] subplot(3,1,1) stem(n,y); ylabel(‘振幅’); subplot(3,1,2) stem(n,yt); ylabel(‘振幅’); subplot(3,1,3) stem(n,d); ylabel(‘振幅’); title(‘差信号’) (1)假定另一个系统为y[n]=x[n]+3.2x[n-2]，修改以上程序，通过绘出的图形判断该系统是否线性系统。 1.2 时变与时不变系统 根据(2.9)的系统，产生两个不同的输入序列x[n]和x[n-D]，根据输出判断是否时不变系统。 clear all; n=0:40; a=2;b=-3; D=10; x=cos(2*pi*0.1*n); xd=[zeros(1,D) x]; num=[2.24 2.49]; den=[1 -0.4]; y=filter(num.den,x); yd=filter(num.den,xd); d=y-yd(1+D:41+D);%计算差值输d[n] subp lot(3,1,1) stem(n,y); ylabel(‘振幅’); title(‘输出y[n]’);grid; subplot(3,1,2)
数字信号处理实验——维纳滤波器设计..
实验一 维纳滤波 1. 实验内容 设计一个维纳滤波器： （1） 产生三组观测数据，首先根据()(1)()s n as n w n =-+产生信号()s n ，将其加噪，（信噪比分别为20,10,6dB dB dB ），得到观测数据123(),(),()x n x n x n 。 （2） 估计()i x n ，1,2,3i =的AR 模型参数。假设信号长度为L ，AR 模型阶数为N ，分析实验结果，并讨论改变L ，N 对实验结果的影响。 2. 实验原理 滤波目的是从被噪声污染的信号中分离出有用的信号来，最大限度地抑制噪声。对信号进行滤波的实质就是对信号进行估计。滤波问题就是设计一个线性滤波器，使得滤波器的输出信号()y n 是期望响应()s n 的一个估计值。下图就是观测信号的组成和信号滤波的一般模型。 观测信号()()()x n s n v n =+ 信号滤波的一般模型 维纳滤波解决从噪声中提取信号的滤波问题，并以估计的结果与真值之间的误差均方值最小作为最佳准则。它根据()()(),1, ,x n x n x n m --估计信号的当前 值，它的解以系统的系统函数()H z 或单位脉冲()h n 形式给出，这种系统常称为最佳线性滤波器。 维纳滤波器设计的任务就是选择()h n ，使其输出信号()y n 与期望信号()d n 误差的均方值最小。
假设滤波系统()h n 是一个线性时不变系统，它的()h n 和输入信号都是复函数，设 ()()()h n a n jb n =+ 0,1, n = 考虑系统的因果性，可得到滤波器的输出 ()()()()()0 *m y n h n x n h m x n m +∞ ===-∑ 0,1, n = 设期望信号()d n ，误差信号()e n 及其均方误差()2 E e n ???? 分别为 ()()()()()e n d n y n s n y n =-=- ()()()()()()22 2 0m E e n E d n y n E d n h m x n m ∞=?? ????=-=--????? ????? ∑ 要使均方误差为最小，需满足： ()() 2 0E e n h j ?????=? 整理得()()0E x n j e n *??-=??，等价于()()0E x n j e n * ??-=?? 上式说明，均方误差达到最小值的充要条件使误差信号与任一进入估计的输入信号正交，这就是正交性原理。 将()()0E x n j e n * ??-=??展开，得 ()()()()00m E x n k d n h m x m +∞ *** =????--=?? ???? ?∑ 整理得 ()()()0 dx xx m r k h m r m k +∞ *=-=-∑ 0,1,2, k = 等价于()()()()()0 dx xx xx m r k h m r k m h k r k +∞ ==-=*∑ 0,1,2, k = 此式称为维纳-霍夫（Wiener-Holf ）方程。解此方程可得到最优权系数 012,,, h h h ，此式是Wiener 滤波器的一般方程。 定义
《数字信号处理》课程实验题目
计电学院《数字信号处理》课程实验 适用专业：电子通信工程专业；实验学时：9 学时 一、实验的性质、任务和基本要求 (一)本实验课的性质、任务 数字信号处理课程实验是数字信号处理课程的有效的补充部分，通过实验，使学生巩固和加深数字信号处理的理论知识的理解和掌握，在实验过程中了解简单但是完整的数字信号处理的工程实现方法和流程。通过实践进一步加强学生独立分析问题和解决问题的能力、实际动手能力、综合设计及创新能力的培养。 （二）基本要求 掌握数字信号处理基本理论知识和滤波器设计及应用。 （三）实验选项
二、实验教学内容 实验一 1、实验目的和要求 1）加深理解时域采样定理、体会使用MATLAB的离散FT函数fft( )来解决涉及模拟信号的问题； 2）加深理解对带通信号的采样特性，学会采用MATLAB解决该问题； 3）加深理解在频率采样法中，过渡点对所设计滤波器特性的影响。 2、实验要求 1）提供MATLAB程序，画出每个步骤的曲线图； 2）写实验报告，包含有对所得结果进行分析和说明。 第一组：张毅雷凌峰白法聪覃昱滔刘强何新文 第二组：邓志强林盛勇李日胜黎少锋梁聪杨晨 实验二 1、实验目的和要求 （1）加深理解采用数字信号处理方法对模拟信号处理的过程、掌握使用MATLAB处理的方法；对一段音乐信号进行处理和输出；要求画出滤波前后语音信号时域波形、信号和滤波器的幅度频率特性曲线、相位频率特性曲线； （2）加深对截断效应的理解； （3）掌握使用MATLAB设计滤波器，并对语音信号处理的方法。对一段音乐信号进行处理和输出；要求画出滤波前后语音信号时域波形、信号和滤波器的幅度频率特性曲线、相位频率特性曲线。 2、实验要求 1）提供MATLAB程序，画出每个步骤的曲线图； 2）写实验报告，包含有对所得结果进行分析和说明。 第九组:汪涛张汉毅巫金敏张经中柳泽举 第六组：罗涛梁乐杰黄乃生 实验三 1、实验目的和要求 掌握采用MATLAB数字滤波器设计软件编制方法。软件要求在界面内有不同类型（高通低通带通带阻）滤波器的选择、或者只对低通滤波器采用不同方法设
信号与系统实验指导全部实验答案
00.51 1.52 2.53-2-1.5-1-0.50 0.5 1 1.52正弦信号 -20-15-10-505101520 -0.4 -0.200.20.40.60.81抽样信 号 -0.500.51 1.52 2.53 00.5 11.5 2矩形脉冲信 号 -1 012345 -0.5 0.5 1 1.5 单位跃阶信号 实验一 连续时间信号的MATLAB 表示 实验目的 1．掌握MATLAB 语言的基本操作，学习基本的编程功能； 2．掌握MATLAB 产生常用连续时间信号的编程方法； 3．观察并熟悉常用连续时间信号的波形和特性。 实验原理： 1. 连续信号MA TLAB 实现原理 从严格意义上讲，MATLAB 数值计算的方法并不能处理连续时间信号。然而，可用连续信号在等时间间隔点的取样值来近似表示连续信号，即当取样时间间隔足够小时，这些离散样值能够被MATLAB 处理，并且能较好地近似表示连续信号。 MATLAB 提供了大量生成基本信号的函数。比如常用的指数信号、正余弦信号等都是MATLAB 的内部函数。为了表示连续时间信号，需定义某一时间或自变量的范围和取样时间间隔，然后调用该函数计算这些点的函数值，最后画出其波形图。 实验内容：
0123-2-1 1 2实 部012 3 -1 12 虚 部0123012取 模0123 -50 5 相角00.2 0.40.6 0.81 -1-0.50 0.51 方波信号 实验编程： （1）t=0:0.01:3; K=2;a=-1.5;w=10; ft=K*exp((a+i*w)*t); A=real(ft); B=imag(ft); C=abs(ft); D=angle(ft); subplot(2,2,1),plot(t,A),grid on;title('实部'); subplot(2,2,2),plot(t,B),grid on;title('虚部'); subplot(2,2,3),plot(t,C),grid on;title('取模'); subplot(2,2,4),plot(t,D),grid on;title('相角'); （2） t=0:0.001:3; y=square(2*pi*10*t,30); plot(t,y); axis([0,1,-1,1]); title('方波信号');