- c - 在位数组中找到第一个零
- linux - Unix 显示有关匹配两种模式之一的文件的信息
- 正则表达式替换多个文件
- linux - 隐藏来自 xtrace 的命令
我编写了这段 C 代码以使用二分法比较两个函数的根。我的第一个函数 (g(x)) 执行正确,但第二个函数 (h(x)) 在屏幕上输出“#1QO”。我在代码中找不到我做错了什么。
你能解释一下吗?任何形式的帮助将不胜感激。谢谢!
#include <stdio.h>
#include <math.h>
typedef double (*DFD) (double);
double g (double x)
{
double y;
y = pow (x,3) - pow (x,2)-1;
return y;
}
double h (double x)
{
double k;
k = 1.0 + (1.0/x) + (1.0 /pow (x,2));
return k;
}
double bisection (DFD f, double x0, double x1,double tol)
{
int i;
double middle;
for (i=1;i<=50;) {
middle = (x0+x1)/2.0;
if (fabs (middle - x0) <tol) return middle;
if (f(middle)* f(x0) <0.0) x1 = middle;
else x0 = middle;
}
}
int main () {
double root_gx, root_hx = 0.0;
root_gx = bisection (g,0,2,0.0005);
printf ("Root found using g(x) = %.3lf\n",root_gx);
root_hx = bisection (h,1,2,0.0005);
printf ("Root found using h(x) = %.3lf\n",root_hx);
printf ("Difference between the two roots = %.3lf\n", (fabs (root_gx- root_hx)));
return 0;
}
编辑:在二分法中初始化 i=1 并将 bisection (h,0,2,0.0005)
更改为 bisection (h,1, 2,0.0005)
成功了 谢谢大家!
最佳答案
值未初始化导致未定义的行为。设置为某个值。
// int i;
int i = 0;
...
for (;i<=50;) {
i
永远不会递增。50
是任意的。
建议使用 double
的二进制精度作为 bisection()
将每次迭代减半。
for (i=0; i <= DBL_MANT_DIG; i++) {
...
}
return middle;
还建议更改算法以允许更改 |difference| 0.0 以满足 0.0 的公差。
// if (fabs (middle - x0) <tol) return middle;
if (fabs (middle - x0) <= tol) return middle;
关于c - 在 C 中使用二分法的两个根之间的区别?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/34374709/
我正在尝试编写一个程序来确定给定函数 (f(x) := ln((sin(x**(1/2))**3) + 2) - 1 的零点,使用二分法。值 a 和 b,它们是二分法中使用的初始值,已经插入到程序中。
二分法是搜索算法中极其典型的方法,其要求输入序列有序并可随机访问。算法思想为 。 输入:有序数组nums,目的数值target 要求输出:如果target存在在数组中,则输出其index,否则输出-
一旦找到正确的根(r3 的值 == 0),退出 while 循环就会遇到问题。我应该在某个地方休息吗?如果是的话——在哪里? 我已经尝试过使用 if、else if 和 else 语句的条件。尝试放置
我在做 this problem来自 lightoj 法官(抱歉提供链接,我不知道如何添加图片)。这是纯基于几何的问题,我的方法是这个导致接受的解决方案。 代码 #include using
是否有一个库函数对列表/元组执行二进制搜索,如果找到则返回项目的位置,如果没有则返回“False”(-1、None 等)? 我在 bisect module 中找到了函数 bisect_left/ri
我正在为二分法编写代码。 我的代码在下面,不知何故循环似乎没有开始。 没有特殊的编译问题。 我认为变量声明/函数原型(prototype)没有任何问题。 谁能帮我找到真正的问题所在? 最佳答案 标准b
用二分查找很容易找到一个整数even if it can be arbitrarily large :先猜数量级,再继续划分区间。 This answer描述了如何找到任意有理数。 设置场景后,我的问
昨天我在另一个 boost 功能上遇到了问题,但幸运的是你们帮助我解决了这些问题。今天我需要知道如何正确使用二分函数。 所以这就是我认为它应该如何工作,但似乎我也弄错了。好的,我想使用: templa
本文实例总结了PHP字符串逆序排列实现方法。分享给大家供大家参考,具体如下: 关于字符串的逆序排列,最简单的使用PHP函数strrev()的测试代码如下: ?
如下所示: ? 1
一. lamda匿名函数 为了解决一些简单的需求而设计的一句话函数 ?
我是一名优秀的程序员,十分优秀!