电子科技大学成都学院数字信号处理复习重点（共五则）

第一篇：电子科技大学成都学院数字信号处理复习重点

数字信号处理复习重点

第1章

（1）求序列的周期、Z变换和逆Z变换，会求简单序列的DTFT。掌握DTFT与Z变换的关系。

（2）判断系统是否为线性时不变系统，判断系统的稳定性与因果性。

第2章

（1）计算简单序列的DFT。掌握DFT的一些重要性质及应用（圆周共轭对称性，时域、频域循环移位性质，循环卷积）。

（2）DFT对连续信号进行频谱分析时可能出现的几个问题。

（3）掌握基-2 DFT—FFT的基本思想及特点（算法思想，运算量，运算流图，结构规则等），会计算FFT算法的运算量（加法和乘法次数），会画运算流图。

（4）计算两个有限长序列的线性卷积和循环卷积，掌握线性卷积与循环卷积的关系，线性卷积等于循环卷积的条件。

第3章

（1）映射变换的两个基本原则。

（2）脉冲响应不变法和双线性变换法的原理与特点。（优点与缺点，适合设计哪些滤波器）

（3）掌握由模拟低通原型到数字各型滤波器的设计步骤。（重点是低通、高通）

第4章

（1）FIR数字滤波器的线性相位特性。（线性相位满足的条件）

（2）窗口设计法设计FIR滤波器时怎样减少振荡，怎样使过渡带变窄。熟悉四种窗函数的特点，掌握窗长对频谱的影响。

（3）理解频率采样设计法的概念及理论依据，掌握设计步骤及要点。

（4）IIR与FIR数字滤器的比较。（掌握两类滤波器的特点）

第5章

（1）IIR滤波器的各种结构以及各种结构的特点（直接型、级联、并联）

（2）已知系统的差分方程或系统函数，用各种结构来实现滤波器。（掌握IIR直接型、级联型和并联型）

第二篇：电子科技大学 数字信号处理 第二次编程作业

第二次编程作业

加载信号(EEGdata.txt)；分析信号的幅度谱,确定低通FIR数字滤波器的指标；分别利用各种窗函数（Rectangular, Hamming, Kaiser）设计此FIR滤波器；对信号加随机噪声，并用设计的FIR滤波器对含噪声信号进行滤波。要求：

1）给出程序，画出滤波器幅度谱及其增益图(dB)；分析设计的滤波器是否达到指标。

2）利用设计的三种滤波器对加载的信号分别进行滤波，对比分析滤波结果。程序：

data = load('EEGdata.txt');n=1:length(data);

xn=0.1*rand(length(data),1)-0.05;x=xn+data;

[h0,w0]=freqz(data);[h1,w1]=freqz(x);figure(1)

plot(w0/pi,abs(h0),w1/pi,abs(h1),'r');legend('Ô-dataÐźÅƵÆ×ͼ','¼ÓÔëÐźÅƵÆ×ͼ');

wp = 0.15*pi;ws = 0.2*pi;DB = ws-wp;

N = ceil(2*3.11*pi/DB+1);wc =(wp+ws)/2;wc = wc/pi;

win = rectwin(N+1);

b = fir1(N, wc, 'low', win);figure(2)freqz(b,1);

wp = 0.15*pi;ws = 0.2*pi;DB = ws-wp;

N = ceil(2*3.11*pi/DB+1);wc =(wp+ws)/2;wc = wc/pi;

win = hamming(N+1);

b1 = fir1(N, wc, 'low', win);figure(3)freqz(b1,1);

fpts = [0.15, 0.2];mag = [1, 0];dev = [0.01, 0.01];

[N, Wn, beta, ftype] = kaiserord(fpts, mag, dev)Kwin = kaiser(N+1, beta);b2 = fir1(N, Wn, Kwin);[h,omega] = freqz(b2,1,512);figure(4)freqz(b2,1);

y=filter(b, 1, x);figure(5)subplot(2,1,1);

plot(n,data,n,x,'g',n,y,'r');title('¾ØÐδ°Â˲¨·ù¶ÈÆ×Ч¹û¶Ô±Èͼ');legend('Ô-ÐźÅdata','¼ÓÔëÐźÅx','Â˲¨ºóÐźÅy');[H,W]=freqz(y);subplot(2,1,2);

plot(w0/pi,abs(h0),W/pi,abs(H),'r');title('¾ØÐδ°Â˲¨ÆµÆ×Ч¹û¶Ô±Èͼ');legend('Ô-ÐźÅdata','Â˲¨ºóÐźÅy');

y1=filter(b1, 1, x);figure(6)subplot(2,1,1);

plot(n,data,n,x,'g',n,y1,'r');title('Hamming´°Â˲¨·ù¶ÈÆ×Ч¹û¶Ô±Èͼ');legend('Ô-ÐźÅdata','¼ÓÔëÐźÅx','Â˲¨ºóÐźÅy1');[H,W]=freqz(y1);subplot(2,1,2);

plot(w0/pi,abs(h0),W/pi,abs(H),'r');title('Hamming´°Â˲¨ÆµÆ×Ч¹û¶Ô±Èͼ');legend('Ô-ÐźÅdata','Â˲¨ºóÐźÅy1');

y2=filtfilt(b2, 1, x);figure(7)subplot(2,1,1);

plot(n,data,n,x,'g',n,y2,'r');title('Kaiser´°Â˲¨·ù¶ÈÆ×Ч¹û¶Ô±Èͼ');legend('Ô-ÐźÅdata','¼ÓÔëÐźÅx','Â˲¨ºóÐźÅy2');[H,W]=freqz(y2);subplot(2,1,2);

plot(w0/pi,abs(h0),W/pi,abs(H),'r');title('Kaiser´°Â˲¨ÆµÆ×Ч¹û¶Ô±Èͼ');legend('Ô-ÐźÅdata','Â˲¨ºóÐźÅy2');

figure(8)

plot(n,data,n,y,'c',n,y1,'r',n,y2,'g');title('ÈýÖÖ´°º¯ÊýÂ˲¨Ð§¹û¶Ô±Èͼ');legend('Ô-dataÐźÅ','¾ØÐδ°Â˲¨','Hamming´°Â˲¨','Kaiser´°Â˲¨');

***6420 0原data信号频谱图加噪信号频谱图 0.10.20.30.40.50.60.70.80.91 50Magnitude(dB)0-50-10000.10.20.30.40.50.60.70.8Normalized Frequency( rad/sample)0.910Phase(degrees)-1000-2000-300000.10.20.30.40.50.60.70.8Normalized Frequency( rad/sample)0.91

50Magnitude(dB)0-50-100-15000.10.20.30.40.50.60.70.8Normalized Frequency( rad/sample)0.910Phase(degrees)-1000-2000-300000.10.20.30.40.50.60.70.8Normalized Frequency( rad/sample)0.91 50Magnitude(dB)0-50-100-15000.10.20.30.40.50.60.70.8Normalized Frequency( rad/sample)0.910Phase(degrees)-500-1000-1500-200000.10.20.30.40.50.60.70.8Normalized Frequency( rad/sample)矩形窗滤波幅度谱效果对比图0.20.10-0.1-0.2 ***原信号data加噪信号x滤波后信号y300350 0.91矩形窗滤波频谱效果对比图1510原信号data滤波后信号y 50 00.10.20.30.40.50.60.70.80.91Hamming窗滤波幅度谱效果对比图0.20.10-0.1-0.2 050100150200原信号data加噪信号x滤波后信号y1250300350 Hamming窗滤波频谱效果对比图1510原信号data滤波后信号y1 50 00.10.20.30.40.50.60.70.80.91Kaiser窗滤波幅度谱效果对比图0.20.10-0.1-0.2 050100150200原信号data加噪信号x滤波后信号y2250300350 Kaiser窗滤波频谱效果对比图20151050 00.10.20.30.40.50.60.70.80.91原信号data滤波后信号y2 三种窗函数滤波效果对比图0.2原data信号矩形窗滤波Hamming窗滤波Kaiser窗滤波 0.150.10.050-0.05-0.1-0.***0200250300350

结果分析：滤波效果皆不尽如人意，原因是加载的噪声信号是随机信号，各频率皆有，滤波器通带部分的噪声无法滤掉。尤其是矩形窗函数和Hamming窗函数构造的滤波器长度大于data的三分之一，无法使用filtfilt函数滤波，只能用filter函数，致使滤波结果有相位平移。

第三篇：数字信号处理复习总结

数字信号处理复习要点

数字信号处理主要包括如下几个部分

