原创力文档- 1、本文档共8页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 5、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 6、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
数值课程设计典型数值算法C语言程序设计
数值计算课程设计说明书
题目： 典型数值算法的C++语言程序设计
院 （系）： 理学院
专业班级：
学 号：
学生姓名：
指导教师：
2012 年 6 月 21 日
目录
1、经典四阶龙格库塔法解一阶微分方程组
2、高斯列主元法解线性方程组
3、牛顿法解非线性方程组
4、龙贝格求积分算法
5、三次样条插值算法
6、M次多项式曲线拟合
经典四阶龙格库塔法解一阶微分方程组
例：x(t)=x+2y,y(t)=3x+2y,x(0)=6,y(0)=4的龙格库塔解法（数值方法课本p400 例9.15）
C++程序为：
#includeiostream
#includecmath
using namespace std;
float f(float t,float x,float y) //定义f函数
{float f;
f=x+2*y;
return(f);
}
float g(float t,float x,float y) //定义g函数
{float g;
g=3*x+2*y;
return(g);
}
int main()
{float h,n,k,t1,t2,x1,x2,y1,y2,f1,f2,f3,f4,g1,g2,g3,g4;
coutplease input t1,x1,y1,h,n;
cint1x1y1hn;
for(k=0;kn;k++)
{f1=f(t1,x1,y1);
g1=g(t1,x1,y1);
f2=f(t1+h/2,x1+h*f1/2,y1+h*g1/2);
g2=g(t1+h/2,x1+h*f1/2,y1+h*g1/2);
f3=f(t1+h/2,x1+h*f2/2,y1+h*g2/2);
g3=g(t1+h/2,x1+h*f2/2,y1+h*g2/2);
f4=f(t1+h,x1+h*f3,y1+h*g3);
g4=g(t1+h,x1+h*f3,y1+h*g3);
x2=x1+h*(f1+2*f2+2*f3+f4)/6;
y2=y1+h*(g1+2*g2+2*g3+g4)/6;
t2=t1+h;
cout(t2,x2,y2)endl;
t1=t2;
x1=x2;
y1=y2;
}
return 0;
}
运行结果为：
please input t1,x1,y1,h,n 0 6 4 0.02 10
(0.02,6.29355,4.53932)
(0.04,6.61562,5.11949)
(0.06,6.96853,5.74397)
(0.08,7.35474,6.41653)
(0.1,7.77697,7.14127)
(0.12,8.23814,7.9226)
(0.14,8.74141,8.76532)
(0.16,9.29021,9.6746)
(0.18,9.88827,10.6561)
(0.2,10.5396,11.7158)
Press any key to continue
若是求解其他相似的微分方程组，可修改f函数和g函数，主函数不变。
2、高斯列主元法解线性方程组
试验原理：
高斯消去法是一个古老的求解线性方程组的方法，但它的改进、变形得到的主元素消去法仍然是计算机上常用的计算方法。
列主元消去法
构造增广矩阵
设系数矩阵为
A=
设增广矩阵为
B=
b.找到每列的绝对值的最大的数
c.行变换(重复b-c)
d.消去
e. 回代
实验过程：
#includestdio.h
#includeconio.h
#includemath.h
#define N 100
#define epsilon 1e-6
float a[N][N+1];
void main( )
{
int i,j,k,n;
float t,s=0;
printf(Input juzhen jieshu:);
scanf(%d,n);
printf(\n);
for(i=0;in;i++)
{
printf(\n\n请输入a%d1--a%d%d系数和向量b%d:,i+1,i+1,n,i+1);
/*实现将每一行中的系数和向量一次性输入，数之间用空格格开，输完后回车确定*/
for(j=0;jn+1;j++) /*将刚才输入的数存入数组*/
scanf(%f,a[i
文档评论（0）