【MATLAB】【数字信号处理】基本信号的仿真与实现_阶跃信号matlab仿真

目的

1、用MATLAB软件实现冲激序列

2、用MATLAB软件实现阶跃序列

3、用MATLAB软件实现指数序列

4、用MATLAB软件实现正弦序列

---

内容与测试结果 

常用奇异信号

• 阶跃信号 heaviside(t+k)
• 冲激信号 dirac(t+k)

1、用MATLAB软件实现冲激序列

程序如下：

% 1 冲激序列
clc;    clear all;
n0 = -10; nf = 50; ns = 1; A = 1;%起点为-1，终点为5，幅值为1 在3出有单位冲激
n = n0 : nf;  %生成离散信号的时间序列
y = dirac(n - ns);  % n = ns处有δ函数，即δ(n-ns)
y = A * sign(y);  %改变幅度
subplot(2, 2, 1); stem(n, y, 'r'); axis([n0, nf, -0.1, 1.1]);  %绘制离散序列
title('δ(k-3)');  %加标题

运行结果如下：

2、用MATLAB软件实现阶跃序列

代码如下：

% 2 阶跃序列
clc;
clear all;
n0 = -10; nf = 50; ns = 2;   %起点为-1，终点为5，在3出有单位阶跃序列
n = n0 : nf;  %生成离散信号的时间序列
f = [zeros(1,ns-n0), ones(1,nf-ns+1)]; %生成离散信号f(n)
%也可用逻辑运算方法产生，f=[(n-ns)>=0]
subplot(2, 2, 2);
stem(n, f, 'filled');
axis([n0, nf, -0.1, 1.1]);
title('u(k-3)');  %加标题

运行结果如下：

3、用MATLAB软件实现指数序列

代码如下：

% 3 指数序列
clc;
clear all;
n0 = -10; nf = 50; A = 0.5; nx = 3/4;
k = n0 : nf;  %生成离散信号的时间序列
xk = A * nx.^k;
subplot(2, 2, 3); 
stem(k, xk, 'g');
title('x(k) = 0.5*(3/4)^k');  %加标题

运行结果如下：

4、用MATLAB软件实现正弦序列

代码如下：

% 4 正弦序列
clc;    clear all;
k = -10 : 0.1 : 20; w = pi/3; phi = pi/5;
y = sin(w * k + phi);
subplot(2, 2, 4);
stem(k, y, 'b');
title('sin(π*K/3 + π/5)');

运行结果如下：

将所有图像显示到一个窗口里

% ones(1,5)
% zeros(2,5)
% x = [zeros(1,2),1,zeros(1,3)]
 
 
 
% 1 冲激序列
clc;    clear all;
n0 = -10; nf = 50; ns = 1; A = 1;%起点为-1，终点为5，幅值为1 在3出有单位冲激序列
n = n0 : nf;  %生成离散信号的时间序列
y = dirac(n - ns);  % n = ns处有δ函数，即δ(n-ns)
y = A * sign(y);  %改变幅度
subplot(2, 2, 1); stem(n, y, 'r'); axis([n0, nf, -0.1, 1.1]);  %绘制离散序列
title('δ(k-3)');  %加标题
 
 
 
% 2 阶跃序列
clc;    clear all;
n0 = -10; nf = 50; ns = 2;   %起点为-1，终点为5，在3出有单位阶跃序列
n = n0 : nf;  %生成离散信号的时间序列
f = [zeros(1,ns-n0), ones(1,nf-ns+1)]; %生成离散信号f(n)
%也可用逻辑运算方法产生，f=[(n-ns)>=0]
subplot(2, 2, 2); stem(n, f, 'filled');
axis([n0, nf, -0.1, 1.1]); title('u(k-3)');  %加标题
 
 
 
% 3 指数序列
clc;    clear all;
n0 = -10; nf = 50; A = 0.5; nx = 3/4;
k = n0 : nf;  %生成离散信号的时间序列
xk = A * nx.^k;
subplot(2, 2, 3); 
stem(k, xk, 'g');
title('x(k) = 0.5*(3/4)^k');  %加标题
 
 
 
% 4 正弦序列
clc;    clear all;
k = -10 : 0.1 : 20; w = pi/3; phi = pi/5;
y = sin(w * k + phi);
subplot(2, 2, 4);   stem(k, y, 'b');    title('sin(π*K/3 + π/5)');
 

注：

• clc：清除命令窗口的内容，对工作环境中的全部变量无任何影响
• clear：清除工作空间的所有变量
• stem(n,x,‘r’);函数的第三个入口参数为绘制线的颜色，“r”表示红色，“k”为黑色，默认浅蓝色。

参考文章：

常用序列的MATLAB代码（一）_matlab冲激序列-CSDN博客

网址：http://t.csdnimg.cn/TEhQS