数值分析-课程设计大作业.docVIP

课程设计 2013年07月20日 设计 题目成绩 课 程 设 计 主 要 内 容 指 导 教 师 评 语 建议：从学生的工作态度、工作量、设计（论文）的创造性、学术性、实用性及书面表达能力等方面给出评价。 签名： 20 年,第一至第五次猴子在夜里藏椰子后，椰子的数目分别为，，，，，再设最后每个人分得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,ex,bx]; B=[ay,cy,dy,ey,by]; %把由函数sub_koch1生成的五点横、纵坐标A,B顺次划分为四组,分别对应四条折线段中 %每条线段两端点的坐标,并依次分别存储在4*4阶矩阵k中,k中第i(i=1,2,3,4)行数字代表第%i条线段两端点的坐标 function w=sub_koch2(A,B) a11=A(1);b11=B(1); a12=A(2);b12=B(2); a21=A(2);b21=B(2); a22=A(3);b22=B(3); a31=A(3);b31=B(3); a32=A(4);b32=B(4); a41=A(4);b41=B(4); a42=A(5);b42=B(5); w=[a11,b11,a12,b12;a21,b21,a22,b22;a31,b31,a32,b32;a41,b41,a42,b42]; 调用函数得到下图 n=5; i=0;whi