湖南大学实验三连续时间傅里叶级数.docVIP

实验三 连续时间傅里叶级数 §3.1 连续时间傅立叶级数的性质 目的 本练习要检验连续时间傅立叶级数（CTFS）的性质。 实验内容： 1．满足的最小周期T是多少？利用这个Ｔ值，用解析法求的CTFS系数。 2．考虑信号，利用CTFS的时间倒置和共轭性质求的CTFS系数。 3．在上画出信号。能预计出什么样的对称性？能够利用CTFS的对称性说明它吗？ 4．考虑信号。利用CTFS的时间倒置和共轭性质求的CTFS系数。 代码： 1. x=sym(cos(2*pi*t)+sin(4*pi*t)); t=linspace(-1,1,1000); subplot(311); ezplot(x,[-1,1]); grid on; k=[-5:5]; syms t; f=x*exp(-i*2*pi*t*k); Fn=int(f,t,0,1); F=abs(Fn); subplot(312); subplot(313); stem(F); 2. x1=sym(cos(2*pi*t)+sin(4*pi*t)); x2=subs(x1,-t,t); y=x1+x2; t=linspace(-1,1,1000); k=[-5:5]; syms t; f=y*exp(-i*2*pi*t*k); Fn=int(f,t,0,1); F=abs(Fn)； subplot(211); stem(Fn); subplot(212); stem(F); 3. x1=sym(cos(2*pi*t)+sin(4*pi*t)) x2=subs(x1,(-1*t),t) y=x1+x2; ezplot(y,[-1,1]); 4. x1=sym(cos(2*pi*t)+sin(4*pi*t)) x2=conj(x1); z=x1-x2; ezplot(z,[-1,1]); §3.2 连续时间傅立叶级数中的能量关系 目的 分别在时域和频域求信号能量，验证帕斯瓦尔定理。 相关知识 一个硬限幅器是一种器件，其输出是即时输入信号符号的函数，具体说就是当输入信号是正时，输出信号等于1；而当是负时，输出信号等于－1。调频（FM）的某些部分或雷达系统中都常用硬限幅器处理某一即时输入信号的相位，而不管任何可能的幅度失真。在本练习中要考虑将信号通过硬限幅器的输出的问题。 中等题 1．求信号的CTFS表示。提示：利用CTFS性质，并根据周期为T的对称方波 具有CTFS系数为 的知识。 2．一个周期信号的基波分量的能量可以定义为，其中是该信号的CTFS。试计算输出和输入中的基波分量的能量，能量有增益或损失吗？能说明能量变化的原因吗？ 3．利用帕斯瓦尔定理求该信号一个周期内的总能量，利用前100个频率，即近似这个和式，这个和式收敛到何值？ 4．为了观察该能量估计值收敛得有多快，试画出该信号能量估计值作为在和式中所用项数个数的函数图。 深入题（可选） 5．利用该信号能量的解析式求下面和式的闭式表达式： 提示：利用帕斯瓦尔定理将总能量的时域和频域表达式联系起来。 代码： 1. syms t; k=[-20:20]; ak=(sin((k+eps)*pi/2)./((k+eps)*pi)) fk=ak.*(1-exp(-i.*(k+eps)*pi)); stem(k,fk); 2. x=sym(cos(pi*t)^2); ex=int(x,0,2) ey=((sin(pi/2)./(pi)).*(1-exp(-i*pi)))^2+((sin(pi/2)./(pi)).*(1-exp(i*pi)))^2 运行结果：ex =1 ey = 0.8106 系统传递的过程中有能量的损失。 3. syms t; k=[-100:100]; ak=(sin((k+eps)*pi/2)./((k+eps)*pi)); fk=ak.*(1-exp(-i.*(k+eps)*pi)); figure(1) stem(k,fk); Ey=sum(abs(fk).^2); x=sym(cos(pi*t)); f=x*exp(-i*2*pi*t*k); fx=int(x.^2,t,0,2); e=2*Ey; Ex=abs(sum(x.^2)); 4. syms t; k=[-30:30]; ak=(sin((k+eps)*pi/2)./((k+eps)*pi)); fk=ak.*(1-exp(-i.*(k+eps)*pi)); y=abs(fk); Py=(y.^2); cumsum(Py,2) stem(k,2*Py) 5. syms t