matlab運(yùn)行牛頓迭代總提示輸入?yún)?shù)不足怎么回事啊
matlab運(yùn)行牛頓迭代總提示輸入?yún)?shù)不足怎么回事啊,下面是程序,
迭代程序
function x=newton(fnew,dfnew,x0,delta,maxx)
for k=1:maxx
x1=x0-fnew(x0)/dfnew(x0);
err= x1-x0;
if err<delta;
break
end
end
函數(shù):
function y=fnew(x)
y=x^2-0.388*x-0.026-10^6*x^3/(8314*473.15)
end
導(dǎo)函數(shù):
function dy=dfnew(x)
dy=2*x-0.388-3*10^6*x^2/(8314*473.15)
end

1592890520(1).png
返回小木蟲查看更多
今日熱帖
京公網(wǎng)安備 11010802022153號(hào)
fnew和dfnew作為參數(shù)傳遞需要先@一下的,可以這樣調(diào)用:
f1=@fnew;
f2=@dfnew;
newton(f1,f2,x0,delta,maxx);,
幫幫忙
提醒一下,程序缺少了更新步