I am facing a problem in solving a system of three equations for three unknown parameters. In order to have a clear understanding of the issue, I have given a detailed description of the origin of the said equations. The equations arise from a differential equation and its boundary conditions from one of the research articles as:
$\omega^2 L C \frac{d^2v(x)}{dx^2} + v(x) - V_{cc}+ V_R \sin(x+\phi) = 0$
Whose General solution is given in the normalized form as
$\frac{v(x)}{V_{cc}}=C_1\cos qx+C_2\sin qx+1+ \frac{q^2}{1-q^2}.\frac{V_R}{V_{cc}}\sin (x+\phi)$
Where
$q=\frac{1}{\omega\sqrt{LC}}$
and
$\omega C\frac{d v(x)}{d x} =i_{c}(x)$
And the boundary conditions are:
$i_{c}(\pi)=\frac{\pi V_{cc}-2V_R \cos\phi}{\omega L} $
and
$v(\pi)=0$
So once I apply boundary conditions I have these two following equations
$0=C_1\cos q\pi+C_2\sin q\pi+1- \frac{q^2}{1-q^2}.\frac{V_R}{V_{cc}}\sin (\phi)\\ i_{c}(\pi)=\frac{\pi V_{cc}-2V_R \cos\phi}{\omega L} = \omega C \begin{bmatrix}-C_1q \sin(q\pi)+C_2q \cos(q\pi)+\frac{q^2}{1-q^2}V_R\cos(\phi) \end{bmatrix}_{x = \pi}$
After substituting the boundary conditions I am solving the above two equations for C1 and C2 constants using "Solve " in Wolfram Mathematica. and I obtain the following
$C_1= -\cos(q\pi)-q \pi \sin(q\pi)+(\frac{q}{1-q^2})\frac{V_R}{V_{cc}} \times\begin{bmatrix} q \cos(q\pi) \sin \phi +(1-2q^2) \sin(q\pi) \cos \phi \end{bmatrix}$
$C_2= -\sin(q\pi)+q \pi \cos(q\pi)+(\frac{q}{1-q^2})\frac{V_R}{V_{cc}} \times \begin{bmatrix} q \sin(q\pi) \sin \phi -(1-2q^2) \cos(q\pi) \cos \phi \end{bmatrix}$
The value of C1 and C2 is back substituted in v(x) of the general solution of the differential equation obtained above. Then I have set of three equations as follows
$V_R= -\frac{1}{\pi} \int_{0}^{2\pi} v(x)\sin(x+\phi)dx$
$\begin{bmatrix}v(x)\end{bmatrix}_{x=2\pi}=0$ $$\begin{bmatrix}\frac{d v(x)}{d x}\end{bmatrix}_{x=2\pi}=0$$
I am trying to solve the above three equations for three unknown parameters $$\phi,q, V_R$$ whereas Vcc is a constant in my equations. The three unknown parameters are reported by the authors of the research article as
$\phi = -41.614\\ q=1.607\\ V_R= 0.9253 \times V_{cc}$
I am solving the above set of three equations using Solve in Wolfram Mathematica but I get the error message as Solve::nsmet: This system cannot be solved with the methods available to Solve. >>. The complete notebook code that I use in Mathematica to solve it is given below
Clear[x, y, A, q, b, c, ysol, ysolsimp]
A = \[Omega]^2 *L*C
b = Vcc
c = VR
eqn = A* y''[x] + y[x] - b + c *Sin[x + phi] == 0
Solut = DSolve[eqn, y, x]
ysol = y /. First@DSolve[eqn, y, x]
ysolsimp = ysol[x] /. {C[1] -> c1, C[2] -> c2 , (x/(Sqrt[C]* Sqrt[L]* \[Omega])) -> q*x } //
FullSimplify
diff = D[ysolsimp, x]
ic = \[Omega]*C*diff
icp = \[Omega]*C*diff /. x -> Pi
icL = (Pi*Vcc - 2*VR*Cos[phi])/(\[Omega]*L)
icp == icL
eq1 = (icp == icL)
eq2 = (ysol[x] /. {C[1] -> c1, C[2] -> c2, x -> \[Pi]} /. (\[Pi]/( Sqrt[C]* Sqrt[L]* \[Omega])) -> q*\[Pi] // FullSimplify) == 0
soleq1 = Solve[ eq1 && eq2, {c1, c2}] /. { (\[Pi]/(Sqrt[C]* Sqrt[L]* \[Omega])) -> q*\[Pi] } // FullSimplify
c1 = c1 /. soleq1
c2 = c2 /. soleq1
inty = \!\( \*SubsuperscriptBox[\(\[Integral]\), \(0\), \(2 \[Pi]\)]\(ysolsimp*\ Sin[x + phi] \DifferentialD]x\)\) // FullSimplify
fourintegral = (-1/\[Pi])*\!\( \*SubsuperscriptBox[\(\[Integral]\), \(0\), \(2 \[Pi]\)]\(ysolsimp*\ Sin[x + phi] \[DifferentialD]x\)\) // FullSimplify
eq4 = (fourintegral == VR)
eq5 = ysol[x] /. {C[1] -> c1, C[2] -> c2 , (x/(Sqrt[C]* Sqrt[L]* \[Omega])) -> q*x , x -> 2 \[Pi]} // FullSimplify
eq5b = (eq5 == 0)
eq6 = D[ysolsimp, x]
eq6b = eq6 /. x -> 2*\[Pi]
eq6c = eq6b == 0
Solve[{eq4 , eq5b , eq6c}, {phi, q, VR}]
Solve::nsmet: This system cannot be solved with the methods available to Solve. >>
The complete Wolfram Mathematica Notebook is also attached to this post. Kindly help me to debug the error. How can I get the same answer as that of the research article? Thanking you.
Attachments: