计算方法
web小奶精
这个作者很懒,什么都没留下…
展开
-
高斯消去法
#include main(){ FILE *f; double a[15][15]; double b[15],s; int i,j,k,n; f=fopen("Gauss.txt","r"); fscanf(f,"%d",&n); for(i=1;i<=n;i++) { for(j=1;j<=n;j++)原创 2017-11-28 19:47:49 · 407 阅读 · 0 评论 -
龙格-库塔法
#include #include #includefloat f(float x,float y){ return 2*x*y; ///一阶方程}float RK(float x,float y,float h) ///龙格-库塔函数{ float k1,k2,k3,k4; k1=f(x,y); k2=f(x+h/原创 2018-01-01 10:08:12 · 689 阅读 · 0 评论 -
龙贝格求积法
#includeusing namespace std;#include"math.h"#define e 1e-5double f(double x){ double y; if (x==0) return y=1.0;//f(0)=1 else y=sin(x)/x; return y;}void romberg(double原创 2017-12-25 20:04:17 · 934 阅读 · 0 评论 -
Romberg(龙贝格求积)
#include float f(float x){ return(exp(x)/(4+x*x));}main(){ float a=0,b=1;//积分上下界 float h=b-a,T1,T2,s,x;//h是区间长度,T1是分为n等分的近似值,T2是分为2n等分的近似值 int i;//次数 T1=h/2*(f(a)+f(b));//公式原创 2017-12-24 13:18:52 · 2055 阅读 · 1 评论 -
追赶法
文件名:zgf.txt104 -1 7 -1 4 -1 5 -1 4 -1 -13-1 4 -1 2 -1 4 -1 6-1 4 -1 -12 -1 4 -1 14 -1 4 -1 -4 -1 4 -1 5 -1 4 -5#include main(){ FILE *f; double a[15],b[15],c[原创 2017-11-26 20:28:33 · 8402 阅读 · 0 评论 -
超松弛迭代法(SOR)
c.txt4 -2 -4 10-2 17 10 3-4 10 9 -7#include#include#define M 3main(){ FILE *f; double a[M][M],b[M]; double x[M]={0,0,0}; double y[M]; int i,j,k,n=100; f=fopen("c.txt",原创 2017-12-04 20:00:47 · 15370 阅读 · 1 评论 -
迭代法
f(x)=x*x*x-3*x-1 迭代公式:x=(3*x+1)/x*x#include #include #define maxrept 20 //迭代最高次数double f(double x){ return (3*x+1)/x*x; //直接写迭代公式}void main(){ int k=0;//迭代次数 double d原创 2017-11-20 22:44:04 · 342 阅读 · 0 评论 -
Stenffensen加速迭代法
f(x)=x*x*x-3*x-1迭代公式:x=(3*x+1)/(x*x)#include #include #define maxrept 1000double f1(double x){ return (3*x+1)/(x*x);}double f(double t){ double y,z; y=f1(t); z=f1原创 2017-11-21 22:55:17 · 2822 阅读 · 0 评论 -
牛顿迭代法
f(x)=x*x*x-3*x-1 迭代公式:x=x-(x*x*x-3*x-1)/(3*x*x-3)#include #include #define maxrept 1000double f(double x){ return (x-(x*x*x-3*x-1)/(3*x*x-3)); //迭代公式}void main(){ int k=0原创 2017-11-20 22:57:40 · 356 阅读 · 0 评论 -
高斯-赛德尔迭代法
b.txt-10.01 9.05 0.12 1.431.22 -4.33 2.67 3.221.25 -3.69 -12.37 0.58#include #define M 3main(){ FILE *f; double a[M][M],b[M]; double x[M]={0,0,0},y[M]; int i,j,k,n=100;原创 2017-12-03 22:13:47 · 4278 阅读 · 0 评论 -
Jacobi迭代法
#include #define M 3main(){ FILE *f; double a[M][M],b[M]; double x[M]={0,0,0},y[M]; int i,j,k,n=100; f=fopen("a.txt","r"); for(i=0;i<M;i++) { for(j=0;j<M;j++)原创 2017-12-03 21:50:11 · 1978 阅读 · 0 评论 -
Lagrange插值多项式
#include #define n 7double x[10],y[10];double lag(double xx){ double l,s=0; int i,k; for(i=0;i<n;i++) { l=1; for(k=0;k<n;k++) if(k-i) l=l*(xx-原创 2017-12-10 17:42:30 · 1940 阅读 · 0 评论 -
平方根法
#include "stdio.h"#includemain(){ FILE *f; double a[15][15],l[15][15],y[15],b[15],x[15]; double s; int i,j,k,n; f=fopen("sqrt.txt","r"); fscanf(f,"%d",&n); for(i=1;i<=n;原创 2017-11-28 20:11:21 · 2629 阅读 · 1 评论 -
最小二乘法(高斯)
# include# include# define N 12 ///12个值# define M 4# include# includeusing namespace std;///这是数据初始化。double x[N]= {0,5,10,15,20,25,30,35,40,45,50,55};double y[N]= {0,1.27,2.16,2.86,3.44,3.87原创 2017-12-19 12:17:23 · 2325 阅读 · 3 评论