实验一 离散时间信号产生,图形显示及文件读存实验
%% 数字信号处理实验1 内容1
n1=-1;
n2=8;
n0=0;
n3=5;
if((n0<n1)||(n0>n2)||(n1>n2))
error('参数必须满足n1<=n0<=n2')
end
n=n1:n2;
x=[(n-n0)==0];
y=[(n-n3)==0];
figure(1);
subplot(211);
stem(n,x);
xlabel('N');
ylabel('theta(n)');
title('班级 学号 姓名');
legend('单位抽样序列');
subplot(212);
stem(n,y);
xlabel('N');
ylabel('目标序列');
title('班级 学号 姓名');
legend('延时的单位抽样序列');
%% 数字信号处理实验1 内容2
n1=-5;
n2=5;
n0=3;
if((n0<n1)||(n0>n2)||(n1>n2))
error('参数必须满足n1<=n0<=n2')
end
n=n1:n2;
x=[(n+n0)>=0];
figure(2)
stem(n,x);
xlabel('N');
ylabel('目标序列');
title('班级 学号 姓名');
legend('前移3个单位的阶跃序列');
axis([-10 10 -0.5 1.5]);
%% 数字信号处理实验1 内容3
n0=-30;
n2=30;
n=n0:n2;
x=5*exp(3*i*n);
figure(3);
subplot(211);
stem(n,real(x),'.');
title('复指数序列 班级 学号 姓名');
ylabel('实部');
grid;
subplot(212);
stem(n,imag(x),'filled');
title('复指数序列 班级 学号 姓名');
ylabel('虚部');
grid;
k=0:1000;
w=(pi/500)*k;
X=x*exp(-i*pi/500).^(n'*k);
magX=abs(X);
angX=angle(X);
figure(4);
subplot(211);
plot(w/pi,magX);
grid;
xlabel('以\pi为单位的频率');
ylabel('模值');
title('幅频特性 班级 学号 姓名');
subplot(212);
plot(w/pi,angX);
grid;
xlabel('以\pi为单位的频率');
ylabel('弧度');
title('相频特性 班级 学号 姓名');
上图为实验1内容1
实验1内容2
实验1内容3(1)
实验1内容3(2)