信号与系统实验报告
实验四 连续时间傅立叶变换
§4.1连续时间傅立叶变换的数字近似
1．求CTFT的解析表达式。可将看作，。
g=sym('exp(-2*t)*Heaviside(t)');
g2=subs(g,'-t','t');
x=g+g2;
fx=fourier(x);
形容茶的优美句子
2．创建一个向量，它包含了在区间t=[0:tau:T-tau] 上（其中和），信号炖羊腿的做法的样本。
clc;
t=0:0.01:10-0.01;
春节有几天>品牌大全y=exp(-2*(t-5)).*(0.5+0.5*sign(t-5))+exp(2*(t-5)).*(0.5+0.5*sign(-t+5));
plot(t,y);
3．键入y=fftshift(tau*fft(y))计算样本。因为对于基本上为零，就能近似用个样本分析中计算出信号的CTFT。
clc;
t=0:0.01:10-0.01;
y=exp(-2*(t-5)).*(0.5+0.5*sign(t-5))+exp(2*(t-5)).*(0.5+0.5*sign(-t+5));
y=fftshift(0.01.*fft(y));我和我的祖国串词
y=abs(y);
plot(t,y);
axis([4,6,-0.1,1.2]);
4.构造一个频率样本向量w，它按照
>> w=-(pi/tau)+(0:N-1)*(2*pi/(N*tau));
与存在向量Y中的值相对应。
5.因为是通过时移与相联系的，所以CTFT就以线性相移项与相联系。利用频率向量w直接由Y计算的样本，并将结果存入x中。去粉刺的方法
clc;
靡靡之音的意思t=0:0.01:10-0.01;
tau=0.01;
N=10/0.01;
y=exp(-2*(t-5)).*(0.5+0.5*sign(t-5))+exp(2*(t-5)).*(0.5+0.5*sign(-t+5));殚精竭虑造句
y=fftshift(0.01.*fft(y));
w=-(pi/tau)+(0:N-1)*(2*pi/(N*tau));
x=exp(j*5*w).*y;
6.利用abs和angle画出在w标定的频率范围内X的幅值和相位。对于相同的值，也画出在1中所导出的解析式表达式的幅值和相位。CTFT的近似值与解析导得的相符吗?若想在一张对数坐标上画出幅值，可以用milogy，这是会注意到，在较高的频率上近似不如在较低的频率上好。因为用了样本近似，所以在时间段长度内，信号变化不大的那些信号的频率分量近似程度会更好一些。
clc;
tau=0.01;
T=10;
t=0:tau:T-tau;
N=T/tau;
y=exp(-2*(t-5)).*(0.5+0.5*sign(t-5))+exp(2*(t-5)).*(0.5+0.5*sign(-t+5));
y=fftshift(0.01.*fft(y));
w=-(pi/tau)+(0:N-1)*(2*pi/(N*tau));
x=exp(i*5*w).*y;
xp=abs(x);
xf=angle(x);
subplot(211)
plot(t,xp);
subplot(212)
plot(t,xf);
7.利用abs和angle画出Y的幅值和相位，它们与X的图比较后怎样？能估计到这一结果吗？
clc
syms s w;
g=sym('exp(-2*s)*Heaviside(s)');
g2=subs(g,'-s','s');
y=g+g2;
fw=fourier(y,s,w);
ff=atan(imag(fw)/real(fw));
fp=abs(fw);
tau=0.01;
T=10;
t=0:tau:T-tau;
N=T/tau;
y=exp(-2*(t-5)).*(0.5+0.5*sign(t-5))+exp(2*(t-5)).*(0.5+0.5*sign(-t+5));
y=fftshift(0.01.*fft(y));
w=-(pi/tau)+(0:N-1)*(2*pi/(N*tau));
x=exp(i*5*w).*y;
