Question

用MATLAB求解复杂的方程式中的未知数？求编程。

函数y=1-(a*(0.00412/4.84)^b)*(x/570)^c，其中a，b，c为未知数。
x 分别取50 100 150 200 260
y 分别取0.988 0.984 0.978 0.973 0.968

Answer 1

可以用数据拟合的方法来求解复杂的方程式中的未知数。解决方法：

a0=[0.4341 0.33535 0.43285];

func1=@(a,x)1-(a(1)*(0.00412/4.84).^a(2))*(x/570).^a(3);

[a,r,J] = nlinfit(x,y,func1,a0);

y1=func1(a,x);

A=[y1 y];

disp('        y1            y')

disp(A)

a0=double(a(1));b0=double(a(2));c0=double(a(3));

disp('模型: y=1-(a*(0.00412/4.84)^b)*(x/570)^c')

str=['a=',num2str(a0),'，b=',num2str(b0),'，c=',num2str(c0)];

fprintf('%s\n',str)

运行结果

a=0.49409，b=0.31618，c=0.6484