24小时热门版块排行榜    

Znn3bq.jpeg
查看: 935  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +9 小聂爱学习 2026-04-16 11/550 2026-04-17 22:34 by chixmc
[考研] 0854求调剂 +21 门路摸摸 2026-04-15 25/1250 2026-04-17 15:45 by qzxyhcsy
[考研] 收到复试调剂但是去不了 +6 小蜗牛* 2026-04-16 6/300 2026-04-17 10:05 by 涵竹刘
[考研] 307中医考研调剂 +6 于以采蘩 2026-04-14 6/300 2026-04-16 16:20 by qingfeng258
[考博] 申博自荐 +3 Linxia林夏 2026-04-13 3/150 2026-04-16 12:55 by 墨荷之露
[考研] 药学求调剂 +14 喽哈加油 2026-04-14 16/800 2026-04-16 10:15 by beilsong20
[考研] 求调剂 +11 小聂爱学习 2026-04-11 15/750 2026-04-15 21:57 by noqvsozv
[考研] 通信工程求调剂!!! +6 zlb770521 2026-04-14 6/300 2026-04-15 20:00 by 学员JpLReM
[考研] 271求调剂 +35 2261744733 2026-04-11 41/2050 2026-04-14 15:36 by zs92450
[教师之家] 转长聘了 +7 简单化xn 2026-04-13 7/350 2026-04-14 14:50 by xindong
[考研] 考研调剂 +13 长弓傲 2026-04-13 14/700 2026-04-14 14:44 by zs92450
[考研] 食品与营养(0955)271求调剂 +15 升格阿达 2026-04-12 16/800 2026-04-14 13:18 by 浮若_安生
[考研] 300分求调剂 (085501机械专硕,本科扬大) +9 xu@841019 2026-04-11 10/500 2026-04-14 08:48 by 木木mumu~
[考研] 290求调剂 +18 柯淮然 2026-04-12 20/1000 2026-04-13 12:56 by cyh—315
[考研] 339求调剂 +4 hanwudada 2026-04-12 4/200 2026-04-13 12:03 by 蓝云思雨
[考研] 一志愿085802 323分求调剂 +13 drizzle_9 2026-04-12 14/700 2026-04-13 10:26 by Faiz5552
[考研] 344 材料专业 求调剂211 无地域要求 +8 hualkop 2026-04-11 8/400 2026-04-12 22:24 by fqwang
[考研] 346分,工科0854求调剂,专硕 +6 moser233 2026-04-12 7/350 2026-04-12 22:11 by fqwang
[考研] 调剂求助 +6 果然有我 2026-04-11 7/350 2026-04-11 16:22 by 明月此时有
[考研] 一志愿985机械学硕380求调剂 +5 关关雎鸠10 2026-04-11 5/250 2026-04-11 10:10 by 知念。A
信息提示
请填处理意见