很简单的方程求解程序,调用mulDNewton函数求解,之前在Matlab 2011b版本上运行成功,现在在Matlab 2018a版本上总是出错,程序代码和出错的提示如下,mulDNewton函数代码也如下。
请教大神该程序应该如何修改,请给出正确的代码,必有重谢!@所有人
mulDNewton迭代函数代码如下:
function [r,m]=mulDNewton(F,x0,eps)
%非线性方程组:F
%初始解:x0
%解的精度:eps
%求得的一组解:r
%迭代步数:m
if nargin==2
eps=1.0e-4;
end
x0 = transpose(x0);
dF = jacobian(F);
m=1;
tol=1;
while tol>eps
tol=1;
w=1;
Fx = subs(F,findsym(F),x0);
dFx = subs(dF,findsym(dF),x0);
F1=norm(Fx);
while tol>=0 %下面的循环是选取下山因子w的过程
r=x0-w*inv(dFx)*Fx; %核心的迭代公式
Fr = subs(F,findsym(F),r);
tol=norm(Fr)-F1;
w=w/2;
end
tol=norm(r-x0);
m=m+1;
x0=r;
if(m>10000) %迭代步数控制
disp('迭代步数太多,可能不收敛!');
return;
end
end
方程程序代码如下:
clc
clear
syms x y z
f1=x^2+y+z^2-12
f2=x^3+y^2+z-8
f3=x+y+z-6
x0=[1 1 1]% 初值
f=[f1;f2;f3]%方程组
[r,b]=mulDNewton(f,x0,1.0e-6)%调用迭代函数
%求解方程,并将结果赋予N01,N02,N03
N01=r(1);
N02=r(2);
N03=r(3);
% 求解的正确结果应该是X=1,y=2, z=3
错误提示代码如下:
警告: FINDSYM will be removed in a future release. Use SYMVAR instead.
> In sym/findsym (line 32)
In mulDNewton (line 19)
In equation (line 11)
错误使用 sym>convertChar (line 1448)
Character vectors and strings in the first argument can only specify a variable or number. To evaluate character vectors and strings representing symbolic
expressions, use 'str2sym'.
出错 sym>tomupad (line 1214)
S = convertChar(x);
出错 sym (line 211)
S.s = tomupad(x);
出错 sym/subs>normalize (line 168)
X = {sym(X)};
出错 sym/subs>mupadsubs (line 157)
[X2,Y2,symX,symY] = normalize(X,Y); %#ok
出错 sym/subs (line 145)
G = mupadsubs(F,X,Y);
出错 mulDNewton (line 19)
Fx = subs(F,findsym(F),x0);
出错 equation (line 11)
[r,b]=mulDNewton(f,x0,1.0e-6)%调用迭代函数