中国传媒大学理工学部
软件课程设计实验报告
系：__通信工程_____ 班级：_____
姓名：_____ 学号：_ ________
项目名称：___ 基于Matlab的通信系统建模与仿真__ ___
? 单人项目 ? 多人项目 多人项目请填写：
小组合作者：______________________________________
实验起止日期 指导老师
2015年 7 月 13 日 —— 7 月 27 日 辅导老师 得 分 第1页共12 页
软件综合课程设计实验报告
模块1名称：_matlab基本功能_____________________
1.1、 模块主要功能： 1.MATLAB的基本元素： 1）数组及数组的初始化； 2）2）一些特殊的矩阵，如zeros(m,n) 、ones(m,n) 、diag(v)、eye(n)。 2.MATLAB的基本运算： 1）数组运算：./ 与 / ； .*与 *； 2）矩阵运算：逆矩阵inv(A)； 3）多项式计算：求值y=ployval(c,x);求根y=roots(p);求系数poly(r)。多项式乘法为对应向量卷积和conv（p1,p2），除法为解卷积deconv(p1,p2);微分polyder(p)，积分polyint(p,k); 4）多项式的插值：y1=interp1(x,y,xi,’method’)； 5）多项式的拟合：p=poly(x,y,n)，y=polyval(p,x)。 3）MATLAB的绘图功能: 1）二维基本绘图：plot(x,y,’属性’);同一窗口多个子图subplot(m,n,k)； 2）图形修饰控制函数：axis,legend,title,xlabel,ylabel,text,grid on/off； 4）三维绘图：plot3(x,y,z,’s’); meshgrid(x,y);mesh(x,y,z,c);surf(x,y,z,c)； 5）其它图形：条形图bar(x);饼状图pie(x);等高线图contour(x,y,z)。 4.MATLAB的编辑功能：1）M文件；2）脚本文件；3）函数文件 1.2、具体任务： 正弦衰减数据：y=sinx*exp(-x/10),其中x=0:pi/5:4*pi,用三次样条法进行插值，并绘制插值前后的曲线。 第2页共12 页
软件综合课程设计实验报告
1.2、 主要内容 插值的意义是：在已知数据中，用较简单的插值函数通过所有样本点，并对临近数据进行估值计算，分为一维插值和二维插值。此处为三次样条法进行插值。使用形式为yi=interp1(x,y,xi,’spline’)其中， x,y为行向量，xi为插值范围内任意点的x坐标，yi为插值运算后的对应y坐标。最后用plot函数将其用不同标记画出来。 1.3、 实现方法： 新建一个m文件，输入代码如下： x0=0:pi/5:4*pi; y0=sin(x0).*exp((-x0)./10); x1=0:pi/30:4*pi; y=sin(x1).*exp((-x1)./10); y1=interp1(x0,y0,x1,'spline'); plot(x1,y1,'b:') hold on plot(x0,y0,'m-') hold on plot(x1,y,'y--') 1.4、 仿真结果
第3页共12 页
软件综合课程设计实验报告
模块2名称： matlab信号的处理________________________
2.1、模块主要功能： 1) 表示常用连续信号； 2）进行信号基本运算； 3）表示常用离散信号； 4）实现离散信号卷积； 5）求解离散系统的单位响应； 6）求解离散系统的零状态响应； 7）求解连续系统的零状态响应； 8) 求解冲激响应及阶跃响应； 9) 实现连续信号卷积； 10)实现连续时间信号傅里叶变换； 11) 分析LTI系统的频域特性； 12) 分析LTI系统的输出响应。 2.2、具体任务： 1、已知差分方程y(k)-y(k-2)=f(k-1) (1)当f(n)=sin(k)u(k)是，求零状态响应y(n) (2)求单位阶跃响应 2、y’’(t)+y(t)=f’(t)求该系统的频率响应，并用matlab绘出其幅频特性和相频特性图。 第4页共12 页
软件综合课程设计实验报告
2.3、主要内容 1、MATLAB 提供求离散系统零状态响应数值解函数 filter，调用格式为 filter(b,a,x),式中 b 和 a 是表示离散系统的向量， x 是输入序列非零样值点行向量，输出向量序号同 x 一样。 stepz(b,a,x0) %求离散系统的单位阶跃响应 2、当系统的频率响应 H（ jw）是 jw 的有理多项式时， MATLAB 信号处理工具箱提供的 freqs 函数可直接计算系统的频率响应的数值解。其调 用格式如下 H=freqs(b,a,w) b=[bM, bM-1,…b0] a=[aN, aN-1,…a0] 其中， a 和 b 分别是 H(jw)的分母和分子多项式的系数向量， w 为形如 w1:p:w2 的向量， 定义系统频率响应的频率范围， w1 为频率起始值， w2 为频率终止值,p 为频率取样间隔。 H 返回 w 所定义的频率点上，系统频率响应的样值。 第5页共12 页