《二次插值法无约束最优化》由会员分享，可在线阅读，更多相关《二次插值法无约束最优化（6页珍藏版）》请在金锄头文库上搜索。
1、Matlab 实践1、 二次插值法无约束最优化 算法说明：在包含f （x）极小值X0的区间【a b】，给定三点xl、x2、x3,其对应的函数值分别为 fl、f2、f3，且满足x1x2x2，转步骤二；否则转步骤三；步骤二：如果f0f2，则令x1v=x2, x2v=x0，转步骤四；否则令x3v=x0，转步骤四; 步骤三：如果f0xOO(2)?fOvfOO (是是是否fOvfOO (2) ?叵1x00 (1) -xOO (3) Ivxw 或 者IfO-fOO (2) Ivfw?输出x0, fOIxOO (1) -xOO (3) Ivxw 或 者IfOO (1) -fOO (3) Ivfw?/输入原函
2、数f,搜索区 间x,允许误差xw,fw。xOO (3) =xO计算xO, fO输出“搜索区 间太小”xOO (1) =xOO (2)xOO (2) =xO在x内确定初始的三个值xOO,fOO=f(xOO)xOO (3) =xOO (2) xOO (2) =xO4fOO=f (xOO), 计算xO, fOxOO (1) =xO结束算法举例：f (x) = (x2-2) 2/2-1, xU0,5CuMSJid 習 f=inline(k.2-2).2/2-1? j ? k?); a=0; b=5; sw=le-5; fw=le-8; kO f O=ercha(fj a b j fw)kO =1.41
3、42fO =-1.00002、拉压杆系的静不定问题。求各杆的轴力Ni及节点C的位移，已知桁架结构如图所示，各杆 横截面积分别为Ai,材料的弹性模量为E。算法说明：假设各杆均受拉力，C点因各杆变形而引起的x方向位移Ax, y方向位移厶y,由几何关系，的变形方程：N LA L = i匸=A x cos a + A y sin a i=1, ,ni EA i iiL 令 K= L,i EAi故 K N - Ax cos a - Ay sin a =i iii再加上平面共点力系的两个平衡方程L N cos a = P cos aiiiL N sin a = P sin aiii共有n+2个方程，其中包
4、括n个轴力和两个待求位移厶X0,方程组可解。 线性方程组，可用矩阵除法直接解出。算法举例：P=1000，a=pi/2，l1=0.5，a1=pi/4，a2=pi/2，a3=3pi/4，s1=s2=s3=1e-3 ，E=100e9，求 Ni,x,Ay0容易算出P cosa1N = N = 2.92893218813452e0021231 + 2 cos a 1PN = 5.85786437626905 e002331 +2 cos a 1dx =0N l1cos a 1dy= 3= 2.071067811865475e-006EsiCouand Tindov p=1000; a=pi/2; 11=0.5; e=100e9; ni dx dy=caili (p11)请依次输入各根杆的斜角值：pi/4 pi/2 3*pi/4 请依次输入各根杆的横截面积：le-dni =1. 0e+002 *2.928932188134525.857864376269052.92893218813452dx 二dy =2.071067811865475e-006
