MATLAB常用信号实现（抽样序列，阶跃序列，矩形序列，指数序列，余弦，虚指数）

单位抽样序列

% 单位抽样序列和延时的单位抽样序列
n=0:10;
x1=[1 zeros(1,10)];x2=[zeros(1,5) 1 zeros(1,5)];%列值实现
subplot(1,2,1);%第一个页面
stem(n,x1);%画离散信号的函数
xlabel ('时间序列n');
ylabel('振幅');
title('单位抽样序列x1');
subplot(1,2,2);%第二个页面显示
stem(n,x2); 
xlabel('时间序列n');
ylabel('振幅');
title('延时了5的单位抽样序列');

%先在matlab中定义单位抽样序列：
function [x,n]=Deltafun(n1,n2,n0)
n=(n1:n2);
x=(n==n0);
stem(n,x);
xlabel('n');
title('单位抽样序列x');

[x,n]= Deltafun (-5,5,3)

单位阶跃序列

k=0:10;
u=[ones(1,11)];
stem(k,u);
xlabel ('时间序列n');
ylabel('振幅');
title('单位阶跃序列');

%先在matlab中定义单位抽样序列：
function [x,k]= Stepfun (n1,n2,n0)
k=(n1:n2);
x=(k>=n0);
stem(k,x);
xlabel('k');
title('单位阶跃序列x');

[x,n]=Stepfun(-5,5,3)

矩形序列实现

function [x,k]= Box (n1,n2)
k=[n1:n2];
x=(k>=n1); % x=(and(k>=n1,k<=n2));
stem(k,x);
xlabel('k');
title('矩形序列x');

[x,n]= Box (-5,5)

指数序列

k=1:30;
x=1.2.^k;
stem(k,x); 
xlabel('时间序列k');
ylabel('振幅');
title('指数序列x=1.2.^k');

余弦序列

k=1:30;
x=2*cos(pi*k/6+pi/3);
stem(k,x); 
xlabel('时间序列k');
ylabel('振幅');
title('正弦函数序列x=2*cos(pi*k/6+pi/3)');

正弦序列改符号即可。

虚指数信号

k=0:30;
x=1.5*exp(j*pi/6*k);
subplot(2,2,1);
stem(k,real(x));%实部
grid on;
title('虚指数序列实部');
axis([0,30,-2,2]);%设置坐标范围
xlabel('k');

subplot(2,2,2);
stem(k,imag(x));
grid on;
title('虚指数序列虚部');
axis([0,30,-2,2]);
xlabel('k');

subplot(2,2,3);
stem(k,abs(x));
grid on;
title('虚指数序列的模');
axis([0,30,0,2]);
xlabel('k');

subplot(2,2,4);
stem(k,angle(x));
grid on;
title('虚指数序列的相角');
axis([0,30,-4,4]);
xlabel('k');

其他帮助函数

rectpuls(t,w)：在时间段 t 内产生高度为 1，宽度为 w，关于 t=0 对称的矩形脉冲信号。 
square(w0,duty)：产生基本频率为 w0（以 2pi 为周期），占空比为 duty（百分比）的周期方波信号。
randn(m,n)：产生一个均值为 0，方差为 1 的随机矩阵，其中 m 为行数，n 为列数。