1、离散时间信号与系统的基本理论、信号的频谱分析

2、离散傅立叶变换、快速傅立叶变换

3、数字滤波器的设计

一、离散时间信号与系统的基本理论、信号的频谱分析

1、离散时间信号：

1）离散时间信号。时间是离散变量的信号，即独立变量时间被量化了。信号的幅值可以是连续数值，也可以是离散数值。2）数字信号。时间和幅值都离散化的信号。

（本课程主要讲解的实际上是离散时间信号的处理）3）离散时间信号可用序列来描述 4）序列的卷积和（线性卷积）

y(n)mx(m)h(nm)x(n)*h(n)

5）几种常用序列

1,n0a)单位抽样序列（也称单位冲激序列）(n),(n)

0,n01,n0b)单位阶跃序列u(n),u(n)

0,n01,0nN1c)矩形序列，RN(n)

0,n其它d)实指数序列,x(n)anu(n)

6）序列的周期性

所有n存在一个最小的正整数N,满足：x(n)x(nN),则称序列x(n)是周期序列，周期为N。(注意：按此定义，模拟信号是周期信号，采用后的离散信号未必是周期的)

7）时域抽样定理：

一个限带模拟信号xa(t)，若其频谱的最高频率为F0，对它进行等间隔抽样而得x(n)，抽样周期为T，或抽样频率为Fs1/T；

只有在抽样频率Fs2F0时，才可由xa(t)准确恢复x(n)。

2、离散时间信号的频域表示（信号的傅立叶变换）

X(j)nx(n)ejn，X(j(2))X(j)

1x(n)X(j)ejnd 2

3、序列的Z变换

X(z)Z[x(n)]nx(n)zn

1）Z变换与傅立叶变换的关系，X(j)X(z)zej

2）Z变换的收敛域

收敛区域要依据序列的性质而定。同时，也只有Z变换的收敛区域确定之后，才能由Z变换唯一地确定序列。

一般来来说，序列的Z变换的收敛域在Z平面上的一环状区域：Rx|z|Rx

x(n)N1nN23）有限长序列：x(n)，0|z|

0其它x(n)N1n右序列：x(n)，|Z|>Rx-

其它0x(n)nN2左序列：x(n)，0其它（|z|0时：0≤|Z|< Rx+；N2≤0时： 0<|Z|< Rx+）双边序列：x(n),n，Rx|z|Rx

常用序列的Z变换：

Z[(n)]1,|z|01,|z|111z

1Z[anu(n)],|z||a|1az11Z[bnu(n1)],|z||b|1bz1Z[u(n)] 逆变换

x(n)12jn1X(z)zdzx，C：收敛域内绕原点逆时针的一条闭合曲线 c1）留数定理：x(n)[X(z)zn1在C内极点留数之和] 2）留数辅助定理：x(n)[X(z)zn1在C外极点留数之和] 3）利用部分分式展开：X(z)Z变换求解。

4、离散时间系统：

T[x(n)]y(n)系统函数：H(j)Y(j)Y(z)，H(z) X(j)X(z)Ak，然后利用定义域及常用序列的1akz1冲激响应：h(n)T[(n)]

5、线性系统：满足叠加原理的系统。T[ax(n)by(n)]aT[x(n)]bT[y(n)]

6、移不变系统：若T[x(n)]Y(n)，则T[x(nk)]Y(nk)

7、线性移不变系统

可由冲激响应来描述（系统的输出相应是输入与单位冲激响应的线性卷积）

y(n)x(n)*h(n)，Y(j)X(j)H(j)，Y(z)X(z)H(z)

8、系统的频率特性可由其零点及极点确定

X(z)bziMiak0i0NA(1zziM1)Akzk(1zk1i1N(zz)ziMMkz1)(zzk1i1N

k)zN（式中，zk是极点，zi是零点；在极点处，序列x(n)的Z变换是不收敛的，因此收敛区域内不应包括极点。）

9、稳定系统：有界的输入产生的输出也有界的系统，即：若|x(n)|，则|y(n)|

线性移不变系统是稳定系统的充要条件：

n|h(n)|

或：其系统函数H(z)的收敛域包含单位园 |z|=1

10、因果系统：n0时刻的输出y(n0)只由n0时刻之前的输入x(n),nn0决定

线性移不变系统是因果系统的充要条件：h(n)0,n0 或：其系统函数H(z)的收敛域在某园外部：即：|z|>Rx

11、稳定因果系统：同时满足上述两个条件的系统。

h(n)0,n0 线性移不变系统是因果稳定系统的充要条件：|h(n)|，n或：H(z)的极点在单位园内 H(z)的收敛域满足：|z|Rx,Rx1

12、差分方程

线性移不变系统可用线性常系数差分方程表示（差分方程的初始条件应满足松弛条件）

aynkbxni

kik0i0NM13、差分方程的解法 1）直接法：递推法 2）经典法

3）由Z变换求解

二、离散傅立叶变换、快速傅立叶变换

1、周期序列的离散傅立叶级数（DFS）

Xp(k)DFS[xp(n)]xp(n)en0N1j2knNkn xp(n)WNn0N11xp(n)IDFS[Xp(k)]N其中：WN=ej2/N

KON1XPke2jknN1NKON1XPkWNkn

2、有限长序列的离散傅立叶变换(DFT)

knX(k)DFT[x(n)]{DFS[x(nN)]}RN(k)x(n)WN，0≤k≤N1

n0N11N1kn x(n)IDFT[X(k)]{IDFS[X(kN)]}RN(n)X(k)WN，0≤n≤N1

Nk0应当注意，虽然x(n)和X(k)都是长度为N得有限长序列，但他们分别是由周期序列xp(n)和Xp(k)截取其主周期得到的，本质上是做DFS或IDFS，所以不能忘记它们的隐含周期性。尤其是涉及其位移特性时更要注意。

3、离散傅立叶变换与Z变换的关系 X(k)X(j)|2X(z)|j2k

NkzeN

4、频域抽样定理

对有限长序列x(n)的Z变换X(z)在单位圆上等间隔抽样，抽样点数为N，或抽样间隔为2/N，当N≥M时，才可由X(k)不失真恢复X(j)。

1zN内插公式：X(z)N

5、周期卷积、循环卷积

周期卷积：xp3(n)xp1(m)xp2(nm)

m0N1X(k)k1k01WNzN1循环卷积：x3(n)x1(n)N1x2(n)xp3(n)RN(n)xp1(m)xp2(nm)RN(n)

m0

6、用周期（周期）卷积计算有限长序列的线性卷积

对周期要求：NN1N21（N1、N2分别为两个序列的长度）

7、基2 FFT算法 1）数据要求：N2M 2）计算效率（乘法运算次数：NM，加法计算次数：NM）（复数运算）（DFT运算：乘法运算次数：N2，加法计算次数：N2）（复数运算）

8、快速卷积（采用FFT计算）

9、分辨率

三、数字滤波器的设计

（一）FIR滤波器的设计

1、特点：可实现严格的线性相位特性、系统是稳定的、因果的、阶数较高

2、实现线性相位的条件（1）h(n)为实数（2）h(n)=h(N-1-n)做一般意义下的FIR滤波器，N是偶数，不适合做高通滤波器 或 h(n)=-h(N-1-n)对称中心：(N-1)/2 适于做希尔伯特变换器，微分器和正交网络。

3、主要设计方法 1）窗函数法

2）频率抽样设计

频率抽样内插公式设计。特点：

频率特性可直接控制。

若滤波器是窄带的，则能够简化系统

若无过渡带样本，则起伏较大。改进办法是增加过渡带样本，采用过渡带的自由变量法，通常使用优化方法求解。可得到较好的起伏特性，但是会导致过渡带宽度加大，改进办法是增加抽样点数。

抽样点的获得采取两种办法：I型抽样及II型抽样。

若要满足线性相位特性，则相位要满足一定要求。

（二）IIR滤波器的设计

1、特点

• 阶数少、运算次数及存储单元都较少 • 适合应用于要求相位特性不严格的场合。

• 有现成的模拟滤波器可以利用，设计方法比较成熟。• 是递归系统，存在稳定性问题。

2、主要设计方法

先设计模拟滤波器，然后转换成数字滤波器。设计过程：

1）先设计模拟低通滤波器Ha(s)：butterworth滤波器设计法等，有封闭公式利用

2）将模拟原型滤波器变换成数字滤波器（1）模拟低通原型先转换成数字低通原型，然后再用变量代换变换成所需的数字滤波器；  模拟低通原型先转换成数字低通原型：HaL(s)HL(z)，主要有冲激不变法、阶跃不变法、双线性变换法等。

 将数字低通原型滤波器通过变量代换变换成所需的数字滤波器。，z1G(Z1)HL(z)HD(Z)

