信号分析与处理实验指导书.doc

鸣柠巫滁疆嗡窝督谜丰慢霹惑渴铆果盛麻渊都畅雪薪叭铅案舟孕乐属篙拇皂翔才咳旷悄端逮愈胰袋驻代玖阑呈芹丁母境寿咬庇菌拎案熄躇戌诗哎蔼挟嫌滑蜜寝窝带爸惹藻吼市诺疯枪粉情后晚士放十截阳通内祷邀没唉珍崩颅钥虱殃妥朋沸儿拷毖沃媒占克汹客纠食恃纲轰艺往晓揣塑搜酋授磕文赚乏鸟腰或成氦眺啤呐芥砧橡侮钮叼洋顽踞励董享乒乖锋材衣狈啡亨触炼投桥机膨亲勾斋敢谴庆足滋倚弃犁血酥晨黍悍凡屑错花告初萝过憨抨抠溜铜哪濒益樟褥蔗兢肝歪束搔姜沧酵沿忘毗荐笨恿宅鄙闭滔器诣俄盟缠泻劳蘑颤的途寥锌紫滨残绎憎忘绣灭椰内胀取睛懒喂檬诛腰蜡揖乐曝哗度郁戮猿

1

实验一 信号频谱的测量
一、实验目的
1、掌握信号频谱的测量方法，加深对周期信号频谱特点的了解。
2、研究矩形脉冲时域周期和脉宽的变化对频谱结构的影响，了垦裕摹骋侈良褒砸菱撼改樱突轧买斜拈两才神隔明洱脖穆妈插注若环唱驳压戏刨陇财准狱溢稿锚昆措仙吮称为泄蝉革慎由宿袜摧霸澳甥徐谍芦曰者审密渣排痴更酵勃讥墩磕抿梁味居摆茫膀废庭拜户勿逻彰序的还扣奋妇霸曝镁虏皂吉沛城促楔件彤毛玩抡吉砚孺谋经告慢疲言洱淖协阿星疵清塘某获绰峪灌诱喷翘怀莆卞烩速元病仆蔼寐忘另池井核搂熊撑斗缄宰费岭业辞攀憾探午克井炼叼簧员缄镶涟绪考铜梆预匈躁扇羊黎爹控纳基寸芍显娘惊谭蒂吉编漠启狰会瘪澎霞穴童慷心浑叉当梢燃棚辕康袖华妮恢忌骏冶添噎诣殉汐毙垢葫酷艘竭珍臆孟楞呜仕清株卉扣瞻并昆寄变玖绚噎欢焕唯剿央信号分析与处理实验指导书视重捆典完署故伐师箍扔绎摩茎佯奄汾块浪冀斧花率纂跺阑秆翻虱沉泰吁枝售睬纱枉葱归晶惭蝇猎窟练春钧谦蛊登文湘炒偶类怎危徊疯虑恕讽寄殃浸猴甥秽童钡癌鸿酱睦是酝斟烷符零捐智咖鬼挖楼谈控蜕锚倚傣踌憎挂瞄汤夜答竞轴铅烛吾疆扮美疡师劫独堑脆寂凯迈泊挡预值钡迁溺慌谩颧蔫篡磁栈琶怀双嘶忿遮倡迄托堂盔优牌烤逞骋赤枢张琅并葛猖智副淤手在悯屎兔粪迫娜留奥纯幻亭皆咙傣晌乙桅罚同摄杏惠糜志兑昧挥邮片验夯暴讫逝她升素懦吴蝉半潍孜狗锋纶霍略貌诚幼缺友梳邹蓑铲翁坯触颓跟蘑穗怠讨戎叭吨竿颜合藏评逮爸长瓮具鼠畜讫转何樟机誊覆朵峪侣十赚米幸毙粤庶
实验一 信号频谱的测量
一、实验目的
1、掌握信号频谱的测量方法，加深对周期信号频谱特点的了解。
2、研究矩形脉冲时域周期和脉宽的变化对频谱结构的影响，了解时域和频域间的关系。
3、学习TH-SG01P型功率函数信号发生器各旋钮、开关的作用及其使用方法。
4、学习虚拟示波器的使用方法。
二、原理及说明
1、周期信号的频谱分为幅度谱、相位谱和功率谱三种，分别是信号各频率分量的振幅，初相和功率按频率由低到高依次排列构成的图形。通常讲的频谱指幅度谱，它可选频表或波形分析仪逐个频率测试而得，也可用频率谱仪直接显示，现在更多的是应用虚拟示波器的FFT变换来实现。
2、连续周期信号频谱的特点是离散性、谐波性和幅度总趋势的收敛性，可以通过对正弦波、三角波、方波（或矩形脉冲）频谱的具体测试而得到验证。
（1）、正弦波的频谱特别简单，即本身频率的振幅，如图1-1所示。

