第三章_DS定义2016S.pptVIP

* function [Xk] = dfs(xn,N) % Computes Discrete Fourier Series Coefficients % --------------------------------------------- % [Xk] = dfs(xn,N) % Xk = DFS coeff. array over 0 = k = N-1 % xn = One period of periodic signal over 0 = n = N-1 % N = Fundamental period of xn % n = [0:1:N-1]; % row vector for n k = [0:1:N-1]; % row vecor for k WN = exp(-j*2*pi/N); % Wn factor nk = n*k; % creates a N by N matrix of nk values WNnk = WN .^ nk; % DFS matrix Xk = xn * WNnk; % row vector for DFS coefficients 3.2.1 DFS 的Matlab 实现 * function [Xk] = dfs(xn,N) n = [0:1:N-1]; % row vector for n k = [0:1:N-1]; % row vecor for k WN = exp(-j*2*pi/N); % Wn factor nk = n*k; % creates a N by N matrix of nk values WNnk = WN .^ nk; % DFS matrix Xk = xn * WNnk; % row vector for DFS coefficients ? function [xn] = idfs(Xk,N) n = [0:1:N-1]; % row vector for n k = [0:1:N-1]; % row vecor for k WN = exp(-j*2*pi/N); % Wn factor nk = n*k; % creates a N by N matrix of nk values WNnk = WN .^ (-nk); % IDFS matrix xn = (Xk * WNnk)/N; % row vector for IDFS values 3.2.1 DFS 的Matlab 实现 * 用 Matlab 计算，程序和结果如下： ? xn=[0,1,2,3]; N=4; xk=dfs(xn,N) ? xk = 6.0000 -2.0000 + 2.0000i -2.0000 - 0.0000i -2.0000 - 2.0000i * 例 下面给出一周期“方波”序列： 其中， m=0, ±1, ±2,…，N 是基本周期， L/N 是占空比 a）确定一种用 L 与 N 描述的 的表达式 b）分别画出当 L =5, N = 20；L=5，N=40；L=5，N=60；L=7，N=60 时 表达式 c）对所得结果进行讨论 * 解：a）由 DFS 定义可得 而： 的幅值可表示为： * b. Matlab 程序如下： % Chapter 3: Example 3.03 L = 5; N = 20;(可变参数） x = [ones(1,L), zeros(1,N-L)]; xn = x * ones(1,3); xn = (xn(:)); n = -N:1:2*N-1; subplot(1,1,1);subplot(2,1,2); stem(n,xn); xlabel(n); ylabel(xtilde(n)) title(Three periods of xtilde(n)) axis([-N,2*N-1,-0.5,1.5]) -20 -10 0 10 20 30 -0.5 0 0.5 1 1.5 n xtilde(n) Three periods of xtilde(n) * % Part (b) L = 5; N = 20;(可变参数） xn = [ones(1,L), zeros(