%%
%該文件用來描述目標函數，不可運行
function?[ffd]?=?mbhs（Yk）
%%
%賦值給段數n、繩索長度L：
n=483;?L=200.2132;
%%????
%描述小車位置，確定繩索線密度m（i）的情況（這里假設有小車的地方線密度m為無車繩索的100倍）?
for?i=1:1:n
????if?i==k
????????m（i）=2.7959*100;?%?車是繩線密度的100倍?
????else
????????m（i）=2.7959;
????end
end
%%
%描述目標函數f：
for?i=1:1:n
????S1=0;
????S2=-1/2*Y（i）;
????for?j=1:1:i
????????S1=S1+Y（j）;
????end
????S（i）=（S1+S2）*L/n*m（i）;
end
f=sum（S）;
%%
%描述目標函數梯度向量（偏導數）fd：
if?nargout>1
????for?i=1:1:n
????????g2=0;
????????g1=1/2*L/n*m（i）;
????????for?j=i+1:1:n
????????????g2=g2+L/n*m（j）;
????????end
????????????g（i）=g2+g1;
????end
?%將g（i）轉化為列向量賦給glie
????glie=[];
????for?i=1:1:n
????????glie=[glie;g（i）];
????end
????fd=glie;
end
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????823??2022-10-17?14:22??小車運行軌跡數值模擬程序\mbhs.m

?????文件???????1118??2018-07-31?21:10??小車運行軌跡數值模擬程序\xlx1.m

?????文件???????1118??2022-10-17?05:55??小車運行軌跡數值模擬程序\xlx2.m

?????文件????????679??2022-10-31?14:17??小車運行軌跡數值模擬程序\yshs.m

?????目錄??????????0??2023-01-18?01:18??小車運行軌跡數值模擬程序

-----------?---------??----------?-----??----

?????????????????3738????????????????????5