计算机网络课程讲到物理层,布置作业的第一题是求f(t)=t (0≤t≤1)的傅里叶系数。

  我们知道任何一个周期函数都可以被傅里叶级数逼近。如果是实值函数,则可以用正弦分量,余弦分量,直流分量来近似。公式如下:$$f(x)=\frac{a_0}{2}+\sum_{n=1}^{\infty}(a_n cosnx+b_n sinnx)$$

  利用正余弦函数的正交性反变换可得到f(x)的系数。本题中认为函数周期是1,进行延拓可以看成是无数个f(t)沿着x轴平移。

  简单计算得$a_0=1$, $a_n=0$, $b_n=-\frac{1}{\pi n}$ 至此题目结束。

  但是我想练习一下最近自己刚刚接触到的numpy和matplotlib,于是动手写了一个python小程序想看看拟合的效果是什么样的。效果图和代码如下。  

import numpy as np
import matplotlib.pyplot as plt
x = np.linspace(0, 1, 100)
y_ = np.zeros_like(x)
y = x
y_ = y_ + 0.5
for n in range(1, 11):
y_ = y_ - 1/np.pi/n * np.sin(2*np.pi*x*n)
plt.plot(x, y, label='y=x')
plt.plot(x, y_, label='approximation')
plt.xlabel('x')
plt.ylabel('y')
plt.title("CFS")
plt.legend()
plt.show()

  代码一直拟合到十阶分量。可看到效果还是不错的,增加x采样点可以增加图像平滑

f(t) = t的傅里叶系数相关推荐

  1. python画方波_今天学会傅里叶画画,明天就是初音未来 - 如何用Python和Blender画任意图形...

    前言: 最近在[编程德鲁伊]的[傅里叶级数可视化]章节中,练习了傅里叶级数基础原理的可视化,演示了周期方波曲线,可以分解为一系列正弦波曲线: 这次拿学到的知识做个实验,试试喜闻乐见的[用傅里叶级数画画 ...

  2. 深入浅出通信原理2021-03-07

    系列文章目录 第一章 深入浅出通信原理之卷积 第二章 深入浅出通信原理之傅里叶变换 第三章 深入浅出通信原理之复指数信号 文章目录 系列文章目录 复指数信号 一.复指数信号的物理意义 二.余弦信号和正 ...

  3. 傅里叶级数、傅里叶变换 (FT)

    目录 Fourier Series (傅里叶级数) 三角函数的正交性 周期为 2 π 2\pi 2π 的函数展开为傅里叶级数 周期为 2 L 2L 2L 的函数展开为傅里叶级数 傅里叶级数的复数形式 ...

  4. 用python画叶子_今天学会傅里叶画画,明天就是初音未来 - 如何用Python和Blender画任意图形...

    前言: 最近在[编程德鲁伊]的[傅里叶级数可视化]章节中,练习了傅里叶级数基础原理的可视化,演示了周期方波曲线,可以分解为一系列正弦波曲线: 这次拿学到的知识做个实验,试试喜闻乐见的[用傅里叶级数画画 ...

  5. 数字信号处理-5-傅里叶分析

    傅里叶分析的主要思想就是将一段信号分解为不同频率正/余弦函数的叠加,不同频率正/余弦函数会有各自的系数和相位.通过傅里叶分析我们就能知道一段信号的频率组成. 1 傅里叶系数 傅里叶级数用公式如下: a ...

  6. 关于傅里叶系数和傅里叶变换中微分性质的思考

    关于傅里叶系数和傅里叶变换中微分性质的思考 Jzn原创,未经允许不可转载 0.写在前面 1.基本知识点的简单概述 2.傅里叶级数中微分性质的应用 3.傅里叶变换中微分性质的应用 4.小结 0. 写在前 ...

  7. 数字信号处理①之信号的合成、分解、相关性及傅里叶系数的计算

    数字信号处理系列博客: 数字信号处理之信号的合成.分解.相关性及傅里叶系数的计算 数字信号处理之从傅里叶系数到DFT 文章目录 一.信号的合成和分解 (一)傅里叶级数 (二)将正弦波合成方波 二.信号 ...

  8. 数字信号处理②之从傅里叶系数到DFT

    数字信号处理系列博客: 数字信号处理之信号的合成.分解.相关性及傅里叶系数的计算 数字信号处理之从傅里叶系数到DFT 文章目录 一.前言 二.对合成信号进行各次谐波分析 (一)合成信号 (二)使用复指 ...

  9. matlab 求矩形脉冲的傅里叶系数

    本节在 matlab heaviside函数生成矩形脉冲 和 matlab 编写计算傅里叶级数函数 的基础上,求一个矩形脉冲序列的傅里叶级数. 这里只求了一个周期矩形脉冲的傅里叶级数,因为我不知道怎么 ...

最新文章

  1. 51nod 1065 最小正字段和 解决办法:set存前缀和,二分插入和二分查找
  2. 个人贷款的条件、流程和注意事项
  3. Dart Hello World 示例程序【Dart专题1】
  4. line-height的一点见识
  5. css x轴不滚动_泡沫轴滚动放松是什么原理:你并不了解泡沫轴
  6. 简单的Net::FTP模块的简单应用
  7. .Net 指针使用以及在栈上创建引用对象
  8. Android Studio模拟器磁盘空间不足(Not enough disk space to run AVD)
  9. 在线思维导图工具-toolfk程序员在线工具网
  10. 量化交易模型只是武器,能否驾驭,靠的是背后的人
  11. 基于扩展卡尔曼滤波的多传感器融合定位
  12. 没有苹果开发账号,只有p12文件和mobileprovision文件进行打包
  13. 【分析】Ceph数据一致性检查 - Scrub的执行
  14. [C/C++]华为8-19笔试题
  15. 新浪云sae springboot项目部署
  16. dw中html文档命名标准,dw中的css是什么意思?
  17. 计算机学院开学标语,新学期开学宣传标语大全
  18. [指北针分类信息软件 v1.5.2.1] 全自动分类信息软件+高效稳定建立SEO外部链接
  19. PUMA560的运动学建模
  20. 虚拟机通过桥接上网(dhcp)

热门文章

  1. php168整站系统安装时所使用的数据库怎么填写,《PHP168整站系统使用大全》
  2. C语言实现--顺序创建链表(尾插法)
  3. 推荐10个常用的CSS动画库
  4. 国防科技大学计算机学院刘洋,记国防科技大学计算机学院计算机研究所所长肖_全封闭式学校...
  5. Angular-translate笔记
  6. c4c语言编译器,c4编译器源码剖析
  7. 虾米音乐与MQA公司达成国内首家合作 将提供MQA音质音源
  8. windows10正版多少钱_惠安正版OA系统好不好用,高性价比无纸化办公系统大概多少钱_泉州登泽...
  9. 本地使用element ui 图标不显示或显示错乱的问题
  10. bat文件重新启动打印服务器,服务器应用自动重新启动IIS批处理 原创