非线性方程求解！论文急用，大家救命！

herotears · 发表于 2006-1-4 17:19

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

上个月小弟在这里请教过一次，感谢大家给予的指导。现在又遇到一个问题，我用matlab写了个割线法迭代求解非线性方程，但是结果却不令人满意，经常会出来复数解，请问高手如何解决？是方程本身的问题还是算法的问题？或者是初值给的不好，小弟最近就要交论文，而这个问题还没有解决，很着急，请大家不吝赐教！我在这里先行谢过！！！ 非线性方程：y2./D-dp./(2*D)-(b*(1-2*a)*(8-12*y2./(a.^2*D)).^(1/3))./((1-a.^2*(8-12*y2./(a.^2*D)).^(2/3)).^(1/2))＝0，未知数为y2。其他已知，对于这个方程，我想变化D，对于不同的D求出不同的y2，再得出D与y2的关系。 我编了一个函数： function y=fun2(y2,D,dp,b,a) y=y2./D-dp./(2*D)-(b*(1-2*a)*(8-12*y2./(a.^2*D)).^(1/3))./((1-a.^2*(8-12*y2./(a.^2*D)).^(2/3)).^(1/2)); 割线法程序： y20=0; y21=0.0001; y22=y21-(fun2(y21,D,dp,b,a))*(y21-y20)./((fun2(y21,D,dp,b,a))-(fun2(y20,D,dp,b,a))); n=1; while(abs(y21-y20)>=1e-12)&(n<=10000) y20=y21;y21=y22; y22=y21-(fun2(y21,D,dp,b,a))*(y21-y20)./((fun2(y21,D,dp,b,a))-(fun2(y20,D,dp,b,a))); n=n+1; end

NASA · 发表于 2006-1-5 10:33

x0=input('input the initial value: x0=') %x0=10; %x1=2; x1=x0-f(x0)/g(x0); ge=1; flag=0; S=10^3; gmax=100; EPS=0.00001; while 1 y0=f(x0); y1=f(x1); k=1; while y0==y1; x0=x0+k*EPS; y0=f(x0); k=5*k; if k>s flag=1; break; end end if flag==1 disp('can not find the needed initial value x0'); break; end x=x1-(x1-x0)*f(x1)/(f(x1)-f(x0)); if abs(x-x1)<EPS x1=x; break; end x0=x1; x1=x; ge=ge+1; if ge>gmax disp('发散或在给定次数内无法求解') break; end end ge x1 试一下这个割线法程序

herotears · 发表于 2006-1-5 17:43

多谢拉，好像很复杂啊，偶是采鸟，回头我好好研究研究，现在没有时间搞这个了，我找到一个曲线救国的方法，不解它了先！

账号		自动登录	找回密码
密码			我要加入

非线性方程求解！论文急用，大家救命！

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（herotears）非线性方程求解！论文急用，大家救...