ここに私の方程式があります：（たわみと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' undefined. 問題はOB関数とfminsearchにあると思います。
誰か助けてくれませんか？