%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)); %奇序列 %}
例2-24 将单位阶跃序列u(n)分解为偶分量和奇分量
最新推荐文章于 2024-07-01 15:59:21 发布