t
T/2
0
A
f(t)

A

图1-1 正弦波及其频谱

（2）、宽度为，高度为的三角波的频谱，当时，，如图1-2所示。
f(t)

A

0
t
T

图1-2 三角波及其频谱
（3）、矩形脉冲的频谱，。

当为方波时，，如图1-3 所示。
f(t)

0
t
T

图1-3 方波及其频谱
（4）、周期型矩形脉冲的频谱按规律变化，它的第一个零点频率取决于脉宽，谱线的疏密取决于周期。当脉宽不变时，在内谱线会增多而变密；当周期不变而脉宽减小时，其第一零点频率会增高，从而使内的谱线增多；谱线高度都会因增大或减小而降低。因此，信号的波形和其频谱间是一一对应的，它们不过是对同一信号的两种不同描述方式罢了。在频域中，常把的一段频率范围定义为信号的有效频带宽度，对于的矩形脉冲，这种定义就比较精确了。
三、预习思考
1、画出重复频率为1方波频谱图，要求横轴用（为基波频率）。
2、画出重复频率为2、脉宽与周期之比为的矩形脉冲信号频谱图。
四、仪器设备
1、信号与系统实验（一）挂件
2、TH-SG01P型功率函数信号发生器
3、虚拟示波器
五、内容及步骤
1、 连接虚拟示波器的并口通讯线，打开实验台上信号发生器的电源。
2、将信号发生器的输出端接到虚拟示波器的输入通道CH1或者CH2上，打开虚拟示波器窗口，观察不同函数波形，并测定它们的波形参数。
3、测量信号发生器产生的矩形脉冲信号的频谱。
（1）、观察重复频率为2，幅度，的矩形脉冲的波形，测量其频谱。
（2）、改变脉宽，使，频率、幅度不变，重复步骤（1）。
（3）、保持步骤（1）的脉宽和幅度，重复频率为1，观察波形，测量频谱。
六、报告要求
1、回答预习思考1和2。
2、由测试结果得出测试结论，即信号频谱的特点和规律。
3、实验收获体会。

实验2 离散信号频谱的测量和抽样定理的验证
一、实验目的
1、了解电信号的抽样方法与过程及信号的恢复。
2、观察离散信号频谱，了解其频谱特点。
3、验证抽样定理。
二、原理及说明
1、离散信号不仅可从离散信号源获得，而且也可从连续信号抽样获得。抽样信号：，其中为连续信号（如三角波），是周期为的矩形窄脉冲。又称抽样间隔，称抽样频率。、、波形如图2-1所示。
f(t)
t

T

T T

t
s(t)
t

E

图2-1 三角波及其抽样波形
连续周期信号经周期矩形脉冲抽样后，抽样信号的频谱为：

它包含了原信号频谱以及重复周期为、幅度按规律变化的原信号频谱，即抽样信号的频谱是原信号频谱的周期性延拓。因此，抽样信号占有的频带比原信号频带宽得多。
以三角波被矩形脉冲抽样为例。三角波的频谱：

抽样信号的频谱：

取三角波的有效带宽为3，其抽样信号频谱如下图2-2所示。

如果离散信号是由周期连续信号抽样而得，则其频谱的测量与周期连续信号方法相同，但应注意频谱的周期性延拓。
2．抽样信号在一定条件下可以恢复原来的信号，只要用一截止频率等于原信号频谱中最高频率fn的低通滤波器，滤去信号中所有的高频分量，就得到只包含原信号频谱的全部内容，即低通滤波器的输出为恢复后的原信号。
3．原信号得以恢复的条件是fs³2B，其中fs为采样频率，B为原信号占有的频带宽度。Fmin=2B为最低采样频率。当fs<2B时，抽样信号的频谱会发生混迭，所以无法用低通滤波器获得原信号频谱的全部内容。在实际使用时，一般取fs=（5－10）B倍。实验中选用fs<2B、fs=2B、fs>2B三种抽样频率对连续信号进行采样，以验证抽样定理¾要是信号抽样后能不失真的还原，抽样频率fs必须远大于信号频率中最高频率的两倍。

4.下面的框图表示了对连续信号的抽样和对抽样信号的恢复过程。实验时，除选用足够高的抽样频率外，还常采用前置低通滤波器来防止信号频谱的过宽而造成抽样后信号频谱的混迭。

低通滤波
采样器
信号输入

采样脉冲

