雅克比迭代+高斯迭代+SOR迭代法Matlab程序

雅克比迭代+高斯迭代+SOR迭代法Matlab程序，同时支持谱半径计算，利于直接比较三种算法

// a.cpp : 定义控制台应用程序的入口点。
//

#include “stdafx.h“
#include
#include
int main（）
{
	double a[10][10] = {
		{ 4 -1 0 0 0 0 0 0 0 0 }
		{ -1 4 -1 0 0 0 0 0 0 0 }
		{ 0 -1 4 -1 0 0 0 0 0 0 }
		{ 0 0 -1 4 -1 0 0 0 0 0 }
		{ 0 0 0 -1 4 -1 0 0 0 0 }
		{ 0 0 0 0 -1 4 -1 0 0 0 }
		{ 0 0 0 0 0 -1 4 -1 0 0 }
		{ 0 0 0 0 0 0 -1 4 -1 0 }
		{ 0 0 0 0 0 0 0 -1 4 -1 }
		{ 0 0 0 0 0 0 0 0 -1 4 }
	};
	double b[10] = { 75-13 2 6 -1214 -4 5 - 5 };
	double x0[10] = { 0 };
	double x[10] = { 0 };
	int n = 50;//迭代次数 
	double zj = 0;//存储中间变量 
	double r;//存储每一次的迭代误差 
	for （int i = 0; i	{
		for （int m = 0; m<10; m++）
			x0[m] = x[m];

		for （int j = 0; j<10; j++）
		{
			for （int l = 0; l<10; l++）
			zj = zj + a[j][l] * x0[l];

			zj = zj - a[j][j] * x0[j];
			x[j] = （b[j] - zj ）/ a[j][j];
			zj = 0;
		}
	}

	printf（“\nx[]=\n“）;
	for （int i = 0; i<10; i++）
		printf（“\t%f\t\n\n“ x0[i]）;

	system（“pause“）;

	return 0;
}

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1125  2017-04-02 16:47  c语言雅各比.cpp
     文件         533  2017-04-02 21:34  gs.m
     文件         533  2017-04-02 21:33  j.m
     文件         622  2017-04-02 22:55  sor.m