полиномиальных узлов в R

Данный вопрос о возможном пути графически представлять узловые многочлены, как даются выражениями here и here в R.полиномиальных узлов в R

К сожалению, я `t didn очень далеко, так как моя математика знание не так хорошо. Это то, что я попытался

phi <- seq(-30, 30, len=1001) 
fx <- phi^5 - 28*phi #these are the polynomials from the website above 
fy <- phi^7 - 32*phi^3

, который дает, при использовании участка (Fx, Fy) r plot of the above function

, и я предполагаю, что это должно было дать восьмерка-узел. Я могу найти нули вышеуказанных функций за указанный интервал, и все это выглядит нормально, поскольку обе функции имеют пары очень близких нулевых значений, что видно на диаграмме узлов. Есть идеи?

источник

2013-12-08 Ian Stuart

Ваш код верен. Вам просто нужно увеличить:

phi <- seq(-2.6, 2.6, len=1e5) 
fx <- phi^5 - 28*phi 
fy <- phi^7 - 32*phi^3 

plot(fx, fy, type="l")

enter image description here

источник

2013-12-08 16:16:03 Roland

Возможно, ОП должен изменить свой вопрос на «как найти различные значения' phi', для которых пары '(fx1, fy1), (fx2, fy2)' идентичны? Как только они будут найдены, для выбора параметров панорамирования и масштабирования для графика. –

@CarlWitthoft Мне нравится эта цитата из статьи OP, связанной с: «Однако мы решили рассчитать t-значения, которые дают двойные точки с помощью компьютера». И они не объясняют, как это сделать. – Roland

Карл прав, проблема в том, что параметры панорамирования и масштабирования различаются для разных пар полиномов, которые определяют разные узлы. Вопрос заключается в том, как найти универсальные параметры масштабирования или возможную процедуру ? –

Try:

phi <- seq(-3, 3, len=1001) 
fx <- phi^5 - 28*phi 
fy <- phi^7 - 32*phi^3 
plot(fx,fy,xlim=c(-50,50),ylim=c(-150,150))

Вы должны настроить пределы своей оси. HF

источник

2013-12-08 16:21:05 JPS

полиномиальных узлов в R

ответ

Смежные вопросы