目录

【实验目的】

【实验设备】

【实验内容】

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绘制连续非周期信号的频谱图。

【实验设备】

  1. 计算机
  2. 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命令求解函数的傅里叶反变换。

通过实验也加深了对于傅里叶变换的物理意义的理解:任何连续测量的时序或信号,都可以表示为不同频率的正弦波信号的无限叠加。傅里叶变换利用直接测量到的原始信号,以累加方式来计算该信号中不同正弦波信号的频率、振幅和相位。反傅里叶变换算法也是一种累加处理,这样就可以将单独改变的正弦波信号转换成一个信号。因此,可以说,傅里叶变换将原来难以处理的时域信号转换成了易于分析的频域信号(信号的频谱),可以利用一些工具对这些频域信号进行处理、加工。最后还可以利用傅里叶反变换将这些频域信号转换成时域信号。在图像学中,从物理效果看,傅里叶变换是将图像从空间域转换到频率域,其逆变换是将图像从频率域转换到空间域。换句话说,傅里叶变换的物理意义是将图像的灰度分布函数变换为图像的频率分布函数,傅里叶逆变换是将图像的频率分布函数变换为灰度分布函数。

信号与系统实验五 信号的傅里叶分析与频谱相关推荐

  1. 信号与系统实验:信号抽样

    已知一个连续时间信号f(t)=sinc⁡(t)f(t)=sinc⁡(t)f(t)=sinc⁡(t),取最高有限带宽频率fm=1Hzf_m=1Hzfm​=1Hz (1)分别显示原连续信号波形和 fs=f ...

  2. matlab卷积实验原理,信号与系统实验报告——卷积(含程序)

    信号与系统实验报告--卷积 电 子 科 技 大 学 实 验 报 告 学生姓名:苏晓菁 学 号:2804301026 指导教师:张鹰 一.实验室名称:信号与系统实验室 二.实验项目名称:离散系统的冲激响 ...

  3. 2021年春季学期-信号与系统-第五次作业参考答案-第十一移小题—MATLAB

    本文是 2021年春季学期-信号与系统-第五次作业参考答案 中的小题解答. ▌第十一题:MABLAT 在MATLAB中,根据矩形周期脉冲信号傅里叶级数分解也锯齿波傅里叶级数分解的公式,绘制前N项级数叠 ...

  4. 2021年春季学期-信号与系统-第五次作业参考答案-第十小题

    本文是 2021年春季学期-信号与系统-第五次作业参考答案 中的小题解答. ▌第十题 10. 利用傅里叶变换公式计算下面信号的傅里叶变换. (1) e−2(t−1)⋅u(t−1)e^{ - 2\lef ...

  5. 2021年春季学期-信号与系统-第五次作业参考答案-第九小题

    本文是 2021年春季学期-信号与系统-第五次作业参考答案 中的小题解答. ▌第九题 9. 求下图所示的半波余弦脉冲的傅里叶变换, 并画出频谱图: 注:对于该信号的频谱计算,可以直接利用公式进行.其中 ...

  6. 2021年春季学期-信号与系统-第五次作业参考答案-第八小题

    本文是 2021年春季学期-信号与系统-第五次作业参考答案 中的小题解答. ▌第八题 8. 关于周期信号x(t)x\left( t \right)x(t)有如下相应的描述: (1) x(t)x\lef ...

  7. 2021年春季学期-信号与系统-第五次作业参考答案-第七小题

    本文是 2021年春季学期-信号与系统-第五次作业参考答案 中的小题解答. ▌第七题 7. 利用信号的对称性,确定下面各个周期信号的傅里叶级数所含有的频率成分: ※ 求解: (1)第一小题 偶对称.奇 ...

  8. 2021年春季学期-信号与系统-第五次作业参考答案-第六小题

    本文是 2021年春季学期-信号与系统-第五次作业参考答案 中的小题解答. ▌第六题 6. 已知周期信号 f(t)f\left( t \right)f(t)在其四分之一周期内的信号波形如下图所示: 分 ...

  9. 2021年春季学期-信号与系统-第五次作业参考答案-第五小题

    本文是 2021年春季学期-信号与系统-第五次作业参考答案 中的小题解答. ▌第五题 5. 设 f(t)f\left( t \right)f(t)为 [−π,π]\left[ { - \pi ,\pi ...

  10. 2021年春季学期-信号与系统-第五次作业参考答案-第四小题

    本文是 2021年春季学期-信号与系统-第五次作业参考答案 中的小题解答. ▌第四题 4. 下图所示的波形参数为:τ=5μs,T=10μs\tau = 5\mu s,\,\,T = 10\mu sτ= ...

最新文章

  1. 【Linux】Linux简单操作之安装jdk
  2. SharePoint的WebService的应用
  3. 在Asp.net页面中实现数据饼图
  4. mysql查找内容某字符串出现的次数_查找MySQL中特定字符串出现的次数?
  5. school‘s movie demonstration
  6. 服务器驱动停止响应,IE驱动程序的命令行服务器已停止工作
  7. 关于eclipse的indigo版中文注释时字体太小的问题(转)
  8. Tomcat 改BUG之 localhost:8080 404
  9. DAG最长路问题 hdu-1224
  10. 程序员很少加班?得全栈开发者得天下?撕开标签的技术圈真实模样
  11. 排序算法-冒泡排序详细讲解(BubbleSort)
  12. java uuid生成算法_Java生成UUID
  13. 电信光纤猫虚拟服务器设置方法,电信光纤无线路由器(光猫)的设置方法
  14. 利用Python下载网页视频
  15. java解压文件 文件过大,Java压缩/解压缩大文件( 1gb)
  16. Java编写图像浏览器_浏览器下载图像(JAVA代码)
  17. 在浏览器输入url回车后发生了什么
  18. 多线程的三种设计模式的介绍
  19. 突变检测软件 测试数据库,突变测试—通过一个简单的例子快速学习这种有趣的测试技术...
  20. 植物神经紊乱会不会导致失眠

热门文章

  1. 最新版网站推广完全手册(2007年)!
  2. 批量转换word为pdf
  3. [matlab也能用来机器学习!?]保存工具箱模型并使用模型预测结果
  4. 简析通达信股票接口测试过程
  5. python分词代码
  6. OpenOffice实现word文档在线预览
  7. linux缺少libeds.so,OpenSSL鏈接libcrypto.a以靜態方式
  8. mac下安装pyaudio
  9. html5json转换为数组,Json转数组 在线
  10. UG模具设计从入门到精通