解二元一次方程组无解组

来源：蜘蛛抓取(WebSpider) 时间：2011-08-08 10:07 标签： 二元一次方程组无解
您当前的位置：&>&
复制带链接的文本推荐到其它网站→
将本页工具分享到：
在线一元二次方程在线计算器
解二次方程式计算
一元二次方程的解法
&script language=JavaScript>
复制带链接的文本推荐其它网站→
上一篇： 下一篇：
如有问题可到本站-
&看不清楚,点击刷新
学习工程技术的天地
您感兴趣的栏目9887人阅读
本科课程参见：《》
牛顿迭代公式
设已知方程f(x)=0的近似根x0 ,则在x0附近f(x)可用一阶泰勒多项式近似代替.因此, 方程f(x)=0可近似地表示为p(x)=0。用x1表示p(x)=0的根,它与f(x)=0的根差异不大.&
设 ,由于x1满足解得
重复这一过程,得到迭代公式：
这就是著名的牛顿迭代公式,它相应的不动点方程为
Jacobi迭代公式解线性方程组
线性方程组基本解法：
若方程组可同解变形为
Jacobi迭代法的计算公式：
/*简单迭代法的代码实现*/
#include&iostream&
#include&string&
#include&cmath&
double e=2.;
double f(double x){
return pow(e,-1*x);
void SimpleDiedai(double x,double d){
double a=x;
double b=f(a);
int k=0;//记录循环的次数
while(((a-b)&d) || ((a-b)&-1*d)){
if(k&100){
cout&&&迭代失败！（可能是函数不收敛）&&&
int main(){
cout&&&请输入初始值x0和要求得结果的精度：&;
double x,d;
cin&&x&&d;
SimpleDiedai (x,d);
/*牛顿迭代法的代码实现*/
#include&iostream&
#include&string&
#include&cmath&
double e=2.;
double f(double x){
double a=pow(e,-1*x);
return x-(x-a)/(1+a);
void NewtonDiedai(double x,double d){
double a=x;
double b=f(a);
int k=0; //记录循环的次数
while(((a-b)&d) || ((a-b)&-1*d)){
if(k&100){
cout&&&迭代失败！（可能是函数不收敛）&&&
int main(){
cout&&&请输入初始值x0和要求得结果的精度：&;
double x,d;
cin&&x&&d;
NewtonDiedai(x,d);
/*雅可比算法的代码实现*/
#include&iostream&
#include&iomanip&
#include&string&
#include&vector&
//函数求数组中的最大值
double MaxOfList(vector&double&x){
double max=x[0];
int n=x.size();
for(int i=0;i&n;i++)
if(x[i]&max) max=x[i];
//雅可比迭代公式
void Jacobi(vector&vector&double& & A,vector&double& B,int n){
vector&double& X(n,0);
vector&double& Y(n,0);
vector&double& D(n,0);
int k=0; //记录循环次数
for(int i=0;i&n;i++){
double tem=0;
for(int j=0;j&n;j++){
if(i!=j) tem += A[i][j]*X[j];
Y[i]=(B[i]-tem)/A[i][i];
cout&&left&&setw(8)&&Y[i]&&& &;
if(k&100){
cout&&&迭代失败！（可能是函数不收敛）&&&
for(int a=0;a&n;a++){
D[a]=X[a]-Y[a];
}while( MaxOfList(D)&0.00001 || MaxOfList(D)&-0.00001);
int main(){
cout&&&请输入方程组未知数的个数n：&;
vector&vector&double& &A(n,vector&double&(n,0));
vector&double&B(n,0);
cout&&&请输入方程组的系数矩阵：&&&
for(int i=0;i&n;i++){
for(int j=0;j&n;j++){
cin&&A[i][j];
cout&&&请输入方程组的值向量：&&&
for(int k=0;k&n;k++){
cin&&B[k];
cout&&&您输入的方程组为：&&&
for(int a=0;a&n;a++){
for(int b=0;b&n;b++){
cout&&A[a][b]&&& &;
cout&&&由雅可比迭代公式求的方程组的解为：&&&
Jacobi(A,B,n);
实验过程原始记录
（1）分别用简单迭代法和牛顿迭代法求方程在x=0.5附近的一个根x*，要求精度为0.00001
（输入0.5 0.000001）简单迭代法得到结果：
（输入0.5 0.000001）牛顿迭代法得到结果：
X0=0.5 &x1=0.=0.567143
（2）用雅可比迭代法求解方程组 &
运行程序，根据提示输入 （3） （10 -1 -2 -1 10 -2 -1 -2 5） & &（7.2 &8.3 &4.2）
实验结果及分析
1、迭代法是一种逐次逼近法，这种方法使用某个固定公式－所谓迭代公式反复校正根的近似值，使之逐步精确化，直至满足精度要求的结果。迭代法是一种求解函数方程f（x）=0的解的方法，他解决了二分法无法求解复根级偶重根的问题，而其提高了收敛速度。迭代的思想是计算方法中基础的求解问题的思想。
2、简单迭代法的求根过程分成两步，第一步先提供根的某个猜测值，即所谓迭代初值，然后将迭代初值逐步加工成满足精度要求的根。迭代法的设计思想是：f (x) = 0等价变换成 然后由迭代公式 逐步球的满足精度的解。实际迭代中不同迭代函数的求解可能影响求的精确解的运算量，甚至可能因为函数发散而无法求解。解题时可通过对导函数的判断而判断函数是否发散，而编写代码时可以通过判断循环次数——即循环过多次而不能从循环中出来时就判断为死循环，无法求得正解
3、简单迭代法往往只是线性收敛，为得出超线性收敛的迭代格式，通常采用近似替代法，&即牛顿公式。迭代函数为 &- / 牛顿法是一种逐步线性化方法。由实验结果可以看到，虽然选取近似公式，但牛顿迭代法仍能得到精度很高的解，而且牛顿迭代法大大提高了收敛速度。
4、由迭代法求解线性方程组的基本思想是将联立方程组的求解归结为重复计算一组彼此独立的线性表达式，这就使问题得到了简化，类似简单迭代法转换方程组中每个方程式可得到雅可比迭代式
迭代法求解方程组有一定的局限性，例如要求方程组的系数矩阵具有某种特殊性质，以保证迭代过程的收敛性，但迭代法同时有十分明显的优点——算法简单，因而编制程序比较容易，所以在实际求解问题中仍有非常大利用价值。
（转载请注明作者和出处：&未经允许请勿用于商业用途）
* 以上用户言论只代表其个人观点，不代表CSDN网站的观点或立场
访问：1009663次
积分：14518
积分：14518
排名：第259名
原创：141篇
评论：1903条
文章：20篇
阅读：76846
文章：14篇
阅读：83856
文章：20篇
阅读：101979
文章：30篇
阅读：446979
CS小硕一枚，兴趣领域：计算机视觉、媒体分析、图像处理；认知心理学，生理心理学。
喜欢C++、Matlab、OpenCV、Qt、UML、VS、XMind等各种优秀的、提高码农生产力的工具。
我的邮箱weilan.非齐次方程组解的结构及其判定_百度文库
两大类热门资源免费畅读
续费一年阅读会员，立省24元！
文档贡献者贡献于
评价文档：
46页免费26页免费17页免费17页免费14页免费55页免费31页免费38页免费30页免费29页免费
喜欢此文档的还喜欢6页免费26页免费60页免费29页1下载券7页免费
非齐次方程组解的结构及其判定|数​学​课​件
把文档贴到Blog、BBS或个人站等：
普通尺寸(450*500pix)
较大尺寸(630*500pix)
大小：308.50KB
登录百度文库，专享文档复制特权，财富值每天免费拿！
你可能喜欢您当前的位置：&>&
复制带链接的文本推荐到其它网站→
将本页工具分享到：
&解方程计算器
例如，输入 a=1, b=8, c=16 和 d=10.点击解三次方程
三次方程公式:
判别(&D) = q3 + r2
q = (3c- b2)/9
r = -27d + b(9c-2b2)
s = r + math.sqrt(判别)
t = r - math.sqrt(判别)
term1 = math.sqrt(3.0)*((-t + s)/2)
r13= 2 * math.sqrt(q)
x1=(- term1 + r13*math.cos(q3/3) )
x2=(- term1 + r13*math.cos(q3+(2*math.PI)/3) )
x3=(- term1 + r13*math.cos(q3+(4*math.PI)/3) )
这个计算器可以帮助你动态地计算三次方程的根。易于理解和使用。
复制带链接的文本推荐其它网站→
上一篇： 下一篇：
如有问题可到本站-
&看不清楚,点击刷新
学习工程技术的天地
您感兴趣的栏目

我要回帖

更多关于 二元一次方程组无解 的文章

 

随机推荐