x=[0.07 0.13 0.28 0.44];
T=[283 296 312 320];
xcal=ones(1,4);
t=find(abs(xcal)==1);
>> l=length(t);
>> while (l>0)
for i=1:4
m=x(i)^2+2*x(i);
n=2*x(i)^2+3*T(i);
o=4*x(i)^2+5*T(i)^2;
p=m+n*o;
q=3/(T(i)^2);
y(i)=exp(q-p);
k=abs(y(i)-x(i));
if k<=0.01
xcal(i)=y(i);
else
x(i)=y(i);
end
end
t=find(abs(xcal)==1);
l=length(t);
end
x
y
xcal
回复此楼