비선형 solve 방법도 알려줬는데,
#코드 일부 nonlinsolution = sp.nonlinsolve((real_part,imaginary_part),(xx,rr)) print("Non-linear solutions:") for solution in nonlinsolution: solution_xx, solution_rr = solution print(f" solution of xx: {solution_xx}, \n solution of rr: {solution_rr}") print('-' * 50) print(f'Time taken to find non-linear solution: {end_time - start_time} seconds') #결과 Non-linear solutions: solution of xx: 0, solution of rr: 0 solution of xx: es*(-er*ps*sin(delta) + er*qs*cos(delta) - es*qs)*(-er*ps*cos(delta) - er*qs*sin(delta) + es*ps)/((ps**2 + qs**2)*(er*ps*cos(delta) + er*qs*sin(delta) - es*ps)), solution of rr: es*(-er*ps*cos(delta) - er*qs*sin(delta) + es*ps)/(ps**2 + qs**2) -------------------------------------------------- Time taken to find non-linear solution: 2.574453115463257 seconds
이번엔 결과가 2세트로 나왔는데
xx=0 and rr=0 한세트와
xx=es*(-er*ps*sin(δ)+er*qs*cos(δ)-es*qs)*(-er*ps*cos(δ)-er*qs*sin(δ)+es*ps)(ps2+qs2)*(er*ps*cos(δ)+er*qs*sin(δ)-es*ps) xx = \frac{es*(-er*ps*sin(\delta) + er*qs*cos(\delta) - es*qs)*(-er*ps*cos(\delta) - er*qs*sin(\delta) + es*ps)}{(ps^2 + qs^2)*(er*ps*cos(\delta) + er*qs*sin(\delta) - es*ps)}
rr=es*(-er*ps*cos(δ)-er*qs*sin(δ)+es*ps)(ps2+qs2) rr = \frac{es*(-er*ps*cos(\delta) - er*qs*sin(\delta) + es*ps)}{(ps^2 + qs^2)}
아 분모 분자 뒷부분 약분되고 -1 남으니, solve 와 같네요. 그러면 여기서는 분모≠0 이 한번에 결과에 표현됐다고 봐야할 듯.
여기에 파일을 끌어 놓거나 왼쪽의 버튼을 클릭하세요.
파일 용량 제한 : 0MB (허용 확장자 : *.*)
업로드 중... (0%)
비선형 solve 방법도 알려줬는데,
이번엔 결과가 2세트로 나왔는데
xx=0 and rr=0 한세트와
답이 다르네요? 다른 건가? 한번 계산을 해봐야 할 듯?
아 분모 분자 뒷부분 약분되고 -1 남으니,
solve 와 같네요.
그러면 여기서는 분모≠0 이 한번에 결과에 표현됐다고 봐야할 듯.