汇编语言实验程序.doc

syms t w
Gt=sym('Heaviside(t+1)-Heaviside(t-1)');
Fw=fourier(Gt,t,w);
FFP=abs(Fw);
ezplot(FFP,[-10*pi 10*pi]);grid;
axis([-10*pi 10*pi 0 2.2])

例一图
syms t w
Fw=sym('1/(1+w^2)');
ft=ifourier(Fw,w,t);
FFP=abs(ft);
ezplot(FFP,[-10*pi 10*pi]);grid;
axis([-10*pi 10*pi 0 0.6])

例二图

R=0.02;
t=-2:R:2;
ft=[zeros(1,50),ones(1,101),zeros(1,50)];
W1=10*pi;
M=500;k=0:M;w=k*W1/M;
Fw=ft*exp(-j*t'*w)*R;
FRw=abs(Fw);
W=[-fliplr(w),w(2:501)];
FW=[fliplr(FRw),FRw(2:501)];
Subplot(2,1,1);plot(t,ft);grid;
xlabel('t');ylabel('f(t)');
title('f(t)=u(t+1)-u(t-1)');
subplot(2,1,2);plot(W,FW);grid on;
xlabel('W');ylabel('F(W)');
title('f(t)的振幅频谱图');

例三图

1、（1）：
R=0.02;
t=-2:R:2;
ft=[zeros(1,75),ones(1,51),zeros(1,75)];
W1=10*pi;
M=500;k=0:M;w=k*W1/M;
Fw=ft*exp(-j*t'*w)*R;
FRw=abs(Fw);
W=[-fliplr(w),w(2:501)];
FW=[fliplr(FRw),FRw(2:501)];
Subplot(2,1,1);plot(t,ft);grid;
xlabel('t');ylabel('f(t)');
title('f(t)=u(2t+1)-u(2t-1)');
subplot(2,1,2);plot(W,FW);grid on;
xlabel('W');ylabel('F(W)');
title('f(t)的振幅频谱图');

（2）：
syms t w
Gt1=sym('Heaviside(t)-Heaviside(t-1)');
Gt2=sym('Heaviside(t+1)-Heaviside(t)');
ft=Gt1.*(1-t)+Gt2.*(1+t);
Fw=fourier(ft,t,w);
%FFw=maple('convert',Fw,'piecewise');
FFP=abs(Fw);
ezplot(FFP,[-20 20]);grid;

（3）：
syms t w
gt1=exp(-1*t);
gt2=sym('heaviside(t)');
gt=gt1*gt2;
fw=fourier(gt,t,w);
ffw=maple('convert',fw,'piecewise');
ffp=abs(ffw);
ezplot(ffp,[-10 10]);grid;

2、（1）：
syms t w
Fw=sym('-2*j*w/(16+w^2)');
ft=ifourier(Fw,w,t);
ezplot(ft);grid;