Can an equation of the form $sin(omega t) + at^2+bt=N$, with $aneq 0$, be solved analytically?
$begingroup$
Can the equation of the form
$$sin(omega t) + at^2+bt=N$$
where $a,b,N in mathbb{R}$ be solved analytically for $a ne 0$?
It, for example, arises when the rotating object falls down subject to gravity. Solving such equation can allow to compute the time of the object touchdown.
trigonometry
$endgroup$
|
show 5 more comments
$begingroup$
Can the equation of the form
$$sin(omega t) + at^2+bt=N$$
where $a,b,N in mathbb{R}$ be solved analytically for $a ne 0$?
It, for example, arises when the rotating object falls down subject to gravity. Solving such equation can allow to compute the time of the object touchdown.
trigonometry
$endgroup$
$begingroup$
Not that I know of. Are you sure you got your equation right? The $sin omega t$ term implies an oscillation due to elastic force, not rotation.
$endgroup$
– Andrei
Dec 17 '18 at 5:53
$begingroup$
$ sin omega t$ is a harmonic oscillation due to the object rotation. It will be moving according to the quadratic law, but will also be spinning.
$endgroup$
– Jennifer M.
Dec 17 '18 at 5:54
$begingroup$
I am just curious how you got that equation.
$endgroup$
– Andrei
Dec 17 '18 at 5:57
$begingroup$
Imagine a rotating object in a freefall. Its Z coordinate is decreasing as $-frac{1}{2} g t^2$ and its rotation would add $sin omega t$ for each surface point.
$endgroup$
– Jennifer M.
Dec 17 '18 at 6:04
$begingroup$
If the domain of $sin omega t$ be $[0,π]$ then it is possible to get approximate solutions.
$endgroup$
– Awe Kumar Jha
Dec 17 '18 at 6:05
|
show 5 more comments
$begingroup$
Can the equation of the form
$$sin(omega t) + at^2+bt=N$$
where $a,b,N in mathbb{R}$ be solved analytically for $a ne 0$?
It, for example, arises when the rotating object falls down subject to gravity. Solving such equation can allow to compute the time of the object touchdown.
trigonometry
$endgroup$
Can the equation of the form
$$sin(omega t) + at^2+bt=N$$
where $a,b,N in mathbb{R}$ be solved analytically for $a ne 0$?
It, for example, arises when the rotating object falls down subject to gravity. Solving such equation can allow to compute the time of the object touchdown.
trigonometry
trigonometry
edited Dec 18 '18 at 4:15
Blue
48.6k870154
48.6k870154
asked Dec 17 '18 at 5:47
Jennifer M.Jennifer M.
1112
1112
$begingroup$
Not that I know of. Are you sure you got your equation right? The $sin omega t$ term implies an oscillation due to elastic force, not rotation.
$endgroup$
– Andrei
Dec 17 '18 at 5:53
$begingroup$
$ sin omega t$ is a harmonic oscillation due to the object rotation. It will be moving according to the quadratic law, but will also be spinning.
$endgroup$
– Jennifer M.
Dec 17 '18 at 5:54
$begingroup$
I am just curious how you got that equation.
$endgroup$
– Andrei
Dec 17 '18 at 5:57
$begingroup$
Imagine a rotating object in a freefall. Its Z coordinate is decreasing as $-frac{1}{2} g t^2$ and its rotation would add $sin omega t$ for each surface point.
$endgroup$
– Jennifer M.
Dec 17 '18 at 6:04
$begingroup$
If the domain of $sin omega t$ be $[0,π]$ then it is possible to get approximate solutions.
$endgroup$
– Awe Kumar Jha
Dec 17 '18 at 6:05
|
show 5 more comments
$begingroup$
Not that I know of. Are you sure you got your equation right? The $sin omega t$ term implies an oscillation due to elastic force, not rotation.
$endgroup$
– Andrei
Dec 17 '18 at 5:53
$begingroup$
$ sin omega t$ is a harmonic oscillation due to the object rotation. It will be moving according to the quadratic law, but will also be spinning.
$endgroup$
– Jennifer M.
Dec 17 '18 at 5:54
$begingroup$
I am just curious how you got that equation.
$endgroup$
– Andrei
Dec 17 '18 at 5:57
$begingroup$
Imagine a rotating object in a freefall. Its Z coordinate is decreasing as $-frac{1}{2} g t^2$ and its rotation would add $sin omega t$ for each surface point.
$endgroup$
– Jennifer M.
Dec 17 '18 at 6:04
$begingroup$
If the domain of $sin omega t$ be $[0,π]$ then it is possible to get approximate solutions.
$endgroup$
– Awe Kumar Jha
Dec 17 '18 at 6:05
$begingroup$
Not that I know of. Are you sure you got your equation right? The $sin omega t$ term implies an oscillation due to elastic force, not rotation.
$endgroup$
– Andrei
Dec 17 '18 at 5:53
$begingroup$
Not that I know of. Are you sure you got your equation right? The $sin omega t$ term implies an oscillation due to elastic force, not rotation.
$endgroup$
– Andrei
Dec 17 '18 at 5:53
$begingroup$
$ sin omega t$ is a harmonic oscillation due to the object rotation. It will be moving according to the quadratic law, but will also be spinning.
$endgroup$
– Jennifer M.
Dec 17 '18 at 5:54
$begingroup$
$ sin omega t$ is a harmonic oscillation due to the object rotation. It will be moving according to the quadratic law, but will also be spinning.
$endgroup$
– Jennifer M.
Dec 17 '18 at 5:54
$begingroup$
I am just curious how you got that equation.
$endgroup$
– Andrei
Dec 17 '18 at 5:57
$begingroup$
I am just curious how you got that equation.
$endgroup$
– Andrei
Dec 17 '18 at 5:57
$begingroup$
Imagine a rotating object in a freefall. Its Z coordinate is decreasing as $-frac{1}{2} g t^2$ and its rotation would add $sin omega t$ for each surface point.
$endgroup$
– Jennifer M.
Dec 17 '18 at 6:04
$begingroup$
Imagine a rotating object in a freefall. Its Z coordinate is decreasing as $-frac{1}{2} g t^2$ and its rotation would add $sin omega t$ for each surface point.
$endgroup$
– Jennifer M.
Dec 17 '18 at 6:04
$begingroup$
If the domain of $sin omega t$ be $[0,π]$ then it is possible to get approximate solutions.
$endgroup$
– Awe Kumar Jha
Dec 17 '18 at 6:05
$begingroup$
If the domain of $sin omega t$ be $[0,π]$ then it is possible to get approximate solutions.
$endgroup$
– Awe Kumar Jha
Dec 17 '18 at 6:05
|
show 5 more comments
2 Answers
2
active
oldest
votes
$begingroup$
Your equation is a transcendental equation. It can not be exactly solved but you can find the roots using numerical methods like Newton's method. However if you are still interested in having a closed form solution, you may use Bhaskara's approximation for domain $[-frac {π}{2},frac {π}{2}]$. For example, consider another related equation,
$$cos omega t + t^2 + bt - N = 0$$
From Bhaskara's approximation formula,
$$cos omega t ≈ frac {π^2 - 4omega^2 t^2}{π^2 + omega^2 t^2}$$
The equation is reduced to a general quartic,
$$t^4+bt^3+(T^2 - (4+N))t^2+T^2bt+T^2(1-N)=0$$
Where $T=frac {π}{omega}$This can be solved using Ferrari's method (https://proofwiki.org/wiki/Ferrari%27s_Method). Although the computational cost is high and the solution may turn really ugly, yet it will give you much required insight into the dependence of the time of flight of the object on various parameters.
$endgroup$
1
$begingroup$
I totally agree with your approach but $$cos(x) simeqfrac{pi ^2-4x^2}{pi ^2+x^2}qquad (-frac pi 2 leq xleqfrac pi 2)$$
$endgroup$
– Claude Leibovici
Dec 17 '18 at 9:02
$begingroup$
Now I have corrected it. Well , thanks for pointing my error.
$endgroup$
– Awe Kumar Jha
Dec 17 '18 at 10:53
$begingroup$
You are welcome ! Interesting problem from a numerical point of veiw. Cheers.
$endgroup$
– Claude Leibovici
Dec 17 '18 at 10:57
add a comment |
$begingroup$
In the same spirit as Awe Kumar Jha, considering $$sin (omega t) + at^2+bt=N$$ let
$$x=omega tqquad alpha=frac a{omega^2}qquadbeta=frac b{omega}$$ to make the equation
$$sin(x)+alpha x^2+beta x=N$$ and use Bhaskara I's sine approximation formula
$$sin(x) simeq frac{16 (pi -x) x}{5 pi ^2-4 (pi -x) x}qquad (0leq xleqpi)$$
$$4 alpha x^4+ 4( beta - pi alpha )x^3+ left(5 pi ^2 alpha -4 pi beta -4 N-16right)x^2+pi (5 pi beta +4
N+16)x-5 pi ^2 N=0tag 1$$ which, as said,can be solved using Ferrari's method (this would be quite tedious).
You can use another approximation (from far, not as good asBhaskara I's one ) minimizing with respect to $k$
$$int_0^pi left[sin(x)-k (pi-x)xright]^2,dx=frac{pi ^5 }{30}k^2-8 k+frac{pi }{2}$$ to get $color{blue}{k=frac{120}{pi ^5}}$ and then $sin(x)simeq k(pi-x)x$. This leads to
$$ (alpha -k)x^2+ (beta +pi k)x-N=0tag 2$$ which is simpler for sure. Using its solution as $x_0$ you can start Newton iterations to polish the root accorging to
$$x_{n+1}=x_n-frac{sin (x_n)+alpha x_n^2+beta x_n-N}{ cos (x_n)+2 alpha x_n+beta}$$
Let us try with $alpha=0.4$, $beta=0.6$, $N=2.345$. Solving $(1)$ gives $x^{(1)}_0 =1.26126$ while solving $(2)$ gives $x^{(2)}_0=1.27312$. Now, the iterates of Newton method starting with $x_0=x^{(2)}_0$
$$left(
begin{array}{cc}
n & x_n \
0 & 1.273117520 \
1 & 1.260970996 \
2 & 1.260976964
end{array}
right)$$ which is the solution for ten significant figures.
Edit
To stay with a quadratic equation, we could be even more empirical (not respecting the values at end points) and write
$$sin(x) simeq a_0+a_1x+a_2x^2$$ The same minimization method would lead to
$$a_0=frac{12 left(pi ^2-10right)}{pi ^3}qquad a_1=-frac{60 left(pi ^2-12right)}{pi ^4}qquad a_2=frac{60 left(pi ^2-12right)}{pi ^5}$$ which would lead to $x^{(3)}_0=1.26757$.
$endgroup$
add a comment |
Your Answer
StackExchange.ifUsing("editor", function () {
return StackExchange.using("mathjaxEditing", function () {
StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix) {
StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
});
});
}, "mathjax-editing");
StackExchange.ready(function() {
var channelOptions = {
tags: "".split(" "),
id: "69"
};
initTagRenderer("".split(" "), "".split(" "), channelOptions);
StackExchange.using("externalEditor", function() {
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled) {
StackExchange.using("snippets", function() {
createEditor();
});
}
else {
createEditor();
}
});
function createEditor() {
StackExchange.prepareEditor({
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: true,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: 10,
bindNavPrevention: true,
postfix: "",
imageUploader: {
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
},
noCode: true, onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
});
}
});
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3043577%2fcan-an-equation-of-the-form-sin-omega-t-at2bt-n-with-a-neq-0-be-sol%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
2 Answers
2
active
oldest
votes
2 Answers
2
active
oldest
votes
active
oldest
votes
active
oldest
votes
$begingroup$
Your equation is a transcendental equation. It can not be exactly solved but you can find the roots using numerical methods like Newton's method. However if you are still interested in having a closed form solution, you may use Bhaskara's approximation for domain $[-frac {π}{2},frac {π}{2}]$. For example, consider another related equation,
$$cos omega t + t^2 + bt - N = 0$$
From Bhaskara's approximation formula,
$$cos omega t ≈ frac {π^2 - 4omega^2 t^2}{π^2 + omega^2 t^2}$$
The equation is reduced to a general quartic,
$$t^4+bt^3+(T^2 - (4+N))t^2+T^2bt+T^2(1-N)=0$$
Where $T=frac {π}{omega}$This can be solved using Ferrari's method (https://proofwiki.org/wiki/Ferrari%27s_Method). Although the computational cost is high and the solution may turn really ugly, yet it will give you much required insight into the dependence of the time of flight of the object on various parameters.
$endgroup$
1
$begingroup$
I totally agree with your approach but $$cos(x) simeqfrac{pi ^2-4x^2}{pi ^2+x^2}qquad (-frac pi 2 leq xleqfrac pi 2)$$
$endgroup$
– Claude Leibovici
Dec 17 '18 at 9:02
$begingroup$
Now I have corrected it. Well , thanks for pointing my error.
$endgroup$
– Awe Kumar Jha
Dec 17 '18 at 10:53
$begingroup$
You are welcome ! Interesting problem from a numerical point of veiw. Cheers.
$endgroup$
– Claude Leibovici
Dec 17 '18 at 10:57
add a comment |
$begingroup$
Your equation is a transcendental equation. It can not be exactly solved but you can find the roots using numerical methods like Newton's method. However if you are still interested in having a closed form solution, you may use Bhaskara's approximation for domain $[-frac {π}{2},frac {π}{2}]$. For example, consider another related equation,
$$cos omega t + t^2 + bt - N = 0$$
From Bhaskara's approximation formula,
$$cos omega t ≈ frac {π^2 - 4omega^2 t^2}{π^2 + omega^2 t^2}$$
The equation is reduced to a general quartic,
$$t^4+bt^3+(T^2 - (4+N))t^2+T^2bt+T^2(1-N)=0$$
Where $T=frac {π}{omega}$This can be solved using Ferrari's method (https://proofwiki.org/wiki/Ferrari%27s_Method). Although the computational cost is high and the solution may turn really ugly, yet it will give you much required insight into the dependence of the time of flight of the object on various parameters.
$endgroup$
1
$begingroup$
I totally agree with your approach but $$cos(x) simeqfrac{pi ^2-4x^2}{pi ^2+x^2}qquad (-frac pi 2 leq xleqfrac pi 2)$$
$endgroup$
– Claude Leibovici
Dec 17 '18 at 9:02
$begingroup$
Now I have corrected it. Well , thanks for pointing my error.
$endgroup$
– Awe Kumar Jha
Dec 17 '18 at 10:53
$begingroup$
You are welcome ! Interesting problem from a numerical point of veiw. Cheers.
$endgroup$
– Claude Leibovici
Dec 17 '18 at 10:57
add a comment |
$begingroup$
Your equation is a transcendental equation. It can not be exactly solved but you can find the roots using numerical methods like Newton's method. However if you are still interested in having a closed form solution, you may use Bhaskara's approximation for domain $[-frac {π}{2},frac {π}{2}]$. For example, consider another related equation,
$$cos omega t + t^2 + bt - N = 0$$
From Bhaskara's approximation formula,
$$cos omega t ≈ frac {π^2 - 4omega^2 t^2}{π^2 + omega^2 t^2}$$
The equation is reduced to a general quartic,
$$t^4+bt^3+(T^2 - (4+N))t^2+T^2bt+T^2(1-N)=0$$
Where $T=frac {π}{omega}$This can be solved using Ferrari's method (https://proofwiki.org/wiki/Ferrari%27s_Method). Although the computational cost is high and the solution may turn really ugly, yet it will give you much required insight into the dependence of the time of flight of the object on various parameters.
$endgroup$
Your equation is a transcendental equation. It can not be exactly solved but you can find the roots using numerical methods like Newton's method. However if you are still interested in having a closed form solution, you may use Bhaskara's approximation for domain $[-frac {π}{2},frac {π}{2}]$. For example, consider another related equation,
$$cos omega t + t^2 + bt - N = 0$$
From Bhaskara's approximation formula,
$$cos omega t ≈ frac {π^2 - 4omega^2 t^2}{π^2 + omega^2 t^2}$$
The equation is reduced to a general quartic,
$$t^4+bt^3+(T^2 - (4+N))t^2+T^2bt+T^2(1-N)=0$$
Where $T=frac {π}{omega}$This can be solved using Ferrari's method (https://proofwiki.org/wiki/Ferrari%27s_Method). Although the computational cost is high and the solution may turn really ugly, yet it will give you much required insight into the dependence of the time of flight of the object on various parameters.
edited Dec 17 '18 at 15:56
answered Dec 17 '18 at 6:38
Awe Kumar JhaAwe Kumar Jha
43813
43813
1
$begingroup$
I totally agree with your approach but $$cos(x) simeqfrac{pi ^2-4x^2}{pi ^2+x^2}qquad (-frac pi 2 leq xleqfrac pi 2)$$
$endgroup$
– Claude Leibovici
Dec 17 '18 at 9:02
$begingroup$
Now I have corrected it. Well , thanks for pointing my error.
$endgroup$
– Awe Kumar Jha
Dec 17 '18 at 10:53
$begingroup$
You are welcome ! Interesting problem from a numerical point of veiw. Cheers.
$endgroup$
– Claude Leibovici
Dec 17 '18 at 10:57
add a comment |
1
$begingroup$
I totally agree with your approach but $$cos(x) simeqfrac{pi ^2-4x^2}{pi ^2+x^2}qquad (-frac pi 2 leq xleqfrac pi 2)$$
$endgroup$
– Claude Leibovici
Dec 17 '18 at 9:02
$begingroup$
Now I have corrected it. Well , thanks for pointing my error.
$endgroup$
– Awe Kumar Jha
Dec 17 '18 at 10:53
$begingroup$
You are welcome ! Interesting problem from a numerical point of veiw. Cheers.
$endgroup$
– Claude Leibovici
Dec 17 '18 at 10:57
1
1
$begingroup$
I totally agree with your approach but $$cos(x) simeqfrac{pi ^2-4x^2}{pi ^2+x^2}qquad (-frac pi 2 leq xleqfrac pi 2)$$
$endgroup$
– Claude Leibovici
Dec 17 '18 at 9:02
$begingroup$
I totally agree with your approach but $$cos(x) simeqfrac{pi ^2-4x^2}{pi ^2+x^2}qquad (-frac pi 2 leq xleqfrac pi 2)$$
$endgroup$
– Claude Leibovici
Dec 17 '18 at 9:02
$begingroup$
Now I have corrected it. Well , thanks for pointing my error.
$endgroup$
– Awe Kumar Jha
Dec 17 '18 at 10:53
$begingroup$
Now I have corrected it. Well , thanks for pointing my error.
$endgroup$
– Awe Kumar Jha
Dec 17 '18 at 10:53
$begingroup$
You are welcome ! Interesting problem from a numerical point of veiw. Cheers.
$endgroup$
– Claude Leibovici
Dec 17 '18 at 10:57
$begingroup$
You are welcome ! Interesting problem from a numerical point of veiw. Cheers.
$endgroup$
– Claude Leibovici
Dec 17 '18 at 10:57
add a comment |
$begingroup$
In the same spirit as Awe Kumar Jha, considering $$sin (omega t) + at^2+bt=N$$ let
$$x=omega tqquad alpha=frac a{omega^2}qquadbeta=frac b{omega}$$ to make the equation
$$sin(x)+alpha x^2+beta x=N$$ and use Bhaskara I's sine approximation formula
$$sin(x) simeq frac{16 (pi -x) x}{5 pi ^2-4 (pi -x) x}qquad (0leq xleqpi)$$
$$4 alpha x^4+ 4( beta - pi alpha )x^3+ left(5 pi ^2 alpha -4 pi beta -4 N-16right)x^2+pi (5 pi beta +4
N+16)x-5 pi ^2 N=0tag 1$$ which, as said,can be solved using Ferrari's method (this would be quite tedious).
You can use another approximation (from far, not as good asBhaskara I's one ) minimizing with respect to $k$
$$int_0^pi left[sin(x)-k (pi-x)xright]^2,dx=frac{pi ^5 }{30}k^2-8 k+frac{pi }{2}$$ to get $color{blue}{k=frac{120}{pi ^5}}$ and then $sin(x)simeq k(pi-x)x$. This leads to
$$ (alpha -k)x^2+ (beta +pi k)x-N=0tag 2$$ which is simpler for sure. Using its solution as $x_0$ you can start Newton iterations to polish the root accorging to
$$x_{n+1}=x_n-frac{sin (x_n)+alpha x_n^2+beta x_n-N}{ cos (x_n)+2 alpha x_n+beta}$$
Let us try with $alpha=0.4$, $beta=0.6$, $N=2.345$. Solving $(1)$ gives $x^{(1)}_0 =1.26126$ while solving $(2)$ gives $x^{(2)}_0=1.27312$. Now, the iterates of Newton method starting with $x_0=x^{(2)}_0$
$$left(
begin{array}{cc}
n & x_n \
0 & 1.273117520 \
1 & 1.260970996 \
2 & 1.260976964
end{array}
right)$$ which is the solution for ten significant figures.
Edit
To stay with a quadratic equation, we could be even more empirical (not respecting the values at end points) and write
$$sin(x) simeq a_0+a_1x+a_2x^2$$ The same minimization method would lead to
$$a_0=frac{12 left(pi ^2-10right)}{pi ^3}qquad a_1=-frac{60 left(pi ^2-12right)}{pi ^4}qquad a_2=frac{60 left(pi ^2-12right)}{pi ^5}$$ which would lead to $x^{(3)}_0=1.26757$.
$endgroup$
add a comment |
$begingroup$
In the same spirit as Awe Kumar Jha, considering $$sin (omega t) + at^2+bt=N$$ let
$$x=omega tqquad alpha=frac a{omega^2}qquadbeta=frac b{omega}$$ to make the equation
$$sin(x)+alpha x^2+beta x=N$$ and use Bhaskara I's sine approximation formula
$$sin(x) simeq frac{16 (pi -x) x}{5 pi ^2-4 (pi -x) x}qquad (0leq xleqpi)$$
$$4 alpha x^4+ 4( beta - pi alpha )x^3+ left(5 pi ^2 alpha -4 pi beta -4 N-16right)x^2+pi (5 pi beta +4
N+16)x-5 pi ^2 N=0tag 1$$ which, as said,can be solved using Ferrari's method (this would be quite tedious).
You can use another approximation (from far, not as good asBhaskara I's one ) minimizing with respect to $k$
$$int_0^pi left[sin(x)-k (pi-x)xright]^2,dx=frac{pi ^5 }{30}k^2-8 k+frac{pi }{2}$$ to get $color{blue}{k=frac{120}{pi ^5}}$ and then $sin(x)simeq k(pi-x)x$. This leads to
$$ (alpha -k)x^2+ (beta +pi k)x-N=0tag 2$$ which is simpler for sure. Using its solution as $x_0$ you can start Newton iterations to polish the root accorging to
$$x_{n+1}=x_n-frac{sin (x_n)+alpha x_n^2+beta x_n-N}{ cos (x_n)+2 alpha x_n+beta}$$
Let us try with $alpha=0.4$, $beta=0.6$, $N=2.345$. Solving $(1)$ gives $x^{(1)}_0 =1.26126$ while solving $(2)$ gives $x^{(2)}_0=1.27312$. Now, the iterates of Newton method starting with $x_0=x^{(2)}_0$
$$left(
begin{array}{cc}
n & x_n \
0 & 1.273117520 \
1 & 1.260970996 \
2 & 1.260976964
end{array}
right)$$ which is the solution for ten significant figures.
Edit
To stay with a quadratic equation, we could be even more empirical (not respecting the values at end points) and write
$$sin(x) simeq a_0+a_1x+a_2x^2$$ The same minimization method would lead to
$$a_0=frac{12 left(pi ^2-10right)}{pi ^3}qquad a_1=-frac{60 left(pi ^2-12right)}{pi ^4}qquad a_2=frac{60 left(pi ^2-12right)}{pi ^5}$$ which would lead to $x^{(3)}_0=1.26757$.
$endgroup$
add a comment |
$begingroup$
In the same spirit as Awe Kumar Jha, considering $$sin (omega t) + at^2+bt=N$$ let
$$x=omega tqquad alpha=frac a{omega^2}qquadbeta=frac b{omega}$$ to make the equation
$$sin(x)+alpha x^2+beta x=N$$ and use Bhaskara I's sine approximation formula
$$sin(x) simeq frac{16 (pi -x) x}{5 pi ^2-4 (pi -x) x}qquad (0leq xleqpi)$$
$$4 alpha x^4+ 4( beta - pi alpha )x^3+ left(5 pi ^2 alpha -4 pi beta -4 N-16right)x^2+pi (5 pi beta +4
N+16)x-5 pi ^2 N=0tag 1$$ which, as said,can be solved using Ferrari's method (this would be quite tedious).
You can use another approximation (from far, not as good asBhaskara I's one ) minimizing with respect to $k$
$$int_0^pi left[sin(x)-k (pi-x)xright]^2,dx=frac{pi ^5 }{30}k^2-8 k+frac{pi }{2}$$ to get $color{blue}{k=frac{120}{pi ^5}}$ and then $sin(x)simeq k(pi-x)x$. This leads to
$$ (alpha -k)x^2+ (beta +pi k)x-N=0tag 2$$ which is simpler for sure. Using its solution as $x_0$ you can start Newton iterations to polish the root accorging to
$$x_{n+1}=x_n-frac{sin (x_n)+alpha x_n^2+beta x_n-N}{ cos (x_n)+2 alpha x_n+beta}$$
Let us try with $alpha=0.4$, $beta=0.6$, $N=2.345$. Solving $(1)$ gives $x^{(1)}_0 =1.26126$ while solving $(2)$ gives $x^{(2)}_0=1.27312$. Now, the iterates of Newton method starting with $x_0=x^{(2)}_0$
$$left(
begin{array}{cc}
n & x_n \
0 & 1.273117520 \
1 & 1.260970996 \
2 & 1.260976964
end{array}
right)$$ which is the solution for ten significant figures.
Edit
To stay with a quadratic equation, we could be even more empirical (not respecting the values at end points) and write
$$sin(x) simeq a_0+a_1x+a_2x^2$$ The same minimization method would lead to
$$a_0=frac{12 left(pi ^2-10right)}{pi ^3}qquad a_1=-frac{60 left(pi ^2-12right)}{pi ^4}qquad a_2=frac{60 left(pi ^2-12right)}{pi ^5}$$ which would lead to $x^{(3)}_0=1.26757$.
$endgroup$
In the same spirit as Awe Kumar Jha, considering $$sin (omega t) + at^2+bt=N$$ let
$$x=omega tqquad alpha=frac a{omega^2}qquadbeta=frac b{omega}$$ to make the equation
$$sin(x)+alpha x^2+beta x=N$$ and use Bhaskara I's sine approximation formula
$$sin(x) simeq frac{16 (pi -x) x}{5 pi ^2-4 (pi -x) x}qquad (0leq xleqpi)$$
$$4 alpha x^4+ 4( beta - pi alpha )x^3+ left(5 pi ^2 alpha -4 pi beta -4 N-16right)x^2+pi (5 pi beta +4
N+16)x-5 pi ^2 N=0tag 1$$ which, as said,can be solved using Ferrari's method (this would be quite tedious).
You can use another approximation (from far, not as good asBhaskara I's one ) minimizing with respect to $k$
$$int_0^pi left[sin(x)-k (pi-x)xright]^2,dx=frac{pi ^5 }{30}k^2-8 k+frac{pi }{2}$$ to get $color{blue}{k=frac{120}{pi ^5}}$ and then $sin(x)simeq k(pi-x)x$. This leads to
$$ (alpha -k)x^2+ (beta +pi k)x-N=0tag 2$$ which is simpler for sure. Using its solution as $x_0$ you can start Newton iterations to polish the root accorging to
$$x_{n+1}=x_n-frac{sin (x_n)+alpha x_n^2+beta x_n-N}{ cos (x_n)+2 alpha x_n+beta}$$
Let us try with $alpha=0.4$, $beta=0.6$, $N=2.345$. Solving $(1)$ gives $x^{(1)}_0 =1.26126$ while solving $(2)$ gives $x^{(2)}_0=1.27312$. Now, the iterates of Newton method starting with $x_0=x^{(2)}_0$
$$left(
begin{array}{cc}
n & x_n \
0 & 1.273117520 \
1 & 1.260970996 \
2 & 1.260976964
end{array}
right)$$ which is the solution for ten significant figures.
Edit
To stay with a quadratic equation, we could be even more empirical (not respecting the values at end points) and write
$$sin(x) simeq a_0+a_1x+a_2x^2$$ The same minimization method would lead to
$$a_0=frac{12 left(pi ^2-10right)}{pi ^3}qquad a_1=-frac{60 left(pi ^2-12right)}{pi ^4}qquad a_2=frac{60 left(pi ^2-12right)}{pi ^5}$$ which would lead to $x^{(3)}_0=1.26757$.
edited Dec 18 '18 at 3:52
answered Dec 17 '18 at 10:42
Claude LeiboviciClaude Leibovici
122k1157134
122k1157134
add a comment |
add a comment |
Thanks for contributing an answer to Mathematics Stack Exchange!
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
Use MathJax to format equations. MathJax reference.
To learn more, see our tips on writing great answers.
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3043577%2fcan-an-equation-of-the-form-sin-omega-t-at2bt-n-with-a-neq-0-be-sol%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
$begingroup$
Not that I know of. Are you sure you got your equation right? The $sin omega t$ term implies an oscillation due to elastic force, not rotation.
$endgroup$
– Andrei
Dec 17 '18 at 5:53
$begingroup$
$ sin omega t$ is a harmonic oscillation due to the object rotation. It will be moving according to the quadratic law, but will also be spinning.
$endgroup$
– Jennifer M.
Dec 17 '18 at 5:54
$begingroup$
I am just curious how you got that equation.
$endgroup$
– Andrei
Dec 17 '18 at 5:57
$begingroup$
Imagine a rotating object in a freefall. Its Z coordinate is decreasing as $-frac{1}{2} g t^2$ and its rotation would add $sin omega t$ for each surface point.
$endgroup$
– Jennifer M.
Dec 17 '18 at 6:04
$begingroup$
If the domain of $sin omega t$ be $[0,π]$ then it is possible to get approximate solutions.
$endgroup$
– Awe Kumar Jha
Dec 17 '18 at 6:05