（2）由模拟原型变成所需型式的模拟滤波器，然后再把它转换成数字滤波器；

 将模拟低通原型滤波器通过变量代换变换成所需的模拟滤波器。HaL(s)HaD(S1)，sF(S1)

 模拟滤波器转换成数字数字滤波器：HaD(s)HD(z)，主要有冲激不变法、阶跃不变法、双线性变换法等

（3）由模拟原型直接转换成所需的数字滤波器

直接建立变换公式：HaL(s)HD(z)，sG(z1)

3、模拟数字转换法（1）冲激不变法

H(z)ZL1[Ha(s)]|tnT

单阶极点情况

NAkAk'skT' H(z)，Ha(s)AApekkk11pzssk1k1kkN

（2）阶跃不变法

H(z)z1ZL1[Ha(s)/s]|tnT z

冲激不变法和阶跃不变法的特点： • 有混叠失真

• 只适于限带滤波器

• 不适合高通或带阻数字滤波器的设计

1z1（3）双线性变换法 sC 11z常数C的计算：1）Cccot(c2)2）C=2/T 特点：

(i)稳定性不变（ii）无混叠

（iii）频率非线性变换，会产生畸变，设计时，频率要做预畸变处理

4、直接法设计IIR数字滤波器 • z平面的简单零极点法

（三）滤波器的网络结构

第四篇：数字信号处理实验教学-电子教案

数字信号处理实验 电子教案

讲义1 Matlab简介及其安装使用说明..................................................2 讲义2 Matlab基本语句..........................................................................8 讲义3 Matlab基本数值运算................................................................13 讲义4 Matlab函数、及其调用方法....................................................16 实验1 常见离散信号产生和实现.........................................................20 实验2 离散系统的时域分析.................................................................22 实验3 FFT算法的应用..........................................................................24 实验4 离散系统的变换域分析.............................................................27 实验5 有限冲激响应数字滤波器设计.................................................32 实验6 无限冲激响应数字滤波器设计.................................................36 实验7 设计性和研究性实验.................................................................41讲义1 Matlab简介及其安装使用说明

一.MATLAB程序设计语言简介

MATLAB，Matrix Laboratory的缩写，是由Mathworks公司开发的一套用于科学工程计算的可视化高性能语言，具有强大的矩阵运算能力。与大家常用的Fortran和C等高级语言相比，MATLAB的语法规则更简单，更贴近人的思维方式，被称之为“草稿纸式的语言”。截至目前，MATLAB已经发展到7.x版，适用于所有32位的Windows操作系统，按NTFS(NT文件系统)格式下完全安装约需 850 MB。MATLAB软件主要由主包、仿真系统和工具箱三大部分组成。

二.MATLAB应用入门

1.MATLAB的安装与卸载

MATLAB软件在用户接口设计上具有较强的亲和力，其安装过程比较典型，直接运行光盘中的安装向导支撑程序SETUP.exe，按其提示一步步选择即可。MATLAB自身带有卸载程序，在其安装目录下有uninstall子目录，运行该目录下的uninstall.exe即可； 也可以通过Windows系统的安装卸载程序进行卸载。2.MATLAB的启动与退出

MATLAB安装完成后，会自动在Windows桌面上生成一个快捷方式，它是指向安装目录下binwin32matlab.exe的链接，双击它即可来到MATLAB集成环境的基本窗口，通常称之为命令窗口。MATLAB的退出与普通WIN32的程序一样，值得一提的是它有一个自身专有的快捷键Ctrl+Q。3.MATLAB界面简介

图 1 MATLAB基本界面——命令窗口

图2

图3

图4

图5

图6 指令历史

图7 1)菜单栏

菜单栏中包括File、Edit、View、Web、Window和Help六个菜单项。这里着重介绍File、help项。

File项：数据输入/输出的接口，包括10个子项，这里重点介绍其中的5个子项： New：新建文件项。有四个选择：M File(*.M，文本格式的MATLAB程序文件，可以直接通过文件名的方式在MATLAB环境下解释运行；Figure(图形);Model(仿真模型文件)和GUI(可视化界面文件)。

Open：打开所有MATLAB支持的文件格式，系统将自动识别并采用相应的程序对文件进行处理。例如, 打开一个.m文件，系统将自动打开M文件编辑器对它进行编辑。

Set Path...：设置工作路径。可以打开路径设置(Set Path)对话框(图2)，将用户自己建立的目录加入MATLAB的目录系统中，以便所编制的文件能够在MATLAB环境中直接调用。

图8 路径设置对话框

单击Add Folder...按钮可以将你的一个文件夹加入到系统路径中；Add with Subfolders...允许把一个文件夹包括其所有的子文件夹加入到系统路径中。这两种操作均可以直观地在右侧的路径栏内看到结果。

选中一个加入的文件夹，你可以利用Move to Top(移至所有路径的最前面),Move Up(上移一个),Move Down(下移一个),Move to Bottom(移至所有路径的最后面)等四个按钮将改变文件在系统路径中的排列位置以利于对文件的搜索使用，也可以利用Remove按钮将其删除。

对路径操作完毕后，按Save按钮予以保存;Help 项：

Matlab Help：打开以html超文本形式存储的的帮助文件主页； Demos：打开matlab演示窗主页；

About Matlab：Matlab注册图标、版本、制造商和用户信息；

图9 Help选项

图10 Help窗口

2)命令行区

进行通用操作，数值计算，编程和数据类型，输入输出，绘图，用户界面等命令，例如，命令：help函数名（*.m文件）； edit编辑函数、文件

对输入命令的解释MATLAB按以下顺序进行：

① 检查它是否是工作空间中的变量，是则显示变量内容。② 检查它是否是嵌入函数，是则运行之。③ 检查它是否是子函数。④ 检查它是否是私有函数。

⑤ 检查它是否是位于MATLAB搜索路径范围内的函数文件或脚本文件。

请注意，如果有两个以上的方案与输入的命令相匹配，MATLAB将只执行第一个匹配。

讲义2 Matlab基本语句

一.程序控制语句

1.循环语句

MATLAB的循环语句包括for循环和while循环两种类型。1)for循环 语法格式：

for 循环变量 = 起始值：步长：终止值 循环体 end 起始值和终止值为一整形数，步长可以为整数或小数，省略步长时，默认步长为1。执行for循环时，判定循环变量的值是否大于(步长为负时则判定是否小于)终止值，不大于(步长为负时则小于)则执行循环体，执行完毕后加上步长，大于(步长为负时则小于)终止值后退出循环。例1 给矩阵A、B赋值。MATLAB 语句及运行结果如下： k=5;a=zeros(k, k)%矩阵赋零初值 for m=1 : k for n=1: k a(m,n)=1/(m+n-1);end end for i=m :-1 : 1 b(i)=i;end 运行结果：  a= 1.0000 0.5000 0.3333 0.2500?0.2000 0.5000 0.3333 0.2500 0.2000 0.1667 0.3333 0.2500 0.2000 0.1667 0.1429 0.2500 0.2000 0.1667 0.1429 0.1250 0.2000 0.1667 0.1429 0.1250 0.1111 b= 1 2?3 4 5 2)while循环 语法格式：

while 表达式 循环体 end 其执行方式为：若表达式为真(运算值非0)，则执行循环体；结果为0)，则退出循环体，执行end后的语句。

例2 a=3;while a a=a-1 end 输出: a=2 a=1 a=0 2.条件转移语句

条件转移语句有if和switch两种。

若表达式为假(运算 1)if语句

MATLAB中if语句的用法与其他高级语言相类似，其基本语法格式有以下几种： 格式一： if 逻辑表达式 执行语句 end 格式二： if 逻辑表达式 执行语句1 else 执行语句2 end 格式三：

if 逻辑表达式1 执行语句1

else? if 逻辑表达式2 执行语句2 end 2)switch语句

switch语句的用法与其他高级语言相类似，其基本语法格式为：

switch表达式(标量或字符串)case 值1 语句1 case 值2 语句2 … otherwise 语句n end

二.绘图语句 常用的MATLAB绘图语句有figure、plot、subplot、stem等，图形修饰语句有title、axis、text等。

