南通大学电子信息学院信号与系统课程设计报告
班级：
姓名：
学号：
指导老师：
设计时间：
2014—2015学年第一学期
一、连续信号的时域分析
二、
1. 信号的产生 (1)阶跃函数
function [t,y1]=jieyue(t1,t2,t0) dt=0.01; ttt=t1:dt:t0-dt; tt=t0:dt:t2; t=t1:dt:t2; n=length(ttt); nn=length(tt); u=zeros(1,n); uu=ones(1,nn); y1=[u,uu]; return 冲激函数
function [t,y2]=chongji(t1,t2,t0) dt=0.01; t=t1:dt:t2; n=length(t); y2(1:n)=0;
y2(1,(t0-t1)/dt+1)=1/dt;
（2）调用上述函数产生信号）2-t ε（，）（4-t δ，-t e ）（t ε，-6s ≤t ≤6s,并画出波形。
Command Window
subplot(3,1,1);
[t1,y1]=jieyue(-6,6,2);
stairs(t1,y1);
axis([-6 6 0 1.5]);
subplot(3,1,2);
[t2,y2]=chongji(-6,6,4);plot(t2,y2); subplot(3,1,3);
[t3,y3]=jieyue(-6,6,0);
y3=exp(-(t3)).*y3;plot(t3,y3);
波形如下图所示：
（3）根据f(t)画出f(2t)和f(1-0.5t)的波形
t=-3:0.01:3;
y=tripuls(t,4,0.6);
subplot(3,1,1);
plot(t,y);
title('f(t)');
xlabel('(a)');
y1=tripuls(2*t,4,0.6); subplot(3,1,2);
plot(t,y1);
title('f(2t)');
xlabel('(b)');
t1=2-2*t;
y2=tripuls(1-0.5*t1,4,0.6); subplot(3,1,3);
plot(t1,y2);
title('f(1-0.5*t)');
xlabel('(c)');
得到波形如下图所示：
已知信号f（t）=（1+t/2）*（u(t+2)-u(t-2)）,用matlab求f（t+2），f(t-2),f(-t).f(2t),-f(t),并绘出时域波形。
Command Window
syms t;
f=sym(‘(t/2+1)*(Heaviside(t+2)- Heaviside(t-2))’);
subplot(2,3,1);ezplot(f,[-3,3]);
y1=subs(f,t,t+2);
subplot(2,3,2);ezplot(y1,[-5,1]);
y2=subs(f,t,t-2);
subplot(2,3,3);ezplot(y2,[-1,5]);
y3=subs(f,t,-t);
subplot(2,3,4);ezplot(y3,[-3,3]);
y4=subs(f,t,2*t);
subplot(2,3,5);ezplot(y4,[-2,2]);
y5=-f;
subplot(2,3,6);ezplot(y5,[-3,3]);
仿真结果如下图所示：
2、信号的卷积
（1）f （t ）=f1(t)*f2(t)
function [k,f]=myconv(f1,f2,k1,k2,p) f=conv(f1,f2)*p; k0=k1(1)+k2(1);
k3=length(f1)+length(f2)-2; k=k0:p:k0+k3*p;
（2）求x1=[]）（）（2-t -t e t εε-和x2=δ(t+3)+δ(t-3)的卷积x1(t)*x2(t),并验证卷积的性质。
Command Window
[t1,f11]=jieyue(-8,8,0);[t1,f12]=jieyue(-8,8,2); f1=f11-f12; x1=exp(-t1).*f1;
[t2,f21]=chongji(-8,8,-3);[t2,f22]=chongji(-8,8,3);
x2=f21+f22;
subplot(3,1,1);plot(t1,x1); subplot(3,1,2);plot(t2,x2); [t3,f]=myconv(x1,x2,t1,t2,0.01) subplot(3,1,3);plot(t3,f); 仿真结果：
(3)已知???≤≤-<<=32,12
0,1t f 1
t t ）（,t 22e t f -=）（,0<=t<=10的卷积f(t)=f1(t)*f2(t)的时域波形图。
Command Window
t11=0; t12=3; t21=0; t22=10;
t1=t11:0.001:t12;
ft1=-sign(t1-2);
t2=t21:0.001:t22;
ft2=exp(-2*t2);
t=t11+t21:0.001:t12+t22; ft=conv(ft1,ft2);
ft=ft*0.001;
subplot(2,2,1);
plot(t1,ft1);
title('f1(t)');
subplot(2,2,2);
plot(t2,ft2);
title('f2(t)');
subplot(2,2,3);
plot(t,ft);
title('f1(t)*f2(t)');
仿真结果：
例：已知两个信号f1(t)=u(t-1)-u(t-2),f2(t)=u(t)-u(t-1),求f(t)=f1(t)*f2(t)的时域波形图。
Command Window
t1=1;
t2=2;
t3=0;
t4=1;
t=0:T:t2+t4;
x1=ones(size(t)).*((t>t1)-(t>t2));
x2=ones(size(t)).*((t>t3)-(t>t4));
y=conv(x1,x2)*T;
subplot(3,1,1),plot(t,x1);
ylabel('x1(t)');
subplot(3,1,2),plot(t,x2);
ylabel('x2(t)');
subplot(3,1,3),plot(t,y(1:(t2+t4)/T+1));
ylabel('y(t)=x1*x2');
xlable('---t/s');
仿真结果：
3.连续系统的响应
（1）已知系统的微分方程为y＇＇(t)+5y＇(t)+6y(t)=f(t),求系统的单位冲激响应h（t）和单位阶跃响应g（t）；
CommandWindow
b=[1];a=[1,5,6];
subplot(2,1,1);impulse(b,a);
subplot(2,1,2);step(b,a);
仿真结果：
（2）对于上述系统，请画出激励f(t)分别为）（t e -t ε、）（t e -2t ε、）（）（t t 2cos ε、
）（t t 2ε时系统的零状态响应的波形，分析与理论计算的结果是否相符。
Command Window
b=[1]; a=[1,5,6]; t=0:0.1:10; f1=exp(-t); f2=cos(2*t); f3=t.^2; f4=(exp(-2*t));
subplot(2,2,1);lsim(b,a,f1,t);grid;title('f(t)=exp(-t)'); subplot(2,2,2);lsim(b,a,f2,t);grid;title('f(t)=cos(2*t)'); subplot(2,2,3);lsim(b,a,f3,t);grid;title('f(t)=t.^2'); subplot(2,2,4);lsim(b,a,f4,t);grid;title('f(t)=exp(-2*t)');
仿真结果：
（3）已知系统的微分方程为y ’’(t)+3y ’(t)+2y(t)=3f(t)+f ’(t),初始条件：y(0+)=1,y ’(0-)=2,求：
1）系统的零输入响应）
（t y x ; 2)激励为f(t)=）（t e -3t 时，系统的零状态响应）（t y f 和全响应y(t)，分析与理论计算的结果是否相符。
二．离散系统的时域分析
（1）已知离散系统的差分方程为：
y(k)+1/3y(k-2)=1/6f(k)+1/2f(k-1)+1/2f(k-2)+1/6f(k-3),
画出单位序列响应、单位阶跃响应的波形。
Command Window
k=0:1:32;
a=[1,0,1/3,0];
b=[1/6,1/2,1/2,1/6];
subplot(2,1,1);hk=impz(b,a,k);stem(k,hk,'k');
subplot(2,1,2);gk=dstep(b,a,k);stem(k,gk,'r');
仿真结果：
（2）已知离散系统的差分方程为：
y(k)-1.5y(k-1)+0.5y(k-2)=f(k)满足初始条件y（-1）=4,y(-2)=10,
用filtic和filter子函数求系统的激励为f(k)=(0.25) ^k*u(k)时的零输入、零状态以及完全响应。
Command Window
a=[1,-1.5,0.5];
b=[1];
N=20;
k=0:N-1;
f=0.25.^k;
f0=zeros(1,N);
y01=[4,10];
fi=filtic(b,a,y01);
y0=filter(b,a,f0,fi);
fi0=filtic(b,a,0);
y1=filter(b,a,f,fi0);
y=filter(b,a,f,fi);
y2=((1/2).^k+1/3*(1/4).^k+2/3).*ones(1,N);
subplot(2,3,1);stem(k,f);title('输入信号f(k)');
subplot(2,3,2);stem(k,y0);title('零输入响应');
subplot(2,3,3);stem(k,y1);title('零状态响应');
subplot(2,3,4);stem(k,y);title('用filter求完全响应');
subplot(2,3,5);stem(k,y2);title('用公式求完全响应');
仿真结果：
（3）已知离散系统的差分方程为y （k ）+3y(k-1)+2y(k-2)=f(k),
（a ）画出单位阶跃响应、单位序列响应的波形；
（b ）画出激励）（）（k 2k f k ε=时的系统零状态响应波形。
Command Window
n=0:1:100; a=[1,3,2]; b=[1,0,0];
subplot(3,1,1);hn=impz(b,a,n);stem(n,hn,'k'); subplot(3,1,2);gn=dstep(b,a,n);stem(n,gn,'r'); f2=2.^n;
subplot(3,1,3);y2=filter(b,a,f2);plot(n,y2,'.'); 仿真结果：
三．信号与系统的频域分析
1、门函数的频谱
，画出τ=10秒时门函数在-2(rad/s）<=w<=2（1）产生宽度为τ的门函数）
g
（t
τ
π（rad/s）频率范围内频谱，记录最大值，观察第一过零点位置：Command Window
dt=0.1;
N=500;
door_width=10;
tao=door_width/2;
t1=-(N-dt):dt:-tao+dt;
t2=-tao:dt:tao;
t3=tao+dt:dt:N-dt;
t=[t1,t2,t3];
f=[zeros(1,length(t1)),ones(1,length(t2)),zeros(1,length(t3))];
w=-2*pi:0.1:2*pi;
F1=f*exp(-j*t'*w)*dt;
plot(w,real(F1));grid;
仿真结果：
（2）改变τ的值为5秒和20秒，重复（1）步骤；
5=τ 20=τ
2、傅里叶变换及其性质的验证 （1）求f(t)=t
2e
-的傅里叶变换，并绘出f(t)及其傅里叶变换的波形图。
Command Window
syms t ;
ft=exp(-2*abs(t)); Fw=fourier(ft);
subplot(2,1,1);ezplot(ft); subplot(2,1,2);ezplot(Fw); 仿真结果：
（2）以门函数）（t g 为分析对象，验证傅里叶变换的时移性质、频移性质； 时移性质：
Command Window
dt=0.1; N=500; door__width=10; tao=door__width/2; t1=-(N-1):dt:-tao+1; t2=-tao:dt:tao; t3=tao-1:dt:N-1;
课程设计报告 课程名称信号与系统课程设计指导教师 设计起止日期 学院信息与通信工程 专业电子信息工程 学生 班级/学号 成绩 指导老师签字
目录 1、课程设计目的 (1) 2、课程设计要求 (1) 3、课程设计任务 (1) 4、课程设计容 (1) 5、总结 (11) 参考文献 (12) 附录 (12)
1、课程设计目的 “信号与系统”是一门重要的专业基础课，MATLAB作为信号处理强有力的计算和分析工具是电子信息工程技术人员常用的重要工具之一。本课程设计基于MATLAB完成信号与系统综合设计实验，以提高学生的综合应用知识能力为目标，是“信号与系统”课程在实践教学环节上的必要补充。通过课设综合设计实验，激发学生理论课程学习兴趣，提高分析问题和解决问题的能力。 2、课程设计要求 （1）运用MATLAB编程得到简单信号、简单信号运算、复杂信号的频域响应图； （2）通过对线性时不变系统的输入、输出信号的时域和频域的分析，了解线性时不变系统的特性，同时加深对信号频谱的理解。 3、课程设计任务 （1）根据设计题目的要求，熟悉相关容的理论基础，理清程序设计的措施和步骤； （2）根据设计题目的要求，提出各目标的实施思路、方法和步骤； （3）根据相关步骤完成MATLAB程序设计，所编程序应能完整实现设计题目的要求； （4）调试程序，分析相关理论； （5）编写设计报告。 4、课程设计容 （一）基本部分 （1）信号的时频分析 任意给定单频周期信号的振幅、频率和初相，要求准确计算出其幅度谱，并准确画出时域和频域波形，正确显示时间和频率。 设计思路： 首先给出横坐标，即时间，根据设定的信号的振幅、频率和初相，写出时域波形的表达式；然后对时域波形信号进行傅里叶变化，得到频域波形；最后使用plot函数绘制各个响应图。 源程序： clc; clear; close all; Fs =128; % 采样频率 T = 1/Fs; % 采样周期 N = 600; % 采样点数 t = (0:N-1)*T; % 时间，单位：S x=2*cos(5*2*pi*t);
信号分析课程设计 信号系统的时域分析 编程实现的卷积积分或卷积和 一、课程设计题目： 基于 MATLAB 的连续时间LTI 系统的时域分析 二、基本要求： ① 掌握连续时不变信号处理的基本概念、基本理论和基本方法； ② 学会 MATLAB 的使用，掌握 MATLAB 的程序设计方法； ③ 学会用 MATLAB 对信号进行分析和处理； ④ 编程实现卷积积分或卷积和，零输入响应，零状态响应； ⑤ 撰写课程设计论文，用信号处理基本理论分析结果。 三、设计方法与步骤： 一般的连续时间系统分析有以下几个步骤: ①求解系统的零输入响应; ②求解系统的零状态响应; ③求解系统的全响应; ④分析系统的卷积；⑤画出它们的图形. 下面以具体的微分方程为例说明利用MATLAB 软件分析系统的具体方法. 1．连续时间系统的零输入响应 描述n 阶线性时不变（LTI ）连续系统的微分方程为： 已知y 及各阶导数的初始值为y(0),y (1)(0)，… y (n-1)(0)， 求系统的零输入响应。 建模 当LIT 系统的输入为零时，其零输入响应为微分方程的其次解（即令微分方程的等号右端为零），其形式为（设特征根均为单根） 其中p 1,p 2,…,p n 是特征方程a 1λ n +a 2λn-1+…+a n λ+a n =0的根，它们可以 用root(a)语句求得。各系数 由y 及其各阶导数的初始值来确定。对此有 1121111n n m n n m m n n m d y d y dy d u du a a a a y b b b u dt dt dt dt dt -++-++?????++=+????++1212()n p t p t p t n y t C e C e C e =++????+120n C C C y ++????+=11220 n n p C p C p C Dy ++????+=
《数字信号处理》课程设计报告 设计题目： IIR滤波器的设计 专业： 班级： 姓名： 学号： 指导教师： 2010年月日
1、设计目的 1、掌握IIR 滤波器的参数选择及设计方法； 2、掌握IIR 滤波器的应用方法及应用效果； 3、提高Matlab 下的程序设计能力及综合应用能力。 4、了解语音信号的特点。 2、设计任务 1、学习并掌握课程设计实验平台的使用，了解实验平台的程序设计方法； 2、录制并观察一段语音信号的波形及频谱，确定滤波器的技术指标； 3、根据指标设计一个IIR 滤波器，得到该滤波器的系统响应和差分方程，并根据差分方程将所设计的滤波器应用于实验平台，编写相关的Matlab 程序； 4、使用实验平台处理语音信号，记录结果并进行分析。 3、设计内容 3.1设计步骤 1、学习使用实验平台，参见附录1。 2、使用录音机录制一段语音，保存为wav 格式，录音参数为：采样频率8000Hz、16bit、单声道、PCM 编码，如图1 所示。 图1 录音格式设置 在实验平台上打开此录音文件，观察并记录其波形及频谱（可以选择一段较为稳定的语音波形进行记录）。 3、根据信号的频谱确定滤波器的参数：通带截止频率Fp、通带衰减Rp、阻带截止频率Fs、阻带衰减Rs。 4、根据技术指标使用matlab 设计IIR 滤波器，得到系统函数及差分方程，并记录得到系统函数及差分方程，并记录其幅频响应图形和相频响应图形。要求设计 第 1页出的滤波器的阶数小于7，如果不能达到要求，需要调整技术指标。 5、记录滤波器的幅频响应和系统函数。在matlab 中，系统函数的表示公式为：
因此，必须记录系数向量a 和b。系数向量a 和b 的可以在Matlab 的工作空间（WorkSpace）中查看。 6、根据滤波器的系统函数推导出滤波器的差分方程。 7、将设计的滤波器应用到实验平台上。根据设计的滤波器的差分方程在实验平台下编写信号处理程序。根据运行结果记录处理前后的幅频响应的变化情况，并试听处理前后声音的变化，将结果记录，写入设计报告。 3.2实验程序 （1）Rs=40; Fs=1400; Rp=0.7; Fp=450; fs=8000; Wp=2*pi*Fp;Ws=2*pi*Fs; [N,Wn]=buttord(Wp,Ws,Rp,Rs,'s'); [b1,a1]=butter(N,Wn,'s'); [b,a]=bilinear(b1,a1,fs); [H,W]=freqz(b,a); figure; subplot(2,1,1);plot(W*fs/(2*pi),abs(H));grid on;title('频率响应'); xlabel('频率');ylabel('幅值');、 subplot(2,1,2); plot(W,angle(H));grid on;title('频率响应'); xlabel('相位(rad)');ylabel('相频特性'); 3.3实验结果（如图）： N =5 Wn=6.2987e+003 第 2页
抽样定理的应用 摘要 抽样定理表示为若频带宽度有限的，要从抽样信号中无失真地恢复原信号，抽样频率应大于2倍信号最高频率。抽样频率小于2倍频谱最高频率时，信号的频谱有混叠。抽样频率大于2倍频谱最高频率时，信号的频谱无混叠。 语音信号处理是研究用数字信号处理技术和语音学知识对语音 信号进行处理的新兴学科，是目前发展最为迅速的学科之一，通过语音传递信息是人类最重要，最有效，最常用和最方便的交换信息手段，所以对其的研究更显得尤为重要。 Matlab语言是一种数据分析和处理功能十分强大的计算机应用 软件，它可以将声音文件变换成离散的数据文件，然后用起强大的矩阵运算能力处理数据。这为我们的本次设计提供了强大并良好的环境！ 本设计要求通过利用matlab对模拟信号和语音信号进行抽样，通过傅里叶变换转换到频域，观察波形并进行分析。 关键词：抽样Matlab
目录 一、设计目的： (2) 二、设计原理： (2) 1、抽样定理 (2) 2、MATLAB简介 (2) 3、语音信号 (3) 4、Stem函数绘图 (3) 三、设计内容： (4) 1、已知g1(t)=cos(6πt),g2(t)=cos(14πt),g3(t)=cos(26πt),以抽样频率 fsam=10Hz对上述三个信号进行抽样。在同一张图上画出g1(t)，g2(t)，g3(t)及其抽样点,对所得结果进行讨论。 (4) 2、选取三段不同的语音信号，并选取适合的同一抽样频率对其进 行抽样，画出抽样前后的图形，并进行比较，播放抽样前后的语音。 (6) 3、选取合适的点数，对抽样后的三段语音信号分别做DFT，画图 并比较。 (10) 四、总结 (12) 五、参考文献 (13)
课程设计报告 课程名称数字信号处理 课题名称数字滤波器设计及在语音信号分析中的应用 专业通信工程 班级1281 学号201213120101 姓名杨俊 指导教师彭祯韩宁 2014年12月5日
湖南工程学院 课程设计任务书 课程名称数字信号处理 课题数字滤波器设计 及在语音信号分析中的应用专业班级通信工程1281班 学生姓名杨俊 学号201213120101 指导老师彭祯韩宁 审批 任务书下达日期2014 年12月5日 任务完成日期2014 年12月13日
《数字信号处理》课程设计任务书 一、课程设计的性质与目的 《数字信号处理》课程是通信专业的一门重要专业基础课，是信息的数字化处理、存储和应用的基础。通过该课程的课程设计实践，使学生对信号与信息的采集、处理、传输、显示、存储、分析和应用等有一个系统的掌握和理解；巩固和运用在《数字信号处理》课程中所学的理论知识和实验技能，掌握数字信号处理的基础理论和处理方法，提高分析和解决信号与信息处理相关问题的能力，为以后的工作和学习打下基础。 数字滤波器是一种用来过滤时间离散信号的数字系统，通过对抽样数据进行数学处理来达到频域滤波的目的。根据其单位冲激响应函数的时域特性可分为两类：无限冲激响应（IIR）滤波器和有限冲激响应（FIR）滤波器。 二、课程设计题目 题目1：数字滤波器设计及在语音信号分析中的应用。 1、设计步骤： （1）语音信号采集 录制一段课程设计学生的语音信号并保存为文件，要求长度不小于10秒，并对录制的信号进行采样；录制时可以使用Windows自带的录音机，或者使用其它专业的录音软件，录制时需要配备录音硬件（如麦克风），为便于比较，需要在安静、干扰小的环境下录音。 然后在Matlab软件平台下，利用函数wavread对语音信号进行采样，记住采样频率和采样点数。 （2）语音信号分析 使用MATLAB绘出采样后的语音信号的时域波形和频谱图。根据频谱图求出其带宽，并说明语音信号的采样频率不能低于多少赫兹。 （3）含噪语音信号合成 在MATLAB软件平台下，给原始的语音信号叠加上噪声，噪声类型分为如下几种：①白
语音信号的数字滤波 ——利用双线性变换法实现IIR数字滤波器的设计一．课程设计的目的 通过对常用数字滤波器的设计和实现，掌握数字信号处理的工作原理及设计方法；熟悉用双线性变换法设计 IIR 数字滤波器的原理与方法，掌握利用数字滤波器对信号进行滤波的方法，掌握数字滤波器的计算机仿真方法，并能够对设计结果加以分析。 二．设计方案论证 1.IIR数字滤波器设计方法 IIR数字滤波器是一种离散时间系统，其系统函数为 假设M≤N，当M＞N时,系统函数可以看作一个IIR的子系统和一个(M-N)的FIR子系统的级联。IIR数字滤波器的设计实际上是求解滤波器的系数和，它 是数学上的一种逼近问题，即在规定意义上（通常采用最小均方误差准则）去逼近系统的特性。如果在S平面上去逼近，就得到模拟滤波器；如果在z平面上去逼近，就得到数字滤波器。 2.用双线性变换法设计IIR数字滤波器 脉冲响应不变法的主要缺点是产生频率响应的混叠失真。这是因为从S平面到Ｚ平面是多值的映射关系所造成的。为了克服这一缺点，可以采用非线性频率压缩方法，将整个频率轴上的频率范围压缩到-π/T～π/T之间，再用z=e sT转换 平面的-π/T～π到Z平面上。也就是说，第一步先将整个S平面压缩映射到S 1 /T一条横带里；第二步再通过标准变换关系z=e s1T将此横带变换到整个Z平面上去。这样就使S平面与Z平面建立了一一对应的单值关系，消除了多值变换性，也就消除了频谱混叠现象，映射关系如图1所示。 图1双线性变换的映射关系 为了将S平面的整个虚轴jΩ压缩到S1平面jΩ1轴上的-π/T到π/T段上，可以通过以下的正切变换实现
《信号处理与系统分析》课程设计指导书 南通大学电子信息学院 信息工程系 2013年5月
前言 《信号处理与系统》是南通大学杏林学院通信信息类专业的一门专业基础课程，其理论性强，是其它后续专业课程的基础。 开设该课程设计的重要意义在于：首先，从帮助学习《信号处理与系统》课程的角度讲，学生借助于计算机，通过系统仿真，可以对信号以及线性系统的分析方法有一个更深入、更直接的认识，巩固理解一些抽象的知识，从而掌握《信号处理与系统》课程中的主要理论与基本原理；其次，从长远意义讲，学生掌握了数值分析软件Matlab的应用方法，为后续专业课的学习打下了坚实的基础；另外实践环节使学生在综合使用现代电子信息技术和手段进行设计、制作和创新方面的能力有所提高，为以后走上工作岗位从事信号分析和系统分析创造了必备的条件。 本课程设计时间为两周，学生根据课程设计指导书进行练习，考核成绩将根据学生出席情况及学习态度、课程设计报告完成情况、最后检查情况综合给出。 编者：李蕴华 2013年5月
课程设计的要求 一、熟练掌握Matlab语言的编程方法； 二、熟悉用于《信号处理与系统分析》的Matlab主要函数的应用； 三、记录实验结果（包括波形和数据），撰写课程设计报告。 主要内容及步骤 一、连续系统的时域分析 1、信号的产生 （1）编写生成连续阶跃信号u(t-t0)及冲激信号δ(t-t0)的函数：function [x,y]=jieyue(t1,t2,t0) 和 function [x,y] =chongji(t1,t2,t0)，信号的时间变量取值区间为t1～t2，t0为阶跃点或冲激点处的时间，x为信号的时间向量，y为相应的信号值向量。（提示：冲激信号可以用时间宽度为dt、高度为1/dt的矩形脉冲来近似表示。当dt很小时，矩形脉冲信号可近似认为是冲激信号。在对该矩形脉冲信号采点取样后（设取样间隔为dt），信号值y的第1+(t0- t1)/dt个元素的值为1/dt，其余元素的值为0。） 参考程序: function [t,y1]=jieyue(t1,t2,t0) dt=0.01; ttt=t1:dt:t0-dt; tt=t0:dt:t2; t=t1:dt:t2; n=length(ttt); nn=length(tt); u=zeros(1,n); uu=ones(1,nn); y1=[u,uu]; return
目录 一、课题名称 (2) 二、内容摘要 (2) 三、设计目的 (2) 四、设计内容及要求 (2) 五、系统方案设计 (3) 六、电路设计及原理分析 (4) 七、电路仿真结果 (7) 八、硬件设计及焊接测试 (8) 九、故障的原因分析及解决方案 (11) 十、课程设计总结及心得体会 (12)
一、课题名称：函数信号发生器的设计 二、内容摘要： 函数信号发生器作为一种常用的信号源,是现代测试领域内应用最为广泛的通用仪器之一。在研制、生产、测试和维修各种电子元件、部件以及整机设备时，都要有信号源，由它产生不同频率不同波形的电压、电流信号并加到被测器件或设备上，用其他仪器观察、测量被测仪器的输出响应，以分析确定它们的性能参数。信号发生器是电子测量领域中最基本、应用最广泛的一类电子仪器。它可以产生多种波形信号,如正弦波,三角波,方波等,因而此次课程设计旨在运用模拟电子技术知识来制作一个能同时输出正弦波、方波、三角波的信号发生器。 三、设计目的： 1、进一步掌握模拟电子技术知识的理论知识，培养工程设计能力和综合分析能力、解决问题的能力。 2、基本掌握常用电子电路的一般设计方法，提高电子电路的设计和实验能力。 3、学会运用Multisim仿真软件对所做出来的理论设计进行仿真测试，并能进一步解决出现的基本问题，不断完善设计。 4、掌握常用元器件的识别和测试，熟悉万用表等常用仪表，了解电路调试的基本方法，提高实际电路的分析操作能力。 5、在仿真结果的基础上，实现实际电路。 四、设计内容及要求： 1、要求完成原理设计并通过Multisim软件仿真部分 （1）RC桥式正弦波产生电路，频率分别为300Hz、1KHz、10KHz、500KHz，输出幅值300mV~5V可调、负载1KΩ。 （2）占空比可调的矩形波电路，频率3KHz，占空比可调范围10%~90%，输出幅值3V、负载1KΩ。 （3）占空比可调的三角波电路，频率1KHz，占空比可调范围10%~90%，输出幅值3V、负载1KΩ。 （4）多用信号源产生电路，分别产生正弦波、方波、三角波，频率范围
数字信号处理课程设计 姓名：刘倩 学号：201014407 专业：信息与计算科学 实验一：常见离散信号产生和实现 一、实验目的： 1、加深对常用离散信号的理解； 2、掌握matlab 中一些基本函数的建立方法。 二、实验原理： 1.单位抽样序列 在MATLAB 中可以利用zeros()函数实现。 如果)(n δ在时间轴上延迟了k 个单位，得到)(k n -δ即： 2．单位阶越序列 在MATLAB 中可以利用ones()函数实现。 3．正弦序列 在MATLAB 中 4．复指数序列 在MATLAB 中 5．指数序列 在MATLAB 中
实验内容：由周期为10的正弦函数生成周期为20的余弦函数。 实验代码： n=0:30; y=sin(0.2*pi*n+pi/2); y1=sin(0.1*pi*n+pi/2); subplot(121) stem(n,y); xlabel ('时间序列n');ylabel('振幅');title('正弦函数序列y=sin(0.2*pi*n+pi/2)'); subplot(122) stem(n,y1); xlabel ('时间序列n');ylabel('振幅'); title('正弦函数序列y=sin(0.2*pi*n+pi/2)'); 实验结果： 实验二：离散系统的时域分析 实验目的：加深对离散系统的差分方程、冲激响应和卷积分析方法的理解。实验原理：离散系统 其输入、输出关系可用以下差分方程描述： 输入信号分解为冲激信号， 记系统单位冲激响应 则系统响应为如下的卷积计算式：
当N k d k ,...2,1,0==时，h[n]是有限长度的（n ：[0，M]），称系统为FIR 系统；反之，称系统为IIR 系统。 在MATLAB 中，可以用函数y=filter(p,d,x)实现差分方程的仿真，也可以用函数 y=conv(x,h)计算卷积，用y=impz(p,d,N)求系统的冲激响应。 实验内容：用MATLAB 计算全解 当n>=0时，求用系数差分方程y[n]+y[n-1]-6y[n-2]=x[n]描述的一个离散时间系统对阶跃输入x[n]=8μ[n]的全解。 实验代码： n=0:7; >> [y,sf]=filter(1,[1 1 -6],8*ones(1,8),[-7 6]); >> y1(n+1)=-1.8*(-3).^n+4.8*(2).^n-2; >> subplot(121) >> stem(n,y); >> title('由fliter 函数计算结果'); >> subplot(122) >> stem(n,y1); >> title('准确结果'); 实验结果： 结果分析：有图可得由fliter 函数得出的结果与计算出的准确结果完全一致。 实验三FFT 算法的应用
信号与线性系统 课程设计 题目 学号姓名学号姓名学号姓名学号姓名 院系 年级 专业 日期
摘要 利用MATLAB对周期为T0的方波信号进行傅里叶级数展开，并绘制离散幅度谱和不同次谐波叠加后的图形。通过观察绘制的各个图像，加深对傅立叶变换和信号的分解与合成的理解。 Abstract Expanded the square wave signal with periodic T0 to Fourier series by MATLAB , and drew the discrete spectrum and plot the patterns after different sub harmonics are superimposed. Through the observation of each image, deepen the understanding of the Fourier transform and signal decomposition and synthesis. 关键词：矩形信号傅里叶级数谐波叠加分解与合成 Keywords: Squarewave signal.Fourier series. Harmonic superposition. Decomposition and synthesis 一、设计目的和要求 本设计主要利用MATLAB绘制信号的离散幅度谱和各次谐波叠加后的波形，通过观察谐波展开次数增加后的波形，进一步掌握信号分解与合成的原理。 培养运用所学知识分析解决问题的能力。 掌握用MATLAB实现通信系统仿真实验的能力。 这里要做一个信号的分解与合成的仿真系统，利用matlab软件的仿真模拟能力来体现信号的分解与合成过程中出现的情况。 MATLAB（矩阵实验室）是MATrix LABoratory的缩写，是一款由美国
随机信号分析课程设计 HEN system office room 【HEN16H-HENS2AHENS8Q8-HENH1688】
随机信号分析课程设计 一、题目： 设计一个抑制载波的复AM 信号，载波为40MHz ，接收带宽10MHz ，调制信号带宽50KHz ，加入高斯白噪声，带内信噪比10dB ： 1．1画出加噪后信号时域波形； 画出功率谱密度； 画出其同相、正交分量的功率谱； 统计方法画出包络概率密度。 二、问题分析： 画出加噪后信号时域波形： 首先，由已知条件先采样产生抑制载波的实AM 离散信号sr ，经过Hilbert 变换求得其解析信号s0，并经过低通滤波器，截止频率fs 限制在接收带宽下，加入噪声v ，得到此时的复AM 信号s ，再画出此时得到的复AM 信号时域波形。 画出功率谱密度： 将信号s 进行fft 变换后求得其功率谱密度，画出图形。 画出其同相、正交分量的功率谱： 将信号s 分解为正交分量)(cos )()(t t A t A s φ=和同相分量 )(sin )(t t A A c φ=，进行fft 变换得到功率谱密度，画出图形。 三、程序代码： f0=4*10^7;%载波信号频率40MHz f1=5*10^4;%调制信号频率50kHz fs=1*10^7;%接收机带宽采样频率10Mhz N=40001;%采样点数 %t=(0:1:N-1)/fs; f =10*f0; %画图范围设置 t0 = 5/f1;
t = 0:1/f:t0; k=1; sr=k*cos(2*pi*f1*t).*cos(2*pi*f0*t);%实am信号 % figure(2) % plot(t,sr) s0=hilbert(sr);%复am信号 h=sin(fs*t)/(pi*t); s0=conv(s0,h); am=max(abs(s0)); % % -------加噪方案（由加噪后信噪比确定高斯白噪声）----- snr=10; %设定加入白噪声后的信噪比为10db(均值为0) Pv=(am/(10^(snr/20)))^2;%噪声方差 % % -------------------------------------------------- % % % ---------加噪声------------- v=rand(1,N); v=v*sqrt(Pv);%白噪声 s=s0+v;%信号加噪声 % % ----------信号画图------------- figure(1) subplot(2,1,1),plot(t,s0); axis([*10^(-4) *10^-4 -10 10]) title('原始信号') subplot(2,1,2),plot(t,s); title(['加噪信号信噪比= ',num2str(snr),' dB. 噪声方差= ',num2str(Pv)]) axis([*10^(-4) *10^-4 -10 10]) %%----------画功率谱--------------- s1=detrend(s);%去趋势 ffs=abs(fft(s1)); theta=angle(s1)-2*pi*f0*t; a=abs(s1); ffs=ffs.*conj(ffs)*2/N;%频谱 %ffs=ffs.^2;%功率谱 figure(2) plot(ffs(1:N/2)); title('加噪信号功率谱') axis([3500 4500 0 4*10^4]) xlabel('*10^4') %%------------画正交同相分量功率谱---------- ac=s.*cos(2*pi*f0*t)-j*(hilbert(s)-s).*sin(2*pi*f0*t); as=-s.*sin(2*pi*f0*t)-j*(hilbert(s)-s).*cos(2*pi*f0*t); as1=detrend(as); ffas=abs(fft(as1)); ffas=(abs(ffas)).^2*2/N; ac1=detrend(ac); ffac=abs(fft(ac1)); ffac=(abs(ffac)).^2*2/N;
《数字信号》课程设计报 告 学院：信息科学与工程 专业班级：通信1201
一、 目的与要求 是使学生通过上机使用Matlab 工具进行数字信号处理技术的仿真练习，加深对《信号分析与处理（自）》课程所学基本理论和概念的理解，培养学生应用Matlab 等工具进行数字信号处理的基本技能和实践能力，为工程应用打下良好基础。 二、 主要内容 1．了解Matlab 基本使用方法，掌握Matlab 数字信号处理的基本编程技术。掌握数字信号的基本概念。 2.用Matlab 生成几种典型数字信号（正弦信号、矩形信号、三角波信号等），并做幅频特性分析 2．Matlab 编程实现典型离散信号（正弦信号、矩形信号、三角信号）的离散傅立叶变换，显示时域信号和频谱图形（幅值谱和相位谱）；以正弦周期信号为例，观察讨论基本概念（混叠、泄漏、整周期截取、频率分辨率等）。 3．设计任意数字滤波器，并对某类型信号进行滤波，并对结果进行显示和分析。 4．利用matlab 求解差分方程，并做时域和频域分析。用matlab 函数求解单位脉冲响应，并利用窗函数分离信号。 5.用matlab 产生窗函数，并做世玉和频域分析。 6.显示图像，理解图像的模型，将图像进行三原色分解和边缘分析。 三．课程设计题目 一、 1） 生成信号发生器：能产生频率（或基频）为10Hz 的周期性正弦波、三角波和方波信号。绘出它们的时域波形 2） 为避免频谱混叠，试确定各信号的采样频率。说明选择理由。 3）对周期信号进行离散傅立叶变换，为了克服频谱泄露现象，试确定截取数据的长度，即信号长度。分析说明选择理由。 4）绘出各信号频域的幅频特性和相频特性 5）以正弦周期信号为例，观察讨论基本概念（频谱混叠、频谱泄漏、整周期截取等）。 二、已知三个信号()i a p n ，经调制产生信号3 1 ()()cos(/4)i i s n a p n i n π==∑，其中i a 为常 数，()p n 为具有窄带特性的Hanning 信号。将此已调信号通过信道传输，描述该信道的差分方程为 得到接收信号()()*()y n s n h n = 1）分析Hanning 信号()p n 的时域与频域特性 2）分析已调信号()s n 的时域与频域特性 () 1.1172(1)0.9841(2)0.4022(3)0.2247(4) 0.2247()0.4022(1)0.9841(2) 1.1172(3)(4)y n y n y n y n y n x n x n x n x n x n --+---+-= --+---+-
数字信号处理课程设计报告《应用Matlab对信号进行频谱分析及滤波》 专业： 班级： 姓名： 指导老师： 二0 0五年一月一日
目录 设计过程步骤（） 2.1 语音信号的采集（） 2.2 语音信号的频谱分析（） 2.3 设计数字滤波器和画出其频谱响应（） 2.4 用滤波器对信号进行滤波（） 2.5滤波器分析后的语音信号的波形及频谱（） ●心得和经验（）
设计过程步骤 2.1 语音信号的采集 我们利用Windows下的录音机，录制了一段开枪发出的声音，时间在1 s内。接着在C盘保存为WAV格式，然后在Matlab软件平台下．利用函数wavread对语音信号进行采样，并记录下了采样频率和采样点数，在这里我们还通过函数sound引入听到采样后自己所录的一段声音。通过wavread函数和sound的使用，我们完成了本次课程设计的第一步。其程序如下： [x,fs,bite]=wavread('c:\alsndmgr.wav',[1000 20000]); sound(x,fs,bite); 2.2 语音信号的频谱分析 首先我们画出语音信号的时域波形；然后对语音信号进行频谱分析，在Matlab中，我们利用函数fft对信号进行快速傅里叶变换，得到信号的频谱特性性。到此，我们完成了课程实际的第二部。 其程序如下： n=1024; subplot(2,1,1); y=plot(x(50:n/4)); grid on ; title('时域信号') X=fft(x,256); subplot(2,1,2); plot(abs(fft(X))); grid on ; title('频域信号'); 运行程序得到的图形：
南通大学电子信息学院信号与系统课程设计报告 班级： 姓名： 学号： 指导老师： 设计时间： 2014—2015学年第一学期
一、连续信号的时域分析 二、 1. 信号的产生 (1)阶跃函数 function [t,y1]=jieyue(t1,t2,t0) dt=0.01; ttt=t1:dt:t0-dt; tt=t0:dt:t2; t=t1:dt:t2; n=length(ttt); nn=length(tt); u=zeros(1,n); uu=ones(1,nn); y1=[u,uu]; return 冲激函数 function [t,y2]=chongji(t1,t2,t0) dt=0.01; t=t1:dt:t2; n=length(t); y2(1:n)=0; y2(1,(t0-t1)/dt+1)=1/dt; （2）调用上述函数产生信号）2-t ε（，）（4-t δ，-t e ）（t ε，-6s ≤t ≤6s,并画出波形。 Command Window subplot(3,1,1); [t1,y1]=jieyue(-6,6,2);
stairs(t1,y1); axis([-6 6 0 1.5]); subplot(3,1,2); [t2,y2]=chongji(-6,6,4);plot(t2,y2); subplot(3,1,3); [t3,y3]=jieyue(-6,6,0); y3=exp(-(t3)).*y3;plot(t3,y3); 波形如下图所示： （3）根据f(t)画出f(2t)和f(1-0.5t)的波形 t=-3:0.01:3; y=tripuls(t,4,0.6); subplot(3,1,1); plot(t,y);
信号与系统课程设计报告信号与系统课程设计题目 信号与系统课程设计报告 ——频分复用通信系统的仿真设计 指导老师：XXX 小组成员： 摘要： 通过对信号与系统这门课程第八章通信系统学习，我们对频分复用（FDMA ）技术产生了浓厚的兴趣，于是决定自己利用MATLAB 强大的仿真功能来对频分复用系统进行仿真。本文首先录制三段不同的语音信号。然后通过推导，确定合适的载波信号的频率，对信号进行调制，调制后整合到一个复用信号上。再在复用信号上加一个随机的高斯白噪声得到在信道中传输的信号。之后根据通过对复用信号的频谱分析，得出切比雪夫滤波器的各项参数，通过设计好的滤波器进行信号分离后分别根据载波信号进行解调，再通过一个低通滤波器，得到原始信号。通过此次对FDMA 的仿真，我们更清楚了解了频分复用的工作原理，以及AM 调制解调方法，和滤波器的设计方法。频分复用技术对与通信系统节省资源有着重要的意义。
关键词： 频分复用 MATLAB 高斯白噪声 引言： 在电话通信系统中，语音信号频谱在300—3400Hz 内，而一条干线的通信资源往往远大于传送一路语音信号所需的带宽。这时，如果用一条干线只传一路语音信号会使资源大大的浪费，所以常用的方法是“复用”，使一条干线上同时传输几路电话信号，提高资源利用率。 本文是基于MATLAB 的简单应用，首先录制三段不同的语音信号。然后选择合适的高频载波，对信号进行调制，调制后整合到一个复用信号上。确定合适的信噪比，在复用信号上加一个随机的高斯白噪声得到在信道中传输的信号。之后根据载波信号设计合适的带通滤波器将三种信号进行分离，信号分离后分别进行同步解调，再通过一个低通滤波器，得到通过频分复用系统传输后得到的各个信号，将得到的信号与原信号对比，要保证信号与原信号吻合较好。 正文：
数字信号处理上机指导 设计一 正余弦信号的谱分析 【一】 设计目的 1. 用DFT 实现对正余弦信号的谱分析； 2. 观察DFT 长度和窗函数长度对频谱的影响； 3. 对DFT 进行谱分析中的误差现象获得感性认识。 【二】 设计原理 一、谱分析原理 数字信号处理方法的一个重要用途是在离散时间域中确定一个连续时间信号的频谱，通常称为频谱分析，更具体地说，它也包括确定能量谱和功率谱。数字频谱分析可以应用在很广阔领域，频谱分析方法是基于以下的观测：如果连续时间信号)(t g a 是带限的，那么它的离散时间等效信号)(n g 的DFT 进行谱分析。然而，在大多数情况下，)(t g a 是在∞<<∞-t 范围内定义的，因此)(n g 也就定义在∞<<∞-n 的无线范围内，要估计一个无限长信号的频谱是不可能的。实用的方法是：先让模拟连续信号)(t g a 通过一个抗混叠的模拟滤波器，然后把它采样成一个离散序列)(n g 。假定反混叠滤波器的设计是正确的，则混叠效应可以忽略，又假设A/D 变换器的字长足够长，则A/D 变换中的量化噪声也可忽略。 假定表征正余弦信号的基本参数，如振幅、频率和相位不随时间改变，则此信号的傅立叶变换)(ω j e G 可以用计算它的DTFT 得到 ∑ ∞ -∞ =-= n n j j e n g e G ωω )()( （1.1） 实际上无限长序列)(n g 首先乘以一个长度为M 的窗函数)(n w ，使它变成一个长为M 的有限长序列，)()()(1n w n g n g =，对)(1n g 求出的DTFT )(1ω j e G 应该可以作为原连续 模拟信号)(t g a 的频谱估计，然后求出)(1ω j e G 在πω20≤≤区间等分为N 点的离散傅立 叶变换DFT 。为保证足够的分辨率，DFT 的长度N 选的比窗长度M 大，其方法是在截断了的序列后面补上N －M 个零。计算采用FFT 算法。 二、MATLAB 函数介绍 1. 输入函数input( ) 格式：R=input(string) 功能：在屏幕上显示input 括号后的’string ’内容，提示用户从键盘输入某值，并将输入的值赋给R 。 例如，在命令窗口输入R=input(‘How many apples ’) 会显示How many apples
数字信号处理 课 程 设 计 院系：电子信息与电气工程学院 专业：电子信息工程专业 班级：电信班 姓名： 学号： 组员：
摘要 滤波器设计在数字信号处理中占有极其重要的地位，FIR数字滤波器和IIR 滤波器是滤波器设计的重要组成部分。利用MATLAB信号处理工具箱可以快速有效地设计各种数字滤波器。课题基于MATLAB有噪音语音信号处理的设计与实现，综合运用数字信号处理的理论知识对加噪声语音信号进行时域、频域分析和滤波。通过理论推导得出相应结论，再利用 MATLAB 作为编程工具进行计算机实现。在设计实现的过程中，使用窗函数法来设计FIR数字滤波器，用巴特沃斯、切比雪夫和双线性变法设计IIR数字滤波器，并利用MATLAB 作为辅助工具完成设计中的计算与图形的绘制。通过对对所设计滤波器的仿真和频率特性分析，可知利用MATLAB信号处理工具箱可以有效快捷地设计FIR和IIR数字滤波器，过程简单方便，结果的各项性能指标均达到指定要求。 关键词数字滤波器 MATLAB 窗函数法巴特沃斯
目录 摘要 (1) 1 引言 (1) 1.1课程设计目的 (1) 1.2 课程设计内容及要求 (1) 1.3课程设计设备及平台 (1) 1.3.1 数字滤波器的简介及发展 (1) 1.3.2 MATLAB软件简介 (2) 2 课程设计原理及流程 (4) 3.课程设计原理过程 (4) 3.1 语音信号的采集 (4) 3.2 语音信号的时频分析 (5) 3.3合成后语音加噪声处理 (7) 3.3.1 噪声信号的时频分析 (7) 3.3.2 混合信号的时频分析 (8) 3.4滤波器设计及消噪处理 (10) 3.4.1 设计IIR和FIR数字滤波器 (10) 3.4.2 合成后语音信号的消噪处理 (13) 3.4.3 比较滤波前后语音信号的波形及频谱 (13) 3.4.4回放语音信号 (15) 3.5结果分析 (15) 4 结束语 (15) 5 参考文献 (16)
实验一 时域采样与频域采样定理的验证实验 1. (1) 时域采样理论与频域采样理论是数字信号处理中的重要理论。要求掌握模拟信号采样前后频谱的变化，以及如何选择采样频率才能使采样后的信号不丢失信息； (2) 要求掌握频率域采样会引起时域周期化的概念，以及频率域采样定理及其对频域采样点数选择的指导作用。 2. 时域采样定理的要点是： ① 对模拟信号()a x t 以T 进行时域等间隔理想采样，形成的采样信号的频谱 会 以采样角频率Ωs （Ωs=2π/T ）为周期进行周期延拓。公式为 ② 采样频率Ωs 必须大于等于模拟信号最高频率的两倍以上，才能使采样信号的频谱不产生频谱混叠。 利用计算机计算上式并不方便，下面我们导出另外一个公式，以便在计算机上进行实验。 理想采样信号 和模拟信号()a x t 之间的关系为： 对上式进行傅里叶变换，得到： 上式中，在数值上x a (nT)＝x(n)，再将ω=ΩT 代入，得到： 上式的右边就是序列的傅里叶变换，即 上式说明采样信号的傅里叶变换可用相应序列的傅里叶变换得到，只要将自变量ω用ΩT 代替即可。 频域采样定理的要点是： ① 对信号x(n)的频谱函数在［0，2π］上等间隔采样N 点，得到: ?(j )a X Ωa a a s 1??(j )FT[()](j j ) k X x t X k T ΩΩΩ∞ =-∞ ==-∑a ?()x t a a ?()()()n x t x t t nT δ∞ =-∞ =-∑j a a ?(j )[()()]e d t n X x t t nT t ΩΩδ∞ ∞--∞=-∞ =-∑?j a ()()e d t n x t t nT t Ωδ∞ ∞ --∞ =-∞ -∑?＝j a a ?(j )()e nT n X x nT ΩΩ∞ -=-∞ =∑j a ?(j )(e )T X X ωωΩΩ==j 2π()(e ) , 0,1,2, ,1 N k N X k X k N ωω= ==-
一、数字信号处理课程设计内容及考核要求 1、课程设计内容： （1）从以下四个题目中任选其中一个题目，根据题目要求完成程序的编制、调试和仿真； （2）按照题目要求撰写课程设计报告，回答题目设定的问题。 2、考核要求： （一）课程设计以（6——8人）小组完成，但不能出现设计报告雷同情况，一经发现，雷同报告均按不合格处理；最终以PPT小组答辩作为考核。
题目二：有限冲激响应滤波器（FIR）的设计1. 设计目的： 1、加深对数字滤波器的常用指标理解。
2、学习数字滤波器的设计方法。 3. 掌握FIR 滤波器的原理。 2. 设计内容： 利用MATLAB 编程，分别用窗函数法和等波纹滤波器法设计两种FIR 数字滤波器，指标要求如下： 通带边缘频率：ππ65.045.021=Ω=ΩP P ，，通带峰值起伏：][1dB P ≤α。 阻带边缘频率：ππ75.03.021=Ω=ΩS S ，，最小阻带衰减：][40dB S ≥α。 3. 设计原理： 图1 一个典型数字低通滤波器的结构 低通滤波器的常用指标： ? ? ?≤Ω≤Ω≤ΩΩ≤Ω+≤Ω≤-πδδδ|||)(|||1|)(|1S S P P P H H ，， （1）通带边缘频率P Ω； （2）阻带边缘频率S Ω； （3）通带起伏P δ； （4）通带峰值起伏])[1(log 2010dB P P δα--=； （5）阻带起伏S δ，最小阻带衰减])[(log 2010dB S S δα-=。 4. 设计步骤： 1．熟悉MATLAB 的开发环境和使用方法。 2．按照实验内容，编写一个.m 脚本文件，利用MA TLAB 函数fir1和窗函数法设计FIR 数字滤波器。具体参数为：b=fir1(N,Wn,’ftype ’,taper)，N 代表滤波器阶数；Wn 代表滤波器的截止频率(归一化频率)，当设计带通和带阻滤波器时，Wn 为双元素相量；ftype 代表滤
信号与系统课程设计 ——利用matlab实现信号的取样与重构 学院: 工业大学城市学院 专业班级：通信工程C131班 姓名：穆永欢 学号：138213 指导老师：安亚军
目录 摘要 (1) 第一章概述 (1) 第二章设计过程 (2) 2.1设计目的 (2) 2.2设计原理 (2) 2.2.1.MATLAB的介绍 (2) 2.2.2连续时间信号 (3) 2.2.3采样定理 (3) 2.2.4信号重构 (4) 2.3设计容 (4) 2.3.1Sa(t)的临界采样及重构 (4) 2.3.2Sa(t)的过采样及重构 (6) 2.3.3Sa(t)的欠采样及重构 (8) 第三章设计结果分析 (10) 第四章心得体会 (11) 参考文献 (12)
摘要： 本次课程设计以信号与系统和数字信号处理这两门理论与实践紧密结合的课程为基础，经过两个学期的理论学习和上机实验后我们已初步掌握MATLAB软件，通过课程设计更加有助于我们进一步理解和巩固所学知识，学习应用MATLAB 软件的仿真技术，初步掌握线性系统的设计方法，提高分析和解决实际问题的能力，培养独立工作能力。 本实验设计是利用MATLAB实现信号的抽样与重构仿真。通过对该连续的Sa 信号进行抽样，在满足采样定理和不满足采样定理即过抽样和欠抽样两种情况下对连续的Sa信号和采样信号进行频谱分析 【关键词】：信号采样 MATLAB 采样周期频谱信号重构 第一章概述： 针对连续信号的采样与重构问题,利用MATLAB仿真软件平台,仿真不同条件下连续信号的采样信号时域波形和采样后信号频谱、重构信号时域波形和重构后误差波形图。通过对采样周期对采样频谱叠加和信号重构精度的影响、以及信号被采样前后在频域的变化对比分析,得出在不同采样频率的条件下,对应采样信号的时域、频域特性以及重构信号与误差信号也随之产生变化,连续信号可以完全恢复过来。本次课程设计应用MATLAB实现连续信号的采样与重构仿真，了解MATLAB软件，学习应用MATLAB软件的仿真技术。它主要侧重于某些理论知识的灵活运用，以及一些关键命令的掌握，理解，分析等。初步掌握线性系统的设计方法，培养独立工作能力。加深理解采样与重构的概念，掌握利用MATLAB分析系统频率响应的方法和掌握利用MATLAB实现连续信号采用与重构的方法。计算在临界采样、过采样、欠采样三种不同条件下重构信号的误差，并由此总结采样频率对信号重构误差的影响。