当前位置：文档库 › C++二元非线性方程组

C++二元非线性方程组

二元非线性方程牛顿迭代法C++代码。

小弟需要解一个二元非线性方程组，虽然说matlab的解法是一两个代码的事，但是我的情况是需要在C++ 里面实现。由于自己不熟悉C++代码和牛顿迭代法的解法，小弟最近在网上找了一些资料。发现很多有问题的例子，纯粹是在误导别人，浪费别人的时间。于是小弟自己在看了解法（牛顿迭代法）的基础上，自己通过C++ 语言实现了一些简单方程组的解法，分享给给大家。

1. 非线性方程组的解法原理

我在网上找了一个资料，应该是中南大学计算机学院的教程，如下：

注意，他这里的例子第二个方程应该是，x1*x2*x2+x1-10*x2+8

经过观察，我发现其实对于2元方程来说：f(x1,x2)和g(x1,x2)。牛顿迭代法里面的迭代量(Δx 1 和Δx 2) 可以如下表示。(我看到网上很多例子用矩阵来表示，我表示看不太懂，而且根据网上贴出来的程序，发现有错误。而我发现对于2元方程来说，用下面的表示方法更直观，更简单。2元方程并不需要再学习矩阵的相关知识。)

()()()2

21121211,,,x x x f x x x f x x f x ??+??-=? ()()()2

21121212,x x x x g x ?+?-=?

2. 代码(运行环境, VS 2010, WIN 7 64)

// Nonlinear.cpp : Defines the entry point for the console application.

#include"stdafx.h"

#include

using namespace std;

int _tmain(int argc, _TCHAR* argv[])

