沈阳航空航天大学研究生数值分析程序设计Word文档下载推荐.docx
- 文档编号:19389328
- 上传时间:2023-01-05
- 格式:DOCX
- 页数:12
- 大小:158.28KB
沈阳航空航天大学研究生数值分析程序设计Word文档下载推荐.docx
《沈阳航空航天大学研究生数值分析程序设计Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《沈阳航空航天大学研究生数值分析程序设计Word文档下载推荐.docx(12页珍藏版)》请在冰豆网上搜索。
课程结束时上交报告;
6.严禁抄袭行为,一旦发现,程序设计成绩为不及格。
一、通识题
1.作出下列函数的图形:
(1)编写M文件对下面函数进行分段曲面绘制:
(2)
.
(3)
及
要求:
a/分别考虑
;
b/
2.用数值方法计算
,
其中n分别取100,200,1000。
3.编写程序求
,其中n,a由键盘输入。
二、求解对称正定线性方程组的Cholesky分解法
1)掌握求解对称正定线性方程组的Cholesky分解法的基本原理;
2)编写求解对称正定线性方程组的Cholesky分解法的通用程序;
3)利用所编写的程序求解方程组
,其中
,
程序设计成绩评定单
指导教师评语:
课程设计成绩
任课教师签字
年月日
一、正文
一、通识题
(1)编写M文件对下面函数进行分段曲面绘制:
(2)
(3)
b/
n=5
n=10
n=20
n=50
结果:
n=100;
y=1.6350
n=200;
y=1.6399
n=1000;
y=1.6439
3.编写程序求
二、求解对称正定线性方程组的Cholesky分解法
答:
1)先求解
然后求解
2)见附录
3)
二、总结
通过对于这次数值分析的程序编写,我学到了很多的编程经验,并且形成了精益求精的数学学习态度,非常感谢老师在课堂上辛苦的讲解知识,这次大作业让我了解并学会应用matlab软件,让我的数学功底又变得扎实了很多,在我学习的过程中与同学的交流让我明白了自己的不足,让我为以后的发展有了明确的方向。
各个学科都离不开数学,只有数学扎实以后的科研和工作才能更顺利。
完稿日期:
2014年12月10日
附录
源程序:
1作出下列函数的图形:
[x,y]=meshgrid(linspace(-1,1));
z1=0.5*exp(-0.75*y.^2-3.75*x.^2-1.5*x).*((x+y)>
1);
z2=0.8*exp(-y.^2-6*x.^2).*(-1<
(x+y)&
(x+y)<
=1);
z3=0.5*exp(-0.75*y.^2-3.75*x.^2+1.5*x).*((x+y)<
=-1);
z=z1+z2+z3;
mesh(x,y,z)
(2)
x=-1:
1/1000:
1;
y=1./(1+25.*x.^2);
plot(x,y)
n=5时
x=-50:
50;
y2=0;
forn=1:
5
y2=y2+(-1)^(n-1)*x.^(2*n-1)/(2*n-1);
end
y1=atan(x);
plot(x,y1,x,y2)
n=10时
10
n=20时
20
n=50时
50
2.用数值方法计算
n=100
y=0;
k=100;
k
y=1/n^2+y;
y
S=1.6350
n=200
k=200;
S=1.6399
n=1000
k=1000;
S=1.6439
a=input('
pleaseinputa:
'
);
n=input('
pleaseinputn:
fori=1:
n;
symsnn;
a_sum(i)=symsum(a*10^nn,0,i-1)
end;
sum(a_sum)
2)functionsolution=cholesky(M,precision)
ifnargin==2
try
digits(precision);
catch
disp('
你输入的精度有误,这里按照缺省的精度(10位有效数字)计算'
)
digits(10);
else
A=vpa(M);
row=size(A,1);
col=size(A,2);
ifndims(A)~=2|col-row~=1
矩阵的大小有误,不能使用cholesky分解'
return
ifdet(M(:
1:
row))==0.0
该方程的系数矩阵行列式为零,无解或有无穷多解,不能使用Cholesky分解法'
B=M(:
row);
ifB==B'
L=chol(double(A(:
row)));
fori=row:
-1:
1
temp=A(i,col);
fork=i+1:
row
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 沈阳 航空航天 大学研究生 数值 分析 程序设计