《数字信号处理实验报告》.docVIP

数字信号处理实验报告一、从给定的程序（文件包Friday.rar）中，选择一个源程序做详细标注。（目的：熟悉Matlab程序） 如下：clear;%清空 close all;%关闭所有图形 figure(1);%创建以1为句柄的图形窗口 delta=1;%赋值delta为1 subplot(3,2,1);%将图像分为六个子区域，在第一个区域内绘制图形 plot([-10 -delta],[0 0],-k);%画一条直线，横坐标范围为-10到-1，纵坐标为0，黑色 hold on;%保留图形 plot([-delta -delta],[0 1],-k);%画第二条条直线，横坐标范围为-1，纵坐标为0到1，黑色 hold on;%保留图形 plot([-delta delta],[1 1],-k);%画第三条直线，横坐标范围为-1到1，纵坐标为1，黑色 hold on;%保留图形 plot([delta delta],[1 0],-k);%画第四条直线，横坐标范围为1，纵坐标为1到0，黑色 hold on;%保留图形 plot([delta 10],[0 0],-k);%画第五条直线，横坐标范围为1到10，纵坐标为0，黑色 axis([-10 10 -0.5 1.5]);%坐标轴上横坐标范围为-10到+10，纵坐标为-0.5到1.5 text(-8,0.7,\delta=1,FontSize,10);%在图形（-8，0.7）处标记delta=1，字体大小为10 subplot(3,2,2);%在第二个区域内绘制图形 Dimension1=41;%定义41维 delta_f1=5/delta/(Dimension1-1);%计算delta_f1等于0.125 S1=zeros(1,Dimension1);%生成1行41列全零矩阵S1 f=5/2/delta;%计算f=2.5 plot([-f f],[0 0],--k);%画一条直线，横坐标范围为-2.5到2.5，纵坐标为0，虚线，黑色 hold on;%保留图形 f1=-f:delta_f1:f;%数组f1取值从-2.5到2.5，间距为0.125 N=-(Dimension1+1)/2;%赋N为-21 for i=1:Dimension1k=i+N;t=k*delta_f1;%循环，当i在1到41中时，k=i-21，t=k*0.125if abs(t)0.0001S1(i)=2*delta;%t的绝对值小于0.0001时，S1(i)=2elseS1(i)=sin(2*pi*delta*t)/pi/t;%否则， S1(i)=sin(2*pi*t)/pi/tend end plot(f1,S1,-k);%画图，横坐标照数组f1取值，纵坐标照s1取值，黑色 axis([-f f -0.5 2.5]);%横坐标范围为-2.5到2.5，纵坐标为-0.5到2.5delta=5;%赋值delta为5 subplot(3,2,3);%将图像分为六个子区域，在第3个区域内绘制图形 plot([-10 -delta],[0 0],-k);%，画一条直线，横坐标范围为-10到-5，纵坐标为0，黑色 hold on;%保留图形 plot([-delta 0],[0 1],-k);%画第二条条直线，横坐标范围为-5到0，纵坐标为0到1，黑色 hold on;%保留图形 plot([0 delta],[1 0],-k);%画第三条条直线，横坐标范围为0到5，纵坐标为1到0，黑色 hold on;%保留图形 plot([delta 10],[0 0],-k);%画第四条直线，横坐标范围为5到10，纵坐标为0，黑色 axis([-10 10 -0.5 1.5]);%横坐标范围为-10到10，纵坐标为-0.5到1.5 text(-8,0.7,\delta=5,FontSize,10);%在图形（-8，0.7）处标记delta=5，字体大小为10 subplot(3,2,4);%将图像分为六个子区域，在第四个区域内绘制图形 Dimension2=41;%赋值 delta_f2=6/delta/(Dimension2-1);%计算得delta_f2=0.03 S2=zeros(1,Dimension2);%生成1行41列全零矩阵S2 f=6/2/delta;%f=0.6 plot([-f f],[0 0],--k);%画一条直线，横坐标范围为-0.6到0.6，纵坐标为0，虚线，黑色 hold on;%保留图形 f2=-f:delta_f2:f;%数组f2取值从-0.6到0.6，间距为0.03 N=-(Dimension2+1)/2;%赋N为-21 for i=1:Dimension2k=i+N;t=k*