
数值分析
文章平均质量分 50
发烧的小龙虾
BUAAer
展开
-
秦九韶算法:计算从1加到100000
#include #include #include using namespace std;int main(){ clock_t start,finish; long x=1,a[100000],i,j,s=0,t; for(i=0;i<100000;i++) a[i]=i+1; start=clock(); for(i=0;原创 2017-10-29 14:20:05 · 1190 阅读 · 0 评论 -
线性方程组的迭代解法:高斯-塞德尔迭代法
#include #include #include #include using namespace std;double max(double x[],int n){ double max_=abs(x[0]); for(int i=1;i<n;i++) if(abs(x[i])>max_) max_=abs(x[i]);原创 2017-10-29 15:23:53 · 1724 阅读 · 0 评论 -
求解方程根的近似解:一般迭代法
#include #include #include using namespace std;double f(double x){ return log(4-x)/2;}int main(){ double x=1,y; int count=1; do { y=f(x); printf("%-3d: ",co原创 2017-10-29 16:08:31 · 1918 阅读 · 0 评论 -
求解方程根的近似解:埃特金加速迭代法
#include #include #include using namespace std;double f(double x){ return log(4-x)/2;}int main(){ double x=1,x1,x2,y; int count=1; do { x1=f(x); x2=f(x1);原创 2017-10-29 16:19:17 · 4861 阅读 · 0 评论 -
求解方程根的近似解:牛顿法
#include #include #include using namespace std;double f(double x){ return exp(2*x)+x-4;}double f1(double x){ return 2*exp(2*x)+1;}int main(){ double x=1,y; int count=1;原创 2017-10-29 16:26:32 · 2529 阅读 · 0 评论 -
求解方程根的近似解:弦截法
#include #include #include using namespace std;double f(double x){ return exp(2*x)+x-4;}int main(){ double a=0,b=1,x; int count=1; do { x=b-f(b)*(b-a)/(f(b)-f(a));原创 2017-10-29 16:34:53 · 1721 阅读 · 0 评论 -
线性方程组的迭代解法:雅可比迭代法
#include #include #include #include using namespace std;double max(double x[],double y[],int n){ double max_=abs(x[0]-y[0]); for(int i=1;i<n;i++) if(abs(x[i]-y[i])>max_)原创 2017-10-29 16:50:08 · 1946 阅读 · 0 评论 -
线性方程组的迭代解法:超松弛迭代法
#include #include #include #include using namespace std;double max(double x[],int n){ double max_=abs(x[0]); for(int i=1;i<n;i++) if(abs(x[i])>max_) max_=abs(x[i]);原创 2017-10-29 16:56:18 · 4956 阅读 · 0 评论 -
利用MATLAB理解最小二乘法
x=[1 2 3 4 6 7 8];y=[2 3 6 7 5 3 2];p=polyfit(x,y,2);plot(x,y,'bo');hold on;x=0:0.01:10;y=p(1)*x.^2+p(2)*x+p(3);plot(x,y,'r');原创 2017-10-29 17:10:37 · 1306 阅读 · 1 评论