信号与系统上机实验(11级试用)

《信号与系统》上机实验

实验指导（11级试用）

说明：

上机16学时，每次2学时，共8次。

共7个实验项目，实验7为综合性实验项目，需4学时。

上机时间：

由任课老师、学习委员和实验室管理员共同决定。

实验1 信号变换与系统非时变性质的波形绘制

1、实验目的

掌握MATLAB画连续信号波形的方法，并能绘制信号变换和系统线性非时变性质的输出波形。

2、预习阅读内容

《信号与系统实验教程》实验1、实验2。《信号与系统分析》第1章。

3、参考程序

《信号与系统实验教程》第2章

例2.1-2程序（如何用MATLAB画波形），

例2.2-3程序（如何进行波形变换），

例2.2-4程序（如何进行波形分解）。

《信号与系统分析》第1章

计算机例题C1.1（如何用MATLAB画波形），

计算机例题C1.2（如何进行波形变换），

4、MATLAB学习要点

●用MA TLAB画波形的三要素

（1）横轴赋值，两种方法

t=0:0.1:30; t = 起点：间隔：终点，强调横轴的间隔。

t=linspace(0,30,300);

t有300点，从起点0到终点30,等间隔自动计算，强调横轴的总点数。

（2）纵轴赋值

f=sin(pi*t.*sign(t)); 根据所给函数写出MA TLAB表达式

（3）波形绘制

figure(1)，plot(t,f,'linewidth',2)

●MATLAB的内部函数和外部函数

（1）内部函数：MA TLAB软件内部就存在，调用就可以了，本次介绍几个常用函数

门函数rectpuls（）

调用格式 y=rectpuls(t) 产生高度为1,宽度为1的门函数

调用格式 y=rectpuls(t,W) 产生高度为1,宽度为W的门函数

三角脉冲函数tripuls （）

调用格式 y=tripuls(t) 产生高度为1,宽度为1的三角脉冲函数 调用格式 y=tripuls(t,w) 产生高度为1,宽度为w 的三角脉冲函数 调用格式 y=tripuls(t,w,s) 产生高度为1,宽度为w 的三角脉冲函数,

-1=0);

画图函数myplot( ),函数名为myplot.m ，可以自动将波形放在图的中间，使其美观。 Matlab 程序为

function myplot(x,y)

% x 为横坐标数组,y 为纵坐标数组. x0=x(1);xe=x(end);

max_y=max(y);min_y=min(y);dy=(max_y-min_y)/10; plot(x,y,'linewidth',2);grid; axis([x0,xe,min_y-dy,max_y+dy]) set(gca,'FontSize',8) ● 建立简单函数的两种方法

