Group Abstract Group Abstract

Message Boards Message Boards

ListPlot3D to show the dependence of a spectrum on a parameter?

Steve Bardwell
Posted 10 years ago

In solving a differential equation, I would like to visualize how the spectrum of the solutions depend on a parameter. The code I am working on is:

maxTime = 100; eqn = {\[CapitalTheta]''[t] + (1/(1 - 0.1 Cos[\[Alpha] *t])) Sin[\[CapitalTheta][t] ] == 0, \[CapitalTheta][0] == 0, \[CapitalTheta]'[0] ==0.1}
sol = ParametricNDSolveValue[eqn, \[CapitalTheta], {t, 0, maxTime}, {\[Alpha] }];
alphaTable = Table[sol[\[Alpha]], {t, 0, maxTime}];
data = Table[Periodogram[alphaTable[\[Alpha]]], {\[Alpha] , 1.5, 2.5, 0.01}];
lowRange = maxTime / 10;
highRange = maxTime / 4;
ListPlot3D[[data], PlotRange -> {{lowRange, highRange}, All}]

My thought is that since 'sol' is a function of t and alpha, I should be able to generate a table for a fixed alpha, and get the spectrum for that using Periodogram. Then, I am trying to build a table of all the spectra and use ListPlot3D to visualize it.

Somewhere I am missing something; while I can get each individual spectrum for a fixed value of alpha, I am failing to merge them into the kind of object that ListPlot3D can work with.
POSTED BY: Steve Bardwell
4 Replies
Sort By:
Steve Bardwell
Steve Bardwell
Posted 10 years ago

Again, thanks for all your help. I have studied your notebook, and I think something is wrong with the results. Your table manipulation logic eludes me, but the results don't make sense (for example, if you do a line plot of the spectrum at a single value of alpha, it doesn't look at all like the corresponding section of the 3D plot). After a day of playing around with it, I think I have a much simpler solution using PeriodogramArray[]:

maxTime = 100; eqn = {\[CapitalTheta]''[
     t] + (1/(1 - 0.1 Cos[\[Alpha]*t])) Sin[\[CapitalTheta][t]] == 
   0, \[CapitalTheta][0] == 0, \[CapitalTheta]'[0] == .1}
sol = ParametricNDSolveValue[
   eqn, \[CapitalTheta], {t, 0, maxTime}, {\[Alpha]}];
alphaTable[\[Alpha]_] := Table[sol[\[Alpha]][t], {t, 0, maxTime}];
spectrum[\[Alpha]_] := PeriodogramArray[alphaTable[\[Alpha]]];
data = Table[spectrum[\[Alpha]], {\[Alpha], 1.9, 2.1, 0.001}];
lowRange = maxTime/10;
highRange = maxTime/4;
ListPlot3D[data, PlotRange -> {{lowRange, highRange}, All, Full}]

enter image description here

The numbers on the axes need to be scaled (alpha goes from 0 to 200 rather than 1.9 to 2.1), but other than that, this result looks right to me.
POSTED BY: Steve Bardwell

Copy and paste works fine for me. However, I am enclosing the file.

Attachments:
POSTED BY: Gianluca Gorni
Steve Bardwell
Steve Bardwell
Posted 10 years ago
POSTED BY: Steve Bardwell

You have to extract the values from the Periodgram plots. Here is one way:

maxTime = 100; eqn = {\[CapitalTheta]''[
     t] + (1/(1 - 0.1 Cos[\[Alpha]*t])) Sin[\[CapitalTheta][t]] == 
   0, \[CapitalTheta][0] == 0, \[CapitalTheta]'[0] == 0.1};
sol = ParametricNDSolveValue[
   eqn, \[CapitalTheta], {t, 0, maxTime}, {\[Alpha]}];
alphaTable[\[Alpha]_] := Table[sol[\[Alpha]][t], {t, 0, maxTime}];
data = Flatten[
   Table[Flatten[
     Cases[Periodogram[alphaTable[\[Alpha]]], 
      Line[pts_] :> Map[Prepend[\[Alpha]], pts], All], 1], {\[Alpha], 
     1.5, 2.5, 0.01}], 1];
lowRange = maxTime/10;
highRange = maxTime/4;
ListPlot3D[data]
POSTED BY: Gianluca Gorni
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Tag limit exceeded
Note: Only the first five people you tag will receive an email notification; the other tagged names will appear as links to their profiles.
Reply Preview
Attachments
Remove
or Discard