MATLAB实二.doc

实验报告格式 院系： 专业 班级： 实验名称： 课程名称： 数字信号处理 实验室： 实验C楼201 成绩 实验人： 学号： 指导教师：赵发勇 实验目的： 1、掌握连续信号的抽样方法，加深对时域采样定理的理解。 2、熟悉时域离散系统的时域特性。 3、掌握常用序列的产生方法。 4、掌握自相关函数的性质和使用方法。 实验内容： xn = 2 4 20 7 n=-1:2 n = -1 0 1 2 需要注意的问题是，在MATLAB中取矩阵或向量的元素的索引号均为正整数，因此，如果需要获得xn=4的值，该怎么取呢？是用xn(0)还是xn(2)呢？请大家试一试。当然如果画此图形当然可以用 stem(n,xn) 来实现。体会两者的区别。 2、序列的加减乘运算 请输入两个长度为5的实序列x=[1 2 3 4 5]和y=[2 4 6 8 10]，对其进行序列运算。 (1) 序列加的MATLAB语句：z=x+y (2) 序列相乘MATLAB语句：z=x.*y (3) 序列元素和的MATALB语句：s=sum(x(n1:n2)) (4) 序列的能量MATALB语句：s=sum(abs(x).^2) (5) 序列的绝对和MATALB语句：s=sum(abs(x(n1:n2))) (6) 序列的右移n0位MATALB语句：y=[zeros(1,n0),x] 对于序列的乘，由于是每个元素对应相乘，因此，在MATLAB中不能用*乘（相量的乘）而必须.*乘（数组的乘）。观察结果。 3、连续信号的抽样，分析抽样序列的特性。例,试选择三种抽样频率进行抽样，观察抽样信号的图形。MATLAB程序如下。 clear all f=50;a=220;n=0:40;% 设置抽样的点数%下面设置抽样周期（对应相应的抽样频率） TS1=0.01;TS2=0.005;TS3=0.0025; %上面共设置了三个抽样周期，下面获得其抽样的离散序列 xn1=a*sin(2*pi*f*n*TS1);xn2=a*sin(2*pi*f*n*TS2);xn3=a*sin(2*pi*f*n*TS3); %用图形的方式显示结果，给同学们一个直观的认识 subplot(311);stem(n,xn1,.); subplot(312);stem(n,xn2,*); subplot(313);stem(n,xn3,-); grid on 类比例，抽样，A,a,f0,fs参数由键盘输入，产生不同的函数和不同的抽样频率，取抽样点为50，并将结果以图形方法显示出来。 采样频率fs分别为100Hz,200HZ和400Hz，观察时间信号的特征。在第三章后可以观察三种采样信号的离散谱。 ．时域离散信号的相关性研究， z=input(随机噪声功率＝);%在MATLAB窗口输入随机信号的功率 a=sqrt(12*z);u=rand(1,100);u=u-mean(u);k=0:99;x=10*cos(0.0008*pi.*k.*k); y=x+u;subplot(311);stem(k,x,b);axis([0,100,-12,12]);title(原信号);xlabel(n);ylabel(x(n)); subplot(312);stem(k,y,b);axis([0,100,-12,12]);title(随机噪声叠加在信号上); xlabel(n);ylabel(y(n)); subplot(313);yc=xcorr(y,100,unbiased);plot(yc); 例练习，其中是白噪声，现要求 (1)、产生均值为0，功率P=0.1的均匀分布白噪声，求的自相关函数；(2)、使的信噪比为10dB，决定A的数值，并画出的图形及其自相关函数的图形。 实验结果分析 序列的加减乘运算 输入两个长度为5的实序列x=[1 2 3 4 5]和y=[2 4 6 8 10]，对其进行序列运算。 x=[1 2 3 4 5] x = 1 2 3 4 5 y=[2 4 6 8 10] y = 2 4 6 8 10 stem(x) 序列默认由1开始 序列加的MATLAB语句：z=x+y z=x+y z = 3 6 9 12 15 (2) 序列相乘MATLAB语句：z=x.*y z=x.*y z = 2 8 18 32 50 序列元素和的MATALB