Welcome Guest! To enable all features please Login. New Registrations are disabled.

Notification

Icon
Error

Login


3 Pages123>
Options
Go to last post Go to first unread
Offline oscampo  
#1 Posted : 29 May 2012 23:14:49(UTC)
oscampo


Rank: Advanced Member

Groups: Registered
Joined: 10/12/2009(UTC)
Posts: 238
Man
Colombia
Location: Cali, Colombia

Was thanked: 79 time(s) in 60 post(s)
Hi, this is a snippet for numerical inverse Laplace transform. You can find time response of a system expresed in F(s) form.

enjoy


numerical inverse Laplace Transform


Download it here

Regards,

Oscar Campo
thanks 3 users thanked oscampo for this useful post.
on 30/05/2012(UTC),  on 30/05/2012(UTC),  on 30/05/2012(UTC)

Wanna join the discussion?! Login to your SMath Studio Forum forum account. New Registrations are disabled.

Offline omorr  
#2 Posted : 30 May 2012 10:09:26(UTC)
omorr


Rank: Administration

Groups: Registered, Advanced Member
Joined: 23/06/2009(UTC)
Posts: 1,740
Man
Serbia

Was thanked: 318 time(s) in 268 post(s)
Hello Oscar,

Thank you for this example Good

By the way, I find very interesting that you defined an expression in this way an then plot it:

Function<->variable

ft() looks like a function, but it is not.

Regards,
Radovan
When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!"
thanks 1 user thanked omorr for this useful post.
on 30/05/2012(UTC)
Offline oscampo  
#3 Posted : 30 May 2012 14:10:43(UTC)
oscampo


Rank: Advanced Member

Groups: Registered
Joined: 10/12/2009(UTC)
Posts: 238
Man
Colombia
Location: Cali, Colombia

Was thanked: 79 time(s) in 60 post(s)
Yes, it's something that I found accidentally. The statement ft(F(s)) appear as a function, but it's a matrix!

Oscar
Offline Basile  
#4 Posted : 30 May 2012 22:18:43(UTC)
Basile


Rank: Advanced Member

Groups: Registered
Joined: 13/10/2010(UTC)
Posts: 171
Location: USA

Was thanked: 8 time(s) in 7 post(s)
Hello,

I tried to open it using the latest SMath.
I get an error message in the last FOR loop: "Function not defined"

Thanks,
Basile
Offline Andrey Ivashov  
#5 Posted : 30 May 2012 22:53:17(UTC)
Andrey Ivashov


Rank: Administration

Groups: Developers, Registered, Knovel Developers, Administrators, Advanced Member
Joined: 11/07/2008(UTC)
Posts: 1,616
Man
Russian Federation

Was thanked: 1978 time(s) in 666 post(s)
Basile Offline , you need to define a function to work with as shown in the first message.

Originally Posted by: oscampo Go to Quoted Post
numerical inverse Laplace Transform


Regards.
Offline Basile  
#6 Posted : 01 June 2012 04:25:59(UTC)
Basile


Rank: Advanced Member

Groups: Registered
Joined: 13/10/2010(UTC)
Posts: 171
Location: USA

Was thanked: 8 time(s) in 7 post(s)
Thank you Andrey,

Basile
Offline Basile  
#7 Posted : 01 June 2012 04:28:48(UTC)
Basile


Rank: Advanced Member

Groups: Registered
Joined: 13/10/2010(UTC)
Posts: 171
Location: USA

Was thanked: 8 time(s) in 7 post(s)
Oscar,

Which Numerical Inverse Laplace Transform do you use?
Any reference?

Thank you,
Basile
Offline oscampo  
#8 Posted : 01 June 2012 16:24:33(UTC)
oscampo


Rank: Advanced Member

Groups: Registered
Joined: 10/12/2009(UTC)
Posts: 238
Man
Colombia
Location: Cali, Colombia

Was thanked: 79 time(s) in 60 post(s)
Hi Basile,
The method is based on work of J. Valsa and L. Brancik: Approximate Formulae for Numerical Inversion of Laplace Transforms, Int. Journal of Numerical Modelling: Electronic Networks, Devices and Fields, Vol. 11, (1998), pp. 153-166

Regards,

Oscar
thanks 1 user thanked oscampo for this useful post.
on 01/06/2012(UTC)
Offline Basile  
#9 Posted : 02 June 2012 00:29:58(UTC)
Basile


Rank: Advanced Member

Groups: Registered
Joined: 13/10/2010(UTC)
Posts: 171
Location: USA

Was thanked: 8 time(s) in 7 post(s)
Oscar,

Thank you,

Basile
Offline oscampo  
#10 Posted : 05 June 2012 02:47:33(UTC)
oscampo


Rank: Advanced Member

Groups: Registered
Joined: 10/12/2009(UTC)
Posts: 238
Man
Colombia
Location: Cali, Colombia

Was thanked: 79 time(s) in 60 post(s)
Bode Diagram Revisited.
This is a function for graphing the Bode plot of a system. Please plot the magnitude and phase on a diferent graph each one, deactivate "grid" and "axis" for a better visualization of logarithmic grid.

Bode Plot

Download it here

Edited by user 25 November 2019 18:33:45(UTC)  | Reason: links updated

Offline oscampo  
#11 Posted : 05 June 2012 22:02:37(UTC)
oscampo


Rank: Advanced Member

Groups: Registered
Joined: 10/12/2009(UTC)
Posts: 238
Man
Colombia
Location: Cali, Colombia

Was thanked: 79 time(s) in 60 post(s)
Root Locus Snippet

This is a snippet for plotting Root Locus of a sistem. The input is the denominator of coefficients of transfer function:

rlocus

Download it here

Oscar

Edited by user 25 November 2019 18:37:55(UTC)  | Reason: links updated

Offline Leandro de Oliveira  
#12 Posted : 25 August 2017 19:22:41(UTC)
Leandro de Oliveira

Rank: Newbie

Groups: Registered
Joined: 25/08/2017(UTC)
Posts: 1
Brazil

The file seems to be locked, I can't put a equation.
Offline Jean Giraud  
#13 Posted : 25 August 2017 23:02:19(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,868
Canada

Was thanked: 980 time(s) in 808 post(s)
.. here is the Bode Plot.

Bode plot.sm (28kb) downloaded 148 time(s).
thanks 2 users thanked Jean Giraud for this useful post.
on 27/08/2017(UTC),  on 16/11/2017(UTC)
Offline Davide Carpi  
#14 Posted : 26 August 2017 00:44:35(UTC)
Davide Carpi


Rank: Advanced Member

Groups: Registered, Advanced Member
Joined: 13/01/2012(UTC)
Posts: 2,639
Man
Italy
Location: Italy

Was thanked: 1323 time(s) in 873 post(s)
Originally Posted by: Leandro de Oliveira Go to Quoted Post
The file seems to be locked, I can't put a equation.


Open the file in SMath Studio and then go to File > Properties > File Attributes and uncheck Read-only. This attribute is required when the worksheet is designed to be a snippet.
If you like my plugins consider to support SMath Studio buying a plan; to offer me a coffee: paypal.me/dcprojects
Offline Jean Giraud  
#15 Posted : 26 August 2017 15:01:49(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,868
Canada

Was thanked: 980 time(s) in 808 post(s)
Originally Posted by: oscampo Go to Quoted Post
Hi, this is a snippet for numerical inverse Laplace transform. You can find time response of a system expressed in F(s) form. enjoy

=================
This 5 years old document does not work 6179.
Here is the Maple companion.

Laplace 5Years.PNG

Offline oscampo  
#16 Posted : 27 August 2017 16:46:55(UTC)
oscampo


Rank: Advanced Member

Groups: Registered
Joined: 10/12/2009(UTC)
Posts: 238
Man
Colombia
Location: Cali, Colombia

Was thanked: 79 time(s) in 60 post(s)
Originally Posted by: Jean Giraud Go to Quoted Post

=================
This 5 years old document does not work 6179.



Only a very little detail to fix this 5 years old document: declare F(s) and change the last line on snippet:

ft( "F(s" ) = augment(t,ft)

for

ft() = augment(t,ft)

InvLaplace update.png




You can download the updated version from here

Edited by moderator 28 August 2017 19:27:40(UTC)  | Reason: edited dropbox link to direct download

thanks 1 user thanked oscampo for this useful post.
on 28/08/2017(UTC)
Offline Jean Giraud  
#17 Posted : 28 August 2017 17:01:53(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,868
Canada

Was thanked: 980 time(s) in 808 post(s)
Originally Posted by: oscampo Go to Quoted Post
You can download the updated version from here


... that's only several pages long of HTML code.
No work sheet. Best is to upload WS in the forum.
For "Laplace stuff", better use maple(,) as show above.

Cheers, Jean

Offline oscampo  
#18 Posted : 28 August 2017 17:32:53(UTC)
oscampo


Rank: Advanced Member

Groups: Registered
Joined: 10/12/2009(UTC)
Posts: 238
Man
Colombia
Location: Cali, Colombia

Was thanked: 79 time(s) in 60 post(s)
Originally Posted by: Jean Giraud Go to Quoted Post

... that's only several pages long of HTML code.
No work sheet. Best is to upload WS in the forum.


Use "save as"

Quote:

For "Laplace stuff", better use maple(,) as show above.


Yes, and also it’s possible use f(t)=Maxima(ilt(F(s),s,t)) but it’s interesting show how you can implement inverse Laplace of a function in SMath.
P.D. The symbolic form of f(t) obtained from ilt() it’s more “elegant” than maple version Wink
thanks 1 user thanked oscampo for this useful post.
on 28/08/2017(UTC)
Offline Jean Giraud  
#19 Posted : 28 August 2017 19:04:34(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,868
Canada

Was thanked: 980 time(s) in 808 post(s)
Win 7 does not have the virtue to save HTML "save as" *.sm
OK: no work sheet => nothing to see. Tanks anyway.

Cheers, Jean
Offline Davide Carpi  
#20 Posted : 28 August 2017 19:24:42(UTC)
Davide Carpi


Rank: Advanced Member

Groups: Registered, Advanced Member
Joined: 13/01/2012(UTC)
Posts: 2,639
Man
Italy
Location: Italy

Was thanked: 1323 time(s) in 873 post(s)
Originally Posted by: Jean Giraud Go to Quoted Post
Win 7 does not have the virtue to save HTML "save as" *.sm


You can if you use firefox (with chrome you have to change the extension manually from .txt to .sm);

However I've edited the link, now it downloads the file, whatever is your favorite browser. Whistle

Edited by user 28 August 2017 19:26:58(UTC)  | Reason: Not specified

If you like my plugins consider to support SMath Studio buying a plan; to offer me a coffee: paypal.me/dcprojects
thanks 1 user thanked Davide Carpi for this useful post.
on 28/08/2017(UTC)
Users browsing this topic
3 Pages123>
Forum Jump  
You cannot post new topics in this forum.
You cannot reply to topics in this forum.
You cannot delete your posts in this forum.
You cannot edit your posts in this forum.
You cannot create polls in this forum.
You cannot vote in polls in this forum.