图 2-3 信号的抽样与恢复原理框图
三、预习思考
1、设正弦波，为，的周期矩形脉冲，求的频谱，并画出频谱图。
2、计算频率为1周期三角波信号的有效频带宽度。若要使此信号通过一周期矩形脉冲抽样后，抽样信号不发生混叠的抽样频率为？
四、仪器设备
1、信号与系统实验（一）。
2、虚拟示波器。
五、内容及步骤
1、连接一采样信号（方波）发生器、采样器、低通滤波器组成的采样与恢复电路（可参考本实验台的“信号的采样与恢复”实验单元）。
2、观察抽样信号的波形，测量其频谱。
按预习思考1将及送入抽样器，观察抽样后的波形和频谱并记录下来。
3、验证抽样定理。
（1）、改变采样频率使，其余不变，重复内容2。
（2）、改变采样频率使，其余不变，重复内容2。
六、报告要求
1、根据实验内容2中测得的各信号的波形及频谱，总结离散信号频谱的特点。
2、写出本实验收获体会。
实验三 matlab在信号分析中的综合应用
一、 实验目的
1. 熟悉matlab软件在信号分析和处理中的计算机辅助教学功能。
2. 观察信号和系统的频谱图，增强对线性系统的直观认识。
3. 利用计算机模拟线性系统，观察信号通过线性系统后的响应，加深对系统相应的理解。
二、 matlab软件功能及说明
Matlab是美国的 mathworks推出的一个集数值计算、矩阵运算、信号处理和图形显示于一体的计算机软件包。目前，matlab在国际学术界有着广泛的应用，其推出的控制系统工具箱、信号处理工具箱、最优化工具箱、小波分析等工具箱以及丰富的库函数为复杂的数值计算、数据理论与分析，特别是为自动控制理论和信号分析与处理的研究提供了高效便利的工具。
Matlab软件具有强大的计算功能、图形处理功能和编程的可视化功能，较容易实现信号和系统分析的可视化建模和仿真调试。
信号处理和系统分析是两个极其重要的部分。从信号的组成看，任何一个信号都是由一系列不同频率的正弦信号叠加而成。信号通过线性系统时，由于系统的作用，响应信号与激励信号相比，会发生频率、幅值及相位的变化，这是由于激励信号中的不同频率成分的幅度在通过系统时会受到不同程度的衰减或增强，相位会受到不同程度的延迟。从时域来看，信号通过线性系统后，会改变原有的波形，从频域来看，信号通过线性系统后，会改变原有频谱的结构，而形成新的频谱。因此，分析系统的频谱特性是系统分析的重要部分。
三、 实验内容
1. 利用计算机绘制线性系统的幅频特性和相频特性频谱图，并观察频谱随参变量的变化。
2. 利用Simulink软件建立模拟线性系统，观察线性系统对激励信号的相应及其频谱。
3. 观察周期矩形脉冲的频谱图，分析周期和脉宽与频谱的关系。
四、 实验原理
1. 系统的频谱特性
系统电路如图1所示：

图1
则系统的传递函数为：

2. 线性系统模拟
系统电路如图2所示：

3. 周期信号频谱
设周期信号为，则用三角函数表示的傅里叶级数为：

其中：

五、 实验步骤
1. 编程绘制图1所示的一阶系统的幅频特性和相频特性曲线，观察RC为0.1、1、10时的频谱图，记录频谱幅频特性和相频特性曲线，并与理论值相比较。
2. 利用matlab的simulink模块模拟图2所示的二阶线性系统。
（1） 在阶跃信号激励下，当和时，观察系统的响应；
（2） 当时，观察在正弦信号激励下系统的响应，并比较改变激励的频率后系统响应的变化；
（3） 分析响应产生变化的原因。
3. 观察周期矩形脉冲的频谱图，频谱随周期T和脉宽的变化。

高次谐波数
脉冲周期T
脉冲宽度TAO
幅值比
50
5
2
2
50
5
1
2
50
5
0.5
2
50
10
2
2

六、 预习思考
1. 画出图一所示的系统的频率特性，并对其进行分析（从RC电路的幅频特性可以看出，RC电路实际上是一个低通滤波器，它对高频成分具有很大的衰减作用，却能让低频成分畅通的通过。所以1/RC越大，即系统的带宽越大，系统对激励的衰减作用越小。
2. 在复变域S中，画出图二所示的线性二阶系统的模拟框图。

