- 1、本文档共4页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
实验三连续信号频谱分析
实验三 连续信号的分析
实验目的
掌握周期信号的频谱—— Fourier 级数的分析方法及其物理意义。
深入理解信号频谱的概念,掌握典型信号的频谱以及 Fourier 变换的主要性质。
实验原理及方法
(3.1)
(3.2)
计算机不能计算无穷多个系数,假设需要计算的谐波次数为N,则总的系数个数为2N+1T和dt3.7)可以近似为:
(3.8)
对于全部的2N+1个系数,上面的计算可以按照矩阵运算实现。
例题:画出占空比矩形脉冲clc;clear;
dt=0.01;
t=-10:dt:10;
T=4;
d1=-10:T:10;
f1=pulstran(t,d1,rectpuls); %产生脉宽为1,周期为4的周期矩形脉冲
N=100;
omega0=2*pi/T;
for n=-N:N
F (100+1+n)=f1*exp(-j*n*omiga0*t)*dt/T; %计算全部2N+1个系数
end
k=-N:N;
w=k*omega0; %o横坐标对应的频率
magF=abs(F);
phaF=angle(F);
subplot(311);
plot(t,f1);
title(周期矩形脉冲);
subplot(312);
stem(k,magF);
title(幅度谱);
subplot(313);
plot(k,phaF);
stem(相位谱);
其中,时间变量的变化步长dt的大小对傅里叶级数系数的计算精度影响非常大,dt越小,精度越高,但计算机计算所花的时间越长。
原信号可以用有限项谐波成分来近似合成,即
(3.3)
例题:
Matlab实现信号合成的程序如下:
f=0;
L=2*N+1;
for n=1:L-1;
f=f+F(n)*exp(j*(n-1-N)*omega0*t); %信号合成
end
周期信号频谱求解另一种方法:先由定义求出信号频谱的数学表达式,再用matlab做出频谱图。函数abs()求幅度,angle()求相位。
3.非周期信号的傅里叶变换
(3.)
Matlab实现傅里叶变换有两种方法,一种利用符号运算的方法计算,另一种是数值计算。
(1)利用符号运算的方法实现。MATLAB的Symbolic Math Toolbox 提供了能直接求解傅里叶变换与反变换的函数fourier( ) 及 ifourier( )。调用格式如下:
( F=fourier(f):它是符号函数f 的傅里叶变换,默认返回函数F是关于的函数;
( F=fourier(f,v):它的返回函数F是关于符号对象v的函数,即;
( F=fourier(f,u,v): 它是对关于u的函数f进行变换,而返回函数F是v的函数,即
( f=ifourier(F): 它是函数F 的傅里叶反变换,默认的独立变量为,默认返回是关于x的函数。如果F=F(x),则ifourier(F)返回关于t的函数
( f=ifourier(F,u): 它的返回函数f是u的函数,而不是默认的x的函数。
( f=ifourier(f,v,u): 它是对关于v的函数F进行变换,而返回关于u的函数f;
这里要注意的是,在调用上述两个函数之前,先要用syms命令对所用到的变量(如t、u、v、)等进行定义,将这些变量定义为符号变量。对于fourier( )中的函数f或ifourier( )中的F,也要用syms将f或F定义为符号表达式。另外,在采用fourier( )及ifourier( )得到的返回函数,仍然是符号表达式。若需要对返回函数作图时,只能用ezplot( )绘图命令,而不能用plot( )命令。如果返回函数中含有等项,用ezplot()也无法作图。
fourier()函数的局限性:用fourier()对某些信号求反变换时,其返回函数可能会包含一些不能直接表达的式子,甚至可能会出现一些屏幕提示为“未被定义的函数或变量”的项;另外,在许多情况下,信号尽管是连续的,但却不可能表示成符号表达式;函数fourier()也不可能对离散信号进行处理。
(2)用数字计算的方法实现。
用数值计算的方法计算连续时间信号的傅里叶变换需要信号是时限信号,也就是当时间|t|大于某个给定时间时其值衰减为零或接近于零,计算机只能处理有限大小和有限数量的数。采用数值计算算法的理论依据是:
(3.11)
若信号为时限信号,当时间间隔T取得足够小时,上式可演变为:
文档评论(0)