24小时热门版块排行榜    

查看: 903  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

gzy_103

金虫 (正式写手)

[求助] 求解方程组 已有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)))'];

这是代码,总是报错,帮忙看一下吧
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

charmantang

禁虫 (正式写手)


感谢参与,应助指数 +1
jjdg: 金币+1, 感谢参与 2014-06-04 00:47:27
本帖内容被屏蔽

4楼2014-06-01 18:18:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

仙木映月

金虫 (正式写手)

myfun2(x)有点奇怪啊,都放在引号里能运算么?
Hardtosay.
2楼2014-05-30 21:19:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gzy_103

金虫 (正式写手)

引用回帖:
2楼: Originally posted by 仙木映月 at 2014-05-30 21:19:01
myfun2(x)有点奇怪啊,都放在引号里能运算么?

你能帮我改一下吗?
3楼2014-05-30 21:29:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师


jjdg: 金币+1, 感谢参与 2014-06-04 00:47: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)-myfun(x0)时出现错误。

» 本帖已获得的红花(最新10朵)

5楼2014-06-02 14:30:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 286分人工智能专业请求调剂愿意跨考! +3 lemonzzn 2026-03-17 3/150 2026-03-19 23:58 by 23Postgrad
[考研] 329求调剂 +5 想上学吖吖 2026-03-19 5/250 2026-03-19 23:56 by 23Postgrad
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +11 yangfz 2026-03-17 11/550 2026-03-19 15:06 by houyaoxu
[考研] 材料与化工求调剂 +7 为学666 2026-03-16 7/350 2026-03-19 14:48 by 尽舜尧1
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[考研] 332求调剂 +3 ydfyh 2026-03-17 3/150 2026-03-19 10:14 by 功夫疯狂
[考研] 274求调剂 +6 S.H1 2026-03-18 6/300 2026-03-19 09:34 by 花店相见
[考研] 085410人工智能专硕317求调剂(0854都可以) +3 xbxudjdn 2026-03-18 3/150 2026-03-18 22:14 by zhq0425
[考研] 311求调剂 +11 冬十三 2026-03-15 12/600 2026-03-18 14:36 by 星空星月
[考研] 297求调剂 +8 戏精丹丹丹 2026-03-17 8/400 2026-03-18 14:30 by laoshidan
[考研] 070300化学319求调剂 +6 锦鲤0909 2026-03-17 6/300 2026-03-18 13:22 by Iveryant
[考研] 302求调剂 +10 呼呼呼。。。。 2026-03-17 10/500 2026-03-18 12:45 by Linda Hu
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[考研] 304求调剂 +3 曼殊2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 266求调剂 +4 学员97LZgn 2026-03-13 4/200 2026-03-14 08:37 by zhukairuo
[考研] 304求调剂 +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
[考研] 290求调剂 +3 ADT 2026-03-13 3/150 2026-03-13 10:19 by peike
信息提示
请填处理意见