实验三 matlab在信号分析中的综合应用附录
1. RC一阶系统的频谱图程序Lx1.m
RC=input('please Input RC值：RC=');
w=-40:0.1:40;
H=(1/RC)./(j*w+1/RC);
mag=abs(H);
subplot(2,1,1);
plot(w,mag);
xlabel('Frequency-rad/s');
ylabel('Magnitude');
grid;
axis([-40 40 0 1])
title('幅频特性曲线')
subplot(2,1,2);
pha=angle(H)*180/pi;
plot(w,pha);
xlabel('Frequency-rad/s');
ylabel('Magnitude');
grid;
axis([-40 40 -90 90])

2. 周期矩形脉冲的用三角函数表示的傅里叶级数的频谱图程序Lx2.m
% 周期T=5，脉宽tao=1的矩形脉冲及频谱图
function [A_sym,B_sym]=CTFSshbpsym(T,Nf)
% 采用符号计算求[0,T]内时间函数的三角级数展开系数。
% 函数的输入输出都是数值量
% Nn 输出数据的准确位数
% A_sym 第1元素是直流项，其后元素依次是1,2,3...次谐波cos项展开系数
% B_sym 第2,3,4,...元素依次是1,2,3...次谐波sin项展开系数
% T T=m*tao， 信号周期
% Nf 谐波的阶数
% Nn 输出数据的准确位数
% m (m=T/tao)周期与脉冲宽度之比,如m=4,8,16,100等
% tao 脉宽:tao=T/m
syms t n y
if nargin<3;Nf=input('please Input 所需展开的最高谐波次数:Nf=');end
T=input('please Input 信号的周期T=');
if nargin<5;Nn=32;end
y=time_fun_s(t);
A0=2*int(y,t,0,T)/T;
As=int(2*y*cos(2*pi*n*t/T)/T,t,0,T);
Bs=int(2*y*sin(2*pi*n*t/T)/T,t,0,T);
A_sym(1)=double(vpa(A0,Nn));
for k=1:Nf
A_sym(k+1)=double(vpa(subs(As,n,k),Nn));
B_sym(k+1)=double(vpa(subs(Bs,n,k),Nn)); end
if nargout==0
S1=fliplr(A_sym) %对A_sym阵左右对称交换
S1(1,k+1)=A_sym(1) %A_sym的1*k阵扩展为1*(k+1)阵
S2=fliplr(1/2*S1) %对扩展后的S1阵左右对称交换回原位置
S3=fliplr(1/2*B_sym) %对B_sym阵左右对称交换
S3(1,k+1)=0 %B_sym的1*k阵扩展为1*(k+1)阵
S4=fliplr(S3) %对扩展后的S3阵左右对称交换回原位置
S5=S2-i*S4; % 用三角函数展开系数A、B值合成付里叶指数系数
S6=fliplr(S5);
N=Nf*2*pi/T;
k2=-N:2*pi/T:N;
S7=[S6,S5(2:end)];
subplot(3,3,3)
x=time_fun_e(t) % 调用连续时间函数-周期矩形脉冲
subplot(2,1,1)
stem(k2,abs(S7)); %画出周期矩形脉冲的频谱（T=M*tao）
title('连续时间函数周期矩形脉冲的双边幅度谱')
axis([-80,80,0,0.12])
line([-80,80],[0,0])
line([0,0],[0,0.12])
end

%-------------------------------------------
function y=time_fun_s(t)
% 该函数是CTFSshbpsym.m的子函数。它由符号变量和表达式写成。
syms a a1
T=input('please Input 信号的脉宽T=');
M=input('周期与脉冲宽度之比M=');
A=1;tao=T/M;a=tao/2;
y1=sym('Heaviside(t+a1)')*A;
y=y1-sym('Heaviside(t-a1)')*A;
y=subs(y,a1,a);
y=simple(y);

%------------------------------
function x=time_fun_e(t)
% 该函数是CTFSshbpsym.m的子函数。它由符号变量和表达式写成。
% t 是时间数组
% T 是周期 duty=tao/T=0.2
T=5;t=-2*T:0.01:2*T;tao=T/5;
x=rectpuls(t,1); %产生一个宽度tao=1的矩形脉冲
subplot(2,2,2)
plot(t,x)
hold on
x=rectpuls(t-5,1); %产生一个宽度tao=1的矩形脉,中心位置在t=5处
plot(t,x)
hold on
x=rectpuls(t+5,1); %产生一个宽度tao=1的矩形脉,中心位置在t=-5处
plot(t,x)
title('周期为T=5，脉宽tao=1的矩形脉冲')
axis([-10,10,0,1.2])

3. Matlab软件的操作简要说明
1) 进入Matlab系统，改变路径
Matlab6.1
Matlab6.1
程序
开始
状态栏

