描述下问题吧:
Octave解多元一次方程(不是方程组),希望得出某一个解。
比如我打算解 x + y = 100,希望得到整数结果,比如(0,100)或者(50,50),
有一个整数解就好。我玩儿了一下,代码如下:
function y = f (x)
y(1) = x(1) + x(2) - 100;
end function
[x, info] = fsolve ("f", [1; 2])
现在得到的结果是:
49.500
50.500
看,都是double类型的解,请问可以有方法让它是整数解吗?
其实fsolve ("f", [X; X])中两个XX是初始值,初始值设的好,就能解出整数解,例如我这个问题,初始值设为20,80,那直接就解出来了。
这也是Octave解方程的局限所在,必须要有初始值。自结了。。。