1．figure figure有两种用法，只用一句figure命令，会创建一个新的图形窗口，并返回一个整数型的窗口编号。figure(n)表示将第n号图形窗口作为当前的图形窗口，并将其显示在所有窗口的最前面;如果该图形窗口不存在，则新建一个窗口，并赋以编号n。2．plot 线型绘图函数。用法为plot(x,y,'s')。参数x为横轴变量，y为纵轴变量，s用以控制图形的基本特征如颜色、粗细等，通常可以省略，常用方法如表1所示。

表1

3．Stem 绘制离散序列图，常用格式stem(y)和stem(x,y)分别和相应的plot函数的绘图规则相同，只是用stem命令绘制的是离散序列图。4．subplot subplot(m,n,i)图形显示时分割窗口命令，把一个图形窗口分为m行，n列，m×n个小窗口，并指定第i个小窗口为当前窗口。5．绘图修饰命令

在绘制图形时，我们通常需要为图形添加各种注记以增加可读性。在plot语句后使用title('标题')可以在图形上方添加标题，使用xlabel('标记')或ylabel('标记')为X轴或Y轴添加说明，使用text(X值、Y值、'想加的标示')可以在图形中任意位置添加标示。例3 画图基本语句如图1所示。

MATLAB 语句及运行结果如下： x=0:0.1*pi:2*pi;%定义x向量

figure(1);%创建一个新的图形窗口，编号为1 subplot(2,2,1);%将窗口划分为2行，2列，在第1个窗口中作图 plot(x,sin(x));%画图

title('正弦线');%给图形加标题 subplot(2,2,2);%在第2个窗口中作图 plot(x,sin(x),'r');%画一正弦波，红色 xlabel('X');%给x轴加说明 ylabel('SIN(X)');%给y轴加说明 subplot(2,2,3);%在第3个窗口中作图 plot(x,sin(x),'--');%画一正弦波，破折线 subplot(2,2,4);%在第4个窗口中作图 plot(x,sin(x),'r+');%画一正弦波，红色+线 text(4,0,'注记');

图1 讲义3 Matlab基本数值运算

一.MATLAB内部特殊变量和常数

 MATLAB内部有很多变量和常数，用以表达特殊含义。常用的有：

1.变量ans： 指示当前未定义变量名的答案；

2.常数eps：表示浮点相对精度，其值是从1.0到下一个最大浮点数之间的差值。该变量值作为一些MATLAB函数计算的相对浮点精度，按IEEE标准，如：，近似为2.2204e-016；

3.常数Inf：表示无穷大。当输入或计算中有除以0时产生Inf；

4.虚数单位i,j：表示复数虚部单位，相当于5.NaN：表示不定型值，是由0/0运算产生的。

6.常数pi：表示圆周率π，其值为3.14***…；

；

 MATLAB中可表示的数字的近似范围从，1.有效数字表示的典型例子如下：1234.56789,123456.789E-2,1.23456789e3（format指令可以控制显示格式）

2.复数形式：3.5+4*j,-2.1-7.4*j（i 也可以）

取绝对值：abs()语法格式：abs(x)。当x为实数时计算x的绝对值；x为复数时得到的是复数的模值；x为字符串时得到各字符的ASCII码。取相角：angle()语法格式：angle(z)。求复矢量或复矩阵的相角，结果为一个以弧度为单位介于-π和+π之间的值。

二.变量

1.变量命名规则

MATLAB中对变量的命名应遵循以下规则：

（1）变量名可以由字母、数字和下划线混合组成，但必须以字母开头；（2）字符长度不能大于31；（3）变量命名区分大小写。2.局部变量和全局变量 局部变量是指那些每个函数体内自己定义的，不能从其他函数和MATLAB工作空间访问的变量。全局变量是指用关键字“global”声明的变量。全局变量名应尽量大写，并能反映它本身的含义。如果需要在工作空间和几个函数中都能访问一个全局变量，必须在工作空间和这几个函数中都声明该变量是全局的。

三.矩阵及其运算

MATLAB具有强大的矩阵运算和数据处理功能，对矩阵的处理必须遵从代数规则。1．矩阵生成1）一般矩阵的生成

对于一般的矩阵MATLAB的生成方法有多种。最简单的方法是从键盘直接输入矩阵元素。直接输入矩阵元素时应注意：各元素之间用空格或逗号隔开，用分号或回车结束矩阵行，用中括号把矩阵所有元素括起来。

例1 在工作空间产生一个3×3矩阵A可用MATLAB语言描述如下： A=［1 2 3;4 5 6;7 8 9］ 或 A=［1 2 3 4 5 6 7 8 9］ 运行结果： A= 1 2 3 4 5 6 7 8 9 Size(A)得到矩阵的大小，ans = 3 3 2）特殊矩阵的生成

对于特殊的矩阵可直接调用MATLAB的函数生成。

用函数zeros生成全0矩阵：格式 B=zeros(m,n)生成m×n的全0阵。用函数ones生成全1矩阵：格式 B=ones(m,n)生成m×n的全1阵。

用函数eye生成单位阵：格式 B=eye(m,n)生成m×n矩阵，其中对角线元素全为1，其他元素为0。2．矩阵的运算

矩阵的运算有基本运算和函数运算两种类型。基本运算包括矩阵的加、减、乘、除、乘方、求转置、求逆等，其主要特点是通过MATLAB提供的基本运算符＋、－、*、/（）、^等即可完成。函数运算主要是通过调用MATLAB系统内置的运算函数来求取矩阵的行列式（det(A)）,求秩（rank(A)）,求逆(inv(A))，求特征值和特征向量（［V，D］＝eig(A)）, 求Jordan标准形（jordan(A)）和矩阵分解等。需要用时可以参阅联机帮助和相关参考书。（.*,.,表示逐个元素的乘积和相除；矩阵X/Y 相当于X*inv(Y), XY相当于inv(Y)*X）

例2 矩阵的基本运算。A=［1, 2, 3;4, 5, 6］;B =［6, 5, 4;3, 2, 1］;C =A+B %计算两个矩阵的和 D =B' %计算矩阵B的转置

E=A*D %做矩阵乘法，必须要满足矩阵乘法的基本要求E应该是2阶方阵 F=det(E)%求E的行列式值 G=E^(-1)%求E的逆

输出结果： C= 7 7 7 7 7 7 D= 6 3 5 2 4 1 E= 28 10 73 28 F=54 G= 0.5185-0.1852-1.3519 0.5185 讲义4 Matlab函数、及其调用方法

在MATLAB语言中，M文件有两种形式：脚本和函数。

脚本没有输入/输出参数，只是一些函数和命令的组合。它可以在MATLAB环境下直接执行，也可以访问存在于整个工作空间内的数据。由脚本建立的变量在脚本执行完后仍将保留在工作空间中可以继续对其进行操作，直到使用clear命令对其清除为止。

函数是MATLAB语言的重要组成部分。MATLAB提供的各种工具箱中的M文件几乎都是以函数的形式给出的。函数接收输入参数，返回输出参数，且只能访问该函数本身工作空间中的变量，从命令窗或其他函数中不能对其工作空间的变量进行访问。

1．函数结构

MATLAB语言中提供的函数通常由以下五个部分组成：

（1）函数定义行: 以function开头，函数名(必须与文件名相同)及函数输入输出参数在此定义；（2）H1行：第一注释行，供lookfor和help在线帮助使用；

（3）函数帮助文件；通常包括函数输入输出参数的含义，调用格式说明；

（4）函数体：它包括进行运算和赋值的所有MATLAB程序代码。函数体中可以包括流程控制、输入/输出、计算、赋值、注释以及函数调用和脚本文件调用等。在函数体中完成对输出参数的计算；

（5）注释。

这五个部分中最重要的是函数定义行和函数体。

函数定义行是一个MATLAB函数所必需的，其他各部分的内容可以没有，这种函数称为空函数。

例1 function sa=circle(r,s)% Circle plot a circle of radii r in the line specified by s % r raddi % s line color % sa area of the circle % % circle(r)use blue line to draw a circle of radii r % circle(r,s)use 's' to draw circle % sa=circle(r)compute the area of the circle and draw it in blue % sa=circle(r,s)compute the area of the circle and draw it in color 's' if nargin>2 error('to many input ');end clf;t=0:pi/100:2*pi;x=r*exp(i*t);if nargout==0 plot(x,s);else sa=pi*r*r;fill(real(x),imag(x),s)end axis('square');%makes the current axis box square in size.例2 function y = imp_fun(n,n0)% IMP_FUN Unit impulse function.% % IMP_FUN(N,N0), where N is a vector of sequential integers, % returns a vector the same length as N with zeros everywhere except % N = N0.[m1,n1] = size(n);if ~(m1 == 1 | n1 == 1)error('The sample vector must be one-dimensional.');end y = zeros(m1,n1);i = find(n >= n0);if isempty(i)return end y(i(1))= 1;% where n = n0, set output to 1 2．函数调用 函数调用的过程实际上就是参数传递的过程。例如，在一个脚本文件里调用函数“max”可采用如下方式： n=1:20;

