MATLAB下的数字信号处理实现示例.pdf

MATLAB 下的数字信号处理实现示例

附录一信号、系统和系统响应

1、理想采样信号序列
(1)首先产生信号 x(n),0<=n<=50
n=0:50; %定义序列的长度是 50
A=; %设置信号有关的参数
a=50*sqrt()*pi;
T=; %采样率
w0=50*sqrt()*pi;
x=A*exp(-a*n*T).*sin(w0*n*T); %pi 是 MATLAB 定义的π,信号乘可采用“.*”
close all %清除已经绘制的 x(n)图形
subplot(3,1,1);stem(x); %绘制 x(n)的图形
title(‘理想采样信号序列’);

(2)绘制信号 x(n)的幅度谱和相位谱
k=-25:25;
W=(pi/)*k;
X=x*(exp(-j*pi/)).^(n’*k);
magX=abs(X); %绘制 x(n)的幅度谱
subplot(3,1,2);stem(magX);title(‘理想采样信号序列的幅度谱’);
angX=angle(X); %绘制 x(n)的相位谱
subplot(3,1,3);stem(angX) ; title (‘理想采样信号序列的相位谱’)

(3)改变参数为: A = 1,α= ,Ω0 = ,T = 1
n=0:50; %定义序列的长度是 50
A=1; %设置信号有关的参数
a=;
T=1; %采样率
w0=;
x=A*exp(-a*n*T).*sin(w0*n*T); %pi 是 MATLAB 定义的π,信号乘可采用“.*”
close all %清除已经绘制的 x(n)图形
subplot(3,1,1);stem(x); %绘制 x(n)的图形
title(‘理想采样信号序列’);
k=-25:25;
W=(pi/)*k;
X=x*(exp(-j*pi/)).^(n’*k);
magX=abs(X); %绘制 x(n)的幅度谱
subplot(3,1,2);stem(magX);title(‘理想采样信号序列的幅度谱’);
angX=angle(X); %绘制 x(n)的相位谱
subplot(3,1,3);stem(angX) ; title (‘理想采样信号序列的相位谱’)

2、单位脉冲序列
在 MatLab 中,这一函数可以用 zeros 函数实现:
n=1:50; %定义序列的长度是 50
x=zeros(1,50); %注意:MATLAB 中数组下标从 1 开始
x(1)=1;
close all;
subplot(3,1,1);stem(x);title(‘单位冲击信号序列’);
k=-25:25;
X=x*(exp(-j*pi/)).^(n’*k);
magX=abs(X); %绘制 x(n)的幅度谱
subplot(3,1,2);stem(magX);title(‘单位冲击信号的幅度谱’);
angX=angle(X); %绘制 x(n)的相位谱
subplot(3,1,3);stem(angX) ; title (‘单位冲击信号的相位谱’)

3、矩形序列
n=1:50
x=sign(sign(10-n)+1);
close all;
subplot(3,1,1);stem(x);title(‘单位冲击信号序列’);
k=-25:25;
X=x*(exp(-j*pi/25)).^(n’*k);
magX=abs(X); %绘制 x(n)的幅度谱
subplot(3,1,2);stem(magX);title(‘单位冲击信号的幅度谱’);
angX=angle(X); %绘制 x(n)的相位谱
subplot(3,1,3);stem(angX) ; title (‘单位冲击信号的相位谱’)

4、特定冲击串: x(n) = δ(n) + (n −1) + (n − 2) + δ(n − 3)
n=1:50; %定义序列的长度是 50
x=zeros(1,50); %注意:MATLAB 中数组下标从 1 开始
x(1)=1;x(2)=;x(3)=;x(4)=1;
close all;
subplot(3,1,1);stem(x);title(‘单位冲击信号序列’);
k=-25:25;
X=x*(exp(-j*pi/)).^(n’*k);