信号与系统实验室 V 信号的傅立叶分析和频谱
目录
【实验目的】
【实验设备】
【实验内容】
1.参考例5-1,实现教材p125,例3-4中傅里叶级数表达式(p126第二行)。分别采用前4、40、400项,画出周期矩形脉冲信号的近似图。
2.参考例5-2,画出上题频谱图。
3. 编程实现以下功能:将以下信号分解成余弦形式的傅立叶级数,再由信号f的直流分量C0开始,依次相加其基波分量、二次谐波、三次谐波等直至N次谐波(N不超过信号f的采样点数量的一半),观察随着谐波次数的增加,前N次谐波之和与原信号f的相似关系,讨论傅立叶系数的物理意义。
4.求以下信号的傅里叶变换F(w),并画出F(w)的幅频、相频图。
5.求以下信号的傅里叶反变换,并画出时域图。
【实验感悟】
【实验目的】
1.学会运用MATLAB分析傅里叶级数展开,深入理解傅里叶级数的物理意义。
2.学会运用MATLAB分析周期信号的频谱特性。
3.学会运用MATLAB求连续非周期信号的傅里叶变换。
4.学会运用MATLAB绘制连续非周期信号的频谱图。
【实验设备】
- 计算机
- MATLAB软件
【实验内容】
1.参考例5-1,实现教材p125,例3-4中傅里叶级数表达式(p126第二行)。分别采用前4、40、400项,画出周期矩形脉冲信号的近似图。
t=-1:0.001:1;%设置起点、步长和终点
ft=0.5*square(2*pi*(t+0.25),50)+0.5;%函数表达式
subplot(221);%绘制子图命令
plot(t,ft);%绘图命令
grid on;
axis([-1,1,-2,2]);%设置坐标轴
title('周期方波信号');
n_max=[4 40 400];
N=length(n_max);
for k=1:N%用循环实现傅里叶级数求和
n=1:2:n_max(k);
b=2./(pi*n);
x=0.5+sin(n*pi/2).*b*cos(2*pi*n'*t);
subplot(2,2,k+1);
plot(t,x);
grid on;
axis([-1,1,-1.5,1.5]);
title(['最大谐波次数=',num2str(n_max(k))]);
end
2.参考例5-2,画出上题频谱图。
n=-30:30;
tau=0.5;
T=1;
w=2*pi;
x=n*tau/T;
fn=tau/T*sinc(x);
stem(n*w,fn);
grid on;
title('\tau=0.5,T=1')
3. 编程实现以下功能:将以下信号分解成余弦形式的傅立叶级数,再由信号f的直流分量C0开始,依次相加其基波分量、二次谐波、三次谐波等直至N次谐波(N不超过信号f的采样点数量的一半),观察随着谐波次数的增加,前N次谐波之和与原信号f的相似关系,讨论傅立叶系数的物理意义。
Um=1;
T=0.02;
w=2*pi/T;
num_points=200;
t=linspace(-T/2,T/2,num_points);
y=Um*abs(sin(2*w*t)).*(t>0);%分解的波形
subplot(311);
plot(t,y);
grid on;
n=randn(38,1);
n=[0;n;0];
n=interp1([0:39],n,linspace(0,39,num_points),'linear');
y=y+0.4*n;
subplot(312)
plot(t,y);
grid on;
axis([-T/2,T/2,-1,2]);
c=50*trapz(t,y); %利用积分函数
i=sqrt(-1);
for k=1:40
y1=y.*cos(k*w*t);
y2=y.*sin(k*w*t);
an=100*trapz(t,y1);
bn=100*trapz(t,y2);
p=-atan(bn/an);
cn=sqrt(an^2+bn^2);
c=c+an*cos(k*w*t)+bn*sin(k*w*t);
subplot(313)
plot(t,c);
grid on;
axis([-T/2,T/2,-1,2]);
title('叠加后')
pause(0.1)
end
4.求以下信号的傅里叶变换F(w),并画出F(w)的幅频、相频图。
syms t ft;
ft=exp(-abs(t))+3*dirac(t);
Fw=fourier(ft);%Fourier函数实现傅里叶变换
subplot(211);grid on;
fplot(abs(Fw));%利用绝对值实现求幅频
title('幅频特征图');
phase=atan(imag(Fw)/real(Fw));%反正切函数实现求相频
subplot(212);grid on;
fplot(phase);
title('相频特征图')
5.求以下信号的傅里叶反变换,并画出时域图。
syms t w;
Fw=3/(3*1i*w+2-w^2);%定义函数
ft=ifourier(Fw,t);%利用ifourier函数实现傅里叶逆变换
fplot(ft);grid on;
title('函数时域图')
【实验感悟】
通过本次实验,我学会了MATLAB中对于傅里叶级数展开的求解方法, 通过循环实现傅里叶级数求和,通过多次滤波实现周期矩形脉冲信号的近似图趋于方波图。也学会了画出信号的频谱图,掌握了函数linspace命令,linspace是Matlab中的均分计算指令,用于产生x1,x2之间的N点行线性的矢量,了解了基于梯形法则的数值积分公式trapz,掌握了通过fourier命令求解函数的傅里叶变换,其中对于求解傅里叶变换后的函数的幅频和相频的方法让我印象深刻,通过abs绝对值函数求幅频,通过atan求解反正切,imag、real分别求解虚部和实部。还有通过ifourier命令求解函数的傅里叶反变换。
通过实验也加深了对于傅里叶变换的物理意义的理解:任何连续测量的时序或信号,都可以表示为不同频率的正弦波信号的无限叠加。傅里叶变换利用直接测量到的原始信号,以累加方式来计算该信号中不同正弦波信号的频率、振幅和相位。反傅里叶变换算法也是一种累加处理,这样就可以将单独改变的正弦波信号转换成一个信号。因此,可以说,傅里叶变换将原来难以处理的时域信号转换成了易于分析的频域信号(信号的频谱),可以利用一些工具对这些频域信号进行处理、加工。最后还可以利用傅里叶反变换将这些频域信号转换成时域信号。在图像学中,从物理效果看,傅里叶变换是将图像从空间域转换到频率域,其逆变换是将图像从频率域转换到空间域。换句话说,傅里叶变换的物理意义是将图像的灰度分布函数变换为图像的频率分布函数,傅里叶逆变换是将图像的频率分布函数变换为灰度分布函数。
下一篇: 傅立叶变换的物理意义
推荐阅读
-
windows下进程间通信的(13种方法)-摘 要 本文讨论了进程间通信与应用程序间通信的含义及相应的实现技术,并对这些技术的原理、特性等进行了深入的分析和比较。 ---- 关键词 信号 管道 消息队列 共享存储段 信号灯 远程过程调用 Socket套接字 MQSeries 1 引言 ---- 进程间通信的主要目的是实现同一计算机系统内部的相互协作的进程之间的数据共享与信息交换,由于这些进程处于同一软件和硬件环境下,利用操作系统提供的的编程接口,用户可以方便地在程序中实现这种通信;应用程序间通信的主要目的是实现不同计算机系统中的相互协作的应用程序之间的数据共享与信息交换,由于应用程序分别运行在不同计算机系统中,它们之间要通过网络之间的协议才能实现数据共享与信息交换。进程间通信和应用程序间通信及相应的实现技术有许多相同之处,也各有自己的特色。即使是同一类型的通信也有多种的实现方法,以适应不同情况的需要。 ---- 为了充分认识和掌握这两种通信及相应的实现技术,本文将就以下几个方面对这两种通信进行深入的讨论:问题的由来、解决问题的策略和方法、每种方法的工作原理和实现、每种实现方法的特点和适用的范围等。 2 进程间的通信及其实现技术 ---- 用户提交给计算机的任务最终都是通过一个个的进程来完成的。在一组并发进程中的任何两个进程之间,如果都不存在公共变量,则称该组进程为不相交的。在不相交的进程组中,每个进程都独立于其它进程,它的运行环境与顺序程序一样,而且它的运行环境也不为别的进程所改变。运行的结果是确定的,不会发生与时间相关的错误。 ---- 但是,在实际中,并发进程的各个进程之间并不是完全互相独立的,它们之间往往存在着相互制约的关系。进程之间的相互制约关系表现为两种方式: ---- (1) 间接相互制约:共享CPU ---- (2) 直接相互制约:竞争和协作 ---- 竞争——进程对共享资源的竞争。为保证进程互斥地访问共享资源,各进程必须互斥地进入各自的临界段。 ---- 协作——进程之间交换数据。为完成一个共同任务而同时运行的一组进程称为同组进程,它们之间必须交换数据,以达到协作完成任务的目的,交换数据可以通知对方可以做某事或者委托对方做某事。 ---- 共享CPU问题由操作系统的进程调度来实现,进程间的竞争和协作由进程间的通信来完成。进程间的通信一般由操作系统提供编程接口,由程序员在程序中实现。UNIX在这个方面可以说最具特色,它提供了一整套进程间的数据共享与信息交换的处理方法——进程通信机制(IPC)。因此,我们就以UNIX为例来分析进程间通信的各种实现技术。 ---- 在UNIX中,文件(File)、信号(Signal)、无名管道(Unnamed Pipes)、有名管道(FIFOs)是传统IPC功能;新的IPC功能包括消息队列(Message queues)、共享存储段(Shared memory segment)和信号灯(Semapores)。 ---- (1) 信号 ---- 信号机制是UNIX为进程中断处理而设置的。它只是一组预定义的值,因此不能用于信息交换,仅用于进程中断控制。例如在发生浮点错、非法内存访问、执行无效指令、某些按键(如ctrl-c、del等)等都会产生一个信号,操作系统就会调用有关的系统调用或用户定义的处理过程来处理。 ---- 信号处理的系统调用是signal,调用形式是: ---- signal(signalno,action) ---- 其中,signalno是规定信号编号的值,action指明当特定的信号发生时所执行的动作。 ---- (2) 无名管道和有名管道 ---- 无名管道实际上是内存中的一个临时存储区,它由系统安全控制,并且独立于创建它的进程的内存区。管道对数据采用先进先出方式管理,并严格按顺序操作,例如不能对管道进行搜索,管道中的信息只能读一次。 ---- 无名管道只能用于两个相互协作的进程之间的通信,并且访问无名管道的进程必须有共同的祖先。 ---- 系统提供了许多标准管道库函数,如: pipe——打开一个可以读写的管道; close——关闭相应的管道; read——从管道中读取字符; write——向管道中写入字符; ---- 有名管道的操作和无名管道类似,不同的地方在于使用有名管道的进程不需要具有共同的祖先,其它进程,只要知道该管道的名字,就可以访问它。管道非常适合进程之间快速交换信息。 ---- (3) 消息队列(MQ) ---- 消息队列是内存中独立于生成它的进程的一段存储区,一旦创建消息队列,任何进程,只要具有正确的的访问权限,都可以访问消息队列,消息队列非常适合于在进程间交换短信息。 ---- 消息队列的每条消息由类型编号来分类,这样接收进程可以选择读取特定的消息类型——这一点与管道不同。消息队列在创建后将一直存在,直到使用msgctl系统调用或iqcrm -q命令删除它为止。 ---- 系统提供了许多有关创建、使用和管理消息队列的系统调用,如: ---- int msgget(key,flag)——创建一个具有flag权限的MQ及其相应的结构,并返回一个唯一的正整数msqid(MQ的标识符); ---- int msgsnd(msqid,msgp,msgsz,msgtyp,flag)——向队列中发送信息; ---- int msgrcv(msqid,cmd,buf)——从队列中接收信息; ---- int msgctl(msqid,cmd,buf)——对MQ的控制操作; ---- (4) 共享存储段(SM) ---- 共享存储段是主存的一部分,它由一个或多个独立的进程共享。各进程的数据段与共享存储段相关联,对每个进程来说,共享存储段有不同的虚拟地址。系统提供的有关SM的系统调用有: ---- int shmget(key,size,flag)——创建大小为size的SM段,其相应的数据结构名为key,并返回共享内存区的标识符shmid; ---- char shmat(shmid,address,flag)——将当前进程数据段的地址赋给shmget所返回的名为shmid的SM段; ---- int shmdr(address)——从进程地址空间删除SM段; ---- int shmctl (shmid,cmd,buf)——对SM的控制操作; ---- SM的大小只受主存限制,SM段的访问及进程间的信息交换可以通过同步读写来完成。同步通常由信号灯来实现。SM非常适合进程之间大量数据的共享。 ---- (5) 信号灯 ---- 在UNIX中,信号灯是一组进程共享的数据结构,当几个进程竞争同一资源时(文件、共享内存或消息队列等),它们的操作便由信号灯来同步,以防止互相干扰。 ---- 信号灯保证了某一时刻只有一个进程访问某一临界资源,所有请求该资源的其它进程都将被挂起,一旦该资源得到释放,系统才允许其它进程访问该资源。信号灯通常配对使用,以便实现资源的加锁和解锁。 ---- 进程间通信的实现技术的特点是:操作系统提供实现机制和编程接口,由用户在程序中实现,保证进程间可以进行快速的信息交换和大量数据的共享。但是,上述方式主要适合在同一台计算机系统内部的进程之间的通信。 3 应用程序间的通信及其实现技术 ---- 同进程之间的相互制约一样,不同的应用程序之间也存在竞争和协作的关系。UNIX操作系统也提供一些可用于应用程序之间实现数据共享与信息交换的编程接口,程序员可以通过自己编程来实现。如远程过程调用和基于TCP/IP协议的套接字(Socket)编程。但是,相对普通程序员来说,它们涉及的技术比较深,编程也比较复杂,实现起来困难较大。 ---- 于是,一种新的技术应运而生——通过将有关通信的细节完全掩盖在某个独立软件内部,即底层的通讯工作和相应的维护管理工作由该软件内部来实现,用户只需要将通信任务提交给该软件去完成,而不必理会它的具体工作过程——这就是所谓的中间件技术。 ---- 我们在这里分别讨论这三种常用的应用程序间通信的实现技术——远程过程调用、会话编程技术和MQSeries消息队列技术。其中远程过程调用和会话编程属于比较低级的方式,程序员参与的程度较深,而MQSeries消息队列则属于比较高级的方式,即中间件方式,程序员参与的程度较浅。 ---- 4.1 远程过程调用(RPC)
-
信号与系统实验室 V 信号的傅立叶分析和频谱