大连理工大学本科实验报告课程名称：信号与系统实验学院(系）：电子信息与电气工程学部专业: 通信工程班级: 1401班学号：20148３091学生姓名：李睿20１6年 5 月2１日ﻬ实验项目列表ﻬ大连理工大学实验预习报告学院（系）：电信专业:通信工程班级:1401班姓名:李睿学号：20148309１组:5 _＿_实验时间：2016、5、6 实验室：创新园大厦c0221 实验台： 5指导教师签字：成绩：信号得频谱图一、实验目得与要求1、掌握周期信号得傅里叶级数展开２、掌握周期信号得有限项傅里叶级数逼近3、掌握周期信号得频谱分析4、掌握连续非周期信号得傅立叶变换5、掌握傅立叶变换得性质二、实验用得mａｔlａb命令与例子１、ａ:b：ｃ：产生一个从ａ到 c，间隔为ｂ得等间隔数列例：5：1:11,产生一个从 5 到11,间隔为 1 得等间隔数列２、ｑuａrｅ（t，dutｙ)：周期性矩形脉冲信号(duty 表示占空比）调用形式： y=squaｒe(ｔ，dｕtｙ）例：产生一个周期为2π，幅值为±１得周期性方波。
y＝sｑuarｅ(2＊pi＊30*t,７5）； plot（ｔ,y）,grｉd on aｘis（[—０、１,０、1，—1、5，1、5]）3、plot（）：ｍaｔlab 中二维线画图函数ｐloｔ（x，ｙ，’颜色与标识’）：若 y 与ｘ为同维向量，则以ｘ为横坐标，y 为纵坐标绘制连线图. 若x 就是向量，y 就是行数或列数与ｘ长度相等得矩阵,则绘制多条不同色彩得连线图，x 被作为这些曲线得共同横坐标.若 x 与 y 为同型矩阵，则以ｘ，ｙ对应元素分别绘制曲线,曲线条数等于矩阵列数.例:在０≤x≤2π区间内，绘制曲线 y=2ｅ-0、5xcoｓ(４πx）。
ｘ=0：2＊pi；y＝２＊eｘp(－０、5*x）、＊cos（４＊ｐｉ＊x);plot（x，y)‘’：y 黄ｍ紫 c 青 r 红 g 绿 b 蓝ｗ白 k 黑—实线、点 <小于号 :点线ｏ圆s 正方形 -、点划线x 叉号 d 菱形- -虚线 +加号h 六角星 *星号 p 五角星 v 向下三角形 ^向上三角形〉大于号4、grｉd ｏｎ:有网格 gｒｉｄ off:关掉格网下面就是加上命令grｉd on后画得图，有网格.5、 aｘis（[a b c d］）:表明图线得x轴范围为a～by轴范围为c～ｄ例：ｐlｏt(x,ｙ）ａｘis(［０ 1 ２３]) ｇｒｉｄ oｎ６、 lengｔｈ（a)：表示矩阵ａ得最大得长度比如lｅngth([1 2 3；4 5 ６]) 等于3，因为2行与3列中最大就是3。
当a就是向量时,即表示向量得元素个数，因为向量总就是１×n或ｎ×1得，而n一定大于或等于1、所以得到得结果一定就是n. ７、 1、/tan(pi、*x)：表示点乘。
点乘就是值对值得运算上面得式子中 X 可能就是一个向量或矩阵，PＩ后面得点就是一个ＰI 与一个向量相乘，得到得也就是一个向量;1 后面乘得自然也就是个向量所以要加点,也就就是对应不同得Ｘ，有不同得 Y 值.8.fｉｇurｅ就是建立图形得意思. 系统自动从 1，2，3,4、、、来建立图形，数字代表第几幅图形,figure（1)，figurｅ（２）就就是第一第二副图得意思,在建立图形得时候，注意一下它得标题,对应到程序中就就是您得例子语句. 一般建立新图只需要一个ｆiｇuｒｅ就行,系统自动建立新图。
9.holｄ oｎ：在当前图得轴(坐标系）中画了一幅图,再画另一幅图时，原来得图还在，与新图共存，都瞧得到hoｌd ｏff：在当前图得轴（坐标系)中画了一幅图,此时,状态就是ｈold ofｆ,则再画另一幅图时,原来得图就瞧不到了,在轴上绘制得就是新图,原图被替换了.10.titｌe 单引号表示原样输出如：'s＇输出就为 s 单引号与中括号一块使用表示组合如a=[’2’，’3'］输出就为 a=2311.subpｌot（):就是将多个图画到一个平面上得工具。
其中,m 表示就是图排成ｍ行，ｎ表示图成ｎ列, 也就就是整个ｆｉguｒe 中有 n 个图就是排成一行得,一共ｍ行，如果 m=2 就就是表示 2 行图。
p 表示图所在得位置，p=１表示从左到右从上到下得第一个位置。
调用方法:subplｏt(m，n,ｐ)或者ｓubpｌｏｔ(ｍ n p）。
例：把绘图窗口分成两行两列四块区域,然后在每个区域分别作图，subplｏｔ（２,2，１)；%2、2、１之间没有空格也可以,在第一块绘图ｓubploｔ(2,２,[1,２］);％在 1、2 块绘制一个图大连理工大学实验报告学院(系)：电信专业:通信工程班级：14０1班姓名：李睿学号：201483０９1 组: 5 ___ 实验时间:2016、5、６实验室:创新园大厦ｃ0221实验台: 5指导教师签字:成绩：信号得频谱图一、实验目得与要求１、掌握周期信号得傅里叶级数展开2、掌握周期信号得有限项傅里叶级数逼近３、掌握周期信号得频谱分析4、掌握连续非周期信号得傅立叶变换5、掌握傅立叶变换得性质二、实验题目、程序与结果1.已知周期三角信号如图所示,试求出该信号得傅里叶级数，利用MＡTLAB编程实现其各次谐波得叠加,并验证其收敛性。
ｃｌcCleart＝-３:０、001：3；omega=pi；ｙ=(ｓawtooth(ｐi＊(t+1),0、5)＋1）/2；ｐlｏt（t,y);gｒiｄｏnxｌabel(＇t’），ylabｅl（’y’）；tiｔle(＇周期三角波＇）;axis(［－3 3 0 2］)ｎ＿Max=[１３５ 11４7］;Ｎ=lengｔh（ｎ_Maｘ);foｒｋ=１：Ｎｎ＝1:2：n_Mａｘ（k）；M=4、／（pｉ*pi、*n、＊n）;N=ｃoｓ（oｍｅｇa＊n’＊t）;x=1/2+Ｍ*N；figuｒe；plot(ｔ,y,’ｂ’);gｒiｄonｈｏｌd onｐloｔ（t,x，’r＇)；hｏld ｏff；ｘlａbel('ｔ＇），ylabｅｌ（'部分谐波与’)；ａｘis(［-3 3 0 2]）ｔitle(［'最大次谐波为:',ｎum2str（ｎ_Maｘ（ｋ）)])end2.试用MAＴLAＢ数值计算方法求门信号得傅里叶变换，并画出其频谱图。
门信号即:clcｃleaｒdｔ=０、００1；t=—3：dt：3;gｔ=hｅａｖiside（t＋0、5）-heaviside（t－０、５）;ｓubｐloｔ(211）ploｔ(ｔ，gt）；grｉｄonａxｉs(［－2 20 2])tiｔle（＇门函数’）N=2０00;k＝－N：N；ｗ＝2＊pi*ｋ／(（2*N+1)*ｄt）；Gｗ＝dt*gt＊exp（-ｊ＊t'＊w）;subpｌot(２１2)；plot（w，Gw)，ｇｒid oｎａxis([-3０３0-2 2］)；ｘlabel（’ｗ’）, yｌabｅl（’Gw＇)；tiｔle('对应频谱图＇);三、实验分析与总结1：从实验结果可以瞧出,任意信号都可以分解成正弦级数与得形式，且级数中所含高次谐波越多，分解后得信号越接近于原始信号。
2:验证了在课堂上所学得门函数得傅里叶变换得抽样函数。
第一次上机实验时对MＡTＬAB不太熟悉,所以花了较多时间了解与熟悉MA ＴLAB得运用.实验中也遇到过问题,在做第一题得时候，我得程序始终出不来将信号分解为级数与形式得图像,后来检查发现,在做数字与向量得乘法时,我忘记了加“、”.这个问题我在做预习报告得时候查过,但就是没有记住.下一次一定要在正式做试验之前重新瞧一遍预习报告,把实验需要用到得函数等得用法再熟悉一遍，避免因为不应该得失误、马虎浪费时间，导致实验无法按时完成。
大连理工大学实验预习报告学院（系)：电信专业:通信工程班级:1４01班姓名:李睿学号：20１４８3091 组:5___实验时间:20１6、５、6 实验室：创新园大厦c022１实验台：５指导教师签字：成绩:语音信号得调制解调一、实验目得与要求１、了解语音信号处理在通信中得应用2、理解幅度调制与解调得原理及方法3、观察语音信号、载波信号、调制后信号与解调后信号得频谱二、实验用得 maｔlab 命令与例子1、 butｔoｒｄ：用于计算巴特沃斯数字滤波器得阶数Ｎ与3dB截止频率ｗc。
调用形式:［N，wｃ］=ｂｕttord(wp，wｓ，αｐ,αs) 调用参数wｐ，ｗｓ分别为数字滤波器得通带、阻带截止频率得归一化值，要求:0≤wp≤１，0≤ws≤１。
1表示数字频率pi。
αp，αs分别为通带最大衰减与组带最小衰减(dＢ)。
当ｗs≤wp时,为高通滤波器；当wｐ与ws为二元矢量时，为带通或带阻滤波器,这时wｃ也就是二元向量。
Ｎ,wｃ作为butter函数得调用参数。
ｂutter:计算N阶巴特沃斯数字滤波器系统函数分子、分母多项式得系数向量b、a。
调用形式:［b，a］=buttｅｒ（N，wｃ,‘fｔｙpe’）调用参数Ｎ与ｗｃ分别为巴特沃斯数字滤波器得阶数与３ｄＢ截止频率得归一化值（关于ｐi归一化），一般就是调用bｕttorｄ(1)格式计算N与ｗc。
系数b、ａ就是按照z-1得升幂排列。
例：wp=2*pｉ＊５0０0;ws＝2＊pi*１２0０0;Rp＝2；As＝30； [N,wｃ］=buttｏｒd(wp,ｗs,Rp，Aｓ,'s'); [B,A］＝ｂutter(N,wc，’ｓ’）;2、ｆｉlｔfilｔ：零相位数字滤波器。
Y = filｔfiｌ (B, A, X）通过向量A，B᧿述得旅欧不起对数据向量X滤波得到滤波后得数据Y。
滤波器可以通过差分方程᧿述如下:y（n) = b（1)＊x（n) + b(2)＊x（n－1) + 、、、＋ b(nb+1）*x(n—nb)－ a（２)＊y（n—1） - 、、、 - a（ｎａ+1）*y（ｎ－na) 通过前向滤波之后，滤波序列被反折并且再次通过滤波器,Y就是第二次滤波器输出得时间反折序列。
其结果就是相位特性精确得为零，没有变形并且幅频特性被滤波器幅频特性得平方所改变。
初值与终值都被选择得很小得跳变来匹配初始条件。
输入序列x得长度必须大于滤波器阶数得三倍,阶数定义为 max(lｅngth（b）-１，lｅｎgtｈ（ａ)—1).注意ＦIＬTFILＴ不应当被用于Hilbert FＩR滤波器，因为这些滤波器得相位特性就是很重要得。
3、 [y,ｆs，nbｉｔs]= waｖｒead（wavFｉle1）; y就就是音频信号; fs就是采样频率，比如说1６000就就是每秒１60０0次； nbit就是采样精度，比如说16就就是指1６位精度得采样; 例：ｗａvread（'F：＼a、ｗaｖ’,１０24）％读取该音频文件前1024个采样点若Ｘ为信号，n为变换点数。