追帖:我编辑的程序如下:clc
clear
cross_sub=xlsread('jcbt2');
r=cross_sub(:,8);
e=cross_sub(:,4);
z=cross_sub(:,6);
[m, n]=size(z);
for i=1:m;
f=@(x)z(i)-r(i)*x(i)/(e(i)-x(i))^(1-e(i));
x0=-1.8;
w=fzero(f,x0);
end
还有一个
clc
clear
syms x r e z;
cross_sub=xlsread('jcbt2');
r=cross_sub(:,8);
e=cross_sub(:,4);
z=cross_sub(:,6);
[m, n]=size(z);
for i=1:m;
eqs1='r(i)*x(i)/(e(i)-x(i))^(1-e(i))-z(i)';
x=solve(eqs1,'x');
x=subs(x);
end
|