In[]:=
OrthopticCurveCardioid::usage="It is the geometric locus of the intersection point of two tangent lines to a curve and orthogonal to each other. References: mathecurve.com: Orthoptic curve; Personal.";
Development:
Obtaining the parametric equations from the orthoptic curve to the cardioid:
In[]:=
ζ[t_]:={aCos[t](1+Cos[t]),aSin[t](1+Cos[t])};​​x[t_]:=aCos[t](1+Cos[t]);​​y[t_]:=aSin[t](1+Cos[t]);​​
y'[t]
x'[t]
//FullSimplify
Out[]=
-Cot
3t
2

In[]:=
-Cot
3t
2
;
In[]:=
m[t_]:=-Cot
3t
2
;
Perpendicularity condition (be two points of the cycloid
t
1
and
t
2
):
In[]:=
Cot
3
t
1
2
Cot
3
t
2
2
==-1;
This condition is only met if the difference of the angles is
π
2
:
In[]:=
3
t
2
2
-
3
t
1
2
==
π
2
;
t
2
==
t
1
+
π
3
;
In[]:=
​​Solvey-y[t]==m[t](x-x[t]),y-yt+
π
3
==mt+
π
3
x-xt+
π
3
,{x,y}//FullSimplify
Out[]=
x
3
8
(4+9Cos[t]+3Cos[2t]-3
3
(Sin[t]+Sin[2t])),y
9
8
Csc
3t
2
Sec
3t
2
Cos
1
6
(π-6t)+Cos
π
6
+2t+Sin[t]+Sin[2t]Sin[3t]
In[]:=
x
1
4
a2+3Cos[t]+3Cos[2t]+3Sin
1
6
(π-6t)-3Sin
π
6
+2t,y
3
8
a
3
Cos[t]+
3
Cos[2t]+3Sin[t]+Sin[2t];​​
ζ
O
[t_]:=
1
4
a2+3Cos[t]+3Cos[2t]+3Sin
1
6
(π-6t)-3Sin
π
6
+2t,
3
8
a
3
Cos[t]+
3
Cos[2t]+3Sin[t]+Sin[2t];
Animation:
In[]:=
Manipulatea=3;​​ζ[t_]:={aCos[t](1+Cos[t]),aSin[t](1+Cos[t])};​​
ζ
O
[t_]:=
1
4
a2+3Cos[t]+3Cos[2t]+3Sin
1
6
(π-6t)-3Sin
π
6
+2t,
3
8
a
3
Cos[t]+
3
Cos[2t]+3Sin[t]+Sin[2t];​​​​Show​​ParametricPlot[ζ[θ],{θ,0,2π},PlotStyle{Blue,Thick}],​​ParametricPlot[
ζ
O
[t],{t,0,t1},ColorFunction->"BlueGreenYellow"],​​GraphicsBlack,PointSize[0.007],​​Point[ζ[t1]],Text[Style["
P
1
",14,Bold],ζ[t1],{0,-1.2}],​​Pointζt1+
π
3
,TextStyle["
P
2
",14,Bold],ζt1+
π
3
,{0,-1.2},​​Point[
ζ
O
[t1]],Text[Style["
ζ
O
",14,Bold],
ζ
O
[t1],{0,-1.2}],​​InfiniteLine[{ζ[t1],
ζ
O
[t1]}],InfiniteLineζt1+
π
3
,
ζ
O
[t1],AspectRatioAutomatic,​​AxesTrue,AxesOrigin{0,0},AxesLabel{x,y},PlotRange{{-2,8},{-6,6}},ImageSize900,Style["Orthoptic curve of a cardioid",Bold,Large],{{t1,0.00002,"Value (t1)"},0.00001,2π,0.00001},ControlPlacementTop,SaveDefinitions->True
Out[]=
​
Orthoptic curve of a cardioid
Value (t1)