多瑙河
2009-11-02, 17:30
clear all;
close all;
syms t;
f=exp(-2*t)*sym('Heaviside(t)');
F=fourier(f);
subplot(3,1,1);
ezplot(f); %符号函数作二维图像指令
title('输入信号');
xlabel('t');
ylabel('f(t)');
subplot(3,1,2);
ezplot(abs(F)); %求绝对值或模
title('幅频特性');
xlabel('\omega');
ylabel('F(\omega)');
im=imag(F); %求虚部
re=real(F); %求实部
phase=atan(im/re); %求相位
subplot(3,1,3);
ezplot(phase);
title('相频特性');
xlabel('\omega');
ylabel('a(\omega)');
上面这段程序在Matlab 6.5下正常运行;可是在Matlab 2008a下却提示出错
Error in ==> ep35 at 5
f=exp(-2*t)*sym('Heaviside(t)');!
为什么啊?
close all;
syms t;
f=exp(-2*t)*sym('Heaviside(t)');
F=fourier(f);
subplot(3,1,1);
ezplot(f); %符号函数作二维图像指令
title('输入信号');
xlabel('t');
ylabel('f(t)');
subplot(3,1,2);
ezplot(abs(F)); %求绝对值或模
title('幅频特性');
xlabel('\omega');
ylabel('F(\omega)');
im=imag(F); %求虚部
re=real(F); %求实部
phase=atan(im/re); %求相位
subplot(3,1,3);
ezplot(phase);
title('相频特性');
xlabel('\omega');
ylabel('a(\omega)');
上面这段程序在Matlab 6.5下正常运行;可是在Matlab 2008a下却提示出错
Error in ==> ep35 at 5
f=exp(-2*t)*sym('Heaviside(t)');!
为什么啊?