数值分析
麦克波菲柯特
人需要先说很多话,然后保持静默
展开
-
牛顿插值
牛顿插值计算f(x)=√x 完整代码#include<stdio.h> //牛顿插值#include<math.h>int main(){ int i,j,n; double a,x[100],f[100][100],xa,ff=0; printf("数据个数:");scanf("%d",&n); for(i=0;i<=n;i++) {原创 2017-04-29 10:26:36 · 711 阅读 · 0 评论 -
复化梯形公式
复化梯形公式重在理解梯形公式的概念,计算f(X)=#include<stdio.h> //梯形公式#include<math.h>//#define f(x) (sin(x)/x)double f(double x){ return x>0? (sin(x)/x):1;}int paw(int x,int y){ int i,a=1; if(y原创 2017-04-29 10:38:54 · 8466 阅读 · 0 评论 -
龙贝格算法
龙贝格算法是以复化梯形公式,复化希普森,复化牛顿-科特斯为基础采用的是逐渐平衡误差的做法,如下:代码#include<stdio.h> //龙贝格#include<math.h>#define maxlen 100double a;double b;double f(double x){ if(x==0) { return 1; }原创 2017-04-29 10:56:58 · 5358 阅读 · 0 评论 -
迭代法
迭代法又叫做逐次逼近法,这个方法比较的简单,就是不断地利用|x-g(x)|去求解 代码段#include<stdio.h> //迭代法、#include<math.h>#define g(x) (pow((x+1),(1.0/3)))int main(){ int k,n=50; double x0,x; printf("x0=");scanf(原创 2017-04-29 11:01:39 · 1033 阅读 · 0 评论 -
牛顿下山法
牛顿法来源于泰勒展开式,牛顿法的几何意义为斜率不断地接近零点,所以简单的代码如下#include<stdio.h> //牛顿下山#include<math.h>#define f(x) (pow(x,3)-x-1)#define ff(x) (pow((x+1),(1.0/3)))int main(){ int i,k,n=500,m=50; double原创 2017-04-29 11:05:09 · 2780 阅读 · 0 评论 -
高斯消去法
高斯消去法主要是理解其内在的原理,这个思想到底是怎么产生的呢,我看还是和高斯小时候的想法有关,相传高斯小时候,老师给他们出了一道题,题目是:1+2+3+4………..+100的和,老师原以为学生们会计算很久,但是高斯思考片刻就解决了,高斯的解决思路是1+99、2+98、3+97、4+96……..现在看来觉得很简单,但是在当时看来这可就是一项伟大的思想,于是…… 代码#include<stdio.h>原创 2017-04-29 11:21:50 · 1035 阅读 · 0 评论 -
高斯约旦法
高斯约旦法的想法就是高斯消去法的进一步扩展,理解起来不算特别困难,但是就是在构造的过程中有些困难。 代码#include<stdio.h> //高斯约旦法#include<math.h>#define max 10void gaosiyuedan(int n,double a[max][max],double x[max]){ int i,j,k,m; doub原创 2017-04-29 11:25:50 · 1962 阅读 · 0 评论 -
欧拉、改进欧拉、四阶龙格-库塔
取h=0.1时用Euler法,改进Euler法,Rung-Kutta方法求其数值解并与精确解进行比较。 输入:求解区间,初值,数值解个数 输出:数值解 代码#include<stdio.h>#include<math.h>#define f(y,x) (y-((2*x)/y))#define F(y,x) (y-(2*x/y))#define ff(y,x) (y-(2*x/y))#原创 2017-04-29 11:32:11 · 5065 阅读 · 0 评论 -
拉格朗日插值
数值分析是一门十分有用的课程,趁着假期,有幸来整理一下 拉格朗日插值,计算f(x)=1/(1+x²)的插值 代码#include<stdio.h> //拉格朗日插值int main(){ int i,j,a,b,n; double xa,h,x[100],f[100],ff=0,s; printf("n="); scanf("%d",&n); printf原创 2017-04-29 10:15:43 · 1123 阅读 · 0 评论