24小时热门版块排行榜

返回列表

shao0722

金虫 (职业作家)

应助: 0 (幼儿园)
金币: 663.8
散金: 1368
红花: 3
帖子: 3142
在线: 166小时
虫号: 3662352
注册: 2015-01-25
性别: GG
专业: 计算机应用技术

[求助] 关于二分法已有3人参与

请教一下大家：

我的程序大概是这样：
对一个某个变量第一个节点，先估一个值，然后按照这个值使用循环，一直算到最后一个节点，最后一个节点需满足某个条件，
想用二分法做，但是由于有循环，不知道函数该怎么写
这个程序该怎么编啊

回复此楼

» 收录本帖的淘帖专辑推荐

程序

» 猜你喜欢

1楼 2016-08-09 22:14:29

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

FMStation

至尊木虫 (知名作家)

应助: 591 (博士)
贵宾: 0.03
金币: 18174.7
红花: 97
帖子: 8891
在线: 1487.4小时
虫号: 2400059
注册: 2013-04-04
专业: 计算机应用技术

【答案】应助回帖

感谢参与，应助指数 +1

http://edisonshih.pixnet.net/blog/post/30455128-%5bc語言數值分析%5d-方程式求解---bisector%28二分法%29

// ==================================
// filename: BiSector.cpp
// use bisector method to find solution
// author  : Edison.Shih.
// Date : 2010.3.6
// ** all rights resever **
// ==================================

#include <stdio.h>
#include <stdlib.h>
#include <math.h>

#define SMALL (double)(10E-4)
#define NO_FIND_NUM    (double)(-9999.0)

// =======================================
// f(x) = x^3 + 2x^2 -5x +1,
// have sol at (-4, -3), (0,1), (1, 2)
double func(double x){
   return (x*x*x + 2*x*x - 5*x + 1);
}

// =======================================
double bisector(double (*fptr)(double),
                     double low,
                     double high,
                     double small)
{// small: error rate
   double x = (low+high) / 2.0;
   double xvalue = 0.0; // f(x) = f( (h+low)/2 )
   double hvalue = 0.0; // f(h)
   double lvalue = 0.0; // f(low)

   int times = 0;
   //printf("==========================\n"

;
   do{
         xvalue = (*fptr)(x);
         hvalue = (*fptr)(high);
         lvalue = (*fptr)(low);
         /*
         printf("==========================\n"

;
         printf("%d times:\n", times++);

         printf("L :%8.4lf,H :%8.4lf,X :%8.4lf\n",
               low, high, x);
         printf("f(L):%8.4lf,f(H):%8.4lf,f(X) :%8.4lf\n",
               lvalue, hvalue, xvalue);
         */

         if(xvalue>0 && lvalue>0 && hvalue<0) low = x;
         else if(xvalue>0 && lvalue<0 && hvalue>0) high = x;
         else if(xvalue<0 && lvalue>0 && hvalue<0) high = x;
         else if(xvalue<0 && lvalue<0 && hvalue>0) low = x;
         else return (NO_FIND_NUM); // no find ans
         x = (high + low) / 2;
   }while(fabs(xvalue) > small);
   // printf("==========================\n"

;

   return x;
}

// =======================================
int main(int argc, char **argv)
{
   double sol = 0.0;
   for(double x=-5.0; x<=5.0; x=x+1.0)
   {
         sol = bisector(func, x, x+1.0, SMALL);
         printf("sol at (%.0lf,%.0lf):",x,x+1);
         if(sol==NO_FIND_NUM) {
               printf("no find sol\n"

;
         }
         else {
               printf("%.4lf\n",sol);
         }
   }
   return 0;
}