dm10203.m

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

% dm10203
% 直流信号的傅里叶变换
display('Please input the value of a')    %屏幕提示输入a值
a = input('a = ');                    %键盘输入a值
syms t;                           %定义符号变量
f = exp(-a*abs(t));  
subplot(121);
ezplot(f);                         %绘制双边指数信号的波形
axis([-2*pi 2*pi 0 1]);
ylabel('时域波形');
F = fourier(f);                     %求傅里叶变换
subplot(122);
ezplot(abs(F));                        %绘制幅度频谱
axis([-3 3 0 2/a]);
ylabel('幅度频谱');
%End