24小时热门版块排行榜

返回列表

当前主题已经存档。

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

wwl5213384

木虫 (正式写手)

小晨

应助: 0 (幼儿园)
金币: 2412.8
红花: 1
帖子: 667
在线: 68.4小时
虫号: 589412
注册: 2008-08-29
性别: MM
专业: 电气科学与工程

[交流] 【求助】100金币求助matlab程序

是用matlab编程解个方程组，望大家多多帮忙，非常感谢
http://d.namipan.com/d/346972521546d30d71a1911e42a8b3cb151d2ec3005a0000

回复此楼

» 猜你喜欢

天天开心，用心生活。

1楼 2009-12-04 11:24:16

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

bluesine

铁杆木虫 (职业作家)

科苑小木虫

数学EPI: 5
应助: 132 (高中生)
贵宾: 1.991
金币: 9582.7
散金: 89
红花: 19
帖子: 3554
在线: 371.7小时
虫号: 869544
注册: 2009-10-12
性别: GG
专业: 数学物理

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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