a=sin(2*pi*n/20); [Y,I]=max(a);该调用过程把变量“a”传给了函数中的输入参数“x”,然后把函数运算的返回值传给输出参数“Y”和“I”。其中，Y是a序列的最大值，I是最大值Y对应的坐标值。

例3 构造：y[n] = δ[n-3]: 调用函数： n = 0:6;y = imp_fun(n,3);stem(n,y)

图1 例4 构造：y[n] = 5δ[n]imp_fun(n,2);stem(n,y)；

图2 实验1 常见离散信号产生和实现

一、实验目的

1、加深对常用离散信号的理解；

2、熟悉使用MATLAB在时域中产生一些基本的离散时间信号。

二、实验原理

1、单位抽样序列

在MATLAB中可以利用

函数实现。

2、单位阶越序列

在MATLAB中可以利用

函数实现:

3、正弦序列

在MATLAB中实现过程如下：

4、复指数序列

在MATLAB中实现过程如下：

5、指数序列

在MATLAB中实现过程如下：

三、预习要求

1、预先阅读实验讲义（MATLAB基础介绍）；

2、讨论正弦序列、复指数序列的性质。

A．绘出信号，当、时、、时的信号实部和虚部图；当B．绘出信号

时呢？此时信号周期为多少？ 的频率是多少？周期是多少？产生一个数字频率为0.9的正弦序列，并显示该信号，说明其周期。

3、使用帮助功能学习square(方波),sawtooth(锯齿波)和sinc函数，并绘图。

四、实验内容

编制程序产生上述5种信号，长度可输入确定，函数需要的参数可输入确定，并绘出其图形。

实验2 离散系统的时域分析

一、实验目的

1、熟悉并掌握离散系统的差分方程表示法；

2、加深对冲激响应和卷积分析方法的理解。

二、实验原理

在时域中，离散时间系统对输入信号或者延迟信号进行运算处理，生成具有所需特性的输出信号，具体框图如下:

其输入、输出关系可用以下差分方程描述：

输入信号分解为冲激信号，记系统单位冲激响应，则系统响应为如下的卷积计算式：

当称系统为IIR系统。

时，h[n]是有限长度的（），称系统为FIR系统；反之，三、预习要求

1、在MATLAB中，熟悉利用函数

2、在MATLAB中，熟悉用函数 响应的过程。

实现差分方程的仿真；

计算卷积，用求系统冲激

四、实验内容

1、以下程序中分别使用conv和filter函数计算h和x的卷积y和y1，运行程序，并分析y和y1是否有差别，为什么要使用x[n]补零后的x1来产生y1；具体分析当h[n]有i个值，x[n]有j个值，使用filter完成卷积功能，需要如何补零？ % Program P2_7 clf;h = [3 2 1-2 1 0-4 0 3];%impulse response x = [1-2 3-4 3 2 1];%input sequence y = conv(h,x);n = 0:14;subplot(2,1,1);stem(n,y);xlabel('Time index n');ylabel('Amplitude');title('Output Obtained by Convolution');grid;x1 = [x zeros(1,8)];y1 = filter(h,1,x1);subplot(2,1,2);stem(n,y1);xlabel('Time index n');ylabel('Amplitude');title('Output Generated by Filtering');grid;

2、编制程序求解下列两个系统的单位冲激响应和阶跃响应，并绘出其图形。要求分别用 filter、conv、impz三种函数完成。,，给出理论计算结果和程序计算结果并讨论。实验3 FFT算法的应用

一、实验目的

1、加深对离散信号的DFT的理解；

2、在MATLAB中实现FFT算法。

二、实验原理

N点序列的DFT和IDFT变换定义式如下： , , 利用旋转因子具有周期性，可以得到快速算法（FFT）。

在MATLAB中，可以用函数反变换。

和计算N点序列的DFT正、三、预习要求

1、在MATLAB中，熟悉函数fft、ifft的使用；

2、阅读扩展练习中的实例，学习在MATLAB中的实现FFT算法的实现；

3、利用MATLAB编程完成计算，绘出相应图形。并与理论计算相比较，说明实验结果的原因。

四、实验内容 1、2N点实数序列

N=64。用一个64点的复数FFT程序，一次算出，形。

并绘出 的图

2、已知某序列在单位圆上的N=64等分样点的Z变换为:。

用N点IFFT程序计算出和。

五、扩展练习

例1：对连续的单一频率周期信号 按采样频率和N =16，观察其DFT结果的幅度谱。

采样，截取长度N分别选N =20解：此时离散序列算并作图，函数fft用于计算离散傅里叶变换DFT，程序如下： k=8;n1=[0:1:19];xa1=sin(2*pi*n1/k);subplot(2,2,1)plot(n1,xa1)xlabel('t/T');ylabel('x(n)');xk1=fft(xa1);xk1=abs(xk1);subplot(2,2,2)stem(n1,xk1)xlabel('k');ylabel('X(k)');n2=[0:1:15];，即k=8。用MATLAB计xa2=sin(2*pi*n2/k);subplot(2,2,3)plot(n2,xa2)xlabel('t/T');ylabel('x(n)');xk2=fft(xa2);xk2=abs(xk2);subplot(2,2,4)stem(n2,xk2)xlabel('k');ylabel('X(k)');

图 不同的截取长度的正弦信号及其DFT结果

计算结果示于图，(a)和(b)分别是N=20时的截取信号和DFT结果，由于截取了两个半周期，频谱出现泄漏；(c)和(d)分别是N=16时的截取信号和DFT结果，由于截取了两个整周期，得到单一谱线的频谱。上述频谱的误差主要是由于时域中对信号的非整周期截断产生的频谱泄漏。

实验4 离散系统的变换域分析

一、实验目的

1、熟悉对离散系统的频率响应分析方法；

2、加深对零、极点分布的概念理解。

二、实验原理

离散系统的时域方程为

其变换域分析方法如下： 频域:

系统的频率响应为:

Z域:

系统的转移函数为:

分解因式:，其中和称为零、极点。

三、预习要求

1.在MATLAB中，熟悉函数tf2zp、zplane、freqz、residuez、zp2sos的使用，其中：[z，p，K]=tf2zp（num，den）求得有理分式形式的系统转移函数的零、极点；zplane（z，p）绘制零、极点分布图；h=freqz(num,den,w)求系统的单位频率响应；[r，p，k]=residuez（num，den）完成部分分式展开计算；sos=zp2sos（z，p，K）完成将高阶系统分解为2阶系统的串联。

2.阅读扩展练习中的实例，学习频率分析法在MATLAB中的实现；

3.编程实现系统参数输入，绘出幅度频率响应和相位响应曲线和零、极点分布图。

四、实验内容

求系统 的零、极点和幅度频率响应和相位响应。

五、扩展练习

例1: 求下列直接型系统函数的零、极点，并将它转换成二阶节形式

解:用MATLAB计算程序如下：

num=[1-0.1-0.3-0.3-0.2];den=[1 0.1 0.2 0.2 0.5];[z,p,k]=tf2zp(num,den);m=abs(p);disp('零点');disp(z);disp('极点');disp(p);disp('增益系数');disp(k);sos=zp2sos(z,p,k);disp('二阶节');disp(real(sos));zplane(num,den)输入到“num”和“den”的分别为分子和分母多项式的系数。计算求得零、极点增益系数和二阶节的系数： 零点: 0.9615-0.5730-0.1443 + 0.5850i-0.1443-0.5850i 极点: 0.5276+0.6997i 0.5276-0.6997i-0.5776+0.5635i-0.5776-0.5635i 增益系数: 1 二阶节: 1.0000-0.3885-0.5509 1.0000 1.15520 0.6511 1.0000 0.28850 0.36300 1.0000-1.0552 0.7679 系统函数的二阶节形式为：

极点图见图:

图 系统函数的零、极点图

例2: 差分方程

所对应的系统的频率响应。

解:差分方程所对应的系统函数为：

