hệ phương trình phi tuyến với matlab

  • Thread starter lab.cdt
  • Ngày mở chủ đề
L

lab.cdt

Author
Tớ vừa mới khảo sát 1 hệ thống động học và rút ra hệ phương trình như sau:
(2*h*y - 2*x*y)*a + (2*x^2 - h*x)*b - h*a*b = 0 (1)
[(b*x)^2 + a^2*(x - h)^2]*(b*x - a*y)^2 - [(2*x - a - h)^2 + (2*y - b)^2]*[b*x^2 - a*(x - h)*y - b*h*x]^2 (2)
Tớ thao tác trên matlab và cho kết quả như sau:
>> syms a b x y h
>> [a, b] = solve((2*h*y - 2*x*y)*a + (2*x^2 - h*x)*b - h*a*b ,((b*x)^2 + a^2*(x - h)^2)*(b*x - a*y)^2 - ((2*x - a - h)^2 + (2*y - b)^2)*(b*x^2 - a*(x - h)*y - b*h*x)^2)
a =

a
h/2
a + h/2
(a^3*h - a^2*h^2 - a*b^2*h + (a*b*h*(- a^2 + a*b*(2*i) + h*a + b^2 + h*b*i))/(a*i + b))/(a^3 - h*a^2 + a*b^2 + h*b^2)
(a^3*h - a^2*h^2 - a*b^2*h + (a*b*h*(a^2 + a*b*(2*i) - h*a - b^2 + h*b*i))/(a*i - b))/(a^3 - h*a^2 + a*b^2 + h*b^2)


b =

b
b
(2*a*b)/(2*a - h)
(- a^2 + a*b*(2*i) + h*a + b^2 + h*b*i)/(2*(a*i + b))
(a^2 + a*b*(2*i) - h*a - b^2 + h*b*i)/(2*(a*i - b))
Tớ có vài thắc mắc:
1.Tại sao 1 nghiệm của hệ lại chính là a,b????
2.Các tham so x, y đâu hết rồi??
3.Theo mình nghĩ là nó có 3 nghiệm thực??
Mong các Pro matlab giải thich hộ tớ........
 
Top