24小时热门版块排行榜

返回列表

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

jinleisheng

应助: (幼儿园)
在线:
虫号: 991781

[交流] 【求助】求助 matlab 一个循环问题已有2人参与

a=0.2;b=0.2;c=5.7;
s=-c/(2*a)+sqrt(-b/a+(c/(2*a))^2);
>> for i=0:0.01:20
for j=0:0.01:20
f1=sym('(x-a-0.2*(exp(-i*x)-1)-0.2*(exp(-j*x)-1))*(-s+x*(x+a*s+c))+(x+a*s+c)');
f2=sym('(1+i*0.2*exp(-i*x)+j*0.2*exp(-j*x))*(-s+x*(x+a*s+c))+1+(x-a-0.2*(exp(-i1*x)-1)-0.2*(exp(-i2*x)-1))*(2*x+a*s+c)');
fun=inline(f1);
dfun=inline(f2);
x=manewton(fun,dfun,0.1,1e-5);
if real(x)<0
for k=1:1000
d(k)=i
d2(k)=j
end
end
end
end 以上是我的程序，我目的是让f1/f2中的i j 变量循环，从而实现inline的循环，但是这样不行。谁能帮助我一下，解决inline命令中函数某个变量循环的问题。先谢了！！！

回复此楼

» 猜你喜欢

1楼 2010-05-24 11:29:38

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

jinleisheng

应助: (幼儿园)
在线:
虫号: 991781

引用回帖:

Originally posted by jinleisheng at 2010-05-24 11:29:38:
a=0.2;b=0.2;c=5.7;
s=-c/(2*a)+sqrt(-b/a+(c/(2*a))^2);
>> for i=0:0.01:20
for j=0:0.01:20
f1=sym('(x-a-0.2*(exp(-i*x)-1)-0.2*(exp(-j*x)-1))*(-s+x*(x+a*s+c))+(x+a*s+c)');
f2=sym('(1+i*0.2* ...

是运行错误

回复此楼