DSP序列基本运算MATLAB代码

一、序列翻转

%用到函数 fliplr() 维度不超过2
clc;clear;
x = 0:10;
y = 0.5.^x;%指数序列
y1 = flip(y);%将序列值翻转
x1 = -fliplr(x);%序列范围翻转
subplot(121),stem(x,y);
subplot(122),stem(x1,y1);

二、奇偶合成

function [xe,xo,m] = evenodd(x,n)
%xe共轭对称部分，xo共轭反对称部分，位置向量 x = xe + xo
%x待处理序列，n位置向量
if any(imag(x) ~= 0)
    error('x不是实序列');
end
m = -fliplr(n);
m1 = min([m,n]);
m2 = max([m,n]);
m = m1:m2;
nm = n(1) - m(1);
n1 = 1:length(x);
x1 = zeros(1,length(m));
x1(n1+nm) = x;
x = x1;
xe = 0.5 * (x + fliplr(x));%这里用的实序列
xo = 0.5 * (x - fliplr(x));

三、序列加减

function [y,n] = seqadd(x1,n1,x2,n2)
%序列之和
n = min(min(n1),min(n2)):max(max(n1),max(n2)); %位置向量
y1 = zeros(1,length(n));
y2 = y1;
y1(find((n>=min(n1)) & (n<= max(n1)) == 1)) = x1;%把x1拓展到n上形成y1
y2(find((n>=min(n2)) & (n<=max(n2)) == 1)) = x2;%把x2拓展到n上形成y2
y = y1 + y2;

四、序列相乘

function [y,n] = seqmult(x1,n1,x2,n2)
%序列之乘积
%实现y(n) = x1(n) * x2(n);
n = min(min(n1),min(n2)):max(max(n1),max(n2)); %位置向量
y1 = zeros(1,length(n));
y2 = y1;
y1(find((n>=min(n1)) & (n<= max(n1)) == 1)) = x1;%把x1拓展到n上形成y1
y2(find((n>=min(n2)) & (n<=max(n2)) == 1)) = x2;%把x2拓展到n上形成y2
y = y1.*y2;

五、求和（累加）

function [y,n] = seqsum(x,n)
%序列累加
n = n;
y = zeros(1,length(x));
for i = 1:length(x)
    for j=1:i
        y(i) = y(i) + x(j);
    end
end

六、序列移位

%只需要改变位置向量即可
n = [1,2,3,4,5];
y = [2,4,6,7,4];
%y[n+m]
m = 2;
n = n+m;

七、卷积

%只需要知道，卷积后的位置向量怎样计算即可，卷积运算可以用matlab自带的conv函数计算即可
%位置向量n1，和位置向量n2，卷积后的位置向量是
n1(1)+n2(1):n1(end)+n2(end)
%向量长度为length(n1) + length(n2) -1