例2-24 将单位阶跃序列u(n)分解为偶分量和奇分量

%2017/8/22日  星期二
%工程北629
%功能：数字信号的偶奇分解
%思想：使用evenodd函数对阶跃信号进行 偶奇分解

clear all; clc

%输入信号
n0 = 0;
n1 = -10;
n2 = +10;
n = [n1:n2];    %序号序列

%信号处理
x = sin(n);    %使用表达式生成阶跃信号
[xeven,xodd,m] = evenodd(x,n);  %进行偶奇分解
grid on

%输出信号
subplot(2,2,1),stem(n,x); xlabel('n'),ylabel('x(n)'),title('阶跃序列u(n)');  grid on
subplot(2,2,2),stem(m,xeven); xlabel('m'),ylabel('xeven(n)'),title('偶函数部分xeven');grid on   %绘制偶序列
subplot(2,2,3),stem(m,xodd);  xlabel('m'),ylabel('xodd(n)'),title('奇函数xodd');grid on

%总结：可以将信号分解：原信号 == 偶信号 + 奇信号

%{
    %奇偶函数
    %将该函数单独编写一个.m文件，命名为evenodd.m  然后该程序单独编写一个.m文件，放在同一个文件夹下，即可
    function [xe, xo, m] = evenodd(x,n)
    if any(imag(x) ~= 0)
        error('x is not a real sequence')
    end
    m = -fliplr(n);     %将序号反转
    m1 = min([m,n]); m2 = max([m,n]); m = m1:m2;    %创建新序号序列
    nm = n(1) - m(1); n1 = 1:length(n);
    x1 = zeros(1,length(m));    %创建空序列
    x1(n1 + nm) = x; x = x1;    
    xe = 0.5*(x + fliplr(x));   %偶序列
    xo = 0.5*(x - fliplr(x));   %奇序列
%}