用MATLAB计算的程序如下： k=256;num=[0.8-0.44 0.36 0.02];den=[1 0.7-0.45-0.6];w=0:pi/k:pi;h=freqz(num,den,w);subplot(2,2,1);plot(w/pi,real(h));grid title('实部')xlabel('omega/pi');ylabel('幅度')subplot(2,2,2);plot(w/pi,imag(h));grid title('虚部')xlabel('omega/pi');ylabel('Amplitude')subplot(2,2,3);plot(w/pi,abs(h));grid title('幅度谱')xlabel('omega/pi');ylabel('幅值')subplot(2,2,4);plot(w/pi,angle(h));grid title('相位谱')xlabel('omega/pi');ylabel('弧度')

图

实验5 有限冲激响应数字滤波器设计

一、实验目的：

1、加深对数字滤波器的常用指标理解。

2、学习数字滤波器的设计方法。

二、实验原理：

低通滤波器的常用指标：

（1）通带边缘频率；

（2）阻带边缘频率；

（3）通带起伏；

（4）通带峰值起伏，（5）阻带起伏，最小阻带衰减。

三、预习要求

1.在MATLAB中，熟悉函数fir1、kaiserord、remezord、remez的使用；

B = fir1(n,Wn,'high','noscale')设计滤波器；

[n,Wn,beta,ftype] = kaiserord(f,a,dev)估计滤波器阶数；

[n,fo,ao,w] = remezord(f,a,dev,fs)计算等波纹滤波器阶数n和加权函数w(ω)；

B=remez(n,f,a)进行等波纹滤波器的设计。

2.阅读扩展练习中的实例，学习FIR滤波器的设计方法及其在MATLAB中的实现；

3.给出FIR数字滤波器的冲激响应，绘出它们的幅度和相位频响曲线，讨论它们各自的实现形式和特点。

数字滤波器有IIR和FIR两种类型，它们的特点和设计方法不同。

四、实验内容：

利用MATLAB编程,分别用窗函数法和等波纹滤波器法设计两种FIR数字滤波器，指标要求如下：

通带边缘频率：，通带峰值起伏：。

阻带边缘频率：，最小阻带衰减：。

五、扩展练习

例1: 用凯塞窗设计一FIR低通滤波器，通带边界频率，阻带衰减

不小于50dB。，阻带边界频率解: 首先由过渡带宽和阻带衰减

来决定凯塞窗的N和

上图给出了以上设计的频率特性，(a)为N=30直接截取的频率特性(b)为凯塞窗设计的频率特性。凯塞窗设计对应的MATLAB程序为： wn=kaiser(30,4.55);nn=[0:1:29];alfa=(30-1)/2;hd=sin(0.4*pi*(nn-alfa))./(pi*(nn-alfa));h=hd.*wn';[h1,w1]=freqz(h,1);或者：b = fir1(29,0.4,kaiser(30,4.55));[h1,w1]=freqz(b,1);plot(w1/pi,20*log10(abs(h1)));axis([0,1,-80,10]);grid;xlabel('归一化频率/p');ylabel('幅度/dB');还可以使用[n,Wn,beta,ftype]=kaiserord(f,a,dev)函数来估计滤波器阶数等，得到凯塞窗滤波器：

fcuts = [0.3 0.5];%归一化频率omega/pi mags = [1 0];devs = [0.05 10^(-2.5)];[n,Wn,beta,ftype] = kaiserord(fcuts,mags,devs);%计算出凯塞窗N，beta的值 hh = fir1(n,Wn,ftype,kaiser(n+1,beta),'noscale');freqz(hh);实际中，一般调用MATLAB信号处理工具箱函数remezord来计算等波纹滤波器阶数N和加权函数W(ω)，调用函数remez可进行等波纹滤波器的设计，直接求出滤波器系数。函数remezord中的数组fedge为通带和阻带边界 例2：利用雷米兹交替算法设计等波纹滤波器，设计一个线性相位低通FIR数字滤波器，其指标为：通带边界频率fc=800Hz，阻带边界fr=1000Hz，通带波动At=40dB，采样频率fs=4000Hz。

阻带最小衰减解:

在MATLAB中可以用remezord 和remez两个函数设计，其结果如图2，MATLAB程序如下： fedge=[800 1000];mval=[1 0];dev=[0.0559 0.01];fs=4000;[N,fpts,mag,wt]=remezord(fedge,mval,dev,fs);b=remez(N,fpts,mag,wt);[h,w]=freqz(b,1,256);plot(w*2000/pi,20*log10(abs(h)));grid;xlabel('频率/Hz');ylabel('幅度/dB');所得图像如下所示：

实验6 无限冲激响应数字滤波器设计

一、实验目的

1、掌握双线性变换法及脉冲相应不变法设计IIR数字滤波器的具体设计方法;

2、熟悉用双线性变换法及脉冲响应不变法设计低通、高通和带通IIR数字滤波器的计算机编程。

二、实验原理

在MATLAB中，可以用下列函数辅助设计IIR数字滤波器：

1)利用buttord和cheb1ord可以确定低通原型巴特沃斯和切比雪夫滤波器的阶数和截止频率； 2)[num,den]=butter（N,Wn）(巴特沃斯)和[num,den]=cheby1（N,Wn）,[num,den]=cheby2（N,Wn）(切比雪夫1型和2型)可以进行滤波器的设计；

3）lp2hp，lp2bp，lp2bs可以完成低通滤波器到高通、带通、带阻滤波器的转换； 4)使用bilinear可以对模拟滤波器进行双线性变换，求得数字滤波器的传输函数系数； 5)利用impinvar可以完成脉冲响应不变法的模拟滤波器到数字滤波器的转换。

三、预习要求

1.在MATLAB中，熟悉函数butter、cheby1、cheby2的使用，其中：

[num,den]=butter（N,Wn）巴特沃斯滤波器设计； [num,den]=cheby1（N,Wn）切比雪夫1型滤波器设计； [num,den]=cheby2（N,Wn）切比雪夫2型滤波器设计。

2.阅读扩展练习中的实例，学习在MATLAB中进行数字滤波器的设计；

3.给出IIR数字滤波器参数和滤波器的冲激响应，绘出它们的幅度和相位频响曲线，讨论它们各自的实现形式和特点。

四、实验内容

利用MATLAB编程，用脉冲响应不变法和双线性变换法设计一个数字带通滤波器，指标要求如下： 通带边缘频率：阻带边缘频率：，，通带峰值起伏：，最小阻带衰减：

。

五、扩展练习例1：设采样周期T=250μs（采样频率fs =4kHz），用脉冲响应不变法和双线性变换法设计一个三阶巴特沃兹滤波器，其3dB边界频率为fc =1kHz。[B,A]=butter(3,2*pi*1000,'s');[num1,den1]=impinvar(B,A,4000);[h1,w]=freqz(num1,den1);[B,A]=butter(3,2/0.00025,'s');[num2,den2]=bilinear(B,A,4000);[h2,w]=freqz(num2,den2);f=w/pi*2000;plot(f,abs(h1),'-.',f,abs(h2),'-');grid;xlabel('频率/Hz ')ylabel('幅值/dB')程序中第一个butter的边界频率2π×1000，为脉冲响应不变法原型低通滤波器的边界频率；第二个butter的边界频率2/T=2/0.00025，为双线性变换法原型低通滤波器的边界频率.图1给出了这两种设计方法所得到的频响，虚线为脉冲响应不变法的结果；实线为双线性变换法的结果。脉冲响应不变法由于混叠效应，使得过渡带和阻带的衰减特性变差，并且不存在传输零点。同时，也看到双线性变换法，在z=-1即Ω=π或f=2000Hz处有一个三阶传输零点，这个三阶零点正是模拟滤波器在ω=∞处的三阶传输零点通过映射形成的。

下图给出了MATLAB计算的结果。

例2: 设计一数字高通滤波器，它的通带为400～500Hz，通带内容许有0.5dB的波动，阻带内衰减在小于317Hz的频带内至少为19dB，采样频率为1,000Hz。wc=2*1000*tan(2*pi*400/(2*1000));wt=2*1000*tan(2*pi*317/(2*1000));[N,wn]=cheb1ord(wc,wt,0.5,19,'s');[B,A]=cheby1(N,0.5,wn,'high','s');[num,den]=bilinear(B,A,1000);[h,w]=freqz(num,den);f=w/pi*500;plot(f,20*log10(abs(h)));axis([0,500,-80,10]);grid;xlabel('')ylabel('幅度/dB')下图给出了MATLAB计算的结果。

