实验三.docx
- 文档编号:24066308
- 上传时间:2023-05-23
- 格式:DOCX
- 页数:23
- 大小:252.93KB
实验三.docx
《实验三.docx》由会员分享,可在线阅读,更多相关《实验三.docx(23页珍藏版)》请在冰豆网上搜索。
实验三
实验报告
实验日期:
学号:
姓名:
实验名称:
连续时间傅里叶级数
总分:
§3.1连续时间傅立叶级数的性质
目的
本练习要检验连续时间傅立叶级数(CTFS)的性质。
相关知识
考虑信号
,式中
。
为了用MATLAB对该信号求值,利用时间向量
>>t=linspace(-1,1,1000);
它创建了在
范围内1000个时间样本的向量。
中等题
1.满足
的最小周期T是多少?
利用这个T值,用定义求
的CTFS系数。
方法一:
clc;
x=sym('cos(2*pi*t)+sin(4*pi*t)');
symst
n=0:
8;
a0=2*int(x,t,-1/2,1/2)
an=2*int(x*cos(2*pi*n*t),t,-1/2,1/2)
bn=2*int(x*sin(2*pi*n*t),t,-1/2,1/2)
cn=double(sqrt(an.^2+bn.^2))
subplot(121)
ezplot(x,[-3,3])%T=1
subplot(122)
stem(n,cn)
方法二:
clc;
x=sym('cos(2*pi*t)+sin(4*pi*t)');
symst
n=-8:
8;
y=x*exp(-i*2*pi*n*t)
z=double(int(y,t,-1/2,1/2))
stem(n,z)
2.考虑信号
,利用CTFS的时间倒置和共轭性质求
的CTFS系数。
方法一:
clc;
x1=sym('cos(2*pi*t)+sin(4*pi*t)');
x2=sym('cos(-2*pi*t)+sin(-4*pi*t)');
x=x1+x2
symst
n=0:
8
a0=2*int(x,t,-1/2,1/2)
an=2*int(x*cos(2*pi*n*t),t,-1/2,1/2)
bn=2*int(x*sin(2*pi*n*t),t,-1/2,1/2)
cn=double(sqrt(an.^2+bn.^2))
subplot(121)
ezplot(x,[-3,3])%T=1
subplot(122)
stem(n,cn)
方法二:
clc;
x1=sym('cos(2*pi*t)+sin(4*pi*t)');
x2=sym('cos(-2*pi*t)+sin(-4*pi*t)');
x=x1+x2
symst
n=-8:
8;
y=x*exp(-i*2*pi*n*t)
z=double(int(y,t,-1/2,1/2))
stem(n,z)
3.在
上画出信号
。
能预计出什么样的对称性?
能够利用CTFS的对称性说明它吗?
clc;
x1=sym('cos(2*pi*t)+sin(4*pi*t)');
x2=sym('cos(-2*pi*t)+sin(-4*pi*t)');
x=x1+x2;
grid;
ezplot(x,[-1,1]);
4.考虑信号
。
利用CTFS的时间倒置和共轭性质求
的CTFS系数。
方法一:
clc;
x1=sym('cos(2*pi*t)+sin(4*pi*t)');
x2=conj(sym('cos(2*pi*t)+sin(4*pi*t)'));
x=x1-x2
symst
n=0:
8
a0=2*int(x,t,-1/2,1/2)
an=2*int(x*cos(2*pi*n*t),t,-1/2,1/2)
bn=2*int(x*sin(2*pi*n*t),t,-1/2,1/2)
cn=double(sqrt(an.^2+bn.^2))
stem(n,cn)
方法二:
clc;
x1=sym('cos(2*pi*t)+sin(4*pi*t)');
x2=conj(sym('cos(2*pi*t)+sin(4*pi*t)'));
x=x1-x2
symst
n=-8:
8;
y=x*exp(-i*2*pi*n*t)
z=double(int(y,t,-1/2,1/2))
stem(n,z)
§3.2连续时间傅立叶级数中的能量关系
目的
分别在时域和频域求信号能量,验证帕斯瓦尔定理。
相关知识
一个硬限幅器是一种器件,其输出是即时输入信号符号的函数,具体说就是当输入信号
是正时,输出信号
等于1;而当
是负时,输出信号
等于-1。
调频(FM)的某些部分或雷达系统中都常用硬限幅器处理某一即时输入信号的相位,而不管任何可能的幅度失真。
在本练习中要考虑将信号
通过硬限幅器的输出
的问题。
中等题
1.求信号
的CTFS表示。
提示:
利用CTFS性质,并根据周期为T的对称方波
具有CTFS系数
为
的知识。
clc;
symst;
k=[-20:
20];
x=(sin((k+eps)*pi/2)./((k+eps)*pi))
y=x.*(1-exp(-i.*(k+eps)*pi))
subplot(2,1,1)
stem(k,y);
f0=y.*exp(i.*(k+eps)*pi*t);
f=simple(sum(f0))
subplot(2,1,2)
ezplot(t,f)
2.一个周期信号的基波分量的能量可以定义为
,其中
是该信号的CTFS。
试计算输出
和输入
中的基波分量的能量,能量有增益或损失吗?
能说明能量变化的原因吗?
x=sym('cos(pi*t)^2');
ex=int(x,0,2)
y=((sin(pi/2)./(pi)).*(1-exp(-i*pi)))^2+((sin(pi/2)./(pi)).*(1-exp(i*pi)))^2
运行结果:
ex=1;ey=0.8106
原因:
信号传递过程中有能量损失。
3.利用帕斯瓦尔定理求该信号一个周期内的总能量,利用前100个频率,即
近似这个和式,这个和式收敛到何值?
k=[-100:
100];
x=(sin(pi.*(k+eps)/2)./((k+eps)*pi)).*(1-exp(-i.*(k+eps)*pi))
y=abs(x);
z=sum(y.^2)
e=2*z
运行结果:
p=0.9959;e=1.9919
即此信号周期为2,一个周期内的功率为0.9959,这个和式的能量近似收敛到2。
4.为了观察该能量估计值收敛得有多快,试画出该信号能量估计值作为在和式中所用项数个数的函数图。
k=[-20:
20];
x=(sin(pi.*(k+eps)/2)./((k+eps)*pi)).*(1-exp(-i.*(k+eps)*pi))
y=abs(x)
z=cumsum(y.^2)
stem(k,2*z)
深入题(可选)
5.利用该信号能量的解析式求下面和式的闭式表达式:
提示:
利用帕斯瓦尔定理将总能量的时域和频域表达式联系起来。
clc;
symst;
ak=1./(2.*k+1).^2;
ft=ak.*exp(i.*k*pi*t);
f=sum(ft)
ezplot(t,f)
§3.3用傅立叶级数综合连续时间信号
目的
学习CTFS系数的连续时间信号的傅立叶分析与综合。
相关知识
一大类连续时间周期信号可以表示成如下和式
式中
是连续时间傅立叶级数(CTFS)。
本练习要综合具有较少非零系数个数的信号,将考虑具有无限个非零CTFS系数的连续时间信号的傅立叶分析与综合。
基本题
对于这些习题要用少数几个非零的傅立叶级数系数构造周期信号的符号表达式。
3个信号的基波周期和非零的CTFS系数给出如下:
对每一信号创建连续时间信号的符号表达式,并用ezplot画出信号的两个周期。
若已知
的图,如何能由两个信号的傅立叶级数系数预计
的图?
如何本来就能根据傅立叶系数预计出3个信号中每一个都应该是实信号?
结果:
x1=4*cos(6*pi*t)+10*cos(2*pi*t)
x2=1/4*sin(4*pi*t)-1/2*sin(3*pi*t)+sin(2*pi*t)-2*sin(pi*t)
x3=-1/2*sin(3/2*pi*t)-2*sin(1/2*pi*t)-1/4*sin(2*pi*t)-sin(pi*t)
clc;
symst;
x1=sym('5*(exp(i*2*pi*t)+exp(-i*2*pi*t))')+sym('2*(exp(i*6*pi*t)+exp(-i*6*pi*t))')
x2=sym('i*(exp(i*pi*t)-exp(-i*pi*t))')-sym('1/2*i*(exp(i*2*pi*t)-exp(-i*2*pi*t))')+sym('1/4*i*(exp(i*3*pi*t)-exp(-i*3*pi*t))')-sym('1/8*i*(exp(i*4*pi*t)-exp(-i*4*pi*t))')
x3=sym('i*(exp(i*1/2*pi*t)-exp(-i*1/2*pi*t))')+sym('1/2*i*(exp(i*pi*t)-exp(-i*pi*t))')+sym('1/4*i*(exp(i*3/2*pi*t)-exp(-i*3/2*pi*t))')+sym('1/8*i*(exp(i*2*pi*t)-exp(-i*2*pi*t))')
subplot(311)
ezplot(x1,[0,2])
subplot(312)
ezplot(x2,[0,4])
subplot(313)
ezplot(x3,[0,8])
中等题
定义
并考虑具有下面基波周期T和傅立叶级数系数
的信号:
4.
5.
对每一信号K=1,3和9。
对每个K值,创建对
的符号表达式,并用ezplot画出各信号的两个周期。
如果信号为复数,要单独分开画出它们的实部和虚部。
%K=1时:
clc;
symst
k=[-1:
1];
y1=1./(k.^2+1)
f1=y1.*(exp(i*2*pi/5.*k*t))
subplot(211)
ezplot(sum(f1),[-5,5])
y2=sign(k)./(-2).^sqrt(k.*k)
f2=y2.*(exp(i*pi/10.*k*t))
subplot(212)
ezplot(t,sum(f2),[-20,20])
%K=3时:
clc;
symst
k=[-3:
3];
y1=1./(k.^2+1)
f1=y1.*(exp(i*2*pi/5.*k*t))
subplot(211)
ezplot(sum(f1),[-5,5])
y2=sign(k)./(-2).^sqrt(k.*k)
f2=y2.*(exp(i*pi/10.*k*t))
subplot(212)
ezplot(t,sum(f2),[-20,20])
%K=9时:
clc;
symst
k=[-9:
9];
y1=1./(k.^2+1)
f1=y1.*(exp(i*2*pi/5.*k*t))
subplot(211)
ezplot(sum(f1),[-5,5])
y2=sign(k)./(-2).^sqrt(k.*k)
f2=y2.*(exp(i*pi/10.*k*t))
subplot(212)
ezplot(t,sum(f2),[-20,20])
6.怎样本就能由傅立叶级数系数预计到哪个信号应该是实信号?
答:
若傅立叶级数系数本身是一个实数,在自身傅里叶级数虚部不抵消的情况下,一般来讲,则它的信号应该是一个虚数
§3.4方波和三角波的傅立叶表示
目的
这个练习要用傅立叶级数分析周期方波,对信号将研究截断的傅立叶级数重构公式,特别要研究随N的增大,
是如何收敛的?
相关知识
一般来说,傅立叶级数系数可有无限个非零值。
譬如,任何具有间断点的信号都一定有一个无限个非零系数的傅立叶级数表示,而对数值计算来说,这是无法实现的。
有限项和对某个相对小的整数N往往是一个很好的近似。
方波
考虑一个基波周期T=2的周期方波,在
区间内该方波由下式表示,
如下图示:
这个练习将分析方波的傅立叶级数表示,且主要集中在方波的不连续点。
1.利用int创建一个符号表达式a,它包含了该方波每个k值的傅立叶级数系数。
这个符号表达式是k的函数。
利用stem画出
内的傅立叶级数的系数。
clc;
symskt
k=-10:
10
ak=exp(-i*pi*k*t)
a=(1/2)*int(ak,t,-1/2,1/2)
stem(k,a)
2.对N=1,3,5和9,对
创建符号表达式。
利用ezplot画出区间
内的
,用hold将4张图画在同一幅图上。
clc;
symskt
%N=1
k=-1:
1
y=exp(-i*pi*k*t)
a=(1/2)*int(y,t,-1/2,1/2)
x1=sum(a.*exp(i*pi*k*t))
subplot(221)
ezplot(t,x1,[-1,1])
%N=3
k=-3:
3
y=exp(-i*pi*k*t)
a=(1/2)*int(y,t,-1/2,1/2)
x2=sum(a.*exp(i*pi*k*t))
subplot(222)
ezplot(t,x2,[-1,1])
%N=5
k=-5:
5
y=exp(-i*pi*k*t)
a=(1/2)*int(y,t,-1/2,1/2)
x3=sum(a.*exp(i*pi*k*t))
subplot(223)
ezplot(t,x3,[-1,1])
%N=9
k=-9:
9
y=exp(-i*pi*k*t)
a=(1/2)*int(y,t,-1/2,1/2)
x4=sum(a.*exp(i*pi*k*t))
subplot(224)
ezplot(t,x4,[-1,1])
3.在
时,
的值是多少?
这个值随N增加而变化吗?
答:
的值为1/2,这个值不会随N增加而变化
4.不用明确地求出
,对每个N值估计一下超量误差值。
这个超量误差随N增加而减小吗?
随
,如何预期这个值的变化?
答:
这个超量误差会随着N的增加而减少,随
时,这个值趋近于0
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 实验