【Lesson 1】9月14日

【Lesson 2】9月28日

（1）優化設計一般實施步驟

1. 設計要求與目的→定義優化設計問題
2. 建立問題的數學模型
3. 選用合適的優化算法
4. 確定必要數據和設計初始點
5. 編寫程序，利用計算機求解計算最優備選方案
6. 分析結果數據和方案，判斷是否符合工程實際，進行參數分析

（2）優化問題建模：選擇和確定各個獨立參數

• 表示：$X=[x_1,x_2,...,x_n{]}^T$
• 維數：變量數 $n$
• 類型：連續/離散變量
• 設計空間：各獨立變量構成的空間
• 超空間：$n>3$ 時的設計空間

（3）目標函數：設計中預期達到的目標

• 表示：$f(X)=f(x_1,x_2,...,x_n)$
• 分類：單目標函數、多目標函數、綜合目標函數

（4）約束條件：設計變量和響應取值的限制條件

• 邊界約束/區域約束
• 性能約束
• 顯約束/隱約束

（5）數學模型：
$$\begin{array}{rl}& \mathrm{f}\mathrm{i}\mathrm{n}\mathrm{d}:X=[x_1,x_2,...,x_n{]}^T——設計變量\\ & \mathrm{m}\mathrm{i}\mathrm{n}:f(X)——目標函數\\ & \begin{array}{c}\begin{array}{rl}\mathrm{s}.\mathrm{t}.：& g_u(X)\ge 0(u=1,2,...,m)\\ & h_v(X)=0(u=1,2,...,p)\end{array}\end{array}\}——約束條件\end{array}$$
① 約束優化與無約束優化；② 線性優化與非線性優化；③ 單目標優化與多目標優化；④ 可行設計點與可行域；⑤ 局部最優與全局最優。

復習：

1. 泰勒展開

2. 等等

優化設計方法分類;

• 任務：① 找到最優解；② 提高求優過程的效率。
• 對于非線性優化問題，基本是使用迭代算法

【Lesson 3】10月12日

迭代算法：按照一定的邏輯結構進行反復的數值計算，尋求函數值不斷下降的設計點，直到最后獲得足夠精度的近似解時就截斷計算。

迭代種植條件

1. 點距準則
   $$||x^{(k+1)}?x^{(k)}||\le \epsilon$$
2. 目標函數下降量準則
   $$\begin{array}{rl}若|f(x^{(k+1)})|?1& ,|f(x^{(k+1)})?f(x^{(k)})|\le \epsilon \\ 若|f(x^{(k+1)})|?1& ,\frac{|f(x^{(k+1)})?f(x^k)|}{|f(x^{(k)})|}\le \epsilon \end{array}$$
3. 梯度準則
   $$||▽f(x^{(k+1)})||\le \epsilon$$

設計變量的坐標變換

• 原因：變量的數值大小相差極大，無法用統一步長尋優，無法使用統一的終止條件，需要采用坐標變換對設計變量進行標準化處理。
• 標準化處理：將設計變量用變量區間內的相對坐標來表示，將設計變量真實值化為變量區間的相對值，是變化范圍在0~1之間。
  $$\begin{gathered} 對于x_{li}<x_{ti}<x_{ui},設計參數x_{ti} \\ {x}_{ri}=\frac{x_{ti}?x_{li}}{x_{ui}?x_{li}}\in [0,1] \end{gathered}$$

優化設計方法分類

一維優化方法

一維優化方法是多維優化方法的基礎
要求出
$$\begin{gathered} X^{(K+1)}=X^{(k)}+{\alpha }^{(k)}{S}^{(k)} \\ \begin{array}{rl}\mathrm{m}\mathrm{i}\mathrm{n}f(X^{(K+1)})=& \mathrm{m}\mathrm{i}\mathrm{n}f(X^{(k)}+{\alpha }^{(k)}{S}^{(k)})\\ \approx & f(X^{(k)})+[▽f(X^{(k)}){]}^T\alpha S^{(k)}+\frac{1}{2}{\alpha }^2[S^{(k)}{]}^{T}H(X^{(k)})S^{(k)}\end{array} \\ {f}^{\prime }(x)\approx \frac{f(x+\Delta x)?f(x)}{\Delta x}=0?\alpha =\frac{[▽f(X^{(k)}){]}^T{S}^{(k)}}{[S^{(k)}{]}^{T}H(X^{(k)})S^{(k)}} \end{gathered}$$

