作业帮 > 综合 > 作业

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

来源:学生作业帮 编辑:搜狗做题网作业帮 分类:综合作业 时间:2024/05/25 12:30:32
用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文件或程序.
用matlab-最速下降法-解方程组,四个未知数:a,b,c,d,我弄了好久没弄出来.
你得给一个初值,初值不合适的话就解不出来.
程序如下,修改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