C语言实现割线法求零点以及详解割线法
生活随笔
收集整理的這篇文章主要介紹了
C语言实现割线法求零点以及详解割线法
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
問題:求下列函數的零點:
下面式算法的圖解:
以上是割線法圖解的過程:
1.給定初值:x=0,? x1=0.1;
2.計算x,x1對應函數上的點y,y1;
3.通過y,y1點畫出割線交x軸于x2;
4:替換:x=x1,x1=x2;
5:再求出x2對應的函數上的點:y2
6.通過y1,y2再作割線交x軸于x3,
7.替換,循環下去,x最后將逼近零點
問題簡單化:
已經有x,x1,求解出f(x),f(x1):
再用:x1=x2;x=x1;這樣替換,循環去逼近零點。
這個算法和牛頓迭代法的優點在于不用求導函數:
C語言代碼:
#include<stdio.h> #include<math.h> //原函數 double f(double x) {return cos(x)-x*exp(x); }int main() {int count=0;double x=0; //設置迭代初值double x1=0.2;while(fabs(x-x1)>=0.000001) {count++;double tmp=x1;x1=x1-f(x1)*(x1-x)/(f(x1)-f(x)); x=tmp;}printf("牛頓迭代法結果:%lf\n",x);printf("迭代次數:%d\n",count); }結果:
?
總結
以上是生活随笔為你收集整理的C语言实现割线法求零点以及详解割线法的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: fitype拟合多参数函数和遗传算法拟合
- 下一篇: C语言实现最简单的2048小游戏