如何同matlab實現拉普拉斯變換

時間 2021-08-30 11:05:14

1樓:木生子識時務

fs=laplace(ft,t,s)   %功能:求“時域”函式ft的laplace變換fs;

ft=ilaplace(fs,s,t)   %功能:求“頻域”函式fs的laplace變換ft。

2樓:匿名使用者

拉普拉斯(laplace)變換是求解微分方程的一個有力的工具,是將微積分運算轉變為乘除運算,大大提供運算效率。

用matlab的實現拉普拉斯變換的函式是laplace(),其逆變換是ilaplace()。

例1:求函式 y=sin2t 的 laplace 變換。

syms t f f

f=sin(2*t ) %原函式

f=laplace(f) %象函式

f =2/(s^2 + 4)

例2:求函式 1/(s(s²+5)) 的 laplace 逆變換。

syms s f f

f=1/(s*(s^2+5)) %象函式

f=ilaplace(f) %原函式

f =1/5 - cos(5^(1/2)*t)/5

例3:求方程y"+2y'-3y=exp(-t),滿足初始條件y(0)=0,y'(0)=1的解。

解:對方程的兩邊取拉氏變換,並考慮到初始條件,則得

s²y(s)-1+2sy(s)-3y(s)=1/(s+1)

以下用matlab求解。

ys=solve('s^2*y-1+2*s*y-3*y=1/(s+1)','y');

simplify(ilaplace(ys))

ans =

(3*exp(t))/8 - exp(-3*t)/8 - exp(-t)/4

求解得微分方程的解

y(t)=(3*exp(t))/8 - exp(-3*t)/8 - exp(-t)/4

matlab中的拉氏變換和反變換!

3樓:匿名使用者

f(s)=4/s*(s+2)

syms s;

f=4/s*(s+2);

ilaplace(f)

ans =

4*dirac(t)+8

可將一個有引數實數t(t≥ 0)的函式轉換為一個引數為複數s的函式。如果對於實部σ >σc的所有s值上述積分均存在,而對σ ≤σc時積分不存在,便稱 σc為f(t)的收斂係數。

syms函式功能:matlab中,syms函式用於建立符號物件。

擴充套件資料:

拉氏變換和反變換syms函式功能:

>> syms x y z

>> e = sym('e');

>> z = e ^ x * sin(y) + e ^ y * sin(x)

z = e^x*sin(y) + e^y*sin(x)

>> diff(z, 'x')

ans = e^y*cos(x) + e^x*log(e)*sin(y)

>> diff(z, 'y')

ans = e^x*cos(y) + e^y*log(e)*sin(x)

>> y = sin(x)

y = sin(x)

>> int(y)

ans = -cos(x)

syms的功能和sym函式相同,但syms可以同時建立多個符號物件,因此在建立多個符號變數時語法上要比使用sym簡單。

4樓:碎裂的記憶

由於你式子後面換行了,不知道是不是f(s)=4/s*(s+2)反正輸入如下:

syms s;

f=4/s*(s+2);

ilaplace(f)

ans =

4*dirac(t)+8

如果方程不對,改一下就行。

還有不明白的問我

Matlab如何實現樹結構,Matlab 如何實現樹結構

shine戚七七 由於matlab函式只能傳值傳遞引數的限制,操作起來稍有點不同。建立初始化連結串列空間 l,maintain initlist n,a1 a2 a3 例子 現在需要建立一樹,每個節點包含資料域,左指標域,和右指標域,估計需要的空間是100 l,maintain initlist 1...

matlab中如何實現座標軸旋轉

1先繪製個基本圖形 x pi 2 0.01 pi plot x,sin x 2 上面座標軸長了,只要 pi 2到pi的 set gca,xlim pi 2 pi 3座標軸只顯示整數和.5,想制定出 pi 2 pi 4 0 pi 4 pi 2 pi 3 4 pi 座標 set gca,xtick pi...

如何利用多核電腦實現matlab的並行運算

matlab並行運算 目前,新購置的電腦大部分都是多核的了,使用matlab進行大量計算時如何有效利用多核呢?matlab目前版本已經比較好的支援多核並行運算了。是用的matlab版本是r2007b。電腦是雙核的。先簡單試試 matlabpool local 2 submitted parallel...