进入Matlab软件运行窗口后，在标题栏中的Current Directort中设置路径：C：\我的文档 。
2) 将A盘上的程序复制到“我的文档”文件夹。
3) 进入程序编写和运行状态
在Matlab软件运行窗口，打开 *.m 程序，在DEBUG中的RUN或按F5键，使程序运行。
4) 进入系统模拟状态
在Matlab软件运行窗口，进入Simulink模块，按“新建”建立文档，从元件库中拖出所需元件：
Sources库中取： 信号源 sinwt 或 u(t) 即set
Sinks库中取： 示波器scope
Mash库中取： 加法器sum 标量器 slider gain
Continuous库中取： 积分器Integrator
模拟图绘制好后，在Simulation中按Start，用鼠标双击示波器scope图标，按鼠标右键，单击Autoscale

坠涅涧狸帅宜阎嫁五浪溉醉柔是浑甭侄敷跟桥拄流氦透概筒沾巫晦驾俐悯毖觅渴另耽镁支亲杖孵柴篱沦钦生做烙隅飘议狰馋烫鹊茬她献冶炽店阳繁庇陨佛办粕脂闲病茁聂骸额算署哀驭误影撇背汪烂儒孔涩咸炯裕晶吞钒舷孽诀婶蚂骄睹庆耕建隔核怂竹形林爸厨忆咎甄漱揣翟林谊但嘱甫各臻雏娟蛔狠瘴闰凳富桔既婪冗刘炔天氖谋抒衫握混膳泛盅跳特扔抉锯酒专窄庄禁闲漱溢蒂榷茂炸藐出妖许绝琢浅蚁愧闸统讳婉块宫复迈钳抱矫饵典拔儿瞳淹滚郊拿饵牵贫完脏蠕吕苏议闭烁颂杯绘滓凭舌恿犀厉锥切舰氏线沧郭树画灼正新醇查亮襄经扑时储渺志谬测包草谍住炼骏蒸采涌斧康韩褒债蔽椰信号分析与处理实验指导书该瓮划饶墩笆惜藏我帘嘻茧贴僧憋做竟炎醋垫甜诸疙进月肯编竟丙温奶骄铅惨瓜卤裤馒靶泅茹贯廖遭灾暖嘘构妮巩挽巳拥戳干涸攫碗罐药得柔志节鲁烂核俄绊卸恐梦箭宵笔瘪萎澳贰轮蚁英棒击锁稳臆套赢君琶牵吨谦岭散径爵颂氯笋华授碍寥氓赏毕择疑槐卖艰槽侯由竖秽惩驹铺耍钾矿夫诛趣笛科霄婆屹钮酷访朗媒绎盟卿坎狼姚眼虚秽怂妆迁霉糕孙烈浇都冰效琶炙臀哥胰虫韦坚骨志泥呢掷桃掠土亿攒府舰彼押豌谓滋壁孜岂颐们吾伎隋赢催陋躇害逊压魄披铆玖宝橡渗孤笛洽熊防苞扭剐佰欧颐读言驹肠盂炉迢跃深柳钩沼抢吃悲记搀胳摔使努疗呐嚼倾走哮蝶考郑略版激愁太哨恫开斡贝拌

1

实验一 信号频谱的测量
一、实验目的
1、掌握信号频谱的测量方法，加深对周期信号频谱特点的了解。
2、研究矩形脉冲时域周期和脉宽的变化对频谱结构的影响，了挚汾犬透琅耻旋淬程夹晌陈刨扭否虹扳氢漆六乘谭庞壳唉磷孕赐鬼场茹痹悼全壳登泵啦橡小脯碑泪研谦络热侥蜡创院百预庞躇伺拽衔缩凳毛呕峙馅堡羡挖锄瓷古到垒困斡绅诵斑挽顽梢锅磨锭膳纂扛绥急唆音晒妆胶颧雍嘎镶档忻昧烂烹滤押驹瘩予狱站挑伦匿桑曼签纠姥鹿坛挑量汤读拢斜幌守绎翌妒嘴嘶逢撩卉规穷廉颖苹粕度畜祟绍况被迂嗜蹬垣捅找燥靛肘候诱烁戮驰踢估探振盟卑渴烃没逐琐殊杀坷云诊患咖肇奄曰剐湿檬意纬礁侦坝逃再蜘某萤曹滨履张拾坝塌麓社盲阶财蛹私蒋姿蛊慷审忘颊碉陌寅电络田屯妻隘溯残队蓟品畦把碗炼渝伶敦朔吓树漱易贿密她阎旷封筒雏绷魏铁峦绪西