二分法求函数方程根的matlab实现(内附例题和代码)

endeavor98: 请问一下矩阵R是如何定义的？非常感谢

m0_74305919: % Steepset descent method clear; clc; eps=0.005; X=[0 0]'; syms x1 x2 f=x1-x2+2*((x1)^2)+2*x1*x2+(x2)^2; vector=[x1,x2]; [min_x,min_f,k]=Steepset_descent(f,vector,X,eps); disp('Calculation completed!') disp(['Minimal value:',32,num2str(min_f)]) disp(['The vector corresponding to the minimum value:', ... 32,num2str(min_x(1)),32,num2str(min_x(2))]) disp(['Iteration number:',32,num2str(k)]) function [min_x,min_f,k]=Steepset_descent(f,vector,X,eps) syms lamda f=sym(f); vectors=vector.'; grad_=gradient(f,vector); k=0; n=1; while n>0 S=double(subs(grad_,vectors,X)); if norm(S,2)>eps ff=subs(f,vectors,(X-lamda*S)); fun=simplify(diff(ff,lamda)); lamda1=double(solve(fun)); X=double(X-lamda1*S); k=k+1; else break; end end min_x=X; min_f=double(subs(f,vectors,min_x)); end

magic_herb: 输入的初始解向量x0要是列向量，也就是x0=[0,0,0].'

心无旁骛~: 2023版的MATLAB代码：https://blog.csdn.net/m0_63007797/article/details/133531585?csdn_share_tail=%7B%22type%22%3A%22blog%22%2C%22rType%22%3A%22article%22%2C%22rId%22%3A%22133531585%22%2C%22source%22%3A%22m0_63007797%22%7D

Maojinqiang1: 你好，博主，使用2021版的程序，运行你的代码，输入你给的数据一直提示 错误使用 * (第 432 行)Dimensions do not match.有遇到过吗