数字信号处理实验实验一离散时间信号与系统及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);。