dm10201.m

来自「各种周期信号的傅里叶变换」· M 代码 · 共 20 行

% dm10201
% 门信号的傅里叶变换
syms t w ;                                      %定义符号变量
ut = sym('Heaviside(t+0.5)-Heaviside(t-0.5)');  %借助阶跃函数定义门函数
subplot(211);  
ezplot(ut);                                     %绘制门信号
hold on;
axis([-1 1 0 1.1]);                             %给定坐标系的范围
plot([-0.5 -0.5], [0 1]);                       %绘制左跳变沿
plot([0.5 0.5], [0 1]);                         %绘制右跳变沿
title('时域波形')
Fw = fourier(ut,t,w);                           %用fourier函数对门信号求傅里叶变换
FFw = maple('convert',Fw,'piecewise');          %将傅里叶变换的符号表达式转换为样条曲线
FFP = abs(FFw);                                 %求幅度频谱 
subplot(212);
ezplot(FFP,[-10*pi 10*pi]);                     %绘制幅度频谱
axis([-10*pi 10*pi 0 1.1]);                     %给定坐标系的范围
title('幅度频谱')
%End