Digital Signal Processing Using Matlab读书笔记 第三章（四）（P51-P55）

例3.10

>> n=-5:10;x=rand(1,length(n))+j*rand(1,length(n));
>> k=-100:100;
>> w=(pi/100)*k;
>> X=x*(exp(-j*pi/100)).^(n'*k);
>> y=conj(x);
>> Y=y*(exp(-j*pi/100)).^(n'*k);
>> Y_check=conj(fliplr(X));
>> error=max(abs(Y-Y_check))

error =

   1.1025e-13

>> 

例11

>> x=rand(1,length(n));
k=-100:100;w=(pi/100)*k;
X=x*(exp(-j*pi/100)).^(n'*k);
y=fliplr(x);
m=-fliplr(n);
Y=y*(exp(-j*pi/100)).^(m'*k);
Y_check=fliplr(X);
error=max(abs(Y-Y_check))

error =

   1.0474e-15

例12

>> n=-5:10;
>> x=sin(pi*n/2);
>> k=-100:100;w=(pi/100)*k;
>> X=x*(exp(-j*pi/100)).^(n'*k);
>> [xe,xo,m]=evenodd(x,n);
>> XE=xe*(exp(-j*pi/100)).^(m'*k);
>> XO=xo*(exp(-j*pi/100)).^(m'*k);
>> XR=real(X);
>> error1=max(abs(XE-XR))

error1 =

   6.4660e-14

>> XI=imag(X);
>> error2=max(abs(XO-j*XI))

error2 =

   6.4711e-14

LTI系统的频域表示

先前我们陈述了LTI系统中傅里叶变换表示是所有信号表示形式中最有用的，这归因于以下结论

复指数（e^jw0n)响应

定义1 频率响应

冲激信号的离散时间傅里叶变换被称为LTI系统的频域响应，被表示为

由式子3.15，我们可以这样表示系统

因此，输出序列是是输入序列在频率w0上系统响应的调制结果，这证明了H是一个频率响应的定义，因为它被乘一个复指数得到了输出y，这个强大的结论可以通过使用LTI系统的线性扩展到复指数的线性结合

总的来说，H的频率响应是一个关于w的复函数，H的长度称为幅频响应，H的角度称为相频响应

正弦序列响应

随机序列响应

y可以通过对Y进行IDFT变换得到，这需要积分操作，在matlab来说并不简单，在第四章中我们将会看到，有另外一种方法z变换和部分分式展开法来计算输入随机序列的输出，在这一章我们专注于稳态响应