2012级信息与计算科学专业上机试题(学生版).doc

2012级信息与计算科学专业 《数学模型与数学实验》上机试题 第一部分：验证试验 试验要求：对每一行的语句后面加上注释，并将程序运行的结果放在程序后面 第一题： [X,Y,Z]=peaks(30);%返回peaks函数上的三个坐标轴上的函数值 pcolor(X,Y,Z); shading interp %伪彩色图 hold on %保持上面图形 contour(X,Y,Z,19, ‘ k ‘ ) % add 19 contour lines in black xlabel( ‘ X-axis ‘ ),ylabel( ‘ Y-axis ‘ )%对x轴和y轴进行说明 title( ‘ PCOLOR and CONTOUR of PEAKS ‘ ) hold off 第二题： x=linspace(0,3*pi); %对x进行赋值，范围是0到3pi z1=sin(x); %z1代表的函数 z2=sin(2*x); %z2代表的函数 z3=sin(3*x);%z3代表的函数 y1=zeros(size(x)); %求y1的值 y3=zeros(size(x)); %求y3的值 y2=y3/2;%求y2的值 plot3(x,y1,z1,x,y2,z2,x,y3,z3);%画x-y1，x-y2和x-y3的图像 grid,xlabel( ‘ x-axis ‘ ),ylabel( ‘ y-axis ‘) ,abel( ‘ z-axis ‘ ) title( ‘ sin(x),sin(2x),sin(3x) ‘ )%加标题 第三题： x=(0 : .1 : 1)*2*pi; %0为开始值，1为终止值，步长为0.1 y=sin(x);%y的函数表达式 pp=spline(x, y); %（x，y）作为采样点 ppd=spderiv(pp); %寻找微分的样条表示 xi=linspace(0, 2*pi); %xi的范围 yi=ppval(pp, xi); %对采样点（x，y）样条插值 yyd=ppval(ppd, xi); %样条插值 plot(x, y, ‘ o ‘, xi, yi, xi, yyd, ‘ - ‘)%画出x-y的图像并且用圆圈标注以及xi-yi的函数图像用直线。 第四题： x=[0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1];%初始值为0，步长为0.1， y=[-.447 1.978 3.28 6.16 7.08 7.34 7.66 9.56 9.48 9.30 11.2];%建立一个向量 n=2; p=polyfit(x, y, n)%曲线n次拟合 xi=linspace(0, 1, 100); %对xi规定范围（0-100，步长为1） z=polyval(p, xi);%求近似值 plot(x, y, o , x, y, xi, z, : );%画x-y和xi-yi的函数图像 xlabel( x ); %标注x轴 ylabel( y=f(x) ),%标注y轴 title( Second Order Curve Fitting ); 第五题： R0=1; a=12*R0;b=9*R0;T0=2*pi; %对a，b和T0进行赋值，a为椭圆的长半轴，b为短半轴 T=5*T0;dt=pi/100;t=[0:dt:T]; %dt控制卫星运动的快慢 f=sqrt(a^2-b^2); th=12.5*pi/180; %求th的值 E=exp(-t/20); %函数E的表达式 x=E.*(a*cos(t)-f);%x的函数表达式 y=E.*(b*cos(th)*sin(t));%y的函数表达式 z=E.*(b*sin(th)*sin(t));%z的函数表达式 plot3(x,y,z,g) %画x-y-z的图像。 [X,Y,Z]=sphere(30);X=R0*X;Y=R0*Y;Z=R0*Z;%画出以30为经度和纬度的球面grid on,hold on,surf(X,Y,Z),shading interp%处理色彩效果 x1=-18*R0;%x轴下边界 x2=6*R0;%x轴上边界 y1=-12*R0;%y轴下边界 y2=12*R0;%y轴上边界 z1=-6*R0;%z轴下边界 z2=6*R0;%z轴上边界 axis([x1 x2 y1 y2 z1 z2])%调整图幅大小 view([117 37]); %观看视角 comet3(x,y,z,0.02);%卫星运动轨迹函数 hold off