Question

求助！用matlab-最速下降法-解方程组，四个未知数：a,b,c,d,感谢。 我弄了好久没弄出来。悬赏20金~！

方程组：
[77.6*(1/a+(1-exp(-0.5*b/d))/b+0.5/c)=8;
77.6*(1/a+(1-exp(-b/d))/b+1/c)=9;
77.6*(1/a+(1-exp(-1.75*b/d))/b+1.75/c)=10;
77.6*(1/a+(1-exp(-4.25*b/d))/b+4.25/c)=11];
最好有matlab的M文件或程序。谢啦

Answer 1

你得给一个初值，初值不合适的话就解不出来。
程序如下，修改x0就是修改初值了。运行会有个警告，表示迭代次数达到最大。
% 程序开始 test_fsolve.m
function test_fsolve
x0 = [1; 1; 1; 1]*1e-3;
xx = fsolve(@eqns, x0)

function y = eqns(x)
[a,b,c,d]=deal(x(1),x(2),x(3),x(4));
y = zeros(4,1);
y(1) = 77.6*(1/a+(1-exp(-0.5*b/d))/b+0.5/c)-8;
y(2) = 77.6*(1/a+(1-exp(-b/d))/b+1/c)-9;
y(3) = 77.6*(1/a+(1-exp(-1.75*b/d))/b+1.75/c)-10;
y(4) = 77.6*(1/a+(1-exp(-4.25*b/d))/b+4.25/c)-11;
% 程序结束 test_fsolve.m

追问

大神，你这个方法不是matlab-最速下降法。
用你这个方法，求不出来的。
请求帮忙用最速下降法求解下。感谢。
我会追加财富值的！

追答

1. 你怎么知道fsolve里面用的不是最速下降？即便不是，为什么肯定其比最速下降差？ 2. 你用最速下降法做过推导呢？还是用其他语言写的算法仿真过？那你用的迭代起点是什么？

追问

那你用fsolve求出来的结果是什么呢？给出来看看！

追答

1. “你得给一个初值，初值不合适的话就解不出来。”
   

2. 你不会自己运行下？我的意思很明白，我找不到合适的初值，所以解不出来。

追问

逻辑性很强，1,2，3