信号与系统实验报告连续时间傅立叶级数.docVIP

实验三 连续时间傅立叶级数 3.1连续时间傅立叶级数的性质 中等题 1.满足x1(t)=x1(t+T)的最小周期T是多少？利用这个T值用解析法求x1(t)的CTFS系数。 程序： clc; x=sym(cos(2*pi*t)+sin(4*pi*t)); t=linspace(-1,1,1000); ezplot(x,[-1,1]); grid; k=[-5:5]; syms t; f=x*exp(-i*2*pi*t*k); Fn=int(f,t,0,1) F=abs(Fn) Fn=double(Fn) subplot(1,2,1); stem(k,Fn); F=double(F) subplot(1,2,2); stem(k,F) 运行结果： Fn =[ 0, 0, 0, 1/2*i, 1/2, 0, 1/2, -1/2*i, 0, 0, 0] F =[ 0, 0, 0, 1/2, 1/2, 0, 1/2, 1/2, 0, 0, 0] 图： 可以看出周期T=1. 2.考虑信号y（t）=x1（t）+x1（-t），利用CTFS的时间倒置和共轭性质求y（t）的CTFS系数。 程序： clc; x=sym(cos(2*pi*t)+sin(4*pi*t)); t=linspace(-1,1,1000); k=[-5:5]; syms t; f=x*exp(-i*2*pi*t*k); Fn=int(f,t,0,1);%积分 Fn1=fliplr(Fn);%倒置 Fn2=Fn+Fn1; F=abs(Fn2) F=double(F)%转换成double型 stem(k,F) grid; 运行结果： F =[ 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0] 分析：利用到了转置函数fliplr。两个式子相加之后的结果刚好为2cos(2*pi*t)，因而不含有虚数部分。 3.在-1=t=t上画出信号y（t）。能预计出什么样的对称性？能够利用CTFS的对称性说明它吗？ 程序： clc; syms t; x=sym(cos(2*pi*t)+sin(4*pi*t)); y1=subs(x,-t,t);%用t代替-t y=y1+x; ezplot(y,[-1,1]) grid; 分析：可见关于x=0对称，为偶函数，因为相加之后将正弦部分抵消。 4.考虑信号。利用CTFS的时间倒置和共轭性质求z（t）的CTFS系数。 程序： clc; x=sym(cos(2*pi*t)+sin(4*pi*t)); t=linspace(-1,1,1000); k=[-5:5]; syms t; f=x*exp(-i*2*pi*t*k); Fn=int(f,t,0,1); Fn1=fliplr(Fn);%倒置 Fn2=conj(Fn1);%共轭 Fn3=Fn-Fn2 F=abs(Fn3) F=double(F) stem(k,F) grid; Fn3 =[ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0] F =[ 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0] 分析：利用了时间转置函数fliplr和虚数函数conj。两个式子是实数，求共轭相减之后为0. 3.2连续时间傅立叶级数中的能量关系 中等题 求信号y（t）的CTFS表示。提示：利用CTFS性质，并根据周期为T的对称方波 具有CTFS系数为 的知识。 程序： clc; syms t; k=[-20:20]; ak=(sin((k+eps)*pi/2)./((k+eps)*pi)) fk=ak.*(1-exp(-i.*(k+eps)*pi)) subplot(2,1,1) stem(k,fk); grid; f=fk.*exp(i.*(k+eps)*pi*t); f1=simple(sum(f)) subplot(2,1,2) ezplot(t,f1) grid; 图： 程序： clc; x=sym(cos(pi*t)^2); ex=int(x,0,2) subplot(2,1,1) ezplot(ex) grid ey=((sin(pi/2)./(pi)).*(1-exp(-i*pi)))^2+((sin(pi/2)./(pi)).*(1-exp(i*pi)))^2 subplot(2,1,2) ezplot(ey) grid 结果： ex =1 ey =0.8106 分析：系统传递的过程中有能量的损失。 程序： clc; k=[-100:100]; fk=(sin(pi.*(k+eps)/2