24小时热门版块排行榜

>论坛更新日志 (3704)
>导师招生 (304)
>休闲灌水 (163)
>招聘信息布告栏 (135)
>论文道贺祈福 (110)
>硕博家园 (87)
>文献求助 (62)
>虫友互识 (55)
>找工作 (43)
>考研 (41)
>生物科学 (29)
>博后之家 (28)
>考博 (23)
>基金申请 (20)
>医学 (18)
>教师之家 (14)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

gzy_103

金虫 (正式写手)

应助: 7 (幼儿园)
金币: 926.9
散金: 65
红花: 1
帖子: 343
在线: 233.4小时
虫号: 2562392
注册: 2013-07-24
性别: GG
专业: 能源化工

[求助] 求解方程组已有2人参与

function y=broyden(x0)
a=eye(length(x0));
x1=x0-myfun2(x0)/a;
n=1;
while(norm(x1-x0)>=1.0e-6)&(n<=100000000)
x0=x1;
x1=x0-myfun2(x0)/a;
p=x1-x0;
q=myfun2(x1)-myfun(x0);
a=a+(q-p*a)'*p/norm(p);
n=n+1;
end
y=x1;
n
function y= myfun2(x)
a=2.6;
y=['x(1)*0.4-0.21+x(2)*0.6','x(2)-x(1)^(1/(1-a))*(2/3)*((1-x(1))^(a/(a-1))*(0.21/0.79)^(a/(a-1))-x(1)^(1/(1-a)))'];

这是代码，总是报错，帮忙看一下吧

回复此楼

» 猜你喜欢

1楼 2014-05-30 20:35:29

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

gzy_103

金虫 (正式写手)

应助: 7 (幼儿园)
金币: 926.9
散金: 65
红花: 1
帖子: 343
在线: 233.4小时
虫号: 2562392
注册: 2013-07-24
性别: GG
专业: 能源化工

送红花一朵

引用回帖:

5楼: Originally posted by change0618 at 2014-06-02 14:30:13
从这个function y= myfun2(x)函数可以看出x(1)不能大于1也不能小于0，大于1则((1-x(1))^(a/(a-1))为负数，小于0则x(1)^(1/(1-a))为负数，但在while循环中x1=x0-myfun2(x0)/a难免出现这种情况，所以在q=myfun2(x1)-m ...

您看这个方程该怎样解呢？