资源描述：
《实验一 离散序列卷积》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、学号Xb10680230___班级：10通信（2）班姓名：杨嘉雯日期：2012.10.8__地点：3#304__实验一：离散时间序列卷积和MATLAB实现（一）简述实验目的及实验原理1、掌握离散卷积计算方法；2、学会线性常系数差分方程的迭代解法；3、学会针对具体系统设计程序；（二）实验内容1、用MATLAB画出如下表达式的脉冲序列(1)写出程序代码n=-1:6;%n取值-1到6，间隔为1f=[0，8.0,3.4,1.8,5.6,2.9,0.7，0];stem(n,f);%画出离散图形xlabel(
2、'n');%x坐标轴为nylabel('f(n)';%y坐标轴为f（n）(2)仿真图形2、用MATLAB计算序列{-201–13}和序列{120-1}的离散卷积。首先用手工计算，然后用MATLAB函数conv(x,y)编程验证。人工计算序列x(n)、h(n)的卷积结果并与计算机计算结果相比较，给出笔算和计算机计算的结果对照表，画出图形。(1)写出程序代码xn=[0-201-130];hn=[0120-10];yn=conv(xn,hn);%计算xn和hn的卷积subplot(2,2,1);%画在2
3、行2列的第1个位子nx=-1:5;%n取值-1到5，间隔为1stem(nx,xn);%画出离散图像xlabel('n');%x坐标轴为nylabel('x(n)');%y坐标轴为x（n）subplot(2,2,2);%画在2行2列的第2个位子nh=-1:4;%n取值-1到4，间隔为1stem(nh,hn,'-r');%画出离散图像xlabel('n');%x坐标轴为nylabel('h(n)');%y坐标轴为h（n）subplot(2,2,3);%画在2行2列的第3个位子ny=0:(7+6-2);
4、%n取值0到11，间隔为1stem(ny,yn);%画出离散图像xlabel('n');%x坐标轴为nylabel('y(n)');%y坐标轴为y（n）(2)仿真图形及相关分析3、参考教材P14页，使用通用卷积convu函数，计算如下卷积已知h(n)=R6(N-4)，x(n)=2R3(N-3)，计算y(n)=x(n)*h(n)编程实现并在一幅图内分别画出h(n)，x(n)，y(n)，标好横坐标与纵坐标。(1)写出程序代码h=ones(1,5);nh=-3:1;%hn取6个整数，从-3到1x=2*o
5、nes(1,3);nx=3:5;%xn取3个整数，从3到5[y,ny]=convu(h,nh,x,nx);%调用定义subplot(3,1,1);%画在3行1列的第1个位子stem(nh,h,'b');%画出离散图形，颜色为蓝色xlabel('nh');%x坐标轴为nhylabel('h');%y坐标轴为xsubplot(3,1,2);%画在3行1列的第2个位子stem(nx,x,'b');%画出离散图形，颜色为蓝色xlabel('nx');%x坐标轴为nxylabel('x');%y坐标轴为xs
6、ubplot(3,1,3);%画在3行1列的第3个位子stem(ny,y,'b');%画出离散图形，颜色为蓝色xlabel('ny');%x坐标轴为nyylabel('y');%y坐标轴为yfunction[y,ny]=convu(h,nh,x,nx)%定义nys=nh(1)+nx(1);nyf=nh(end)+nx(end);y=conv(h,x);ny=nys:nyf;(2)仿真图形及相关分析4、编写教材P31第15题的MATLAB程序，标出关键语句的含义，并画出输入和输出图形。(1)写出程序
7、代码xn=[1,2,3,4,2,1];a=[10.5];%差分方程系数b=[102];%差分方程系数ys=0;%初始条件y(-1)=0xi=filtic(b,a,ys);%由初始条件计算等效初始条件的输入序列xiyn=filter(b,a,xn,xi);%调用filter解差分方程，求y(n)n=0:length(yn)-1;%n的范围stem(n,yn,'.');%画离散图像xlabel('n');%x坐标轴为nylabel('y(n)');%y坐标轴为y（n）(2)仿真图形及相关分析5、编写教
8、材P32第18题的MATLAB程序，标出关键语句的含义，并画出输入和输出图形。(1)写出程序代码h1n=[1,1/2,1/4,1/8,1/16,1/32];h2n=[1,1,1,1,1,1];h3n=[1/4,1/2,1/4,zeros(1,97)];%创建一个1行97列的0矩阵h12n=conv(h1n,h2n);%计算h1n与h2n的卷积h12n=[h12n,zeros(1,89)];vn=h12n+h3n;a=[1,-0.9,0.81];b=[1,1,1];hn=filter
