数学系10级2班陈志明20101021223.doc

数学系10级2班 陈志明 20101021223 实验一 高射炮的控制区域 曲线簇及其包络曲线的绘制 n=input(input n:); input n:20 alpha=(2:n)*pi/(2*n); for k=1:n-2 a=alpha(k); v1=cos(a);v2=sin(a); t0=v2/4.9; t=(0:16)*t0/16; x(k,:)=v1*t; y(k,:)=v2*t-4.9*t.^2; end plot(x,y) hold on g=9.8; t=1/g:0.001:sqrt(2)/g; x=sqrt(t.^2-1/g^2); y=1/g-0.5*g*t.^2; plot(x,y) Koch分形曲线的绘制 p=[0 0;10 0]; a=[cos(pi/3) -sin(pi/3);sin(pi/3) cos(pi/3)]; for k=1:5 n=max(size(p)); d=diff(p)/3; q=p(1:n-1,:); p(5:4:4*n-3,:)=p(2:n,:); p(2:4:4*n-6,:)=q+d; p(3:4:4*n-5,:)=p(2:4:4*n-6,:)+d*a; p(4:4:4*n-4,:)=q+2*d; end plot(p(:,1),p(:,2)) 实验二 连续计息问题 导函数的值能反映函数的变化。考虑函数在内的图像。 syms x; fun=x*x*cos(x^2+3*x-4); fplot(fun,[-2,2]) grid 观察其导函数 dfun=diff(fun) dfun = 2*x*cos(x^2+3*x-4)-x^2*sin(x^2+3*x-4)*(2*x+3) dfun=char(dfun); hold on fplot(dfun,[-2,2],r) 考虑内函数的变化情况 hold off fplot(dfun,[-0.5,0.5],r) grid