- 1、本文档共7页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
2012级信息与计算科学专业上机试题(学生版)
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
文档评论(0)