安徽建筑工业学院-《数值分析》课程设计.doc

PAGE PAGE 1 课程设计 设计题目 设计题目 《数值分析》课程设计 学生姓名 **** 学 号 #### 专业班级 指导教师 2013年07月20日 设计 题目 《数值分析》课程设计 成绩 课 程 设 计 主 要 内 容 实验一：三题全做；实验二：四题全做；实验三：3.1 3.2已做；实验四 4.1 4.2 4.3 4.4已做；实验五：5.1 5.2 5.3 5.4已做；实验六：6.2 6.3 6.4 已做；实验七：7.3 7.5 7.6 7.8 已做；实验八：8.2 8.8已做。共计26道题。 部分题通过matlab中notebook实现程序，notebook使用非常方便。由于篇幅有限部分题调用的函数源代码没有给出，本作业所使用的m文件均在附录中 指 导 教 师 评 语 建议：从学生的工作态度、工作量、设计（论文）的创造性、学术性、实用性及书面表达能力等方面给出评价。 签名： 20 年 月 日 1.1水手、猴子和椰子问题 算法分析:设椰子起初的数目为,第一至第五次猴子在夜里藏椰子后，椰子的数目分别为，，，，，再设最后每个人分得x个椰子，由题意得： 利用逆向递推方法求解： n=input(n=); for x=1:n p=5*x+1; for k=1:5 p=5*p/4+1; end if p==fix(p) break end end disp([x,p]) 执行代码后得： n= 1023 15621 （输入n=1000000） 即最后每个人分得1023个椰子，椰子总数为15621 1.2当时，选择稳定的算法计算积分 由 得 由上式可知求时，的误差的影响被缩小了。n=100时的近似值为0。 matlab 代码为 fprintf(稳定算法:\n) y0=0; n=100; plot(n,y0,r*); hold on fprintf(y[100]=%10.6f,y0); while(1) y1=1/10*[(1-exp(-n))/n-y0]; fprintf(y[%10.0f]=%10.6f,n-1,y1);plot(n-1,y1,r*) if(n=1) break;end y0=y1;n=n-1; if mod(n,3)==0,fprintf(\n),end,end (具体值已省略) 编程实现得下图。 由图可知，该算法是稳定的。 1.3绘制静态和动态的Koch分形曲线 Koch曲线程序koch.m function koch(a1,b1,a2,b2,n) %koch(0,0,9,0,3) %a1,b1,a2,b2为初始线段两端点坐标,n为迭代次数 %例如a1=0;b1=0;a2=9;b2=0;n=3; %第i-1次迭代时由各条线段产生的新四条线段的五点横、纵坐标存储在数组A、B中 [A,B]=sub_koch1(a1,b1,a2,b2); for i=1:n for j=1:length(A)/5; w=sub_koch2(A(1+5*(j-1):5*j),B(1+5*(j-1):5*j)); for k=1:4 [AA(5*4*(j-1)+5*(k-1)+1:5*4*(j-1)+5*(k-1)+5),BB(5*4*(j-1)+5*(k-1)+1:5*4*(j-1)+5*(k-1)+5)]=sub_koch1(w(k,1),w(k,2),w(k,3),w(k,4)); end end A=AA; B=BB; end plot(A,B) hold on axis equal %由以(ax,ay),(bx,by)为端点的线段生成新的中间三点坐标并把这五点横、纵坐标依次分别存%储在数组A,B中 function [A,B]=sub_koch1(ax,ay,bx,by) cx=ax+(bx-ax)/3; cy=ay+(by-ay)/3; ex=bx-(bx-ax)/3; ey=by-(by-ay)/3; L=sqrt((ex-cx).^2+(ey-cy).^2); alpha=atan((ey-cy)./(ex-cx)); if (ex-cx)0 alpha=alpha+pi; end dx=cx+cos(alpha+pi/3)*L; dy=cy+sin(alpha+pi/3)*L; A=[ax,cx,dx