步驟：

• 確定搜索區間
• 確定步長求解最優

進退法

1. 試探計算
   初選一個初始點 $x_1$ 和初始步長 $h_0$，前進點為 $x_2=x_1+h_0$，計算函數值 $y_1=f(x_1)$ 和 $y_2=f(x_2)$
   （1）當 $y_2>y_1$ 時，
2. $y_2<y_1$ 前進計算
   計算第二個前進點
3. 后退運算

目的是：形成一個“高-低-高”的區間

算例：

【2021.11.2】

變尺度法

DFP變尺度法能夠適用于多參數函數 $f(x)$ 優化，下面以例子來說明。

使用DFP變尺度法對函數 $f(x)=4(x_1?5{)}^2+(x_2?6{)}^2$ 進行優化，初始點為 $[00{]}^T$ ，要求迭代結束條件為梯度 $\xi$ 小于0.01。

該函數的梯度向量為

$$▽f(x)={\left[\begin{array}{cc}8(x_1?5)& 2(x_2?6)\end{array}\right]}^T$$

初始變量值、梯度為

$$\begin{gathered} x^{(0)}={\left[\begin{array}{cc}0& 0\end{array}\right]}^T \\ ▽f(x^{(0)})={\left[\begin{array}{cc}?40& ?12\end{array}\right]}^T \end{gathered}$$

初始海塞矩陣為

$$H^{(0)}=\left[\begin{array}{cc}1& 0\\ 0& 1\end{array}\right]$$

第一次迭代后的變量值為

$$x^{(1)}=x^{(0)}+\alpha S^{(0)}=x^{(0)}?\alpha H^{(0)}▽f(x^{(0)})={\left[\begin{array}{cc}40\alpha & 12\alpha \end{array}\right]}^T$$

則第一次迭代后的函數值為

$$f(x^{(1)})=4(40\alpha ?4{)}^2+(12\alpha ?6{)}^2=6544{\alpha }^2?1744\alpha +136$$

該值取最小時 $\alpha =109/818$，將其帶入 $x^{(1)}$ 和 $▽f(x^{(1)})$。有

$$x^{(1)}={\left[\begin{array}{cc}5.33& 1.60\end{array}\right]}^T,▽f(x^{(1)})={\left[\begin{array}{cc}2.64& ?8.8\end{array}\right]}^T$$

可分別計算變量和梯度的差值

$$\begin{gathered} △x^{(0)}=x^{(1)}?x^{(0)}={\left[\begin{array}{cc}5.33& 1.60\end{array}\right]}^T \\ △g^{(0)}=▽f(x^{(1)})?▽f(x^{(0)})={\left[\begin{array}{cc}42.64& 3.2\end{array}\right]}^T \end{gathered}$$
變尺度法中海塞矩陣的迭代矯正公式為

$$\begin{gathered} E^{(k)}=\frac{△x^{(k)}[△x^{(k)}{]}^T}{[△x^{(k)}{]}^T△g^{(k)}}?\frac{H^{(k)}△g^{(k)}[△g^{(k)}{]}^T{H}^{(k)}}{[△g^{(k)}{]}^T{H}^{(k)}△g^{(k)}} \\ {H}^{(k+1)}=H^{(k)}+E^{(k)} \end{gathered}$$

故可以計算出

$$\begin{array}{rl}{H}^{(1)}& =H^{(0)}+E^{(0)}\\ & \\ & =\left[\begin{array}{cc}1& 0\\ 0& 1\end{array}\right]+\frac{△x^{(0)}[△x^{(0)}{]}^T}{[△x^{(0)}{]}^T△g^{(0)}}?\frac{H^{(0)}△g^{(0)}[△g^{(0)}{]}^T{H}^{(0)}}{[△g^{(0)}{]}^T{H}^{(0)}△g^{(0)}}\\ & \\ & =\left[\begin{array}{cc}0.128& ?0.038\\ ?0.038& 1.104\end{array}\right]\end{array}$$
由此，可以依次邏輯繼續向下迭代，直至最終使得梯度小于0.01。

總結

以上是生活随笔為你收集整理的现代机械优化设计的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。