数字信号处理应用的几个基本序列:
1 单位样本序列
function mainImseq()
clc
clear
disp('生成抽样序列');
y=imseq(,,); %调用样本函数,此时序列下标以1开头(~5之间5个数,下标为1的抽样值为1) %子函数imseq:抽样函数
function [x,n]=imseq(n0,n1,n2)
n=[n1:n2];
x=[(n-n0) == ]
2 单位阶越序列
产生u(n)
function mainImseq()
clc
clear
disp('生成单位阶跃序列');
y=stepseq(,,); %调用样本函数,此时序列下标以1开头(~10之间10个数,下标为4之后的抽样值全为1) %子函数stepseq:阶跃函数
function [x,n]=stepseq(n0,n1,n2)
n=[n1:n2];
x=[(n-n0) >= ]
3 实值指数序列
n=[0:10]; x=(0.9).^n; plot(x,'*')
4复值指数序列
n=[0:10];
x=exp((2+3j)*n);
plot(x,'*')
5正弦序列
n=[0:10];
x=3*cos(0.1*pi*n+pi/3)+2*sin(0.5*pi*n);
plot(x,'*')
6随机序列
两种伪随机序列
rand(1,N)产生一个长度为N,其值在[0,1]之间均匀分布的随机序列;
rand(m,n)生成m行n列的均匀分布的伪随机数 eg:rand(1,10);
rand(m,n,'double')生成指定精度的均匀分布的伪随机数,参数还可以是'single'
rand(RandStream,m,n)利用指定的RandStream(我理解为随机种子)生成伪随机数
randn(1,N)产生一个长度为N,均值为0,方差为1的高斯型随机序列
主要语法:和上面一样
7周期序列
把一个序列x(n)=[1 2 3 4]产生为3个周期:
x=[1 2 3 4];
xtilde=[x,x,x] %因变量符号可以为任意,如y等
方法二:
>> x=[1 2 3 4];
>> y=x'*ones(1,5);
>> y=y(:);
>> y=y'