本文主要是介绍计算方法/数值分析牛顿下山法C/C++实现方法,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!
一、前言(吐槽)
我真的不知道该说些什么好,其实也怪自己。老师将实验发下来,然后我就跟这着老师的实验任务做,实验word上用的什么流程图我就用那个流程图,然而我根据老师的流程图遍完了程序,怎么验算结果都不对,就开始找自己的问题,根本没去想流程图,后来才发现,是自己太年轻了,老师给的流程图错了,其实都怨自己,一味的去抄这个流程图,没有去分析这个流程图,这次狠狠的给我上了一课,浪费了好多时间。
二、算法流程图
具体流程图哪里错了,我就不告诉大家了,下面是错误的流程图,大家自己分析一下吧。挺简单的错误。
三、代码
#include <iostream>
using namespace std;
int main(){bool flag = false;//用来判断是第几层循环double X0,X1,r;cout<<"请输入初值:"<<endl;cin>>X0;double e;cout<<"请输入误差限:"<<endl;cin>>e;int m;cout<<"请输入下山最大次数:"<<endl;cin>>m;int n;cout<<"请输入迭代最大次数:"<<endl;cin>>n;int k=0;int i=0;while(true){//说明下山成功进行迭代if(r !=1){cout<<"第"<<i<<"步下山因子:"<<r<<endl;}r = 1.0;if(3*X0*X0-1 != 0){while (true){if(flag){flag = false;break;}X1 = X0 - r*(X0*X0*X0 - X0 - 1) / (3.0 * X0*X0 - 1);if(abs(X1*X1*X1 - X1 - 1) > abs(X0*X0*X0 - X0 - 1)){cout<<"第"<<i<<"步下山因子:"<<r<<endl;cout<<"--------------"<<endl;i=i+1;r=r*0.5;if(i>=m){cout<<"重新输入x0"<<endl;return 0;}elsecontinue;}else{if (abs(X1 - X0) <e){cout<<"近似根为:"<<X1<<endl;return 0;}else{k=k+1;X0=X1;if(k>=n){cout<<"迭代失败!"<<endl;return 0;}else{flag = true;continue;}}}}} else{cout<<"奇异标志"<<endl;return 0;}}}
这篇关于计算方法/数值分析牛顿下山法C/C++实现方法的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!