[2018年最新整理]DSP第一次实验_12151173.docxVIP

数字信号处理实 验 报 告姓名：李森学号院:宇航学院 数字信号处理实验一离散信号的产生及计李森一．实验目的：复习和巩固数字信号处理中离散信号的产生和运算学习和掌握用MATLAB产生离散信号的方法学习和掌握用MATLAB对离散信号进行运算二．实验原理1． 用MATLAB函数产生离散信号信号是数字信号处理的最基本内容。没有信号，数字信号处理就没了工作对象。MATLAB7.0内部提供了大量的函数，用来产生常用的信号波形。例如，三角函数（sin,cos）,指数函数（exp）,锯齿波函数（sawtooth）, 随机数函数（rand）等。 产生被噪声污染的正弦信号用随机数函数产生污染的正弦信号。 产生单位脉冲序列和单位阶跃序列按定义，单位脉冲序列为单位阶跃序列为。 矩形脉冲信号：在MATLAB 中用rectpuls 函数来表示，其调用形式为：y=rectpuls(t,width)，用以产生一个幅值为1，宽度为width,相对于t=0 点左右对称的矩形波信号，该函数的横坐标范围由向量t 决定，是以t=0 为中心向左右各展开width/2 的范围，width 的默认值为1。例：以t=2T(即t-2×T=0)为对称中心的矩形脉冲信号的MATLAB 源程序如下：（取T=1）t=0:0.001:4;T=1;ft=rectpuls(t-2*T,2*T);plot(t,ft);grid on; axis([0 4 –0.5 1.5]); 周期性矩形波（方波）信号在MATLAB 中用square 函数来表示，其调用形式为：y=square(t,DUTY)，用以产生一个周期为2π、幅值为±1 的周期性方波信号，其中的DUTY参数表示占空比，即在信号的一个周期中正值所占的百分比。例如频率为30Hz 的周期性方波信号的MATLAB 参考程序如下：t=-0.0625:0.0001:0.0625;y=square(2*pi*30*t,75);plot(t,y);axis([-0.0625 0.0625 –1.5 1.5]);grid on ; MATLAB中信号的运算 乘法和加法：离散信号之间的乘法和加法，是指它的同序号的序列值逐项对应相乘和相加。三．实验步骤及内容内容：1．产生噪声污染的正弦信号程序：d=rand(1,50)-0.5;x=0:1:49;y=sin(2*pi*x/50)+rand(size(x));plot(x,y);title(噪声污染的正弦信号);xlabel(x);ylabel(y); 2．产生以下信号，并作图（1）幅度为5,频率为3Hz的正弦信号程序：t=0:0.01:5*pi;y=5*sin(2*pi*3*t);plot(t,y);grid on;axis([0 6 -6 6]);title(正弦信号);xlabel(t);ylabel(y);（2）幅度为13，脉宽为2，时间范围为-6~6s的矩形脉冲信号程序：t=-6:0.001:6;y=13*rectpuls(t,2);plot(t,y);grid on;axis([-6 6 -15 15]);plot(t,y);grid on;axis([-6 6 -15 15]);title(矩形脉冲信号);xlabel(t);ylabel(y);（3）幅度为1，频率为4Hz，相位为π/3的锯齿波信号程序：t=0:0.01:1;y=sawtooth(8*pi*t+pi/3);plot(t,y);grid on;axis([0 1 -1 1]);title(锯齿波信号);xlabel(t);ylabel(y);（4）幅度为1，脉冲宽度占空比为30，周期为6的周期方波信号程序：t=0:0.01:12;y=square(2*pi*t/6,30);plot(t,y);grid on;axis([0 12 -1 1]);plot(t,y);grid on;axis([0 12 -2 2]);title(?ü?úD?·?2¨D?o?);xlabel(t);ylabel(y);3．信号的基本运算：（1）作图：y=sin(π*t/5)+4cos(π*t/4)（2）在同一幅图中画出以下三个函数其中，采样频率fs=1000，k为采样点，从1取到500。代码：t=0:0.01:100;y=sin(pi*t/5)+4*cos(pi*t/4);plot(t,y);title(图一);xlabel(t);ylabel(y); 3(2)代码：k=0:500;fs=1000;x1=sin(2*pi*100*k/fs);subplot(3,1,1);plot(k,x1);title(x1(t));xlabel(t);ylabel(x1);x2=sin(2*pi*10*k/fs);subplo