如)]2()1([)]1()([)(------=t t t t t t f εεεε （1）自定义外部函数 zdyf.m 见例2.2-3程序

function y=zdyf(t)

y=t.*(u(t)-u(t-1))-(u(t-1)-u(t-2)); （2）定义在线函数，见例C1.2程序

zdyf=inline('t.*(u(t)-u(t-1))-(u(t-1)-u(t-2))');

5、上机内容

● 用MA TLAB 画出习题1-8的波形。 ● 用MA TLAB 画出习题1-10的波形。

6、思考题

若已知)3/2-1(t f 的波形如图所示。试画出)(t f 的波形。

实验2 微分方程的符号计算和波形绘制

1、实验目的

掌握MATLAB 解微分方程的方法，并能绘制零输入响应、零状态响应和全响应的波形。

2、预习阅读内容

《信号与系统实验教程》实验4。《信号与系统分析》第2章。

3、参考程序

《信号与系统实验教程》第2章

例2.4-3程序（如何用MA TLAB 的符号计算方法解微分方程并画波形）， 《信号与系统分析》第2章

计算机例题C2.1（如何用MATLAB 的符号计算方法解微分方程），

计算机例题C2.2（如何用MATLAB 的符号计算方法求零输入响应并画波形），

4、MATLAB 学习要点

● 微分方程的符号计算

用Matlab 的符号计算方法dsolve 函数可以计算微分方程的解析式。 其调用格式

r=dsolve('eq1,eq2,...', 'cond1,cond2,...', 'v') r=dsolve('eq1','eq2',...,'cond1','cond2',...,'v')

式中，eq1,eq2,...表示常系数微分方程（组）。cond1,cond2,...表示初始条件。V 表示求解的变量。在缺省情况下，默认变量为t 。

参见参考程序，或《信号与系统实验教程》第40页 ● 符号量与数值量的转换

函数subs （） 参见参考程序，或《信号与系统实验教程》第33页 ● 在一幅图上画多条曲线，并标记图例。

plot(t,y_n,t,yzi_n,'m:',t,yzs_n,'r-.','linewidth',2) legend('全响应','零输入响应','零状态响应',0) 参见参考程序，或《信号与系统实验教程》第20页 ● 微分和积分的符号计算

函数diff （） 《信号与系统实验教程》第33页 函数int （） 《信号与系统实验教程》第34页

5、上机内容

● 用MA TLAB 计算习题2-1，并画出系统响应的波形。 ● 用MA TLAB 计算习题2-3，并画出系统响应的波形。

6、思考题

描述某线性时不变系统的方程为

)(8)(2)(6)(5)(t f t f t y t y t y +'=+'+''

已知初始条件为0)0(',3)0(=-=--y y ，输入)()(t e t f t

ε-=，能否用MA TLAB 的符号计算方

法计算系统全响应的解析解。

实验3 绘制连续系统的四幅图

1、实验目的

用MA TLAB 绘制连续系统的零极点图、冲激响应和阶跃响应图、频率响应和波特图。 四幅图为：零极点图，冲激响应和阶跃响应图（两个响应竖排），频率响应（幅频响应和相频响应竖排），波特图。

2、预习阅读内容

《信号与系统实验教程》实验9。《信号与系统分析》第3章。

3、参考程序

《信号与系统实验教程》实验9

例4.2-1程序（如何用MA TLAB 画零极点图、冲激响应和阶跃响应波形图）， 例4.2-2程序（如何用MA TLAB 画频率响应图和波特图）， 《信号与系统分析》第3章

计算机例题C3.5（如何用MATLAB 画零极点图、冲激响应和阶跃响应波形图），

计算机例题C3.10（如何用MA TLAB 画频率响应图）， 计算机例题C3.11（如何用MA TLAB 画波特图），

4、MATLAB 学习要点

● 冲激响应和阶跃响应的数值计算

Matlab 提供了画系统零极点图的函数。 pzmap( ) 画系统冲激响应的函数。impulse( ) 画系统阶跃响应的函数。step( )

参见参考程序。 《信号与系统实验教程》第107页 ● 频率响应的数值计算

Matlab 提供了专用绘制频率响应的函数。freqs( ) 波特图的函数， bode( ) 参见参考程序。 《信号与系统实验教程》第108页 ● 在一幅图上画多幅子图。subplot( );

参见参考程序，或《信号与系统实验教程》第18页 ● 多项式的乘法。conv( ) , 信号与系统实验教程》第8页

5、上机内容

● 用MA TLAB 计算习题3-28。画出四幅图 ● 用MA TLAB 计算习题3-29。画出四幅图

6、思考题

描述某线性时不变系统的方程为

)(8)(2)(6)(5)(t f t f t y t y t y +'=+'+''

已知初始条件为0)0(',3)0(=-=--y y ，输入)()(t e t f t

ε-=，能否用MA TLAB 的拉氏变换方

法和时域方法计算系统全响应的解析解。

实验4 离散信号波形绘制和迭代法计算

1、实验目的

掌握MATLAB绘制离散信号波形的方法，用迭代法计算差分方程的方法。

2、预习阅读内容

《信号与系统实验教程》实验11、12。《信号与系统分析》第4章。

3、参考程序

《信号与系统实验教程》实验11

例5.1-1程序（如何用MA TLAB画离散信号波形图），

例5.1-2程序（如何用MA TLAB进行波形的变换），

《信号与系统实验教程》实验12

例5.2-1程序（如何用MA TLAB进行迭代法计算），

《信号与系统分析》第4章

计算机例题C4.1（如何用MATLAB画离散信号波形图），

计算机例题C4.2（如何用MATLAB进行波形的变换），

计算机例题C4.3（如何用MATLAB进行迭代法计算），

4、MATLAB学习要点

●画离散信号波形的三要素

与连续信号画图的方法相同，只是画图的函数是 stem( )

参见参考程序。《信号与系统实验教程》第121页

●冲激函数的编写（外部函数）

delta( )

参见参考程序。《信号与系统实验教程》第122页

●循环语句。for….end

参见参考程序，或《信号与系统实验教程》第42页

●数值与字符串的转换。num2str( ) , str2num( ), 《信号与系统实验教程》第27页

●检测数组的长度，length( )。

●特殊矩阵，零矩阵zeros( )等。

5、上机内容

●用MA TLAB计算习题4-4。

●用MA TLAB计算习题4-6。并显示前15项数值，画出其波形。

6、思考题

能否用迭代方法研究离散系统的性质？

用MA TLAB分析习题4-5。

实验5 绘制离散系统的四幅图

1、实验目的

用MA TLAB绘制离散系统的零极点图、冲激响应和阶跃响应图、零输入响应和零状态响应图和频率响应。

四幅图为：零极点图，冲激响应和阶跃响应图（两个响应竖排），零输入响应和零状态响应图（两个响应竖排），频率响应（幅频响应和相频响应竖排）。

2、预习阅读内容

《信号与系统实验教程》实验14。《信号与系统分析》第5章。

3、参考程序

《信号与系统实验教程》实验14

例5.4-1程序（如何用MA TLAB画零极点图、冲激响应和频率响应图），

例5.4-2程序（如何用MA TLAB画全响应图），

例5.4-3程序（如何用MA TLAB画零输入响应和零状态响应和频率响应图），

《信号与系统分析》第5章

计算机例题C5.3（如何用MATLAB画零极点图、冲激响应和阶跃响应波形图），

计算机例题C5.4（如何用MATLAB画频率响应图），

4、MATLAB学习要点

●冲激响应和阶跃响应的数值计算

Matlab提供了画系统零极点图的函数。 zplane( )

画系统冲激响应的函数。impz( )

画系统阶跃响应的函数。stepz( )

参见参考程序。《信号与系统实验教程》第142页

●频率响应的数值计算

Matlab提供了专用绘制频率响应的函数。freqz( )

参见参考程序。《信号与系统实验教程》第143页

●全响应的数值计算。filter( )

参见参考程序，或《信号与系统实验教程》第144页

5、上机内容

●用MA TLAB计算习题5-8。画出四幅图

●用MA TLAB计算习题5-18。画出四幅图

6、思考题

能否用MA TLAB的Z变换方法计算离散系统响应的解析解。如习题5-17，5-22.

实验6 周期信号的分解与合成

1、实验目的

用MA TLAB 绘制周期信号的频谱。理解周期信号的傅里叶级数分解。

根据频谱用MA TLAB 计算前N 项之和，理解周期信号的合成。观察吉布斯现象。

2、预习阅读内容

《信号与系统实验教程》实验5。《信号与系统分析》第6章。

3、参考程序

《信号与系统实验教程》实验5

例3.1-1程序（如何用MA TLAB 实现周期信号的合成）， 例3.1-2程序（如何用MA TLAB 画周期信号的频谱）， 《信号与系统分析》第6章

计算机例题C6.1（画单边频谱图），计算机例题C6.2（用单边频谱合成周期信号）， 计算机例题C6.3（用双边频谱合成周期信号），计算机例题C6.4（画双边频谱图），

计算机例题C6.6（用FFT 画双边频谱图）。

4、MATLAB 学习要点

● 求和的数值计算

用矩阵相乘可实现求和计算。如a=[1 2 3 4 5]; b=[5 4 3 2 1]。

'*5

1

b a b

a i i

i =∑=

参见参考程序。如x=b_n*sin(omega_0*n'*t)

y1=A_n*exp(j*omega_0*n'*t) ● 单边频谱的计算 数值计算：

a(k+1)=2/T*trapz(t,f.*cos(k*w*t)); b(k+1)=2/T*trapz(t,f.*sin(k*w*t)); 符号计算：

an=2/pi*int('exp(-t/2)*cos(2*n*t)',0,pi); bn=2/pi*int('exp(-t/2)*sin(2*n*t)',0,pi);

参见参考程序。 ● 双边频谱的计算

数值计算：FFT 算法 《信号与系统分析》305页 符号计算：

Fn=1/pi*int('exp(-t/2)*exp(-j*2*n*t)',0,pi); 参见参考程序。

5、上机内容

● 用MA TLAB 画出图1、2的频谱和谐波次数N=20时合成的波形。观察N 变化时合成波形和吉布斯现象。其中图1、2是一个周期的波形。

6、思考题

用FFT 函数画出习题6-5、6-6的频谱图。参考《信号与系统分析》402页

图1 图2

实验7 滤波器频域分析与卷积分析

1、实验目的

周期信号输入滤波器的响应可采用频域分析和时域分析方法。通过改变RC 低通滤波器的时间常数或截止频率，观察滤波器的频域和时域性能。先用频域分析、再用卷积分析的方法验证分析的结果。

2、预习阅读内容

《信号与系统实验教程》实验3、7。《信号与系统分析》第2、6章。

3、参考程序

《信号与系统实验教程》实验3、7

例2.3-4程序（如何用MA TLAB 实现卷积积分），

例3.3-1程序（如何用MA TLAB 画滤波器的频率特性）， 例3.3-2程序（如何用MA TLAB 实现频域分析）， 《信号与系统分析》第2、6章 计算机例题C2.4（计算卷积积分），

计算机例题C6.5（用频域方法计算周期信号通过滤波器的响应），

4、上机内容

滤波器为RC 低通滤波器，输入信号为周期锯齿波如图1所示。当RC=0.5、2时，用MA TLAB 分析系统的输入频谱、输出频谱以及系统的时域响应。再用卷积分析的方法画出滤波器输出波形。

5、方法

(1)确定周期信号f(t)的频谱n

F 。基波频率Ω。 (2)确定系统函数)(Ωjn H 和冲激响应)(t h 。 (3)计算输出信号的频谱

n

n F jn H Y )(Ω= (4)系统的时域响应

∑∞

-∞

=Ω=

n t

jn n

e

Y t y )(

MATLAB 计算为

y=Y_n*exp(j*w0*n'*t); (5)系统的时域响应的卷积计算。

)()()(t f t h t y *=

6、要求（画出4幅图）：

(1)在一幅图中画输入信号f(t)和输入信号幅度频谱|F(j ω)|。用多个子图画出。 (2)画出滤波器的幅度频谱|H(j ω)|。

(3)在一幅图中画输出信号y(t)和输出信号幅度频谱|Y(j ω)|。用多个子图画出。 (4)在一幅图中画输入信号f(t)、冲激响应波形h(t)和输出信号y(t)。

信号与系统实验报告1

学生实验报告 （理工类） 课程名称：信号与线性系统专业班级：M11通信工程 学生学号：1121413017 学生姓名：王金龙 所属院部：龙蟠学院指导教师：杨娟

20 11 ——20 12 学年第 1 学期 金陵科技学院教务处制 实验报告书写要求 实验报告原则上要求学生手写，要求书写工整。若因课程特点需打印的，要遵照以下字体、字号、间距等的具体要求。纸张一律采用A4的纸张。 实验报告书写说明 实验报告中一至四项内容为必填项，包括实验目的和要求；实验仪器和设备；实验内容与过程；实验结果与分析。各院部可根据学科特点和实验具体要求增加项目。 填写注意事项 （1）细致观察，及时、准确、如实记录。 （2）准确说明，层次清晰。 （3）尽量采用专用术语来说明事物。 （4）外文、符号、公式要准确，应使用统一规定的名词和符号。 （5）应独立完成实验报告的书写，严禁抄袭、复印，一经发现，以零分论处。 实验报告批改说明 实验报告的批改要及时、认真、仔细，一律用红色笔批改。实验报告的批改成绩采用百分制，具体评分标准由各院部自行制定。 实验报告装订要求

实验批改完毕后，任课老师将每门课程的每个实验项目的实验报告以自然班为单位、按学号升序排列，装订成册，并附上一份该门课程的实验大纲。

实验项目名称：常用连续信号的表示 实验学时： 2学时 同组学生姓名： 无 实验地点： A207 实验日期： 11.12.6 实验成绩： 批改教师： 杨娟 批改时间： 一、实验目的和要求 熟悉MATLAB 软件；利用MATLAB 软件，绘制出常用的连续时间信号。 二、实验仪器和设备 586以上计算机，装有MATLAB7.0软件 三、实验过程 1. 绘制正弦信号)t Asin t (f 0?ω+=（），其中A=1，πω2=，6/π?=； 2. 绘制指数信号at Ae t (f =），其中A=1，0.4a -=； 3. 绘制矩形脉冲信号，脉冲宽度为2； 4. 绘制三角波脉冲信号，脉冲宽度为4；斜度为0.5； 5. 对上题三角波脉冲信号进行尺度变换，分别得出)2t (f ，)2t 2(f -； 6. 绘制抽样函数Sa （t ）,t 取值在-3π到+3π之间； 7. 绘制周期矩形脉冲信号，参数自定； 8. 绘制周期三角脉冲信号，参数自定。 四、实验结果与分析 1．制正弦信号)t Asin t (f 0?ω+=（），其中A=1，πω2=，6/π?= 实验代码： A=1;

信号与系统实验指导书

信号与系统实验指导书 通信教研室 辽宁工业大学 2009年8月

目录 实验一时域卷积积分- 1 - 实验二连续系统时域分析- 10 - 实验三离散系统时域分析- 13 - 实验四二阶低通滤波器的响应- 16 -

1 MATLAB 在信号与系统中的应用 1.1.1 实验目的 (1) 练习连续信号的产生。 (2) 练习傅里叶变换的分析。 (3) 练习连续系统分析。 1.1.2 实验仪器 计算机、MATLAB 软件环境。 1.1.3 实验内容 在下面的实验操作中，认真保存、记录每项操作的作用和目的。 （一） 练习连续信号的产生 已知连续信号()()sin()j t x t e t αω+=+。要求编写程序文件siggen.m ，完成以下功能： (1) 在0≤t ≤5之间，产生该信号。其中0.6,5αω=-=。 (2) 在3个子图上分别画出该信号、信号的实部和虚部，并对图形进行标注。 此外， (3) 将(2)中产生的图形文件以bmp 格式保存到桌面。 （二） 练习连续傅里叶变换的分析 已知信号12()sin(2)2cos(2)s t f t f t ππ=+，其中f 1=47Hz ，f 2=88Hz 。要求： (1) 在0≤t ≤5范围内，步长增量为0.001，求出该信号的傅里叶变换； (2) 在2个子图上，分别绘制该信号的波形和幅频、相频响应曲线图。 （三） 练习连续系统分析 某LTI 系统输入信号为信号1110()0 t u t ≤≤?=? ?其它 ，系统的冲激响应为0.2()t h t e -=，

长度为15。要求： (1) 在时间间隔为0.5前提下，完成系统的卷积计算； (2) 在2个子图上，绘出输入信号和输出信号曲线，并进行标注。 1.1.4 实验预习 (1) 读懂各例题实验程序，了解MA TLAB基本操作方法。 (2) 根据实验内容预先编写实验程序。 1.1.5 实验报告 (1) 列写实验内容和已调试通过的实验程序，并按实验记录完成实验报告，打印实验程序产生的曲线图形。 (2) 思考题 ①若通信信号由有用信号和信道噪声组成，该信号该如何产生？ ②连续系统分析的方法

浙江大学 信号与系统实验-基础实验

本科实验报告 课程名称：信号与系统实验 姓名：Wzh 院系：信电学院 专业：信息工程 学号：xxxxxxx 指导教师：周绮敏、史笑兴、李惠忠 2017年6月 1 日 Copyright As one member of Information Science and Electronic Engineering Institute of Zhejiang University, I sincerely hope this will enable you to acquire more time to do whatever you like instead of struggling on useless homework. All the content you can use as you like. I wish you will have a meaningful journey on your college life. ——W z h

实验报告 课程名称：信号与系统实验指导老师：史笑兴、周绮敏、李惠忠成绩：__________________ 实验名称：实验一MATLAB基本实验实验类型：设计型 一、第一次基本实验 1、利用Matlab自带的sinc函数，在时间区间[-4,4]上产生sinc信号，并画出信号图形。 2、利用./运算符，在时间区间[ -4*pi , 4*pi ]上产生Sa信号，并画出信号图形。 具体要求： （1）将图形窗口分为上下两部分，sinc信号画在上图，Sa信号画在下图。 （2）对两个信号分别设置合适的坐标显示范围。 【思考题】sinc函数与Sa函数二者的关系为何？用表达式表示。 【代码】 【运行结果】

西电数字信号处理上机实验报告

数字信号处理上机实验报告 14020710021 张吉凯 第一次上机 实验一： 设给定模拟信号()1000t a x t e -=，t 的单位是ms 。 (1) 利用MATLAB 绘制出其时域波形和频谱图（傅里叶变换），估计其等效带宽（忽略谱分量降低到峰值的3%以下的频谱）。 (2) 用两个不同的采样频率对给定的()a x t 进行采样。 ○1()()15000s a f x t x n =以样本秒采样得到。 ()()11j x n X e ω画出及其频谱。 ○2()()11000s a f x t x n =以样本秒采样得到。 ()() 11j x n X e ω画出及其频谱。 比较两种采样率下的信号频谱，并解释。 （1）MATLAB 程序： N=10; Fs=5; T s=1/Fs; n=[-N:T s:N]; xn=exp(-abs(n)); w=-4*pi:0.01:4*pi; X=xn*exp(-j*(n'*w)); subplot(211) plot(n,xn); title('x_a(t)时域波形'); xlabel('t/ms');ylabel('x_a(t)'); axis([-10, 10, 0, 1]); subplot(212); plot(w/pi,abs(X)); title('x_a(t)频谱图'); xlabel('\omega/\pi');ylabel('X_a(e^(j\omega))');

ind = find(X >=0.03*max(X))*0.01; eband = (max(ind) -min(ind)); fprintf('等效带宽为%fKHZ\n',eband); 运行结果： 等效带宽为12.110000KHZ

信号系统实验报告

电子工程系 信号与系统课程实验报告 2011-----2012学年第一学期 专业: 电子信息工程技术班级: 学号 : 姓名: 指导教师: 实常用连续时间信号的实现

一、实验目的 （1）了解连续时间信号的特点； （2）掌握连续时间信号表示的向量法和符号法； （3）熟悉MATLAB Plot函数等的应用。 二、实验原理 1、信号的定义 信号是随时间变化的物理量。信号的本质是时间的函数。 2、信号的描述 1）时域法 时域法是将信号表示成时间的函数f(t)来对信号进行描述的方法。信号的时间特性指的是信号的波形出现的先后，持续时间的长短，随时间变化的快慢和大小，周期的长短等。 2）频域（变换域）法 频域法是通过正交变换，将信号表示成其他变量的函数来对信号进行描述的方法。一般常用的是傅立叶变换。信号的频域特性包括频带的宽窄、频谱的分布等。 信号的频域特性与时域特性之间有着密切的关系。 3、信号的分类 按照特性的不同，信号有着不同的分类方法。 （1）确定性信号：可以用一个确定的时间函数来表示的信号。 随机信号：不可以用一个确定的时间函数来表示，只能用统计特性加以描述的信号。 （2）连续信号：除若干不连续的时间点外，每个时间点在t上都有对应的数值信号。离散信号：只在某些不连续的点上有数值，其他时间点上信号没有定义的信号。 （3）周期信号：存在T，使得等式f(t+T)=f(t)对于任意时间t都成立的信号。非周期信号：不存在使得等式f(t+T)=f(t)对于任意时间t都成立的信号。 绝对的周期信号是不存在的，一般只要在很长时间内慢走周期性就可以了。 （4）能量信号：总能量有限的信号。 功率信号:平均功率有限切非零的信号。 （5）奇信号：满足等式f(t)=--f(--t)的信号。偶信号：满足等式f(t)=f(--t)的信号。 三、涉及的MATLAB函数 1、plot函数 功能：在X轴和Y轴方向都按线性比例绘制二维图形。 调用格式： Plot（x,y）:绘出相x对y的函数线性图。 Plot（x1,y1,x2,y2,…..）:会出多组x对y的线性曲线图。 2、ezplot函数 功能：绘制符号函数在一定范围内的二维图形。简易绘制函数曲线。 调用格式： Ezplot (fun):在[-2π，2π]区间内绘制函数。 Ezplot (fun，［min，max］)：在［min,max］区间内绘函数。 Ezplot (funx,funy):定义同一曲面的函数，默认的区间是[0, 2π]。】 3、sym函数 功能：定义信号为符号的变量。 调用格式：sym(fun):fun为所要定义的表达式。 4、subplot函数

信号与系统实验总结及心得体会

信号与系统实验总结及心得体会 2011211204 刘梦颉2011210960 信号与系统是电子信息类专业的一门重要的专业核心基础课程，该课程核心的基本概念、基本理论和分析方法都非常重要，而且系统性、理论性很强，是将学生从电路分析领域引入信号处理与传输领域的关键性课程，为此开设必要的实验对我们加强理解深入掌握基本理论和分析方法，以及对抽象的概念具体化有极大的好处，而且为后续专业课程的学习提供了理论和大量实验知识储备，对以后的学术科研和创新工作都是十分重要的。下面我将从实验总结、心得体会、意见与建议等三方面作以总结。 一．实验总结 本学期我们一共做了四次实验，分别为：信号的分类与观察、非正弦周期信号的频谱分析、信号的抽样与恢复（PAM）和模拟滤波器实验。 1.信号的分类与观察 主要目的是：观察常用信号的波形特点以及产生方法，学会用示波器对常用波形参数进行测量。主要内容是：利用实验箱中的S8模块分别产生正弦信号、指数信号和指数衰减正弦信号，并用示波器观察输出信号的波形，测量信号的各项参数，根据测量值计算信号的表达式，并且与理论值进行比较。 2.非正弦信号的频谱分析 主要目的是：掌握频谱仪的基本工作原理和正确使用方法，掌握非正弦周期信好的测试方法，理解非正弦周期信号频谱的离散性、谐波性欲收敛性。主要内

容是：通过频谱仪观察占空比为50%的方波脉冲的频谱，和占空比为20%的矩形波的频谱，并用坐标纸画图。 3.信号的抽样与恢复 主要目的是：验证抽样定理，观察了解PAM信号的形成过程。主要内容是：通过矩形脉冲对正弦信号进行抽样，再把它恢复还原过来，最后用还原后的图形与原图形进行对比，分析实验并总结。 4.模拟滤波器实验 主要目的是：了解RC无源和有源滤波器的种类、基本结构及其特性，比较无源和有源滤波器的滤波特性，比较不同阶数的滤波器的滤波效果。主要内容：利用点频法通过测试无源低通、高通、带通和有源带阻，以及有源带通滤波器的幅频特性，通过描点画图形象地把它们的特点表现出来。 通过对信号与实验课程的学习，我掌握了一些基本仪器的使用方法，DDS 信号源、实验箱、示波器、频谱仪等四种实验仪器。初步了解了对信号的测试与分析方法对以前在书本上看到的常见信号有了更加具体的认识，使得书本上的知识不再那么抽象。 DDS信号源，也就是函数发生器，可以产生固定波形，如正弦波、方波或三角波，频率和幅度可以调节。实验箱是很多个信号实验装置的集合，可谓集多种功能于一身，其中包括函数发生器、模拟滤波器、函数信号的产生与测量、信号的抽样与恢复等模块。示波器能把抽象的电信号转换成具体的图像，便于人们研究各种电现象的变化过程。利用示波器能观察各种不同的信号幅度随时间变化的波形曲线，还可以用它测试各种不同的电量，如电压、电流、频率、相位差、

数字信号处理上机实验代码

文件名：tstem.m（实验一、二需要） 程序： f unction tstem(xn,yn) %时域序列绘图函数 %xn:被绘图的信号数据序列，yn:绘图信号的纵坐标名称（字符串）n=0:length(xn)-1; stem(n,xn,'.'); xlabel('n');ylabel('yn'); axis([0,n(end),min(xn),1.2*max(xn)]); 文件名：tplot.m（实验一、四需要） 程序： function tplot(xn,T,yn) %时域序列连续曲线绘图函数 %xn:信号数据序列,yn:绘图信号的纵坐标名称（字符串） %T为采样间隔 n=0;length(xn)-1;t=n*T; plot(t,xn); xlabel('t/s');ylabel(yn); axis([0,t(end),min(xn),1.2*max(xn)]); 文件名：myplot.m（实验一、四需要）

%(1)myplot;计算时域离散系统损耗函数并绘制曲线图。function myplot(B,A) %B为系统函数分子多项式系数向量 %A为系统函数分母多项式系数向量 [H,W]=freqz(B,A,1000) m=abs(H); plot(W/pi,20*log10(m/max(m)));grid on; xlabel('\omega/\pi');ylabel('幅度(dB)') axis([0,1,-80,5]);title('损耗函数曲线'); 文件名：mstem.m（实验一、三需要） 程序： function mstem(Xk) %mstem(Xk)绘制频域采样序列向量Xk的幅频特性图 M=length(Xk); k=0:M-1;wk=2*k/M;%产生M点DFT对应的采样点频率（关于pi归一化值） stem(wk,abs(Xk),'.');box on;%绘制M点DFT的幅频特性图xlabel('w/\pi');ylabel('幅度'); axis([0,2,0,1.2*max(abs(Xk))]); 文件名：mpplot.m（实验一需要）

信号与系统实验报告_1(常用信号的分类与观察)

实验一：信号的时域分析 一、实验目的 1.观察常用信号的波形特点及产生方法 2.学会使用示波器对常用波形参数的测量 二、实验仪器 1.信号与系统试验箱一台（型号ZH5004） 2.40MHz双踪示波器一台 3.DDS信号源一台 三、实验原理 对于一个系统特性的研究，其中重要的一个方面是研究它的输入输出关系，即在一特定的输入信号下，系统对应的输出响应信号。因而对信号的研究是对系统研究的出发点，是对系统特性观察的基本手段与方法。在本实验中，将对常用信号和特性进行分析、研究。 信号可以表示为一个或多个变量的函数，在这里仅对一维信号进行研究，自变量为时间。常用信号有：指数信号、正弦信号、指数衰减正弦信号、复指数信号、Sa(t)信号、钟形信号、脉冲信号等。 1、信号：指数信号可表示为f(t)=Ke at。对于不同的a取值，其波形表现为不同的形式，如下图所示： 图1―1 指数信号 2、信号：其表达式为f（t）=Ksin（ωt＋θ），其信号的参数：振幅K、角频率ω、与初始相位θ。其波形如下图所示：

图1－2 正弦信号 3、指数衰减正弦信号：其表达式为其波形如下图： 图1－3 指数衰减正弦信号 4、Sa(t)信号：其表达式为：。Sa(t)是一个偶函数，t= ±π,±2π,…,±nπ时，函数值为零。该函数在很多应用场合具有独特的运用。其信号如下图所示：

图1－4 Sa（t）信号 5、钟形信号（高斯函数）：其表达式为：其信号如下图所示： 图1－5 钟形信号 6、脉冲信号：其表达式为f(t)=u(t)-u(t-T)，其中u(t)为单位阶跃函数。其信号如下图所示： 7、方波信号：信号为周期为T，前T/2期间信号为正电平信号，后T/2期间信号为负电平信号，其信号如下图所示 U(t)

数字信号处理上机实验(第三版)

数字信号处理实验（Matlab） 实验一: 系统响应及系统稳定性 %实验1：系统响应及系统稳定性 close all;clear all %======内容1：调用filter解差分方程，由系统对u(n)的响应判断稳定性====== A=[1,-0.9];B=[0.05,0.05]; %系统差分方程系数向量B和A x1n=[1 1 1 1 1 1 1 1 zeros(1,50)]; %产生信号x1(n)=R8(n) x2n=ones(1,128); %产生信号x2(n)=u(n) hn=impz(B,A,58); %求系统单位脉冲响应h(n) subplot(2,2,1);y='h(n)';tstem(hn,y); %调用函数tstem绘图 title('(a)系统单位脉冲响应h(n)');box on y1n=filter(B,A,x1n); %求系统对x1(n)的响应y1(n) subplot(2,2,2);y='y1(n)';tstem(y1n,y); title('(b)系统对R8(n)的响应y1(n)');box on y2n=filter(B,A,x2n); %求系统对x2(n)的响应y2(n) subplot(2,2,4);y='y2(n)';tstem(y2n,y); title('(c)系统对u(n)的响应y2(n)');box on %===内容2：调用conv函数计算卷积============================ x1n=[1 1 1 1 1 1 1 1 ]; %产生信号x1(n)=R8(n) h1n=[ones(1,10) zeros(1,10)]; h2n=[1 2.5 2.5 1 zeros(1,10)];

信号与系统实验报告

实验三 常见信号的MATLAB 表示及运算 一、实验目的 1．熟悉常见信号的意义、特性及波形 2．学会使用MATLAB 表示信号的方法并绘制信号波形 3. 掌握使用MATLAB 进行信号基本运算的指令 4. 熟悉用MATLAB 实现卷积积分的方法 二、实验原理 根据MATLAB 的数值计算功能和符号运算功能，在MA TLAB 中，信号有两种表示方法，一种是用向量来表示，另一种则是用符号运算的方法。在采用适当的MA TLAB 语句表示出信号后，就可以利用MA TLAB 中的绘图命令绘制出直观的信号波形了。 1.连续时间信号 从严格意义上讲，MATLAB 并不能处理连续信号。在MATLAB 中，是用连续信号在等时间间隔点上的样值来近似表示的，当取样时间间隔足够小时，这些离散的样值就能较好地近似出连续信号。在MATLAB 中连续信号可用向量或符号运算功能来表示。 ⑴ 向量表示法 对于连续时间信号()f t ，可以用两个行向量f 和t 来表示，其中向量t 是用形如12::t t p t =的命令定义的时间范围向量，其中，1t 为信号起始时间，2t 为终止时间，p 为时间间隔。向量f 为连续信号()f t 在向量t 所定义的时间点上的样值。 ⑵ 符号运算表示法 如果一个信号或函数可以用符号表达式来表示，那么我们就可以用前面介绍的符号函数专用绘图命令ezplot()等函数来绘出信号的波形。 ⑶ 常见信号的MATLAB 表示 单位阶跃信号 单位阶跃信号的定义为：10()0 t u t t >?=? 0); %定义函数体，即函数所执行指令

数字信号处理上机实验答案(全)1

第十章 上机实验 数字信号处理是一门理论和实际密切结合的课程，为深入掌握课程内容，最好在学习理论的同时，做习题和上机实验。上机实验不仅可以帮助读者深入的理解和消化基本理论，而且能锻炼初学者的独立解决问题的能力。本章在第二版的基础上编写了六个实验，前五个实验属基础理论实验，第六个属应用综合实验。 实验一 系统响应及系统稳定性。 实验二 时域采样与频域采样。 实验三 用FFT 对信号作频谱分析。 实验四 IIR 数字滤波器设计及软件实现。 实验五 FIR 数字滤波器设计与软件实现 实验六 应用实验——数字信号处理在双音多频拨号系统中的应用 任课教师根据教学进度，安排学生上机进行实验。建议自学的读者在学习完第一章后作实验一；在学习完第三、四章后作实验二和实验三；实验四IIR 数字滤波器设计及软件实现在。学习完第六章进行；实验五在学习完第七章后进行。实验六综合实验在学习完第七章或者再后些进行；实验六为综合实验，在学习完本课程后再进行。 10.1 实验一: 系统响应及系统稳定性 1.实验目的 （1）掌握 求系统响应的方法。 （2）掌握时域离散系统的时域特性。 （3）分析、观察及检验系统的稳定性。 2.实验原理与方法 在时域中，描写系统特性的方法是差分方程和单位脉冲响应，在频域可以用系统函数描述系统特性。已知输入信号可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应，本实验仅在时域求解。在计算机上适合用递推法求差分方程的解，最简单的方法是采用MA TLAB 语言的工具箱函数filter 函数。也可以用MATLAB 语言的工具箱函数conv 函数计算输入信号和系统的单位脉冲响应的线性卷积，求出系统的响应。 系统的时域特性指的是系统的线性时不变性质、因果性和稳定性。重点分析实验系统的稳定性，包括观察系统的暂态响应和稳定响应。 系统的稳定性是指对任意有界的输入信号，系统都能得到有界的系统响应。或者系统的单位脉冲响应满足绝对可和的条件。系统的稳定性由其差分方程的系数决定。 实际中检查系统是否稳定，不可能检查系统对所有有界的输入信号，输出是否都是有界输出，或者检查系统的单位脉冲响应满足绝对可和的条件。可行的方法是在系统的输入端加入单位阶跃序列，如果系统的输出趋近一个常数（包括零），就可以断定系统是稳定的[19]。系统的稳态输出是指当∞→n 时，系统的输出。如果系统稳定，信号加入系统后，系统输出的开始一段称为暂态效应，随n 的加大，幅度趋于稳定，达到稳态输出。 注意在以下实验中均假设系统的初始状态为零。 3．实验内容及步骤 （1）编制程序，包括产生输入信号、单位脉冲响应序列的子程序，用filter 函数或conv 函数求解系统输出响应的主程序。程序中要有绘制信号波形的功能。 （2）给定一个低通滤波器的差分方程为

信号与系统课程设计应用MATLAB实现连续信号的采样与重构仿真

设计题目应用MATLAB实现连续信号的采样与重构仿真 1、设计目的 信号与系统课程设计是学习《信号与系统》课程必要的教学环节。由于该课程是专业基础课，需要通过实践了巩固基础知识，为使学生取得最现代化的设计技能和研究方法，课程设计训练也就成为了一个重要教学环节。通过对信号与系统一书的重新认识，我们将学习如何利用MATLAB软件进行仿真与重构并加深对滤波器的理解，这样的课程设计出了对我们的学习起着只关重要的作用，还可以很好的培养我们自己的动手能力。本次课程设计，我们会引入一个模拟的信号，通过MATLAB软件的防真技术来实现对它的分析、理解与学习。 MATLAB软件是今年来比较长用的一种数学软件，它有很强大的功能，主要侧重于某些理论知识的灵活运用。本次课程设计的目的是：增加对仿真软件MATLAB的感性认识，熟悉MATLAB软件平台的使用和MATLAB编程方法及常用语句；、初步掌握MATLAB的编程方法和特点；加深理解采样与重构的概念，应用MATLAB编程实现对信号的采样与重构；分别计算在临界采样、过采样、欠采样三种不同条件下重构信号的误差，并由此总结采样频率对信号重构误差的影响；学生需要自拟题目，根据自己手中的资料独立思考与分析，明确实习内容，制定实习步骤与方案，独立完成作业。 2、原理说明 2.1.1MATLAB MATLAB是美国Math Works公司产品，MATLAB现已被广泛于数学、通信、信号处理、自动控制、神经网络、图形处理等许多不同学科的研究中。并越来越多的应用到我们的学习生活中来，是目前通信工程上最广泛应用的软件之一。最初的MATLAB 只是一个数学计算工具。但现在的MATLAB已经远不仅仅是一个“矩阵实验室”，它已经成为一个集概念设计、算法开发、建模仿真，实时实现于一体的集成环境，它拥有许多衍生子集工具。 沈阳大学

四川大学信号与系统第一次实验报告(题目二)

周期信号? ??<<-≤<-=21,5.110,5.0)(t t t t t x ，周期T=2. (1)写出x(t)的复指数形式和三角函数形式傅里叶级数表示； ?????=∴-=-+-===-+-=====? ?????---为奇，为偶，k 20)1(-1])5.1()5.0([21)(x 10])5.1()5.0([21)(x 12 22222221102110000π ππππωππωk k a k dt e t dt e t dt e t T a dt t dt t dt t T a T k k t jk t jk T t jk k T ) k cos(2)(x )(x 1jk t a t e a t n k t n k ππ∑∑+∞=+∞-∞=== ， (2)利用“分析公式”或傅里叶级数性质求出谱系数； 由（1）知，?????=为奇，为偶，k 2022πk k a k （3）编程，对x(t)进行频谱分析，具体要求： 程序如下： set(gcf,'color','w') %设置背景颜色为白色 fs=128; %采样频率fs=128Hz tp = 1/fs; %采样时间间隔 N = fs*6*4; %采样点数，总采样时常为4秒 n = -N:N-1; %采样点序列 t = n * tp; %采样点时间序列 x=-0.5*sawtooth(pi*t,0.5); % 产生信号x(t) subplot(311); plot(t,x); %画信号的时域波形

axis([-6,6,-0.5,0.5]); %规定坐标轴的取值范围 xlabel('时间(s)'); %横轴的名称、单位 ylabel('时域'); %纵轴的名称、单位 title('x(t)'); %图的名称 %%%%%%以下对信号进行FFT变换%%%%%% Nf = 512; %做512点的FFT y=fft(x,Nf)/Nf; %进行fft变换——复指数形式谱系数mag=abs(y); %求幅度谱 theta = angle(y)/pi*180; %求相位谱 %%%%%%修正幅度谱（三角函数形式谱系数）%%%%%% delta_1 = [1, 2*ones(1,length(y)-1)]; mag = mag.*delta_1; %修正后的幅度谱 %%%%%%修正相位谱%%%%%% delta_2= (mag>0.01); %判别式，利用逻辑运算实现 % 将幅值为0的频率分量的相位置为0” theta = theta.*delta_2; %修正后的相位谱 f=(0:Nf-1)'*fs/Nf; %进行对应的频率转换 subplot 312 bar(f, mag, 0.1); %画幅度谱 axis([0,5,0,0.5]);xlabel('频率(Hz)');ylabel('幅值'); subplot 313 bar(f, theta, 0.1); %画相位谱 axis([0,20,-200,200]);xlabel('频率(Hz)');ylabel('相位'); %%%%%%以下将主要频谱分量叠加，实现信号的重建%%%%%% i=2; %定义循环变量

数字信号处理上机报告-一

数字信号处理上机报告-一

数字信号处理第一次上机实验报告 实验一： 设给定模拟信号()1000t a x t e -=，的单位是ms 。 (1) 利用MATLAB 绘制出其时域波形和频谱图（傅里叶变换），估计其等效带宽（忽略谱分量降低到峰值的3%以下的频谱）。 (2) 用两个不同的采样频率对给定的进行采样。 ○1 。 ○2 。 比较两种采样率下的信号频谱，并解释。 实验一MATLAB 程序： （1） N=10; Fs=5; Ts=1/Fs; n=[-N:Ts:N]; xn=exp(-abs(n)); w=-4*pi:0.01:4*pi; X=xn*exp(-j*(n'*w)); subplot(211) plot(n,xn); title('x_a(t)时域波形'); xlabel('t/ms');ylabel('x_a(t)'); t ()a x t ()()15000s a f x t x n =以样本秒采样得到。()()11j x n X e ω画出及其频谱()()11000s a f x t x n =以样本秒采样得到。()() 11j x n X e ω画出及其频谱

axis([-10, 10, 0, 1]); subplot(212); plot(w/pi,abs(X)); title('x_a(t)频谱图'); xlabel('\omega/\pi');ylabel('X_a(e ^(j\omega))'); ind = find(X >=0.03*max(X))*0.01; eband = (max(ind) -min(ind)); fprintf('等效带宽为 %fKHZ\n',eband); 运行结果：

《信号与系统》课程教学大纲

《信号与系统》课程教学大纲 课程编码：A0303051 总学时：64 理论学时：64 实验学时：0 学分：4 适用专业：通信工程 先修课程：电路，高等数学，复变函数与积分变换，线性代数 一、课程的性质与任务 《信号与系统》是电类专业的一门重要的专业课程。它的任务是研究信号和线性非时变系统的基本理论和基本分析方法，要求掌握最基本的信号变换理论，并掌握线性非时变系统的分析方法，为学习后续课程，以及从事相关领域的工程技术和科学研究工作奠定坚实的理论基础。通过本课程的学习，学生将理解信号的函数表示与系统分析方法，掌握连续时间系统和离散时间系统的时域分析和频域分析，连续时间系统的S域分析和离散时间系统的Z域分析，以及状态方程与状态变量分析法等相关内容。通过实验，使学生掌握利用计算机进行信号与系统分析的基本方法，加深对信号与线性非时变系统的基本理论的理解，训练学生的实验技能和科学实验方法，提高分析和解决实际问题的能力。

二、课程学时分配 教学章节理论实践 第一章：信号与系统导论6 第二章：连续系统的时域分析8 第三章：信号与系统的频域分析18 第四章：连续系统的复频域分析10 第五章：系统函数的零、极点分析8 第六章：离散系统的时域分析6 第七章：离散系统的Z域分析8 总计64 三、课程的基本教学内容及要求 第一章信号与系统导论（6学时） 1.教学内容 （1）历史的回顾，应用领域，信号的概念 （2）系统的概念，常用的基本信号 （3）信号的简单处理，单位冲激函数 2.重点及难点 教学重点：信号的描述、阶跃信号与冲激信号；信号的运算；线性时不变系统判据；系统定义 教学难点：信号及其分类，信号分析与处理，系统分析 3.课程教学要求

信号与系统实验一

实验一 基本运算单元 一、 实验目的 1．熟悉由运算放大器为核心元件组成的基本运算单元； 2．掌握基本运算单元的测试方法。 二、 实验设备与仪器 1．THKSS-A/B/C/D/E 型信号与系统实验箱； 2．实验模块SS12； 3．双踪示波器。 三、 实验内容 1．设计加法器、比例运算器、积分器、微分器四种基本运算单元电路； 2．测试基本运算单元特性。 四、 实验原理 1．运算放大器 运算放大器实际就是高增益直流放大器，当它与反馈网络连接后，就可实现对输入信号的求和、积分、微分、比例放大等多种数学运算，运算放大器因此而得名。运算放大器的电路符号如图1-1所示： 图1-1 运算放大器的电路符号 由图可见，它具有两个输入端和一个输出端：当信号从“-”端输入时，输出信号与输入信号反相，因此称“-”端为反相输入端；而从“+”端输入时，输出信号与输入信号同相，因此称“+”端为同相输入端。运算放大器有以下的特点： （1）高增益 运算放大器的电压放大倍数用下式表示： )1(0 + --= u u u A 式中，u o 为运放的输出电压；u +为“+”输入端对地电压；u -为“-”输入端对地电压。不加反馈（开环）时，直流电压放大倍数高达104～106。 （2）高输入阻抗 运算放大器的输入阻抗一般在106Ω～1011Ω范围内。 （3）低输出阻抗 运算放大器的输出阻抗一般为几十到一、二百欧姆。当它工作于深度负反馈状态时，其闭环输出阻抗更小。 为使电路的分析简化，人们常把上述的特性理想化，即认为运算放大器的电压放大倍数和输入阻抗均为无穷大，输出阻抗为零。据此得出下面两个结论： 1）由于输入阻抗为无穷大，因而运放的输入电流等于零。

数字信号处理上机实验答案完整版

数字信号处理上机实验 答案 HEN system office room 【HEN16H-HENS2AHENS8Q8-HENH1688】

第十章上机实验 数字信号处理是一门理论和实际密切结合的课程，为深入掌握课程内容，最好在学习理论的同时，做习题和上机实验。上机实验不仅可以帮助读者深入的理解和消化基本理论，而且能锻炼初学者的独立解决问题的能力。本章在第二版的基础上编写了六个实验，前五个实验属基础理论实验，第六个属应用综合实验。 实验一系统响应及系统稳定性。 实验二时域采样与频域采样。 实验三用FFT对信号作频谱分析。 实验四 IIR数字滤波器设计及软件实现。 实验五 FIR数字滤波器设计与软件实现 实验六应用实验——数字信号处理在双音多频拨号系统中的应用 任课教师根据教学进度，安排学生上机进行实验。建议自学的读者在学习完第一章后作实验一；在学习完第三、四章后作实验二和实验三；实验四IIR数字滤波器设计及软件实现在。学习完第六章进行；实验五在学习完第七章后进行。实验六综合实验在学习完第七章或者再后些进行；实验六为综合实验，在学习完本课程后再进行。 functiontstem(xn,yn) %时域序列绘图函数 %xn:信号数据序列，yn:绘图信号的纵坐标名称（字符串） n=0:length(xn)-1; stem(n,xn,'.');boxon xlabel('n');ylabel(yn); axis([0,n(end),min(xn),*max(xn)]) 实验一: 系统响应及系统稳定性 1.实验目的 （1）掌握求系统响应的方法。 （2）掌握时域离散系统的时域特性。 （3）分析、观察及检验系统的稳定性。 2.实验原理与方法 在时域中，描写系统特性的方法是差分方程和单位脉冲响应，在频域可以用系统函数描述系统特性。已知输入信号可以由差分方程、单位脉冲响应或系统函数求出系统对于该输入信号的响应，本实验仅在时域求解。在计算机上适合用递推法求差分方程的解，最简单的方法是采用MATLAB语言的工具箱函数filter函数。也可

MATLAB实验教程

目录实验一、MATLAB 基本应用 实验二、实验二信号的时域表示 实验三、连续信号卷积 实验四、典型周期信号的频谱表示 实验五、傅立叶变换性质研究 实验六、系统的零极点分析 实验七离散信号分析

实验一 MATLAB 基本应用 一、实验目的：学习MATLAB的基本用法，了解 MATLAB 的目录结构和基本功能以及MATLAB在信号与系统中的应用。 二、实验内容： 例一已知x的取值范围，画出y=sin(x)的图型。 参考程序：x=0:0.05:4*pi; y=sin(x); plot(y) 例二计算y=sin(π/5)+4cos(π/4) 例三已知z 取值范围，x=sin(z);y=cos(z);画三维图形。 z=0:pi/50:10*pi; x=sin(z); y=cos(z); plot3(x,y,z)

xlabel('x') ylabel('y') zlabel('z') 例四已知x的取值范围，用subplot函数绘图。 参考程序：x=0:0.05:7; y1=sin(x); y2=1.5*cos(x); y3=sin(2*x); y4=5*cos(2*x); subplot(2,2,1),plot(x,y1),title('sin(x)') subplot(2,2,2),plot(x,y2),title('1.5*cos(x)') subplot(2,2,3),plot(x,y3),title('sin(2*x)') subplot(2,2,4),plot(x,y4),title('5*cos(2*x)')

单位阶跃信号 u(t+3)-2u(t) 其中 u(t)=1/2+(1/2)sign(t) Sigh(t)是符号函数t>0时为1，t<0时为-1； clear all %利用maple中的heaviside函数来实现u(t+3)-2u(t) syms t f=maple('Heaviside(t+3)-2*Heaviside(t)'); t1=-3*pi:0.01:3*pi;%注意精度问题 ff=subs(f,t,t1); figure(1); plot(t1,ff); axis([-5,5,-1.2,1.2]); title('u(t+3)-2u(t)'); %利用自己编写的heaviside函数来实现u(t+3)-2u(t) t=-5:0.01:5; f=Heaviside(t+3)-2*Heaviside(t); figure(2); plot(t,f); axis([-5,5,-1.2,1.2]); title('u(t+3)-2u(t)') %利用符号函数来实现u(t) t=-5:0.05:5; f=sign(t); ff=1/2+1/2*f; figure(3); plot(t,ff); axis([-5 5 -0.1 1.1]); title('u(t)') %利用符号函数来实现u(t+3)-2u(t) t=-5:0.01:5; f=(1/2+1/2*sign(t+3))-2*(1/2+1/2*sign(t)); figure(4) plot(t,f),axis([-5,5,-1.2,1.2]),title('u(t+3)-2u(t)')

数字信号处理上机实验答案(第三版)

实验1：系统响应及系统稳定性 实验程序清单: close all;clear all %======容1：调用filter解差分方程，由系统对u(n)的响应判断稳定性====== A=[1,-0.9];B=[0.05,0.05]; %系统差分方程系数向量B和A x1n=[1 1 1 1 1 1 1 1 zeros(1,50)]; %产生信号x1(n)=R8(n) x2n=ones(1,128); %产生信号x2(n)=u(n) hn=impz(B,A,58); %求系统单位脉冲响应h(n) subplot(2,2,1);y='h(n)';stem(hn, 'y'); %调用函数tstem绘图 title('(a) 系统单位脉冲响应h(n)'); y1n=filter(B,A,x1n); %求系统对x1(n)的响应y1(n) subplot(2,2,2);y='y1(n)';stem(y1n,'y'); title('(b) 系统对R8(n)的响应y1(n)'); y2n=filter(B,A,x2n); %求系统对x2(n)的响应y2(n) subplot(2,2,4);y='y2(n)';stem(y2n,'y'); title('(c) 系统对u(n)的响应y2(n)'); %===容2：调用conv函数计算卷积============================ x1n=[1 1 1 1 1 1 1 1 ]; %产生信号x1(n)=R8(n) h1n=[ones(1,10) zeros(1,10)]; h2n=[1 2.5 2.5 1 zeros(1,10)]; y21n=conv(h1n,x1n); y22n=conv(h2n,x1n); figure(2) subplot(2,2,1);y='h1(n)';stem(h1n,'y'); %调用函数tstem绘图 title('(d) 系统单位脉冲响应h1(n)'); subplot(2,2,2);y='y21(n)'; stem(y21n,'y'); title('(e) h1(n)与R8(n)的卷积y21(n)'); subplot(2,2,3);y='h2(n)'; stem(h2n, 'y'); %调用函数tstem绘图 title('(f) 系统单位脉冲响应h2(n)'); subplot(2,2,4);y='y22(n)';stem(y22n,'y'); title('(g) h2(n)与R8(n)的卷积y22(n)'); %=========容3：谐振器分析======================== un=ones(1,256); %产生信号u(n) n=0:255; xsin=sin(0.014*n)+sin(0.4*n); %产生正弦信号 A=[1,-1.8237,0.9801];B=[1/100.49,0,-1/100.49]; %系统差分方程系数向量B和A y31n=filter(B,A,un); %谐振器对u(n)的响应y31(n) y32n=filter(B,A,xsin); %谐振器对u(n)的响应y31(n) figure(3) subplot(2,1,1);y='y31(n)';stem(y31n,'y'); title('(h) 谐振器对u(n)的响应y31(n)'); subplot(2,1,2);y='y32(n)';stem(y32n,'y'); title('(i) 谐振器对正弦信号的响应y32(n)');