Лучший способ представить это - PullRequest
Новая
       28

Лучший способ представить это

0 голосов
/ 26 ноября 2018

Я пытаюсь построить ДОС бесконечной линейной цепи.По следующему: 

s = 10; (*Number of decimation*) 
t = 1;
step = 0.01;
SubStar[E] = w + I*\[Eta];
\[Eta] = 0.01;
\[Epsilon] = 0;


Subscript[\[Beta], 1] = t*1/(SubStar[E] - \[Epsilon])*t;

Subscript[\[Beta], 0] = 1;



For[j = 2, j < s + 1, j++, 
 Subscript[\[Beta], j] = 
  Subscript[\[Beta], j - 1]*1/(SubStar[E] - \[Epsilon] - 2 \!\(
\*SubsuperscriptBox[\(\[Sum]\), \(i = 1\), \(j - 1\)]
\*SubscriptBox[\(\[Beta]\), \(i\)]\))*Subscript[\[Beta], j - 1]];

Subscript[G, 1, 1] = 1/(SubStar[E] - \[Epsilon] - \!\(
\*SubsuperscriptBox[\(\[Sum]\), \(j = 1\), \(s\)]
\*SubscriptBox[\(\[Beta]\), \(j\)]\));

G22 = (1/(1 - Subscript[G, 1, 1]*Subscript[G, 1, 1]))*Subscript[G, 1, 
   1];

\[Rho] = -1/\[Pi]*Im[G22];

Я хочу построить [Rho] как функцию w.Итак, если я позволю w неопределенным в начале (код выше) и заговор: ListPlot[Table[{w, \[Rho]}, {w, -2.5, 2.5, step}], Joined -> True, Frame -> True]
Это работает хорошо, но очень медленно, возможно, потому что mathematica делает все это алгебраически.Если я изменю w, например, используя: w = Range [-2.5,2.5, step], я могу построить что-то вроде этого: 

ListPlot[{\[Rho]}, Joined -> True, Frame -> True, PlotStyle -> Orange]

Это работает довольно быстро, но горизонтальная ось не в порядкепо сравнению с предыдущим сюжетом.

Затем я подумал сделать что-то вроде, определив 

 \[Rho][w]:= Module[{.....},....]

И поработав с этим как с функцией и графиком в конце.Какие-либо предложения?Что я должен делать?

Edit1: после предложения:

Сравнение между ними

1 Ответ

0 голосов
/ 27 ноября 2018

Почему бы не просто это? 

Plot[\[Rho], {w, -2.5, 2.5}, Frame -> True],

или еще лучше 

Plot[\[Rho], {w, -2.5, 2.5}, Frame -> True, PlotRange -> All]

Возможно, попробуйте это перед построением графика, но это отнимает много времени. 

\[Rho] = Simplify[-1/\[Pi]*Im[G22]]

Я попробовал это, что было интересно ... 

For[j = 2, j < s + 1, j++, 
 Print[j];
  Subscript[\[Beta], j] = 
  Simplify[Subscript[\[Beta], 
    j - 1]*1/(SubStar[E] - \[Epsilon] - 2 \!\(
\*SubsuperscriptBox[\(\[Sum]\), \(i = 1\), \(j - 1\)]
\*SubscriptBox[\(\[Beta]\), \(i\)]\))*Subscript[\[Beta], j - 1]]]

10-я итерация - убийца.Прервано этим сообщением: -

enter image description here

...