关于#MATLABfsolve#的问题，如何解决？

Question

clc;
clear;
close all;

a2=14 ;
d=0.8;

xr1=0.1;              
xg1=1;           
xb1=1;  

for i=1:9
    
    xr1(i+1)=0.1+0.1*i;
x0=[3,3,3,3];     
fun=@(x)[
        (2*a2*a2)^(1/2) - (x(4)+2*d+x(2)),...
        (2*a2*a2)^(1/2) - (x(1)+2*d+x(3)),...
        ((x(1)*x(1))/(2*a2*a2))/((x(2)*x(3))/(a2*a2))  - xr1./xg1,...
        ((x(1)*x(1))/(2*a2*a2))/((x(4)*x(4))/(2*a2*a2))- xr1./xb1,...
        ((x(2)*x(3))/(a2*a2))/((x(4)*x(4))/(2*a2*a2))  - xg1./xb1,...
       ];
x=fsolve(fun,x0);
Xr = roundn(x(1),-2);
Xg = roundn(x(2),-2);
Yg = roundn(x(3),-2);
Xb = roundn(x(4),-2);
xr2=Xr*Xr/2/a2/a2;
xg2=Xg*Yg/a2/a2;
xb2=Xb*Xb/2/a2/a2;
x1=xr2/xg2;
x2=xg2/xg2;
x3=xb2/xg2;
fprintf('xr1 xg1 xb1: %.2f  %.f  %.f',x1,x2,x3);fprintf('\n');
fprintf('\n');

end

输入xr1=[0:0.2:1];
拿解出来的x验证，xr1=[0.15:0.05:0.55]
解出来的x值不对

Answer 1

检查一下公式有没有问题，如果没有问题的话，可能不是这个函数的问题，有可能是程序数值误差方面设计的不够好，这种求x的也可能有多个解，一般可以采用目标函数进行评估，或者可以代入方程进行评估，不一定就是程序的问题