回 帖 发 新 帖 刷新版面

主题:请大侠指点下这个程序

function f=hehe(x)
[v]=[60 50 40];
[t]=[143 130 116];
[s]=[907 722 534];
f1=x(1)*exp(2*x(3)*s(1)+x(2)*t(1))-(x(1)+x(2)*v(1)+x(3)*v(1).^2);
f2=x(1)*exp(2*x(3)*s(2)+x(2)*t(2))-(x(1)+x(2)*v(2)+x(3)*v(2).^2);
f1=x(1)*exp(2*x(3)*s(3)+x(2)*t(3))-(x(1)+x(2)*v(3)+x(3)*v(3).^2);
f=[f1;f2;f3];

x0=[0.0075;0.0002;0.00006];                     % 初始值
% 使用fsolve函数求解非线性方程组
[x,f]=fsolve(@hehe,x0);
disp '            ********  计 算 结 果  ********'  
fprintf (1,'             滚阻系数        a = %3.4f mm \n',x(1))
fprintf (1,'             滚阻常数        b = %3.4f mm \n',x(2))
fprintf (1,'             风阻系数        e = %3.4f mm \n',x(3))
disp '           ********  计算精度检验  ********'  
fprintf (1,'          方程1函数值       f1 = %3.4e \n',f(1))
fprintf (1,'          方程2函数值       f2 = %3.4e \n',f(2))
fprintf (1,'          方程3函数值       f3 = %3.4e \n',f(3))

为什么会出现
[color=FF0000]??? Undefined function or variable 'f3'.

Error in ==> hehe at 8
f=[f1;f2;f3];
Error in ==> fsolve at 180
        fuser = feval(funfcn{3},x,varargin{:});[/color]

回复列表 (共2个回复)

沙发

[quote]
f1=x(1)*exp(2*x(3)*s(1)+x(2)*t(1))-(x(1)+x(2)*v(1)+x(3)*v(1).^2);
f2=x(1)*exp(2*x(3)*s(2)+x(2)*t(2))-(x(1)+x(2)*v(2)+x(3)*v(2).^2);
[quote]f1=x(1)*exp(2*x(3)*s(3)+x(2)*t(3))-(x(1)+x(2)*v(3)+x(3)*v(3).^2);[/quote]
[/quote]

板凳


呵呵,谢谢

我来回复

您尚未登录,请登录后再回复。点此登录或注册