24小时热门版块排行榜    

查看: 920  |  回复: 14
当前主题已经存档。

Works

银虫 (正式写手)

什么程序???????????????
回首向来萧瑟处,归去,也无风雨也无晴
11楼2009-12-21 11:12:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ljtwcx

至尊木虫 (知名作家)

楼主加油,翻下书就有了
12楼2009-12-21 13:16:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hellogod

木虫 (正式写手)

是这样的吗?

★ ★
formleaf(金币+1,VIP+0):谢谢你的参与 12-21 19:43
wwl5213384(金币+1,VIP+0):谢谢,呵呵 12-22 19:22
理论上讲,matlab有程序可以解吧。不用自己编的。我记得是这样的啦。可以随便查一下matlab的书就行了。如果是老师的作业那就要自己写。
13楼2009-12-21 19:18:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bluesine

铁杆木虫 (职业作家)

科苑小木虫

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
formleaf(金币+3,VIP+0):谢谢bluesine的解答 12-22 12:32
wwl5213384(金币+30,VIP+0):谢谢你,我先试试吧 12-22 19:18
程序and结果
结果显示:a=0.75时,迭代次数达到50就收敛了;
结果为x1=x2=0,x3=x4=-1;
ps:初值全都设为1,当然也可以设置成其他的


clc;
a=0.75;
x1=1;
x2=1;
x3=1;x4=1;
for k=1:1000
    x1=13.5/4*(4*x4/(1-4*x4))^(a-1)...
        *exp((-1.47+0.56*(a-1))/(8.6*10^(-2)));
    x2=4.5/4*(4*x4/(1-4*x4))^a...
        *exp((-1.92+0.56*a)/(8.6*10^(-2)));
    x3=3/4*(4*x4/(1-4*x4))^(-1)...
        *exp(-9.89*10^(-2));
    x4=x3+(1-a)*x1-a*x2+(a-3/4);
    if(mod(k,50)==0)
        disp(['k=',int2str(k)]);
        disp(['[x1,x2,x3,x4]=',int2str([x1,x2,x3,x4])]);
    end
end


results:
k=50
[x1,x2,x3,x4]=0  0 -1 -1
k=100
[x1,x2,x3,x4]=0  0 -1 -1
k=150
[x1,x2,x3,x4]=0  0 -1 -1
k=200
[x1,x2,x3,x4]=0  0 -1 -1
k=250
[x1,x2,x3,x4]=0  0 -1 -1
k=300
[x1,x2,x3,x4]=0  0 -1 -1
k=350
[x1,x2,x3,x4]=0  0 -1 -1
k=400
[x1,x2,x3,x4]=0  0 -1 -1
k=450
[x1,x2,x3,x4]=0  0 -1 -1
k=500
[x1,x2,x3,x4]=0  0 -1 -1
k=550
[x1,x2,x3,x4]=0  0 -1 -1
k=600
[x1,x2,x3,x4]=0  0 -1 -1
k=650
[x1,x2,x3,x4]=0  0 -1 -1
k=700
[x1,x2,x3,x4]=0  0 -1 -1
k=750
[x1,x2,x3,x4]=0  0 -1 -1
k=800
[x1,x2,x3,x4]=0  0 -1 -1
k=850
[x1,x2,x3,x4]=0  0 -1 -1
k=900
[x1,x2,x3,x4]=0  0 -1 -1
k=950
[x1,x2,x3,x4]=0  0 -1 -1
k=1000
[x1,x2,x3,x4]=0  0 -1 -1
板凳要做十年冷文章不发一个字
14楼2009-12-22 10:32:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wwl5213384

木虫 (正式写手)

小晨

已经解决谢谢版主提醒
天天开心,用心生活。
15楼2010-03-28 13:54:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wwl5213384 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见