MATLAB ------- 使用 MATLAB 比较占空比不同的周期方波的 DFS 图

ulppknot

6 e9 f. R- s/ Y8 p* A- o周期方波的表达式我就不提炼了，既然你来看文章了，我就不信你周期方波都不知道。

占空比定义为：L/N，N为周期，下面我们对 L = 5，N = 20； L = 5，N = 40；L = 5，N = 60；L = 7，N = 60分别画出它们的DFS图作比较。

' S0 o; d" ^1 l* q我觉得有必要先把对应的周期方波给出来（每幅图都只给出一个周期）：
: m3 \8 J* z6 e; b; r

• clc;clear;close all;
• % Periodic square wave parameters L = 5,N = 20
• L = 5;
• N = 20;
• n = 0:N-1;
• xn = [ones(1,L),zeros(1,N-L)];
• subplot(2,2,1)
• stem(n,xn);
• title('L = 5,N = 20');
• xlabel('n');ylabel('xtilde(n)');
• % Periodic square wave parameters L = 5,N = 40
• L = 5;
• N = 40;
• n = 0:N-1;
• xn = [ones(1,L),zeros(1,N-L)];
• subplot(2,2,2)
• stem(n,xn);
• title('L = 5,N = 20');
• xlabel('n');ylabel('xtilde(n)');
• % Periodic square wave parameters L = 5,N = 60
• L = 5;
• N = 60;
• n = 0:N-1;
• xn = [ones(1,L),zeros(1,N-L)];
• subplot(2,2,3)
• stem(n,xn);
• title('L = 5,N = 20');
• xlabel('n');ylabel('xtilde(n)');
• % Periodic square wave parameters L = 7,N = 60
• L = 7;
• N = 60;
• n = 0:N-1;
• xn = [ones(1,L),zeros(1,N-L)];
• subplot(2,2,4)
• stem(n,xn);
• title('L = 5,N = 20');
• xlabel('n');ylabel('xtilde(n)');
• 
  # M' N1 z" ]& |

            
/ }5 G( }! s0 u
  i2 q# E# c' R
% ]5 Z. T1 {% h+ G% Z上面程序中的DFS函数定义在另一篇博文中：【 MATLAB 】离散傅里叶级数（DFS）及 IDFS 的 MATLAB 实现
3 C( z5 i9 b  _, r) W

1 Z. f- v& I" {! x7 w7 L& Q3 u9 R
3 z! U* t# G+ ]% K, R8 N
直接给出MATLAB脚本：

" f- ~. B' a; H0 I

• clc;clear;close all;
• % Periodic square wave parameters L = 5,N = 20
• L = 5;
• N = 20;
• k = [-N/2:N/2];
• xn = [ones(1,L),zeros(1,N-L)];
• Xk = dfs(xn,N);
• magXk = abs( [ Xk(N/2+1:N),Xk(1:N/2+1) ] );
• subplot(2,2,1)
• stem(k,magXk);
• title('DFS of Square: L = 5,N = 20');
• xlabel('k');ylabel('Xtilde(k)');
• % Periodic square wave parameters L = 5,N = 40
• L = 5;
• N = 40;
• k = [-N/2:N/2];
• xn = [ones(1,L),zeros(1,N-L)];
• Xk = dfs(xn,N);
• magXk = abs( [ Xk(N/2+1:N),Xk(1:N/2+1) ] );
• subplot(2,2,2)
• stem(k,magXk);
• title('DFS of Square: L = 5,N = 40');
• xlabel('k');ylabel('Xtilde(k)');
• % Periodic square wave parameters L = 5,N = 60
• L = 5;
• N = 60;
• k = [-N/2:N/2];
• xn = [ones(1,L),zeros(1,N-L)];
• Xk = dfs(xn,N);
• magXk = abs( [ Xk(N/2+1:N),Xk(1:N/2+1) ] );
• subplot(2,2,3)
• stem(k,magXk);
• title('DFS of Square: L = 5,N = 60');
• xlabel('k');ylabel('Xtilde(k)');
• % Periodic square wave parameters L = 7,N = 60
• L = 7;
• N = 60;
• k = [-N/2:N/2];
• xn = [ones(1,L),zeros(1,N-L)];
• Xk = dfs(xn,N);
• magXk = abs( [ Xk(N/2+1:N),Xk(1:N/2+1) ] );
• subplot(2,2,4)
• stem(k,magXk);
• title('DFS of Square: L = 7,N = 60');
• xlabel('k');ylabel('Xtilde(k)');
• 
  6 t% K3 X5 _: i6 K% z4 `; T; \

              




. |  U% j+ t# @. z7 F
9 t% w2 k0 Q  N5 E7 e

yxlk

谢谢分享

chencol

谢谢分享