matlab中周期延拓,如何周期延拓函数?
在傅里叶级数中我们需要将定义在一个区间[a,b]上的函数f(x)进行周期延拓。
如何写出这个延拓的函数F(x)的表达式?
答:这个延拓后的函数的表达式是:
F(x)=f(x-(b-a)*floor((x-a)/(b-a)) (公式1)
其中floor(u)就是我们所熟悉的取整函数[u]。
如果f(x)定义在区间[-a, a]上,则周期延拓后的函数为
F(x)=f(x-2*a*floor((x+a)/(2*a)) (公式2)
如果f(x)定义在区间[-Pi, Pi]上,则周期延拓后的函数为
F(x)=f(x-2*Pi*floor((x+Pi)/(2*Pi)) (公式3)
例1 将函数f(x)=abs(x) (-Pi<=x<=Pi) 延拓为周期为2*Pi的函数F(x)。(同济大学《高等数学》下册,311页,例5)
解 根据公式3,周期延拓后的函数为
F(x)=abs(x-2*Pi*floor((x+Pi)/(2*Pi))
f(x)的图形
周期延拓后的函数F(x)的图形
作图的Mathematica程序:f[x_] := Piecewise[{{-x, -Pi <= x < 0}, {x, 0 <= x <= Pi}}]
L := Pi
F[x_] := f[x - 2*L*Floor[(x + L)/(2*L)]]
A = Plot[f[x], {x, -Pi, Pi}, PlotStyle -> {Red, AbsoluteThickness[3]},AspectRatio -> Automatic, Ticks -> {Range[-2*L, 2*L, Pi/2], Range[-3, 3, 1]}]
B = Plot[F[x], {x, -5*L, 5*L}, AspectRatio -> Automatic];
Show[A, B, PlotRange -> {{-5*L, 5*L}, {-1, 4}}, Ticks -> {Range[-5*L, 5*L, Pi], Range[-3, 3, 1]}]
例2 将函数f(x)=-1 (-Pi<=x<0),f(x)=1(0<=x
解 根据公式3,周期延拓后的函数为
F(x)=f(x-2*Pi*floor((x+Pi)/(2*Pi))
f(x)的图形
周期延拓后的函数F(x)的图形
作图的Mathematica程序:f[x_] := Piecewise[{{-1, -Pi <= x < 0}, {1, 0 <= x <= Pi}}]
L := Pi
F[x_] := f[x - 2*L*Floor[(x + L)/(2*L)]]
A = Plot[f[x], {x, -Pi, Pi}, PlotStyle -> {Red, AbsoluteThickness[3]},AspectRatio -> 1/2, Ticks -> {Range[-2*L, 2*L, Pi/2], Range[-1/2, 3/2, 1/2]}]
B = Plot[F[x], {x, -4*L, 4*L}, AspectRatio -> 1/2, PlotStyle -> {Blue, AbsoluteThickness[2]}];
Show[B, A, PlotRange -> {{-4*L, 4*L}, {-2, 2}}, Ticks -> {Range[-5*L, 5*L, Pi], Range[-2, 2, 0.5]}]
例3 将函数f(x)=x^2-3*x (2<=x<=5) 延拓为周期为3 的函数F(x)
解 根据公式1,周期延拓后的函数为
F(x)=f(x-3*floor((x-2)/3)
f(x)的图形
周期延拓后的函数F(x)的图形
\
作图的Mathematica程序:f[x_] := Piecewise[{{x*Sin[x/2], 2 <= x < 5}}]
F[x_] := f[x - 3*Floor[(x - 2)/3]]
A = Plot[f[x], {x, -4, 11}, PlotStyle -> {Red, AbsoluteThickness[3]}, Ticks -> {Range[-4, 12, 1], Range[-1, 5, 1]}]
B = Plot[F[x], {x, -4, 11}, PlotStyle -> {Blue, AbsoluteThickness[2]}];
Show[B, A, Ticks -> {Range[-4, 12, 1], Range[-1, 10, 1]}]
matlab中周期延拓,如何周期延拓函数?相关推荐
- matlab得到计算机时间,Matlab中处理日期与时间的函数
Matlab用三种格式来表示日期与时间 (1)双精度型日期数字:一个日期型数字代表从公元0年到某一日期的的天数,例如,2008年8月26日夜时0点钟被表示为733646,而这同一天的中午12点就被表示 ...
- 将函数当做参数 matlab,Matlab中如何将(自定义)函数作为参数传递给另一个函数...
假如我们编写了一个积分通用程序,想使它更具有通用性,那么可以把被积函数也作为一个参数.在c/c++中,可以使用函数指针来实现上边的功能,在matlab中如何实现呢?使用函数句柄--这时类似于函数指针的 ...
- matlab中饼图绘制程序,MATLAB中绘制二维饼图的函数是(? ? ) (1.5分)
MATLAB中绘制二维饼图的函数是(? ? ) (1.5分) 答:pie 计算机的内存储器比外存储器读写速度 ______ 答:快 新民主主义革命时期,提出了反帝反封建这个任务就等于真正懂得了什么叫反 ...
- MatLab中的多项式拟合(lsqcurvefit函数的非线性曲线拟合)
MatLab中的多项式拟合(lsqcurvefit函数的非线性曲线拟合) lsqcurvefit函数 函数原型 x=lsqcurvefit(fun, x0, xdata, ydata) == 其中的参 ...
- MatLab中rand randi randn 三个函数的解析
MatLab中rand randi randn 三个函数的解析 1.rand 概括:生成(0,1)之间的均匀分布的伪随机数 示例1: >>rand %什么参数也没有,生成一个(0,1)的数 ...
- Matlab中一维和二维滤波函数及其典型用法【随时更新】
Matlab中一维和二维滤波函数及其典型用法[随时更新] 一维滤波函数 均值滤波器movmean 中值滤波器movmedian 二维滤波函数 用fspecial预定义各种二维滤波器 主要用于自查,如果 ...
- matlab对数收益直方图,科学网—MATLAB中绘制数据直方图的新函数histogram2 - 王福昌的博文...
MATLAB中有命令hist3() 可以绘制直方图,竖坐标是频数,这与一些教科书中用纵轴表示频率的做法不一致,有些时候不便于使用.当然,使用者可以自己编写定制能够在纵轴绘出频率的直方图.在MATLAB ...
- Matlab中filter.m和filtfilt.m函数C语言实现
一.一些基础知识 filter.m函数是依据z变换的一些知识进行的滤波方法. filtfilt.m则还有另一个名字是零相位滤波,顾名思义,通过filtfilt函数滤波后的信号,幅值会发生变化,但相位不 ...
- matlab中用于小数取整的函数的用法
matlab中小数取整的函数大约有四个:floor.ceil.round.fix 若 A = [-2.0, -1.9, -1.55, -1.45, -1.1, 1.0, 1.1, 1.45, 1.55 ...
- matlab变量区表示函数,MATLAB中的工作区,变量和函数
本文概述 工作空间 工作区包含我们在MATLAB中工作时创建的所有变量. 每当我们为变量分配值时, 它都会自动在工作空间中获取空间. 关闭环境后, 工作空间变量将消失, 因此请将这些变量保存在文件中以 ...
最新文章
- Optional 是个好东西,你真的会用么?
- 2dpsk调制解调实验matlab_【无线通信篇】你知道ASK调制吗?
- 解析Linux操作系统文件目录
- 日志组件logback介绍及配置使用方法
- Nginx的安全控制及SSL加密介绍
- 从零开始学前端:CSS盒子模型属性 --- 今天你学习了吗?(CSS:Day14)
- VMware ESXi 嵌套虚拟化 部署问题总结
- 博客园文章中图片太大显示不全的解决办法
- html3d房子立体图片,如何制作3D立体图片
- 即时通讯开发资料分享
- java动态时钟_Java实现动态模拟时钟
- win10专业版如何快速恢复出厂设置的方法
- 忍者必须死3突然服务器维修,《忍者必须死3》3月12日2:00停服维护公告
- 芯动力——硬件加速设计方法——学习笔记(1)
- 股票做空和做多是什么意思
- 80亿美元侵权诉讼的随想
- Psychopy | 第2期:从Stroop看条件与循环
- VPS究竟有什么用途呢?
- Hexo-fluid主题添加51LA统计
- BuuCTF难题详解| Misc | [HDCTF2019]你能发现什么蛛丝马迹吗
热门文章
- 程序君带你畅聊发送短信验证码
- 各路牛人【找工总结】
- mpc 安全多方计算协议_一文揭秘跨链、密钥管理、合约隐私背后的技术, 你不可不知的安全多方计算(MPC)...
- Vue eslint 报错 eval can be harmful解决办法
- 袁卫统计学第4版笔记和课后习题答案
- 中国文章剽窃严重 —— 腾讯有暗中保护盗版的嫌疑!
- 软件框架-无绪开发5
- “跑路风波”的内在缘由?P2P网络信贷将何去何从?
- 7-2 jmu-Java-03面向对象-06-继承覆盖综合练习-Person、Student、Employee、Company (15分)
- ios快捷指令:一键登录/登出南京大学校园网