数字信号处理实验及参考程序

数字信号处理实验实验一离散时间信号与系统及MA TLAB实现

1.单位冲激信号：

n = -5:5;

x = (n==0);

subplot(122);

stem(n, x);

2.单位阶跃信号：

x=zeros(1,11);

n0=0;

n1=-5;

n2=5;

n = n1:n2;

x(:,n+6) = ((n-n0)>=0);

stem(n,x);

3.正弦序列：

n = 0:1/3200:1/100;

x=3*sin(200*pi*n+1.2);

stem(n,x);

4.指数序列

n = 0:1/2:10;

x1= 3*(0.7.^n);

x2=3*exp((0.7+j*314)*n);

subplot(221);

stem(n,x1);

subplot(222);

stem(n,x2);

5.信号延迟

n=0:20;

Y1=sin(100*n);

Y2=sin(100*(n-3));

subplot(221);

stem(n,Y1);

subplot(222);

stem(n,Y2);

6.信号相加

X1=[2 0.5 0.9 1 0 0 0 0];

X2=[0 0.1 0.2 0.3 0.4 0.5 0.6 0.7];

X=X1+X2;

stem(X);

7.信号翻转

X1=[2 0.5 0.9 1];

n=1:4;

X2=X1(5-n);

subplot(221);

stem(n,X1);

subplot(222);

stem(n,X2);

8.用MATLAB计算序列{-2 0 1 –1 3}和序列{1 2 0 -1}的离散卷积。a=[-2 0 1 -1 3];

b=[1 2 0 -1];

c=conv(a,b);

M=length(c)-1;

n=0:1:M;

stem(n,c);

xlabel('n');

ylabel('幅度');

9．用MA TLAB计算差分方程

当输入序列为时的输出结果。

N=41;

a=[0.8 -0.44 0.36 0.22];

b=[1 0.7 -0.45 -0.6];

x=[1 zeros(1,N-1)];

k=0:1:N-1;

y=filter(a,b,x);

stem(k,y)

xlabel('n');

ylabel('幅度')

10.冲激响应impz

N=64;

a=[0.8 -0.44 0.36 0.22];

b=[1 0.7 -0.45 -0.6];

x=[1 zeros(1,N-1)];

k=0:1:N-1;

y=impz(a,b,N);

stem(k,y)

xlabel('n');

ylabel('幅度')

11.传递函数频率响应

a=[0.8 -0.44 0.36 0.22];%分子的系数数组

b=[1 0.7 -0.45 -0.6];%分母的系数数组

n=(0:500)*pi/500%在pi范围内取501个采样点

[h,f]=freqz(a,b,n);%求系统的频率响应

subplot(2,1,1),plot(n/pi,abs(h));grid

%作系统的幅度频响图

axis([0,1,1.1*min(abs(h)),1.1*max(abs(h))]);

ylabel('幅度');

subplot(2,1,2),plot(n/pi,angle(h));grid %作系统的相位频响图axis([0,1,1.1*min(angle(h)),1.1*max(angle(h))]);

ylabel('相位');

xlabel('以pi为单位的频率');

12.系统零极点图

a=[0.8 -0.44 0.36 0.22];

b=[1 0.7 -0.45 -0.6];

h=zplane(a,b);