🌟 前言
Wassup guys,我是Edison 😎
今天是C语言每日一练,第154天!
Let’s get it!

文章目录
编写用牛顿迭代法求方程根的函数。
方程为 a x 2 + b x 2 + c x + d = 0 ax^2+bx^2+cx+d=0 ax2+bx2+cx+d=0,系数a,b,c,d 由主函数输入。
求 x x x 在 1 1 1 附近的一个实根。求出根后,由主函数输出。
牛顿迭代法的公式是: − f ( x 0 ) f ′ ( x 0 ) -\frac{f(x_0 )}{f'(x_0)} −f′(x0)f(x0) ,设迭代到 ∣ x − x 0 ∣ ≤ 1 0 − 5 |x-x_0|\leq10^{-5} ∣x−x0∣≤10−5 时结束。
牛顿迭代法是取 x 0 x_0 x0 之后,在这个基础上,找到比 x 0 x_0 x0 更接近的方程的根,一步一步迭代,从而找到更接近方程根的近似根。
设 r r r 是 f ( x ) = 0 f(x)=0 f(x)=0 的根,选取 x 0 x_0 x0 作为 r r r 初始近似值。
过点 ( x 0 , f ( x 0 ) ) (x_0, f(x_0)) (x0,f(x0)) 作为曲线 y = f ( x ) y=f(x) y=f(x) 的切线 L L L,
L L L 的方程为 y = f ( x 0 ) + f ′ ( x 0 ) ∗ ( x − x 0 ) y=f(x_0)+f'(x_0)*(x-x_0) y=f(x0)+f′(x0)∗(x−x0),
求出 L 与 x 轴交点的横坐标 x 1 = x 0 − f ( x 0 ) / f ′ ( x 0 ) x_1=x_0-f(x_0)/f'(x_0) x1=x0−f(x0)/f′(x0),称 x x x 为 r r r 的一次近似值,
过点 ( x 1 , f ( x 1 ) ) (x_1,f(x_1)) (x1,f(x1)) 作为曲线 y = f ( x ) y=f(x) y=f(x) 的切线,并求该切线与 x 轴的横坐标 x 2 = x 1 − f ( x 1 ) / f ′ ( x 1 ) x_2=x_1-f(x_1)/f'(x_1) x2=x1−f(x1)/f′(x1),称 x 2 x_2 x2 为 r r r 的二次近似值,
重复以上过程,得 r r r 的近似值 x n x_n xn。
上述过程即为牛顿迭代法的求解过程。
程序流程分析👇
(1) 在 1 1 1 附近找任一实数作为 x 0 x_0 x0 的初值,我们取 1.5 1.5 1.5,即 x 0 = 1.5 x_0=1.5 x0=1.5
(2) 用初值 x 0 x_0 x0 代入方程中计算此时的 f ( x 0 ) f(x_0) f(x0) 及 f ′ ( x 0 ) f'(x_0) f′(x0);程序中用变量 f f f 描述方程的值,用 f d fd fd 描述方程求导之后的值。
(3) 计算增量 h = f / f d h=f/fd h=f/fd。
(4) 计算下一个 x : x = x 0 − h x:x=x_0-h x:x=x0−h。
(5) 用新产生的 x x x 替换原来的 x o x_o xo,为下一次迭代做好准备。
(6) 若 ∣ x − x 0 ∣ > = 1 e − 5 |x-x_0|>=1e-5 ∣x−x0∣>=1e−5,则转到第 (3) 步继续执行,否则转到步骤 (7)。
(7) 所求 x x x 就是方程 a x 3 + b x 2 + c x + d = 0 ax^3+bx^2+cx+d=0 ax3+bx2+cx+d=0 的根,将其输出。
本程序的编写既可用while,也可用do...while,二者得到的结果是一样的,只是在赋初值时稍有不同。
while结构需要先判定条件,即先判断 ∣ x − x 0 ∣ > = 1 e − 5 |x-x_0|>=1e-5 ∣x−x0∣>=1e−5 是否成立,这样对于 x x x, x 0 x_0 x0 我们要在 1 1 1 附近取两个不同的数值作为初值;
do...while结构是先执行一次循环体,得到 x x x 的新值后再进行判定,这样程序开始只需给 x x x 赋初值。
这里我们采用do...while结构来实现。
程序的主体结构如下👇
由于程序中用到了绝对值函数fabs(), 所以在程序的开始要加上头文件#include <math.h>。
流程图如下所示👇
编写程序时要注意的一点是判定 ∣ x − x 0 ∣ > = 1 e − 5 |x-x_0|>=1e-5 ∣x−x0∣>=1e−5。
从牛顿迭代法的原理可以看出:迭代的实质就是越来越接近方程根的精确值,最初给 x 0 x_0 x0 所赋初值与根的精确值是相差很多了,正是因为这个我们才需要不断地进行迭代,也就是程序中循环体的功能。
在经过一番迭代之后所求得的值之间的差别也越来越小,直到求得的某两个值的差的绝对值在某个范围之内时,便可结束迭代。
若我们把判定条件改为 ∣ x − x 0 ∣ < 1 e − 5 |x-x_0|<1e-5 ∣x−x0∣<1e−5,则第一次的判断结果必为假,这样就不能进入循环体再次执行。
定义 solution()函数求方程的根。solution()函数的代码如下👇
完整代码📝
#include <stdio.h>
#include <math.h>
float solution(float a, float b, float c, float d)
{
float x0, f, fd, h;
float x = 1.5;
do
{
x0 = x;
f = a * x0 * x0 * x0 + b * x0 * x0 + c * x0 + d;
fd = 3 * a * x0 * x0 + 2 * b * x0 + c;
h = f / fd;
x = x0 - h;
} while (fabs(x-x0) >= 1e-5);
return x;
}
int main()
{
float a, b, c, d;
float x;
printf("请输入方程的系数:");
scanf("%f %f %f %f", &a, &b, &c, &d);
x = solution(a, b, c, d);
printf("\n");
printf("所求方程的根为:x=%f\n", x);
return 0;
}
运行结果👇
代码解释👇
我正在寻找执行以下操作的正确语法(在Perl、Shell或Ruby中):#variabletoaccessthedatalinesappendedasafileEND_OF_SCRIPT_MARKERrawdatastartshereanditcontinues. 最佳答案 Perl用__DATA__做这个:#!/usr/bin/perlusestrict;usewarnings;while(){print;}__DATA__Texttoprintgoeshere 关于ruby-如何将脚
我在用Ruby执行简单任务时遇到了一件奇怪的事情。我只想用每个方法迭代字母表,但迭代在执行中先进行:alfawit=("a".."z")puts"That'sanalphabet:\n\n#{alfawit.each{|litera|putslitera}}"这段代码的结果是:(缩写)abc⋮xyzThat'sanalphabet:a..z知道为什么它会这样工作或者我做错了什么吗?提前致谢。 最佳答案 因为您的each调用被插入到在固定字符串之前执行的字符串文字中。此外,each返回一个Enumerable,实际上您甚至打印它。试试
几个月前,我读了一篇关于rubygem的博客文章,它可以通过阅读代码本身来确定编程语言。对于我的生活,我不记得博客或gem的名称。谷歌搜索“ruby编程语言猜测”及其变体也无济于事。有人碰巧知道相关gem的名称吗? 最佳答案 是这个吗:http://github.com/chrislo/sourceclassifier/tree/master 关于ruby-寻找通过阅读代码确定编程语言的rubygem?,我们在StackOverflow上找到一个类似的问题:
为什么以下不同?Time.now.end_of_day==Time.now.end_of_day-0.days#falseTime.now.end_of_day.to_s==Time.now.end_of_day-0.days.to_s#true 最佳答案 因为纳秒数不同:ruby-1.9.2-p180:014>(Time.now.end_of_day-0.days).nsec=>999999000ruby-1.9.2-p180:015>Time.now.end_of_day.nsec=>999999998
?博客主页:https://xiaoy.blog.csdn.net?本文由呆呆敲代码的小Y原创,首发于CSDN??学习专栏推荐:Unity系统学习专栏?游戏制作专栏推荐:游戏制作?Unity实战100例专栏推荐:Unity实战100例教程?欢迎点赞?收藏⭐留言?如有错误敬请指正!?未来很长,值得我们全力奔赴更美好的生活✨------------------❤️分割线❤️-------------------------
嗨~大家好,这里是可莉!今天给大家带来的是7个C语言的经典基础代码~那一起往下看下去把【程序一】打印100到200之间的素数#includeintmain(){ inti; for(i=100;i 【程序二】输出乘法口诀表#includeintmain(){inti;for(i=1;i 【程序三】判断1000年---2000年之间的闰年#includeintmain(){intyear;for(year=1000;year 【程序四】给定两个整形变量的值,将两个值的内容进行交换。这里提供两种方法来进行交换,第一种为创建临时变量来进行交换,第二种是不创建临时变量而直接进行交换。1.创建临时变量来
关闭。这个问题是off-topic.它目前不接受答案。想改进这个问题吗?Updatethequestion所以它是on-topic用于堆栈溢出。关闭11年前。Improvethisquestion我不经常使用ruby-通常它加起来相当于每两个月或更长时间编写一次脚本。我的大部分编程都是使用C++进行的,这与ruby有很大不同。由于我与ruby之间的差距如此之大,我总是忘记语言的基本方面(比如解析文本文件和其他简单的东西)。我想每天练习一些基本的东西,我想知道是否有一些我可以订阅的网站,并且会向我发送当天的Ruby问题或类似的东西。有人知道这样的站点/Internet服务吗?
我最近开始编写Ruby代码,但我对block参数有误解。以下面的代码为例:h={#Ahashthatmapsnumbernamestodigits:one=>1,#The"arrows"showmappings:key=>value:two=>2#ThecolonsindicateSymbolliterals}h[:one]#=>1.Accessavaluebykeyh[:three]=3#Addanewkey/valuepairtothehashh.eachdo|key,value|#Iteratethroughthekey/valuepairsprint"#{value}:#{ke
如果特定语言环境中缺少翻译,如何配置i18n以使用en语言环境翻译?当前已插入翻译缺失消息。我正在使用RoR3.1。 最佳答案 找到相似的question这里是答案:#application.rb#railswillfallbacktoconfig.i18n.default_localetranslationconfig.i18n.fallbacks=true#railswillfallbacktoen,nomatterwhatissetasconfig.i18n.default_localeconfig.i18n.fallback
在我的双语Rails4应用程序中,我有一个像这样的LocalesController:classLocalesController用户可以通过此表单更改其语言环境:deflocale_switcherform_tagurl_for(:controller=>'locales',:action=>'change_locale'),:method=>'get',:id=>'locale_switcher'doselect_tag'set_locale',options_for_select(LANGUAGES,I18n.locale.to_s)end这有效。但是,目前用户无法通过URL更改