lab.cdt
18-05-2012, 02:30 AM
Tớ có vấn đề như sau: tớ muốn giải một hệ phương trình phi tuyến trong matlab tùy thuộc các hệ số nhập vào. Tớ dùng lệnh fsolve để giải quyết với mong muốn là chỉ nhận nghiệm thực.Code trong m-file của tớ như sau:
function x = T(a,b,c,d,g,h)
a = input('a:');
b = input('b:');
c = input('c:');
d = input('d:');
g = input('g:');
h = input('h:');
function F = W(m,n,k,p,q,r)
syms x1y1
F = [m*x1+n*y1+k*x1*y1;(x1-p)^2+(y1-q)^2-r];
end
x0 = [1;1];
options = optimset('Display','off');
x = fsolve('W(a,b,c,d,g,h)',x0,options);
end
Khi gọi hàm T thì chương trình báo lỗi ở hàm con W. Tớ biết là khi dùng fsolve để giải thì hàm con W có thông số nhập vào là một vecto tức là khai báo W(x). Nhưng nếu vậy thì làm sao có thể truyền giá trị của các biến ở hàm chính T vô được.
Nhờ các cậu giải thích hộ tớ......
function x = T(a,b,c,d,g,h)
a = input('a:');
b = input('b:');
c = input('c:');
d = input('d:');
g = input('g:');
h = input('h:');
function F = W(m,n,k,p,q,r)
syms x1y1
F = [m*x1+n*y1+k*x1*y1;(x1-p)^2+(y1-q)^2-r];
end
x0 = [1;1];
options = optimset('Display','off');
x = fsolve('W(a,b,c,d,g,h)',x0,options);
end
Khi gọi hàm T thì chương trình báo lỗi ở hàm con W. Tớ biết là khi dùng fsolve để giải thì hàm con W có thông số nhập vào là một vecto tức là khai báo W(x). Nhưng nếu vậy thì làm sao có thể truyền giá trị của các biến ở hàm chính T vô được.
Nhờ các cậu giải thích hộ tớ......