最速下降法解方程组matlab编程,用matlab-最速下降法-解方程组,四个未知数：a,b,c,d,我弄了好久没弄出来....

问题描述：

用matlab-最速下降法-解方程组,四个未知数：a,b,c,d,我弄了好久没弄出来.

方程组：

[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文件或程序.

1个回答

分类：

综合

2014-09-23

问题解答：

我来补答

你得给一个初值,初值不合适的话就解不出来.

程序如下,修改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，3

展开全文阅读