例3： 设计一巴特沃兹带通滤波器，其３dB边界频率分别为f2=110kHz和f1=90kHz，在阻带f3 = 120kHz处的最小衰减大于１０dB，采样频率fs=400kHz。w1=2*400*tan(2*pi*90/(2*400));w2=2*400*tan(2*pi*110/(2*400));wr=2*400*tan(2*pi*120/(2*400));[N,wn]=buttord([w1 w2],[0 wr],3,10,'s');[B,A]=butter(N,wn,'s');[num,den]=bilinear(B,A,400);[h,w]=freqz(num,den);f=w/pi*200;plot(f,20*log10(abs(h)));axis([40,160,-30,10]);grid;xlabel('频率/kHz')ylabel('幅度/dB')下图给出了MATLAB计算的结果。

例4: 一数字滤波器采样频率fs=1kHz，要求滤除100Hz的干扰，其３dB的边界频率为95Hz和105Hz，原型归一化低通滤波器为： w1=95/500;w2=105/500;[B,A]=butter(1,[w1, w2],'stop');[h,w]=freqz(B,A);f=w/pi*500;plot(f,20*log10(abs(h)));axis([50,150,-30,10]);grid;xlabel('频率/Hz')ylabel('幅度/dB')下图为MATLAB的计算结果。

实验7 设计性和研究性实验

设计性实验1 图像信号的抽取与插值

实验目的

1、熟悉图像处理常用函数和方法；

2、培养通过查阅文献解决问题的能力。实验要求

给出一个二维灰度图像，3、编程实现对该图像的任意比例的放大及缩小；

4、编程实现对该图像的任意角度旋转；

5、解决缩放及旋转时产生的锯齿等不图像不平滑问题。实验提示

6、利用上采样、下采样等方法对信号进行缩放变换；

7、观察对图像进行缩放或旋转时，图像是否会出现锯齿等不平滑现象？

8、分析产生锯齿现象的原因；

9、查阅文献了解解决锯齿现象的方法。（例如平滑滤波、双线性插值、双立方插值等处理）

设计性实验2 语音及音乐信号的采样、滤波

实验目的

1、理解采样率和量化级数对语音信号的影响；

2、设计滤波器解决实际问题。实验要求

利用电脑的声卡录一段语音信号及音乐信号，（1）观察使用不同采样率及量化级数所得到的信号的听觉效果，从而确定对不同信号的最佳的采样率；

（2）分析音乐信号的采样率为什么要比语音的采样率高才能得到较好的听觉效果；（3）注意观察信号中的噪声（特别是50hz交流电信号对录音的干扰，设计一个滤波器去除该噪声。实验提示

（1）推荐录音及播放软件：CoolEdit；

（2）分析语音及音乐信号的频谱，根据信号的频率特性理解采样定律对信号数字化的工程指导意义；

（3）可用带阻滤波器对50Hz交流电噪声进行去噪处理；

（4）也可研究设计自适应滤波器对50Hz噪声及其它随机环境噪声进行滤波处理。设计性实验3 双音多频（DTMF）信号的合成和识别

二、实验目的

1、了解电话按键音形成的原理，理解DTMF音频产生软件和DTMF解码算法；

2、利用FFT算法识别按键音；

三、实验要求

（1）设计音频产生函数，音频信号见下图，每个数据信号持续半秒；（2）实现解码函数：接受（1）产生的DTMF信号，识别信号的频率，并生成包含拨号数字的序列；

四、实验提示

（1）DTFT音频可以用两个正弦波按比例叠加产生；

（2）检测算法可以用FFT算法的DFT，或是用一组滤波器实现；

（3）Goertzel算法可以实现调谐滤波器；

设计性实验4 音乐信号处理

五、实验目的

1、了解回声的产生和梳妆滤波器；

2、混音效果的原理和均衡器的设计；

六、实验要求

（1）设计函数实现一段语音或音乐的回声产生；

（2）设计均衡器，使得得不同频率的混合音频信号，通过一个均衡器后，增强或削减某些频率区域，以便修正低频和高频信号之间的关系；

七、实验提示

（1）回声产生可以使用梳妆滤波器，y(n)=x(n)+ax(n-R), a<1(回声

zRH(z),1R1z衰减系数)；或者传输函数为的全通滤波器实现；比较这两种实现方式的区别，分析为什么会有这样的区别；

（2）可以用许多一阶和二阶参数可调的滤波器级联来实现均衡器的功能，滤波器的结构选择结构要求是调整方便，最好调一个参数只影响一个应用指标，且可调参数少；

第五篇：成都电子科技大学章程(征求意见稿).

电子科技大学章程

(征求意见稿)

序

言

电子科技大学原名成都电讯工程学院，是1956年在周恩来总理的亲自部署下，由交通大学（现上海交通大学、西安交通大学）、南京工学院（现东南大学）、华南工学院（现华南理工大学）的电讯工程有关专业合并创建而成，1988年更名为电子科技大学，1997年首批成为国家“211工程”重点建设大学，2001年进入国家“985工程”重点建设大学行列。

电子科技大学以电子信息科学技术为特色，坚持“以人才培养为根本，走内涵式发展道路”的办学理念，大力实施人才强校战略、学科拓展战略、国际化战略，着力提高人才培养质量和科技创新能力，以服务国家、国防和地方经济建设为己任，努力建设电子信息学科领域具有世界先进水平的一流大学。

为促进电子科技大学依法自主办学、自主管理，建立现代大学制度，根据《中华人民共和国教育法》《中华人民共和国高等教育法》等国家相关法律规定，制定本章程。

第一章

总

则

第一条

电子科技大学（以下简称学校），英译文：University of Electronic Science and Technology of China。

学校中文校名字体由郭沫若手稿集字组成（原“成都电讯工程学院”校名由郭沫若手书）。

第二条

学校的校训为：求实求真、大气大为。其字体由毛泽东手稿集字组成。

第三条

学校设有清水河校区、沙河校区、九里堤校区三个校区。学校注册地址为：四川省成都市高新区（西区）西源大道2006号。学校可视需要设立和调整校区及校址。

同创新和自由探索。

第十三条

学校学科专业设置涉及经济学、法学、教育学、文学、理学、工学、医学、军事学、管理学等学科门类。学校可根据经济社会发展、知识创新和科技进步的需要，结合学校战略规划、发展定位，遵循高等教育规律、人才成长规律、科学技术发展趋势和经济社会发展需要，按照国家相关规定，设置、调整学科专业，优化学科专业结构。

第十四条

学校依法颁发学业证书和学位证书。

学校执行国家学位制度，依法授予学士、硕士及博士学位。学校可向国内外卓越的学者或著名的社会活动家授予名誉博士学位。

第十五条

学校徽章由数字“1956”、字母“e”、学校英文简称“UESTC”、中文全称“电子科技大学”字样、深蓝和橘红色交互的两个椭圆等部分构成，整个画面呈圆形。

第十六条

学校校徽为题有校名的长方形证章，教职工校徽为红底白字，本科生校徽为白底红字，研究生校徽为黄底红字。

第十七条

学校确定每年9月29日为校庆日。第十八条

学校网址是http://www.xiexiebang.com。

第二章

学

生

第十九条

学生是指被学校依法录取、取得入学资格，具有学校学籍的受教育者。

第二十条

学生在校期间依法享有下列权利：

（一）参加学校教育教学计划安排的各项活动，使用学校提供的教育教学资源；

（二）参加社会服务、勤工助学，在校内组织、参加学生社团及科技文娱体育等活动；

（三）获得相应的荣誉称号和奖励；

（四）在思想品德、学业成绩等方面获得公正评价，完成学校规定学业后获得相应的学历证书、学位证书；

（五）对学校给予的处分或者处理有异议，向学校、教育行政部

第三章

教职工

第二十九条

学校教职工由教师、其他专业技术人员、管理人员和工勤人员等组成。

学校依法自主聘任各类人员，按照精简、效能的原则设置岗位，依据岗位职责及任职条件公开招聘，择优聘用。

第三十条

学校对教职工定期进行考核，考核结果作为各类人员聘用、晋升、流动、奖惩的依据。

第三十一条

教职工在教育教学、科学研究、管理服务等方面表现优异的，学校予以奖励。

教职工违反学校规章制度的，学校责令其改正或给予相应的纪律处分；考核不符合任职要求的，学校可变更其岗位、解除聘任合同或劳动合同。

第三十二条

学校教职工依法享有下列权利：

（一）进行教育教学活动，开展教育教学改革和实验；

（二）从事科学研究、学术交流，参加专业的学术团队，在学术活动中充分发表意见；

（三）按照工作职责使用学校的公共资源；

（四）在品德、能力和业绩方面获得公平评价，公平获得各类奖励、荣誉称号和自身发展所需的机会和条件；

（五）按时获取工资报酬，享受国家规定的福利待遇以及寒暑假期的带薪休假；

