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值不对