信息奥赛一本通1241:二分法求函数的零点

在这里插入图片描述
代码如下:

//1241:二分法求函数的零点
#include<iostream>
#include<iomanip>
using namespace std;

int i,j;
double f(double x)
{
	return x*x*x*x*x-15*x*x*x*x+85*x*x*x-225*x*x+274*x-121;
}
void qiugen(int l,int r)//在l,r之间求根 
{	
	int mid=(l+r)/2;
	double mid1=mid*1.0/100000000;//得到实际的小数 
	
	if(r-l<10)
	{
		cout<<fixed<<setprecision(6)<<mid1<<endl;
		return;
	}
	else	
	{
		double s=f(mid1);
		if(s==0)
			{
				cout<<fixed<<setprecision(6)<<mid1<<endl;
				return;
			}
		else if(s<0)//更新二分查找右边界 
		{
			qiugen(l,mid);	
		}
		else//更新二分查找左边界
		{
			qiugen(mid,r);
		}
	}	
}
int main()
{
	
	qiugen(150000000,240000000);//放大100000000倍 
	return 0;
}
  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
; end else begin if (pc_src) begin pc <= jump_target; end else begin pc <= pc + 4; end end end // Connect ALU inputs always @ (posedge clk)二分法是一种常用的函数零点的方法。具体实现如下: 1. 确定一个 begin if (reset) begin alu_in1 <= 32'h00000000; alu_in2 <= 区间 [a, b],使得 f(a) 和 f(b) 异号,即 f(a) * f(b) <32'h00000000; end else begin if (opcode == 6'b000010) begin alu 0。 2. 然后取区间的中点 c = (a + b) / 2,计算 f(c)_in1 <= pc + 4; alu_in2 <= 32'h00000000; end else begin 的值。 3. 如果 f(c) 等于 0,则 c 就是函数零点;否则,根据 f(c) 与 f(a) 或 f(b) 的符号确定新的区间: - 若 f(c) 与 f(a alu_in1 <= reg_file[rs]; alu_in2 <= immediate; end end end // Connect data memory inputs always @ (posedge clk) begin if (reset) begin alu_in2 <= 32'h000) 同号,则零点在区间 [c, b] 中,令 a = c; - 若 f(c) 00000; end else begin alu_in2 <= reg_file[rt]; end end // Connect register与 f(b) 同号,则零点在区间 [a, c] 中,令 b = c。 4. 重复 file output to ALU input always @ (posedge clk) begin if (reset) begin alu_in1 <=步骤 2 和 3,直到满足终止条件为止,比如区间长度小于某个阈值 32'h00000000; end else begin alu_in1 <= reg_file[rs]; end end或者迭代次数超过某个值。 需要注意的是,二分法只能用于单峰函数或者 // Connect register file output to data memory input always @ (posedge clk) begin if (reset) begin al单调函数零点解。此外,使用二分法解时需要注意避免浮点数误差问题。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

mengdicfm

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值