24小时热门版块排行榜    

查看: 374  |  回复: 2

jinleisheng

[交流] 【求助】求助 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命令中函数某个变量循环的问题。先谢了!!!
回复此楼

» 猜你喜欢

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

ilovey

金虫 (小有名气)

是结果不对,还是运行错误
2楼2010-05-24 23:52:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jinleisheng

引用回帖:
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* ...

是运行错误
3楼2010-05-25 10:59:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jinleisheng 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见