2012年03月_zhangchao3322218

11月 10月 08月 03月

原创 Romberg积分法计算定积分（C语言）

/*-------------Romberg积分法计算定积分-------参考教材：《数值分析》李乃成，梅立泉，科学出版社《计算方法教程》第二版凌永祥，陈明逵*///误差界eps%被积函数为f(x)=(x^3+sin(x))/x;积分区间为[0.3,0.8]#include#include#includeint main(void){/* int

2012-03-30 18:48:49 6557 1

原创 Romberg积分法计算定积分（Matlab程序）

%Romberg积分法计算定积分%参考教材：《数值分析》李乃成，梅立泉，科学出版社%《计算方法教程》第二版凌永祥，陈明逵clear;clc;close all;format long% %被积函数为f(x)=4/(1+x^2);积分区间为[0,1]% b=1;a=0;h=b-a;eps=10^(-5);%误差界eps%被积函数为f(x)=(x^3+sin(x))/x;

2012-03-30 18:48:05 25045 2

原创函数插值计算（Matlab程序）

%Runge函数插值计算%参考教材：《数值分析》李乃成，梅立泉，科学出版社%《计算方法教程》第二版凌永祥，陈明逵%Runge函数f(x)=(1+25x^2)^(-1),(-1clear;clc;close all;format short;fprintf('输入N：(5,10,15,20)\n');N=input('');Ni=0:N;%format sh

2012-03-30 18:47:22 4397

原创对阵正定矩阵的楚列斯基（Cholesky）分解（C语言）

//定理2.2.3：对阵正定矩阵的楚列斯基（Cholesky）分解//设A为n阶对阵正定矩阵，则存在一个可逆的下三角矩阵G，使得//A=GG’，当限定G的对角元为正时，这种分解是唯一的////--------A=GG’的分解算法-------//参考教材：《数值分析》李乃成，梅立泉，科学出版社// 《计算方法教程》第二版凌永祥，陈明逵#include#inc

2012-03-30 18:46:31 8347

原创对称矩阵的LDLT分解（C语言）

/*对称矩阵的LDLT分解定理2.2.2：对称矩阵的三角分解：设A是n阶对称矩阵，若A的各阶顺序主子式均不等于0，则A可以唯一地分解为 A=LDL’-------------A=LDL’的分解算法-------参考教材：《数值分析》李乃成，梅立泉，科学出版社《计算方法教程》第二版凌永祥，陈明逵*/#include#includeint min(

2012-03-30 18:45:56 6455 1

原创高斯—赛德尔迭代法解线性方程组（C语言）

/*-------------高斯—赛德尔迭代法解线性方程组---------参考教材：《数值分析》李乃成，梅立泉，科学出版社《计算方法教程》第二版凌永祥，陈明逵*/#include#includeint main(void){ double A[5][5] = {{28,-3,0,0,0}, {-3,38,-10,0,-5}, {0

2012-03-30 18:45:09 32072

原创 Jacobi迭代法解线性方程组（C语言）

/*-------------Jacobi迭代法解线性方程组---------参考教材：《数值分析》李乃成，梅立泉，科学出版社《计算方法教程》第二版凌永祥，陈明逵*/#include#includeint main(void){ double A[5][5] = {{28,-3,0,0,0}, {-3,38,-10,0,-5}, {0

2012-03-30 18:44:32 22817

原创求解三对角方程组的追赶法（Matlab程序）

clear all;clc;fprintf('输入n：(10,20,30)\n');n=input('');na=zeros(1,n);b=zeros(1,n);c=zeros(1,n-1);A=zeros(n,n);for i=1:n for j=1:n if j==i A(i,j)=2; elsei

2012-03-30 18:43:49 22999

原创三对角阵的LU分解和三对角方程组的求解（C语言）

/*三对角阵的LU分解和三对角方程组的求解-------------A=LU的分解算法-------参考教材：《数值分析》李乃成，梅立泉，科学出版社《计算方法教程》第二版凌永祥，陈明逵*/#include#include#includeint main(void){ int i,j,n; int N; printf("请输入 N(10,

2012-03-30 18:42:32 7281

原创列主元Gauss消去法解线性方程组（C语言）

//Gauss消去法解线性方程组//参考教材《计算方法教程》第二版，西安交通大学出版社#include#includeint main(void){ float A[7][7]={{3,-5,6,4,-2,-3,8}, {1,1,-9,15,1,-9 ,2}, {2,-1,7,5,-1,6,11}, {-1,1,3,2,7,-1,-2},

2012-03-30 18:41:20 7891

原创平方根法和改进的平方根法解线性方程组（Matlab程序）

%定理2.2.3：对阵正定矩阵的楚列斯基（Cholesky）分解%设A为n阶对阵正定矩阵，则存在一个可逆的下三角矩阵G，使得%A=GG’，当限定G的对角元为正时，这种分解是唯一的%--------A=GG’的分解算法-------%参考教材：《数值分析》李乃成，梅立泉，科学出版社clear;clc;A=[9,18,9,-27; 18,45,0,-45; 9,0,

2012-03-30 18:39:52 23781 2

原创共轭梯度法解线性方程组（Matlab程序）

%-------共轭梯度法解线性方程组-----------%---Conjugate Gradient method-------%参考教材《数值分析》李乃成&梅立泉，科学出版社2011clear;clc;% A=[10,-1,-2;-1,10,-2;-1,-1,5];% b=[72,83,42]';A=[2 0 1;0 1 0;1 0 2];b=[3 1 3]';

2012-03-30 18:38:38 29272 4

原创逐次超松弛迭代法解线性方程组（Matlab程序）

%---逐次超松弛迭代法-----%---successive over-reaxation iteration methodclear;clc;A=[10,-1,-2;-1,10,-2;-1,-1,5];b=[72,83,42]';N=length(b); %解向量的维数fprintf('库函数计算结果：');x=inv(A)*b %库函数计算结果

2012-03-30 18:38:03 18489

原创高斯—赛德尔(Gauss - Seidel)迭代法解线性方程组（Matlab程序）

%---高斯—赛德尔迭代法-----%---Gauss - Seidel iteration methodclear;clc;% A=[10,-1,-2;-1,10,-2;-1,-1,5];% b=[72,83,42]';A=[ 28,-3,0,0,0; -3,38,-10,0,-5; -10,0,25,-15,0; 0,0,-15,45,0;

2012-03-30 18:37:28 34865

原创雅克比(Jacobi)迭代法解线性方程组（Matlab程序）

%---雅克比迭代法-----%---Jacobi iteration methodclear;clc;% A=[10,-1,-2;-1,10,-2;-1,-1,5];% b=[72,83,42]';A=[ 28,-3,0,0,0; -3,38,-10,0,-5; -10,0,25,-15,0; 0,0,-15,45,0; 0,-5,0,

2012-03-30 18:36:30 59238 1

原创矩阵的LU分解（Matlab程序）

%LU分解，Doolittle分解%L为单位下三角阵，U为上三角阵clear;clc;A=[3 -5 6 4 -2 -3 8; 1 1 -9 15 1 -9 2; 2 -1 7 5 -1 6 11; -1 1 3 2 7 -1 -2; 4 3 1 -7 2 1 1; 2 9 -8 11 -1 -4 -1;

2012-03-30 18:34:43 19589