数字信号处理实验指导书1
- 格式：doc
- 大小：35.50 KB
- 文档页数：4
《数字信号处理》
实验指导书
唐山学院信息与控制工程实验教学中心
DSP实验室
2008年9月
实验一常见离散信号的Matlab产生和图形显示
一．实验目的
1.了解Matlab语言的主要特点，学习简单的数组赋值、运算、绘图等程序设计方法。
2.掌握常用序列的Matlab产生方法。
二．实验内容
1.用MATLAB产生信号序列：
（1）采样信号序列
（2）单位脉冲序列
（3）矩形序列
（4）单位阶跃序列
（5）正弦序列
（6）指数序列
2.序列的运算
（1）序列移位
（2）序列翻转
三．程序设计
1.实验内容：用MA TLAB产生信号序列：
（1）产生采样信号序列
参考程序
n=0:1:49;T=0.001;
A=444;
a=50*sqrt(2)*pi;
W=20*sqrt(2)*pi;
xa=A.*exp(-a*n*T).*sin(W*n*T);
stem(xa,’r’);
title(‘序列xa’);
xlabel(‘n’);
ylabel(‘xa’);
(2) 产生单位脉冲序列、矩形序列
参考程序：
xc=ones(1,10);
ha=ones(1,20);
xb=zeros(1,10);xb(1)=1;
hb=[1,2.5,2.5,1];
subplot(2,2,1);stem(xb,'r');
title('序列 xb');
xlabel('n');
ylabel('xb');
%pause;
subplot(2,2,2);stem(xc,'r');
title('序列 xc');
xlabel('n');
ylabel('xc');
%pause;
subplot(2,2,3);stem(ha,'r');
title('序列 ha');
xlabel('n');
ylabel('ha');
%pause;
title('序列 hb');
xlabel('n');
ylabel('hb');
subplot(2,2,4);stem(hb,'r');
%pause;
（3）单位阶跃序列（参照矩形序列程序，产生单位阶跃序列的程序自己设计）
（4）正弦序列 )/2sin()(ϕπ+=s n f f A n x
N=100;
A=2;
f=50;
fs=1000;%可改变采样频率观察，得到的正弦序列有何变化
fai=0;
n=0:N-1;
x=A*sin(2*pi*f*n/fs+fai);
stem(x,'r');
title('正弦序列 x(n)');
xlabel('n');
ylabel('x');
figure(2)
plot(x);
（5）实指数序列)()(n u a n x n
=
N=10;
%deta=1;
%w=1*pi/3;
n=0:N;
%x=exp((deta+j*w)*n);
a=-1.5;%改变a的值观察序列图形有什么特点
x=a.^n;
stem(x,'r');
title('实指数序列x(n)');
xlabel('n');
ylabel('x');
2.序列运算
（1）序列移位
x=[3,11,7,0,-1,4,2];
nx=[-3:3];
[y,ny]=sigshift(x,nx,2);%右移2位，可以左移2位
subplot(2,1,1);
stem(nx,x);
axis([-4 7 0 11]);
ylabel('x');
subplot(2,1,2);
stem(ny,y);
axis([-4 7 0 11]);
ylabel('y(n)');
（2）序列翻转
clear all
x=[3,11,7,0,-1,4,2];
nx=[-3:3];
yf=fliplr(x);%产生序列x的翻转序列
subplot(2,1,1);
stem(nx,x);
axis([-7 8 0 11]);
subplot(2,1,2);
stem(nx,yf);
axis([-7 8 0 11]);
四．实验报告要求
1.记录产生的采样信号序列、单位脉冲序列和矩形序列的波形，比较与理论结果是否一致。
2.记录自行编写的产生单位阶跃序列的程序及运行结果，比较与理论结果是否一致。
3.改变采样频率观察，得到的正弦序列有何变化，记录并分析结果。
4.改变指数序列a的值，分析序列图形与a的关系。
5.记录原序列x(n)及其移位、翻转后的序列。
