fgp211.m

这是数字信号处理——matlab释义与实现一书的源代码

%《数字信号处理教程——MATLAB释义与实现》第二章演示程序fgp211
% 连续信号和它的采样序列和等价序列
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
dt=0.1;tf=10;t=0:0.1:tf;            % 模拟自变量设定
x=1.5+sin(0.5*t);                   % 连续信号生成
subplot(2,1,1)                      % 画图
plot(t,x),grid on,hold on
T=0.5;n=0:round(tf/T);              % 数字自变量设定
xs=1.5+sin(0.5*n*T);                % 离散信号生成
stem(n*T,xs)                        % 画图
title('连续信号在采样点附近的面积')
% 画出n=12处的样本,左右各延伸0.5T，构成小矩形。
n1=12;xf=[(n1-0.5)*T,(n1-0.5)*T,(n1+0.5)*T,(n1+0.5)*T,(n1-0.5)*T];  % 矩形四个顶角回到第一顶角5点的x坐标
yf=[0,1.5+sin(0.5*(n1-0.5)*T),1.5+sin(0.5*(n1+0.5)*T),0,0];         % 矩形四个顶角回到第一顶角5点的y坐标
yf1=[0,1.5+sin(0.5*n1*T),1.5+sin(0.5*n1*T),0,0];                    % 上缘正弦的5点的y坐标
fill(xf,yf,'r')                                 % 画出矩形，填满红色
subplot(2,1,2)
plot(t,x),grid on,hold on                       % 画出连续曲线
fill(xf,yf1,'g')                                % 画出上缘正弦，填满绿色
stem(n*T,xs)                                    % 画出序列样本图
stairs(n*T,xs)                                  % 绘制序列的阶梯图
title('用采样值乘以采样周期所得的面积')
set(gcf,'color','w')                            % 置图形背景色为白