- #include<stdio.h>
- #include<math.h>
- double resistance;
- double func(double x) //original function
- {
- return -5.775e-7*x*x+3.9083e-3*x+1-resistance/100.0;
- }
- double func1(double x) //derivative function
- {
- return -5.775e-7*2.0*x+3.9083e-3;
- }
- int Newton(double *x,double precision,int maxcyc)
- {
- double x1,x0;
- int k;
- x0=*x;
- for(k=0;k<maxcyc;k++)
- {
- if(func1(x0)==0.0)
- {
- printf("derivative =0! on iterative process.\n");
- return 0;
- }
- x1=x0-func(x0)/func1(x0);
- if(fabs(x1-x0)<precision || fabs(func(x1))<precision)
- {
- *x=x1;
- printf("root near the value.%lf\n",*x);
-
- return 1;
- }
- else
- {
- x0=x1;
- }
- }
- printf("iterative repetition overflow.\n");
- return 0;
- }
-
- int main(int argc, char** argv)
- {
- double x,precision;
- int maxcyc;
- printf("resistance of Pt100(85~200):");
- scanf("%lf",&resistance);
- printf("initial iterative value x0:");
- scanf("%lf",&x);
- printf("max iterative repetition：");
- scanf("%d",&maxcyc);
- printf("iterative precision：");
- scanf("%lf",&precision);
- if(Newton(&x,precision,maxcyc)==1)
- {
- }
- else
- {
- printf("iteration failed！\n");
- }
- return 0;
- }
我給你寫了個牛頓迭代代碼，不要害怕這么長的代碼，主要目的還是展示牛頓迭代的原理和所需要的參數，核心代碼就這么幾行，其他的都可以在移植到單片機的時候刪掉。
- for(k=0;k<maxcyc;k++)
- {
- x1=x0-func(x0)/func1(x0);
- if(fabs(x1-x0)<precision || fabs(func(x1))<precision)
- {
- *x=x1;
- return 1;
- }
- else x0=x1;
- }
至于如何從AD讀數來推導出電阻值，我算了下，確實是初中數學問題，就不貼了。 |
[復制鏈接]
QQ好友和群
QQ空間
騰訊微博
騰訊朋友
收藏
淘帖
頂
踩
