（课件）模型实验七.pptVIP

模型实验（七） 铅球投掷的研究 软件调查 MATLAB优化工具箱能求解的优化模型 fsolve：求解非线性方程组 求解方程组：F(X)=0 X=fsolve(FUN,X0) lsqcurvefit：非线性最小二乘法 求解 min sum{(FUN(X,XDATA)-YDATA).^2} X=lsqcurvefit(FUN,X0,XDATA,YDATA) 1.化简模型 2.整理数据 α v h s 李梅素 40.27 13.16 2.20 19.40 隋新梅 39.00 13.95 2.04 21.66 李梅素 38.69 13.51 2.00 20.30 黄志红 37.75 13.58 2.02 20.76 李梅素 37.60 13.75 1.90 20.95 李梅素 35.13 14.08 1.95 21.76 A1=2*9.8*(1-sinα), A2=2*cos α, A3=v2. A=[A1, A2, A3] 非线性方程 (1 1 1)’*((x(1)-g)^2*x(2)^2+x(3)^2)+ +A*[x(1)*x(2)^2 x(1)*x(2)*x(3) -1]’=0 用fsolve求解 或使用更多的数据用lsqcurvefit求解 * * 投掷铅球出手速度的模型 观测数据: α v h s 李梅素 40.27 13.16 2.20 19.40 隋新梅 39.00 13.95 2.04 21.66 李梅素 38.69 13.51 2.00 20.30 黄志红 37.75 13.58 2.02 20.76 李梅素 37.60 13.75 1.90 20.95 李梅素 35.13 14.08 1.95 21.76 如何使用这些数据估计未知的参数 F, t0, v0? 优化工具箱3.0 (MATLAB 7.0 R14) 连续优化 离散优化 无约束优化 非线性 极小 fminunc 非光滑(不可 微)优化 fminsearch 非线性 方程(组) fzero fsolve 全局 优化 暂缺 非线性 最小二乘 lsqnonlin lsqcurvefit 线性规划 linprog 0-1规划 bitprog 一般ＩＰ(暂缺) 非线性规划 fmincon fminimax fgoalattain fseminf 上下界约束 fminbnd fmincon lsqnonlin lsqcurvefit 约束线性 最小二乘 lsqnonneg lsqlin 约束优化 二次规划 quadprog 令F/m=x(1), t0=x(2), v0=x(3), 为所求的未知变量。