Есть одно большое выражение для Maxima. На саму формулу смотреть не обязательно, поскольку она имеет вторичное значение. Главное - это тот факт, что H(ω), как ни странно, зависит от ω.
H(omega):=[(C[0]+C[1]*cos(omega*Delta[t])+C[2]*cos(2*omega*Delta[t]))
-%i*(C[1]*sin(omega*Delta[t])+C[2]*sin(2*omega*Delta[t]))]
/[(1-D[1]*cos(omega*Delta[t])-D[2]*cos(2*omega*Delta[t])
+%i*(D[1]*sin(omega*Delta[t])+D[2]*sin(2*omega*Delta[t])))];Есть вторая формула, для φ(ω), в которой входным параметром является ω, но высчитывается она путём определённых операций над H(ω):
phi(omega):=atan(realpart(H(omega))/imagpart(H(omega)));
Но когда я пытаюсь построить график φ(ω), maxima сообщает мне, что нет такой переменной, по которой он будет строиться:
plot2d([phi(omega)],[omega,0,15],
[plot_format,gnuplot]);
plot-options-parser: unknown plot option: [omega,0,15]
Как мне извернуться и таки построить этот график? Какие входные параметры указать?
Расходимся, парни. Я неправильно ввёл H(ω): там нельзя было использовать квадратные скобки.