（六）对学校教育教学、管理工作和教育行政部门的工作提出意见和建议；通过教职工代表大会或者其他形式，参与学校的民主管理；

（七）参加进修或者其他形式的培训；

（八）就职务聘用、福利待遇、评优评奖、纪律处分等事项表达异议；

（九）法律法规、本章程及聘约规定的其他权利。第三十三条

学校教职工依法履行下列义务：

（一）遵守宪法、法律和职业道德，为人师表；

（二）遵守学术规范，恪守学术道德；

（二）遵守规章制度，履行聘约，完成教育教学、科学研究、管

用。

（五）按照学习型党组织的要求，加强党员理论学习，坚持用中国特色社会主义理论体系武装头脑，坚定走中国特色社会主义道路的信念。组织党员学习党的路线方针政策和决议，学习党的基本知识，学习科学、文化、法律和业务知识。

（六）领导学校的思想政治工作和德育工作，促进和谐校园建设。

（七）领导学校的工会、共青团、学生会等群众组织和教职工代表大会。

（八）做好统一战线工作。对学校内民主党派的基层组织实行政治领导，支持他们按照各自的章程开展活动。支持无党派人士等统一战线成员参加统一战线相关活动，发挥积极作用。

中国共产党电子科技大学委员会由党员代表大会选举产生，每届任期5年；中国共产党电子科技大学委员会常务委员会由委员会全体会议选举产生。

在中国共产党电子科技大学委员会全体会议闭会期间，其职责由其常务委员会履行。

党委会、党委常委会由党委书记主持。

第三十六条

校长是学校行政负责人，在学校党委领导下全面负责学校的教学、科研和其他管理工作。

学校设校长一人，对外代表学校，对内负责全校校务。学校可设副校长和校长助理等若干人，协助校长行使职权。校长的主要职责是：

（一）拟订发展规划，制定具体规章制度和工作计划并组织实施；

（二）组织教学活动、科学研究和思想品德教育；

（三）拟订内部组织机构的设置方案，推荐副校长人选，任免内部组织机构的负责人；

（四）聘任与解聘教师以及内部其他工作人员，对学生进行学籍管理并实施奖励或者处分；

（五）拟订和执行经费预算方案，保护和管理校产，维护学校的合法权益；

门负责人、纪委副书记，根据议题内容需要列席的其他人员。列席人员由书记或书记委托的副书记确定。

学校党委常委会按照《电子科技大学党委常委会议事决策制度》行使决策权。

第四十二条

校长通过《电子科技大学校长办公会制度》行使决策权。

校长办公会由校长或校长授权的校领导主持。校长办公会的议题，在充分发扬民主的基础上，由校长或授权的校领导作出明确决定。如暂不能作出决定，也应有明确说明。

出席校长办公会的成员为：校长、副校长，学校党委书记、副书记，党委常委，校长助理。列席校长办公会的人员为：学校行政有关部门负责人、监察处处长，根据议题内容需要列席的其他人员。列席人员由校长或校长授权的校领导确定。

第三节

内部组织机构

第四十三条

学校围绕人才培养的根本任务，根据学科专业发展和科学研究需要设置教学科研单位，包括学院（教学部）和具有独立建制的科研院所等。

第四十四条

学校根据工作需要，可设置、变更或者撤销学校的职能部门、直属单位，并可根据实际情况合理调整各部门（单位）的职能。各部门（单位）根据学校的授权，履行管理和服务职责，为全校师生员工提供优质服务。

第四十五条 学校实行校、院（教学部、实验室、中心、所）两级管理体制，教学科研单位是学校组织实施办学活动的基本单位。

第四十六条

学校设立学院。学院在学校授权范围内实行自主管理，开展人才培养、科学研究、学科建设等工作。学校鼓励并支持学院根据自身条件和发展阶段，形成办学特色。

具有独立建制的教学部、实验室、研究院（所、中心）以及根据有关规定由学校授权设立的相应级别的教学科研机构，享有与学院同等的权利，履行相应义务。

第四十七条

学校根据教学、科研的需要以及学校发展的具体情形，可增设、变更、合并或撤销学院、教学部或研究院所等。

（四）领导本单位的思想政治工作。

（五）做好本单位党员干部的教育和管理工作。

（六）领导本单位工会、共青团、学生会等群众组织和教职工代表大会。

第五十四条

党政职能部门、直属单位实行部门领导负责制。第五十五条

学校根据需要可设置各种非常设机构。校级非常设机构的设立和撤销，根据其不同职能，由学校党委常委会或校长办公会讨论决定。

第五十六条

学校依据国家有关规定，与其他高等学校、科研机构、国家机关和企事业单位签订协议，联合设立相关机构，开展合作办学、合作研究与技术开发、社会实践等。

第四节

学术组织

第五十七条

学校设置学术委员会，由学术委员会主任委员主持开展工作。

学术委员会是学校最高学术机构，统筹行使对学术事务的咨询、评定、审议和决策权。

学校充分发挥学术委员会在学科建设、学术评价和学术发展等事项上的重要作用，积极探索教授治学的有效途径，尊重并支持学术委员会独立行使职权，为学术委员会正常开展工作提供必要的条件保障。

学术委员会遵循民主原则，按照少数服从多数集体决策。校学术委员会人员组成由校长办公会讨论决定。学术委员会按照其章程开展工作。

第五十八条

学术委员会可以就具体事项设立若干专门委员会；可以根据需要，在教学科研单位设置学术委员会或者委托教授委员会等基层学术组织承担相应职责。

各专门委员会和教学科研单位学术委员会根据校学术委员会的授权或者各自章程开展工作，向校学术委员会报告工作，接受校学术委员会的指导和监督。

第五十九条

学校设置学位评定委员会。

学位评定委员会负责批准授予、撤销学位及处理学位授予争议等

1生参与学校民主管理和监督的重要组织形式，按照各自章程行使职权、履行职责。

第六十五条

学校设置校务委员会，作为学校的顾问咨询机构，依据其工作条例开展学校建设发展战略研究，对关系学校全局的重大决策开展调查研究和顾问咨询，协助科学决策和民主决策。

第六十六条

校内各民主党派、无党派人士依据法律和各自章程开展活动，参与学校民主管理。

第六十七条

学校依法设置工会、团委等群众组织，各群众组织在学校党委的领导下，依法履行各自的职责。

第五章

外部关系

第六十八条

学校依据国家法律法规及本章程自主管理学校内部事务，不受任何组织和个人的非法干涉。

第六十九条

学校附属的具有独立法人资格的单位，依据法律和学校规定实行相对独立的运营与管理。

第七十条

学校主动接受社会监督和评价，依法实行信息公开制度，及时向社会发布办学信息。

第七十一条

学校利用现代化教育手段和多样化办学机制，开展多种形式的高等学历教育和非学历教育培训，为社会提供多样化的优质教育服务。开展非学历教育培训以不影响学历教育为前提。

第七十二条

学校根据国家需要和自身能力，积极开展面向老少边穷地区的对口支援。

第七十三条

学校积极服务社会，以服务求支持，以贡献求发展，推动协同创新，促进科技成果转化。

学校加强与所在地方的沟通与合作，根据自身条件为所在地方的发展提供服务。

第七十四条

电子科技大学校友包括在电子科技大学学习或工作过的学生和教职工，被学校授予各种名誉学位和荣誉职衔的中外各界人士。

第七十五条

学校依法设立电子科技大学校友总会。

电子科技大学校友总会是学校依法注册成立的非营利性组织，按

3加强对各类国有资产的监管，防止国有资产流失，提高资产使用效率。

学校对占有或使用的全部资产实行校、院分级监管，确保资产的安全完整；学校对外投资等资产使用形式形成的资产应依法合理经营，实现资产的保值增值。

第八十三条

学校设立公共服务机构，为教职工和学生提供图书、档案、网络信息、后勤保障、医疗卫生等服务，保障学校教学、科研、管理等各项工作的顺利开展。

第七章

附

则

第八十四条

学校的分立、合并与终止，由举办者决定。学校终止，本章程自动失效；学校分立、合并，根据规定启动章程修订程序。

第八十五条

本章程经学校教代会讨论、校长办公会审议通过后，由学校党委审定，报教育部核准。

第八十六条

本章程由校长办公会提议修改，章程修订案的审核程序依据第八十五条的规定执行。

第八十七条

学校办公室监督本章程的执行，受理对违反本章程的行为、活动的举报和投诉。

第八十八条

本章程由学校党委负责解释。

第八十九条

本章程经教育部核准后，自公布之日起颁布施行。