Solve simple sinusoidal equation

Good day!
Please help me to solve a system of equations:





A simple example. However, I need an algorithmic solution that is suitable for any three points on the graph.
Our desired sine describes three unknown parameters: amplitude, frequency and phase. A system of equations can be solved, the number of unknowns is not less than the number of equations. The conclusion is that we can find the 3 parameters by 3 any points. However, I was not able to write code to do this. Can you help me?
Thank you in advance

- - - Updated - - -

This reminds me of the problem where you find the unique circle which contains all 3 points.

If you can transfer the 3 points on the sinewave, geometrically, to become 3 points on an x-y grid...

Then you might approach it as finding a circle containing all 3 points...

which would give you the sine formula, etc.

Is the variable "frequency" an integer? If so, the first and third equations are not independent; that is, sin(phase)=0 implies that sin(phase+(freq*pi))=0 for all integer values of 'freq', so the system is underdetermined. If you were assuming integer frequency, this system of equations will always be underdetermined.

Also: keep in mind that the periodicity of sine and cosine means that for non-integer frequency, there are an infinite number of solutions corresponding to alias frequencies. So the key to solving this is:

(a) stick to non-integer frequencies if applicable;
(b) decide what half-unit interval of frequencies applies.

In practical terms, (b) means describing the solution with an arc-sine function whose range you know, and offsetting/reflecting it to produce a solution in your acceptable interval. Since each half of the sine cycle is symmetrical about an odd multiple of pi/4, the alias solutions should be arranged in pairs of frequencies spaced <= 0.5 apart, separated by empty intervals of >= 0.5.

Could you describe the original problem you're trying to solve? I might be able to fully constrain the desirable solutions with more information.

If we have the general system of three equations:

A*sin(θ)=x
A*sin(θ+ω/M)=y
A*sin(θ+ω/N)=z

where x,y,z,M and N are known, since sin(a+b)=sin(a)*cos(b)+cos(a)*sin(b), then:

A*sin(θ)=x
A*sin(θ)*cos(ω/M)+A*cos(θ)*sin(ω/M)=y
A*sin(θ)*cos(ω/N)+A*cos(θ)*sin(ω/N)=z

substituting now the first one in the other two:

x*cos(ω/M)+sqrt(1-x²)*sin(ω/M)=y
x*cos(ω/N)+sqrt(1-x²)*sin(ω/N)=z

that are of the form a*sin(α)+b*cos(α) that is rewritable as c*sin(α+φ), so that:

a=c*cos(φ)
b=c*sin(φ)

then:

φ=arctg(b/a) and c=a/cos(φ)

applying this to the first of our two last equations:

φ=arctg(x/sqrt(1-x²)) and c=sqrt(1-x²)/cos(φ)

-----------------------
| ω=M*[arcsin(y/c)-φ] |
-----------------------

dividing now the two equation one each other:

A*sin(θ)*cos(ω/M)+A*cos(θ)*sin(ω/M).....y
---------------------------------------- = --
A*sin(θ)*cos(ω/N)+A*cos(θ)*sin(ω/N).....z

from which:

z*sin(θ)*cos(ω/M)+z*cos(θ)*sin(ω/M) = y*sin(θ)*cos(ω/N)+y*cos(θ)*sin(ω/N)

dividing by cos(θ) and rearranging:

tg(θ)*[z*cos(ω/M)-y*cos(ω/N)]=y*sin(ω/M)-z*sin(ω/N)

then:

----------------------------------------------------------------
| θ = arctg{[y*sin(ω/M)-z*sin(ω/N)]/[z*cos(ω/M)-y*cos(ω/N)]} |
----------------------------------------------------------------

and

A=x/sin(θ)

Please check the calculations are correct since I didn't verified them