HP PRIME -- How to plot a phase Bode plot ?



#2

Hi,

I want to create a program that is able to plot a curve without phase discontinuities on +/- pi rad... So The use of the ARG() command is not enough ... The classic tricks that I use with GNUPLOT can not be used here ...

Here is a description of what I'd like to get:

- first ask for transfer function F(s)(with input command)

F(s)= (1+10*s)^2/(1+s)^3 for example...

- use the subst cmd to replace 's' by 'i*ALOG(X)' where X is the pulsation in rad/s. So I can plot for Xmin <LOG(X)< Xmax

- create the F1 var function ready to plot in Function APP...
F1(X)=phase(X)...
where phase(X) is something like

phase(X)=ARG(F(i*ALOG(X))+2k*PI

The difficulty is to calculate k ...

On Gnuplot I look continuously during plotting process at ARG(f(xold))-ARG(f(x)) and if a jump occurs I increment or decrement k according to the sign of the jump...

Another idea is to use CAS command zero or solve and save all roots of arg(F(X))=+/-Pi in list L1

then create F1(X)=F1(X)+-2*PI*Heaviside(X-L1(n))in a loop for n=1 to size(L1)...

I'm looking for advice and assistance for such a program...

Thank you in advance if someone just spend a little time on this little problem ;o)


Possibly Related Threads…
Thread Author Replies Views Last Post
  HP PRIME - Bode and Black plots - first attempt dg1969 0 1,077 10-20-2013, 01:32 PM
Last Post: dg1969
  HP 48G -- entering more than one equation into the PLOT application Randal B 1 1,120 01-24-2013, 12:25 AM
Last Post: Chris Dreher
  Help with creating residual plot Phillip B. 2 1,251 10-10-2012, 02:35 PM
Last Post: Gilles Carpentier
  OT: Calculator as a plot device Ethan Conner 11 3,198 05-07-2012, 08:41 PM
Last Post: Matt Agajanian
  HP 39gII plot speed Alex Mark (Greece) 5 2,034 02-12-2012, 03:04 PM
Last Post: Gilles Carpentier
  50g won't plot summation function, but emulator will. Hal Bitton in Boise 2 1,228 04-12-2011, 07:14 PM
Last Post: Hal Bitton in Boise
  unable to plot UTPN stroran 3 1,433 04-21-2010, 12:10 PM
Last Post: Kiyoshi Akima
  50G, Plot-Table-Printing to 82240 IR printer George Bailey (Bedford Falls) 5 1,771 08-11-2008, 12:25 PM
Last Post: George Bailey (Bedford Falls)
  Annotated plot axes tick marks on the 50g? Ken Ratkevich 0 748 12-02-2006, 05:33 AM
Last Post: Ken Ratkevich
  CAS vs Plot Mike (Stgt) 4 1,461 11-22-2004, 03:54 AM
Last Post: Mike (Stgt)

Forum Jump: