《数学系10级2班陈志明20101021223》由会员分享,可在线阅读,更多相关《数学系10级2班陈志明20101021223(7页珍藏版)》请在金锄头文库上搜索。
1、数学系 10 级 2 班 陈志明 20101021223实验一 高射炮的控制区域曲线簇及其包络曲线的绘制 n=input(input n:);input n:20 alpha=(2:n)*pi/(2*n); for k=1:n-2a=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/g2); y=1/g-0.5*g*t.2; plot(x,y)K
2、och 分形曲线的绘制 p=0 0;10 0; a=cos(pi/3) -sin(pi/3);sin(pi/3) cos(pi/3); for k=1:5n=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)实验二 连续计息问题导函数的值能反映函数的变化。考虑函数 在)43cos()(2xxf内的图像。2, syms x; fun=x*x*cos(x2+3*x-4); fplot(fun,-2,2) grid观察其导函数 dfun=diff(fun)dfun =2*x*cos(x2+3*x-4)-x2*sin(x2+3*x-4)*(2*x+3) dfun=char(dfun); hold on fplot(dfun,-2,2,r)考虑 内函数的变化情况5.0, hold off fplot(dfun,-0.5,0.5,r) grid