{

double x1,x2; //变量x1,x2;

double x1_0,x2_0; //变量x1,x2迭代前的值

double x1_1,x2_1; //变量x1,x2迭代后的值

double delta_x1,delta_x2; //变量

double function_F,function_G; //方程f(x1,x2) 和g(x1,x2)

double F_x1,F_x2; //方程f的偏导

double G_x1,G_x2; //方程g的偏导

double Error=1e-10; //给定的误差

double Error_cout; //用来计算的误差

int Steps; //运行步骤

int Max_Steps=100000; //最大运行步骤

x1 = 0.5;

x2 = 0.5; //初始化初值

x1_0 = x1;

x2_0 = x2;

for(Steps=1;Steps<=Max_Steps;Steps++)

{

// 方程组, f(x1,x2), g(x1,x2)

function_F = x1_0*x1_0-10.0*x1_0+x2_0*x2_0+8.0; // 10.0的意思是double类型，10 是指整数类型

function_G = x1_0*x2_0*x2_0+x1_0-10.0*x2_0+8.0;

// 偏导,

F_x1 = 2.0*x1_0-10.0;

F_x2 = 2.0*x2_0;

G_x1 = x1_0*x1_0+1.0;

G_x2 = 2.0*x1_0*x2_0-10.0;

// 迭代

delta_x1 = -function_F/(F_x1+F_x2);

delta_x2 = -function_G/(G_x1+G_x2);

x1_1 = x1_0 + delta_x1;

x2_1 = x2_0 + delta_x2;

x1_0 = x1_1; //迭代，将新计算的变量赋值给原始变量

x2_0 = x2_1;

Error_cout = fabs(delta_x1)+fabs(delta_x2);

if(Error_cout

printf(" Run steps is %d\n",Steps);

printf(" x1 is %lf\n",x1_0);

printf(" x2 is %lf\n",x2_0);

return 0;

}

运行结果为：

当然，对于有些方程用这个方法是解不了的，大家可以查一下阻尼牛顿法！

二元一次方程组专项练习及答案

《二元一次方程组》专项练习及答案 §8.1二元一次方程组一、填空题 1、二元一次方程4x-3y=12，当x=0，1，2，3时，y=____ 2、在x+3y=3中，若用x 表示y ，则y= ，用y 表示x ，则x= 3、已知方程(k 2-1)x 2+(k+1)x+(k-7)y=k+2，当k=______时，方程为一元一次方程；当k=______ 时，方程为二元一次方程。 4、对二元一次方程2(5-x)-3(y-2)=10，当x=0时，则y=____；当y=0时，则x=____。 5、方程2x+y=5的正整数解是______。 6、若(4x-3)2+|2y+1|=0，则x+2=。 7、方程组???==+b xy a y x 的一个解为???==3 2y x ，那么这个方程组的另一个解是。 8、若21=x 时，关于y x 、的二元一次方程组? ??=-=-212by x y ax 的解互为倒数，则=-b a 2。二、选择题 1、方程２ｘ－３ｙ＝５，ｘｙ＝3，33=+y x ，３ｘ－ｙ+２ｚ＝０，62=+y x 中是二元一次方程的有（）个。Ａ、１Ｂ、２Ｃ、３Ｄ、４ 2、方程2x+y=9在正整数范围内的解有（） A 、1个 B 、2个 C 、3个 D 、4个 3、与已知二元一次方程5x-y=2组成的方程组有无数多个解的方程是（） A 、10x+2y=4 B 、4x-y=7 C 、20x-4y=3 D 、15x-3y=6

4、若是m y x 25与2214-++n m n y x 同类项，则n m -2的值为（） A 、1 B 、－1 C 、－3 D 、以上答案都不对 5、在方程(k 2-4)x 2+(2-3k)x+(k+1)y+3k=0中，若此方程为二元一次方程，则k 值为（） A 、2 B 、-2 C 、2或-2 D 、以上答案都不对． 6、若???-==1 2y x 是二元一次方程组的解，则这个方程组是（） A 、?? ?=+=-5253y x y x B 、???=--=523x y x y C 、???=+=-152y x y x D 、???+==132y x y x 7、在方程3)(3)(2=--+x y y x 中，用含x 的代数式表示y ，则（） A 、35-=x y B 、3--=x y C 、35+=x y D 、35--=x y 8、已知ｘ＝３－ｋ，ｙ＝ｋ＋２，则ｙ与ｘ的关系是（）Ａ、ｘ＋ｙ＝５Ｂ、ｘ＋ｙ＝１Ｃ、ｘ－ｙ＝１Ｄ、ｙ＝ｘ－１ 9、下列说法正确的是（）Ａ、二元一次方程只有一个解Ｂ、二元一次方程组有无数个解Ｃ、二元一次方程组的解必是它所含的二元一次方程的解Ｄ、三元一次方程组一定由三个三元一次方程组成 10、若方程组???=+=+16 156653y x y x 的解也是方程３ｘ＋ｋｙ=10的解，则ｋ的值是（ =）Ａ、ｋ＝６ = Ｂ、ｋ＝１０Ｃ、ｋ＝９Ｄ、ｋ＝ 10 1 三、解答题 1、解关于x 的方程)1(2)4)(1(+-=--x a x a a

二元一次方程组计算题50道(答案)

.. 中考真题 50 道中考真题之《二元一次方程组计算题》 -----专项练习50题（有答案） 1．（2012?德州）已知，则a+b 等于（） A. 3 B C. 2 D. 1 2．（2012菏泽）已知???==1 2 y x 是二元一次方程组81mx ny nx my +=??-=?的解，则n m -2的算术平方根为（） A ．±2 B ． 2 C ．2 D ． 4 3．（2012临沂）关于x 、y 的方程组3, x y m x my n -=?? +=?的解是1,1,x y =??=? 则m n -的值是（） A ．5 B ．3 C ．2 D ．1 4．（2012?杭州）已知关于x ，y 的方程组，其中﹣3≤a ≤1，给出下列结论： ①是方程组的解； ②当a=﹣2时，x ，y 的值互为相反数； ③当a=1时，方程组的解也是方程x+y=4﹣a 的解； ④若x ≤1，则1≤y ≤4．其中正确的是（） A ．①② B ．②③ C ．②③④ D ．①③④ 5. （2012广东湛江）请写出一个二元一次方程组，使它的解是． 6．（2012广东）若x ，y 为实数，且满足|x ﹣3|+ =0，则（）2012的值是 1 ．

7．（2012安顺）以方程组的解为坐标的点（x ，y ）在第象限． 8．(2012?连云港)方程组的解为． 9.（2012?广州）解方程组． 10．（2012广东）解方程组：． 11.（2012?黔东南州）解方程组． 12、（2012湖南常德）解方程组：???==+1-25y x y x 13. （2011湖南益阳，2，4分）二元一次方程21-=x y 有无数多个解，下列四组值中不是．．该方程的解的是 A ．0 12 x y =???=-?? B ．11x y =??=? C ．1 0x y =??=? D ．11x y =-??=-? 14. （2011四川凉山州，3，4分）下列方程组中是二元一次方程组的是（） A ．12xy x y =??+=? B ． 523 13x y y x -=???+=?? C ． 20 135x z x y +=?? ? -=?? D ．5723 z x y =???+=?? 15. （2011广东肇庆，4，3分）方程组?? ?=+=-4 22 y x y x 的解是 ① ②

二元一次方程组经典练习题+答案解析100道(1)

二元一次方程组练习题100道（卷一）（范围：代数：二元一次方程组）一、判断 1、??? ??-==312y x 是方程组?????? ?=-=-9 1032 6 5 23y x y x 的解 …………（） 2、方程组?? ?=+-=5 231y x x y 的解是方程3x -2y =13的一个解（） 3、由两个二元一次方程组成方程组一定是二元一次方程组（） 4、方程组???????=-++=+++2 5323 473 5 23y x y x ，可以转化为???-=--=+27651223y x y x （） 5、若(a 2-1)x 2 +(a -1)x +(2a -3)y =0是二元一次方程，则a 的值为±1（） 6、若x +y =0，且|x |=2，则y 的值为2 …………（） 7、方程组? ? ?=+-=+81043y x x m my mx 有唯一的解，那么m 的值为m ≠-5 …………（） 8、方程组?? ???=+=+62 3 131 y x y x 有无数多个解 …………（） 9、x +y =5且x ，y 的绝对值都小于5的整数解共有5组 …………（） 10、方程组?? ?=+=-3513y x y x 的解是方程x +5y =3的解，反过来方程x +5y =3的解也是方程组? ??=+=-351 3y x y x 的解 ………（） 11、若|a +5|=5，a +b =1则3 2 -的值为b a ………（） 12、在方程4x -3y =7里，如果用x 的代数式表示y ，则4 37y x += （）二、选择： 13、任何一个二元一次方程都有（）（A ）一个解；（B ）两个解；（C ）三个解；（D ）无数多个解； 14、一个两位数，它的个位数字与十位数字之和为6，那么符合条件的两位数的个数有（）（A ）5个（B ）6个（C ）7个（D ）8个 15、如果? ? ?=+=-423y x a y x 的解都是正数，那么a 的取值范围是（）（A ）a <2；（B ）34- >a ；（C ）3 4 2<<-a ；（D ）3 4 -