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中周期延拓,如何周期延拓函数?的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 11.云计算平台(数据科学概论)
- 下一篇: SQL Server 查询分析器提供的所