牛顿法和最速下降法求解非线性方程_第1页
牛顿法和最速下降法求解非线性方程_第2页
牛顿法和最速下降法求解非线性方程_第3页
牛顿法和最速下降法求解非线性方程_第4页
牛顿法和最速下降法求解非线性方程_第5页
已阅读5页，还剩3页未读， 继续免费阅读

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

1、这是我们的题目要求非线性方程的寻优求解X2课题内容：求解非线性方程f(x)=血x-7=0，分别利用最速下降法和牛顿法求解该非线性方程的非零实根的近似值。课题要求：1、设计良好的人机交互GUI界面。2、程序要求有注释。3、给定误差精度，对比算法结果。因为目标函数有两个解，一个在0附近，一个在2附近。所以这两种方法使用时都要从0附近取一个初值计算一次得出一个解，再从2附近取一个初值计算一次得出另一个解。两种算法都是求得解之后通过while比较误差精度，直到达到所设置的误差精度后才输出解，但这个过程中计算出来的所有解包括最终结果都会画在gui界面上,让计算过程变得直观。两种方法除了核心计算公式不同外

2、，其他程序基本相似，如果大家要借鉴使用，除了改动目标函数之外，还要确定合适的初值，初值的确定大家可以把目标函数解的附近取值，牛顿法解的收敛条件大家可以在网上找一下，因为如果初值给的不收敛，那么计算会出问题。程序代码(1)牛顿法程序clearclai=1；x(i)=0.3;%清除显示屏%迭代次数%给定0附近的一个初值方程不同初值要改%方程不同初值要改whileabs(x(i)_x(i+1)1*10A(_6)plot(x(i+1),f(x(i+1)*)%比较误差精度holdon%保持图形x(i+1)=0.2;i=i+1;x(i+1)=x(i)-(f(x(i)/f1(x(i);%牛顿迭代法计算公式e

3、nda1=x(end)%输出0附近的解j=1;从2附近取一个初值y(j)=0.2;%方程不同初值要改y(j+1)=1.6;%方程不同初值要改whileabs(y(j)_y(j+1)1*10入(_6)plot(y(j+1),f(y(j+1)+)holdonj=j+1；y(j+1)=y(j)_(f(y(j)/f1(y(j);enda2=y(end)%输出2附近的解x0=inline(f(x);x1=fzero(x0,_1,1)x2=fzero(x0,1,2)nwc1=x1_a1nwc2=x2_a2%对比系统函数fzero算法结果，计算求出解的误差(2)最速下降法clearclastep=0.01;

4、x=0.5;i=0;%步长%初值为0附近方程不同初值要改loss=f(x);whileabs(loss)0.01%比较误差精度plot(x,f(x)*)%做出图形holdon%保持图形x=x_step*f1(x);%最速下降法计算公式loss=f(x);i=i+1;end方程不同初值要改xy=1.5;%初值为2附近时los=f(y);j=0；whileabs(los)0.01plot(y,f(y)+)holdony=y-step*f1(y);los=f(y);j=j+1；endyx0=inline(f(x);x1=fzero(x0,-1,1)x2=fzero(x0,1,2)dwc1=x1-x%对比算法结果dwc2=x2-y函数程序functiona=f(x)a=sin(x)-(x.A2)/

人人文库> 全部分类> 图纸下载 > 毕业设计

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论