这是我的方程式:(只有偏转和M(力矩)方程很重要)
M=@(x,F) -(w/2)*(x.^2) + F*x + Fb*(x-L/2).*heaviside(x-L/2);
deflection=@(x,c1,c2,F) (1/EI)*(-w*(x.^4)/24 + F*(x.^3)/6 + (1/6)*Fb*((x-L/2).^3).*heaviside(x-L/2) + c1*x + c2);
我写了一个目标函数来找到c1和c2常数以及F值来优化我的弯矩。
OB1=@(c1,c2,F)(deflection(0,c1,c2,F)).^2 + (deflection(L,c1,c2,F)).^2 + (deflection(L/2,c1,c2,F)).^2;
[c] = fminsearch(@(c,F)OB1(c(1),c(2),c(3)),[-0.5 0 1]);
C1=c(1);
C2=c(2);
f=c(3);
给我错误信息:'f'未定义。我认为问题在于OB功能和fminsearch。
有人可以帮忙吗?