离散时间系统及离散卷积

1. DFT与IDFT

function [XK] = dft(xn,N)
%UNTITLED 此处显示有关此函数的摘要
%   此处显示详细说明
n=0:1:N-1;
k=n;
WN=exp(-1i*2*pi/N);
nk=n\'*k;
WNnk=WN.^nk;
XK=xn*WNnk
end

function [XK] = idft(Xk,N)
%UNTITLED2 此处显示有关此函数的摘要
%   此处显示详细说明
n=0:1:N-1;
k=n;
WN=exp(-1i*2*pi/N);
nk=n\'*k;
WNnk=WN.^(-nk);
XK=(xn*WNnk)/N
end

2 差分方程求序列及DFT

clc;
clear;
b=[1,0];
a=[1,-1];
%y(n)-y(n-1)=x(n)
 y_1=1;  x0=0;
 x=1:10;        %激励序列x(n)
y_0=1;
y=[y_0];
for n=2:10
     y=[y,y(end)+x(n)];    %序列递推
end
subplot(1,2,1)
stem(y);
xlabel(\'n\');
ylabel(\'y(n)\');
title(\'递推后的序列y\');
subplot(1,2,2)
XK=dft(y,10);
stem(XK);
xlabel(\'n\');
ylabel(\'X(k)\');
title(\'DFT:y\');

 

 3 幅频特性与相频特性

%实验一（2）  幅频特性与相频特性计算，画出零极点
clc;             %清除命令行
close all;       %清除窗口
figure 
b=[0,1];
a=[1,-2,1];              %选用系统函数为：H（z）=z/(1-z)^2
[r,p,k]=residue(b,a)     %求z反变换
zplane(b,a);title(\'零极点图\');       %零极点图，z=1是二阶极点
figure
n=0:10;
h=impz(b,a,11);
stem(n,h,\'filled\'); title(\'冲激响应h(n)\');    %冲激响应图,h(n)=n*u(n)
figure
w=pi*freqspace(100);
freqz(b,a,w); title(\'幅频特性与相频特性\');     %系统频率响应的幅频特性和相频特性图

 

 

 

 4 FFT 求卷积

clear;
xn=sin(0.4*[1:15])\';hn=0.9.^(1:20)\';
M=length(xn);N=length(hn);
nx=1:M;nh=1:N;
%L>=N+M+1;
L=pow2(nextpow2(M+N-1));
tic
Xk=fft(xn,L);Hk=fft(hn,L);
Yk=Xk.*Hk;
yn=ifft(Yk,L);
toc
tic,
yn1=conv(xn,hn);
toc
subplot(4,1,1),stem(1:15,xn,\'filled\',\'k\');title(\'序列一\');
subplot(4,1,2),stem(1:20,hn,\'filled\',\'k\');title(\'序列二\');
subplot(4,1,3),stem(1:(M+N-1),yn(1:(M+N-1)),\'filled\',\'k\');title(\'fft ifft\');
subplot(4,1,4),stem(1:(M+N-1),yn1,\'filled\',\'k\');title(\'conv\');

 

 

可以看到内置函数的速度要快很多

5 线性卷积

%实验一（3）     序列的线性卷积计算
clc;                      
clear;
x=[3,11,7,0,-1,4,2];       %序列1 
h=[2,3,0,-1,2,1];          %序列2
y=conv(x,h)                %求线性卷积,并输出结果
subplot(4,1,1)   
stem([0:length(x)-1],x,\'filled\');           %序列1图像
xlabel(\'n\'); ylabel(\'幅度\'); title(\'x\');      
subplot(4,1,2)
stem([0:length(h)-1],h,\'filled\');           %序列2图像
xlabel(\'n\'); ylabel(\'幅度\'); title(\'h\')
subplot(4,1,3)
stem([0:length(y)-1],y,\'filled\');            %线性卷积结果图像
xlabel(\'n\'); ylabel(\'幅度\'); title(\'线性卷积\')
% dft方法
M=size(x,2);N=size(h,2);
L=M+N-1;
x1=zeros(1,L);x1(1:M)=x;
h1=zeros(1,L);h1(1:N)=h;
xk=dft(x1,L);hk=dft(h1,L);
wk=xk.*hk  
y2=idft(wk,L)
subplot(4,1,4)
stem([0:L-1],y2,\'filled\');
xlabel(\'n\'); ylabel(\'幅度\'); title(\'DFT 线性卷积\');