数模实验课报告重邮Word文档下载推荐.docx
《数模实验课报告重邮Word文档下载推荐.docx》由会员分享，可在线阅读，更多相关《数模实验课报告重邮Word文档下载推荐.docx（50页珍藏版）》请在冰点文库上搜索。
ifx<
0&
x~=-3;
y=x^2+x-6;
elseifx>
=0&
x<
5&
x~=2&
x~=3;
y=x^2-5*x+6;
else
y=x^2-x-1;
然后在Matlab主窗体中调用计算：
2、已知：
s=1+2+22+23+24+…+2x（x=0，1，2，…，n，n<
32）
用循环结构及m函数实现。
先建一个名为loop.m文件，代码如下：
functions=loop（n）
s=0;
forx=0:
1:
n;
s=s+2^x;
然后在Matlab主窗体中调用计算（取n=32）：
3、输入一个百分制成绩，要求输出成绩等级优、良、中、及格、不及格。
其中90分—100分为优，80分—89分为良，70分—79分为中，60分—69分为及格，60分以下为不及格。
用swicth实现，用m文件编写。
先建一个名为switch_case.m文件，代码如下：
functionswitch_case（grade）
switchgrade
casenum2cell（90:
100）
grade='
优'
casenum2cell（80:
89）
良'
casenum2cell（70:
79）
中'
casenum2cell（60:
69）
及格'
casenum2cell（0:
59）
不及格'
otherwise
输入有误'
end
实验三matlab绘图
PPT中的练习：
%4
x=linspace（0,2*pi,30）;
y=sin（x）;
z=cos（x）;
plot（x,y,'
r'
x,z,'
go'
title（'
正弦余弦'
xlabel（'
加X轴标记'
ylabel（'
加Y轴标记'
grid
text（1,0.5,'
陈志尧\bullet\pi\rightarrow'
%gtext（'
sinx'
）%自由标注
%gtext（{'
'
cosx'
}）
;
axis（[0503]）%设定最大和最小值
%axis（'
square'
）%将当前图形设置为方形
auto'
）%将坐标系统返回到自动缺省状态
axis（'
equal'
）%两个坐标因子设成相等
off'
）%关闭坐标系统
on'
）%显示坐标系统
%8
figure
x=linspace（0,2*pi,60）;
y1=sin（x）;
y2=cos（x）;
plot（x,y1,x,y2,'
--'
axis（[02*pi-22]）
legend（'
sin（x）'
cos（x）'
%可自由拖动
%9
ezplot（'
[0,pi]）%隐函数，
%实际上画的是cos函数，sin^2+cos^2=1
cos（t）^3'
sin（t）^3'
[0,2*pi]）
axis（[-11-11]）
exp（x）+sin（x*y）'
[-2,0.5,0,2]）
%12
fplot（'
myfun1'
[-1,2]）
）
tanh'
[-2,2]）
[tanh（x）,sin（x）,cos（x）]'
2*pi*[-11-11]）
tanh（x）,sin（x）,cos（x）'
%16
x=logspace（-1,2）;
loglog（x,exp（x）,'
-s'
双对数'
gridon
x=0:
.1:
10;
semilogy（x,10.^x）
y为对数轴'
x=-100:
100;
y=x.^3;
subplot（221）;
plot（x,y）
y=x^3函数图'
subplot（222）;
semilogx（x,y）
y=x^3x轴半对数坐标系'
subplot（223）;
semilogy（x,y）
y=x^3y轴半对数坐标系'
subplot（224）;
loglog（x,y）
y=x^3全对数坐标系'
%18
x=[-2.5:
0.25:
2.5];
y=exp（-x.*x）;
stairs（x,y）;
%绘制阶梯图形命令
stairsplot'
bar（x,y）;
%绘制条形图命令
条形图'
x=[01100];
%正方形顶点坐标向量
y=[00110];
fill（x,y,'
y'
%绘制并以黄色填充正方形图
黄色填充正方形图-1'
x=[0:
0.025:
2*pi];
y=sin（3*x）;
fill（x,y,[000]）;
黑色填充正方形图-y=sin（3*x）'
%24
pi/50:
10*pi;
plot3（sin（t）,cos（t）,t）
sin（t）,cos（t）三维图'
rotate3d%旋转
x=-3:
3;
y=1:
5;
[X,Y]=meshgrid（x,y）;
%产生一个以向量x为行，向量y为列的矩阵
Z=（X+Y）.^2;
plot3（X,Y,Z）
meshgrid'
y1=sin（t）,y2=cos（t）;
plot3（y1,y2,t）;
helix'
）,text（0,0,0,'
origin'
sin（t）'
）,ylabel（'
cos（t）'
）,zlabel（'
t'
gridon;
%28
surf（X,Y,Z）
Z=（X+Y）^2surf'
shadingflat%将当前图形变得平滑，不加时会有格子
mesh（X,Y,Z）
Z=（X+Y）^2mesh'
%39
zoomon
%40
mesh（X,Y,Z）
Z=（X+Y）.^2图1'
mesh（X,Y,Z）;
view（50,-34）
Z=（X+Y）.^2图2'
view（-60,70）
Z=（X+Y）.^2图3'
view（0,1）
Z=（X+Y）.^2图4'
%44
m=moviein（20）;
%建立一个20个列向量组成的矩阵
forj=1:
20
plot（fft（eye（j+10）））%绘制出每一幅眼球图并保存到m矩阵中
title（'
不断变化的眼球'
m（:
j）=getframe;
%截取动画帧，保存到m变量中
pause（.1）%延时语句，越大越慢，1就是延时1秒。
movie（m,1,10）;
%第二个数是播放画面次数，第二个是每秒播放帧数。
%播放时，第一遍是截取动画帧，所以看起来很慢；
%接下来就是按照预设的速度播放了
%45
n=30;
[x,y,z]=sphere;
m=moviein（30）;
surf（x,y,z）
一个直径不断变化的球体'
movie（m,30）;
%程序有错，没搞明白。
[x,y,z]=peaks（30）;
axis（[-33-33-1010]）
m=moviein（15）;
fori=1:
15
view（-37.5+24*（i-1）,30）
动画曲面peaks'
i）=getframe;
movie（m）;
%只有截取动画帧，改了参数还是不会播放。
%48极坐标图
theta=linspace（0,2*pi）,
rho=sin（2*theta）.*cos（2*theta）;
polar（theta,rho,'
g'
sin（2θ）.*cos（2θ）的极坐标图'
%49散点图
loadseamount%seamounts水下山脉纵坐标都是负的
scatter（x,y,5,z）
水下山脉'
%49平面等值线图
[X,Y]=meshgrid（-2:
.2:
2,-2:
3）;
Z=X.*exp（-X.^2-Y.^2）;
[C,h]=contour（X,Y,Z）;
多个等值线二维图'
clabel（C,h）%沿线标识法
colormapcool
%50空间等值线图
[x,y,z]=peaks;
subplot（121）
contour3（x,y,z,16,'
s'
）
grid
x-axis'
y-axis'
z-axis'
山峰的三维等值线图'
subplot（122）
contour（x,y,z,16,'
）,ylabel（'
山峰的二维等值线图'
%51三维散点图
[x,y,z]=sphere（16）;
X=[x（:
）*.5x（:
）*.75x（:
）];
Y=[y（:
）*.5y（:
）*.75y（:
Z=[z（:
）*.5z（:
）*.75z（:
S=repmat（[1.75.5]*10,prod（size（x））,1）;
C=repmat（[123],prod（size（x））,1）;
scatter3（X（:
）,Y（:
）,Z（:
）,S（:
）,C（:
）,'
filled'
）,view（-60,60）
三维散点图'
课内作业：
%
（1）
pi/10:
2*pi;
plot（x,y1,'
-g'
x,y2,'
--r'
LineWidth'
6）;
二条曲线'
%
（2）
0.001:
vt=10*exp（-0.2*t+pi*t.^2）;
plot（t,vt）;
x（θ）=sin（θ）,y=cos（θ-2.5）直角坐标系图'
theta=linspace（0,3*pi）;
x=sin（theta）;
y=cos（theta-2.5）;
z=sqrt（x.^2+y.^2）;
polar（theta,z,'
x（θ）=sin（θ）,y=cos（θ-2.5）极坐标系图'
%（3）
z=exp（2*x+3*y）;
plot3（x,y,z）;
三维曲线z=e.^（2*x+3*y）'
%（4）
y=-3:
z=x+sin（x.*y）;
z=x+sinxy表格图'
gridon
[x,y]=meshgrid（-3:
3,-2:
[c,h]=contour（x,y,z）;
clabel（c,h）
z=x+sinxy平面等值线图'
%（5）
text（pi,0,'
\rightarrow（pi,0）'
text（pi,-1,'
\rightarrow（pi,-1）'
二条曲线-标注版'
%（6）
x=[12822132];
explode=[1000000];
colormaphsv;
labels={'
重庆'
四川'
湖北'
河南'
天津'
福建'
云南'
};
pie（x,explode,labels）;
0101004班籍贯分布'
课外作业：
%1
%心形线不会画，以下是球体
axis（[-11-11-11]）
surf（x,y,z）
x'
ylabel（'
zlabel（'
z'
球体1'
surf（x,y,z）;
球体2'
球体3'
球体4'
球体'
%2
[x,y,r]=sphere;
r=（x.^2+y.^2）./x;
%此处有个警告
surf（x,y,r）
surf（x,y,r）
球面1'
surf（x,y,r）;
球面2'
球面3'
球面4'
%3
x=linspace（0,2*pi）;
y=sin（x）.*cos（5*x）;
axis（[02*pi-0.960.96]）;
y=sin（x）.*cos（5*x）'
y=abs（1000*sin（4*x））+1;
y=abs（1000*sin（4*x））+1'
%5
x=-6:
2;
subplot（211）
plot绘制sin（x）'
subplot（212）
[-6,2]）
fplot绘制sin（x）'
%6
y=linspace（0,2*pi）;
z=sin（x）.*cos（y）;
plot3（x,y,z）
z=sin（x）*cos（y）的三维网格图'
z=sin（x）*cos（y）的三维曲面图'
实验四MATLAB应用
1、确定企业年设备能力与年劳动生产率的关系
某市电子工业公司有14个所属企业，各企业的年设备能力与年劳动生产率统计数据如下表。
试分析企业年设备能力与年劳动生产率的关系。
若该公司计划新建一个设备能力为9.2千瓦/人的企业，估计劳动生产率将为多少？
解：
企业年设备能力与年劳动生产率的关系程序为：
x=[2.82.83.02.93.43.94.04.84.95.25.45.56.27.0];
y=[6.76.97.27.38.48.89.19.810.610.711.111.812.112.4];
X=[ones（size（x'
））,x'
],pause
[c,cint,r,rint,stats]=regress（y'
X,0.05）,pause
rcoplot（r,rint）
运行结果：
c=3.10031.4481
cint=2.25933.9414
1.26551.6308
r=-0.4551-0.2551-0.24480.00010.37600.05190.2071-0.25140.40380.06930.17970.73490.0212-0.8374
rint=
-1.24690.3366
-1.08470.5744
-1.08880.5993
-0.85300.8531
-0.46821.2202
-0.84270.9465
-0.67971.0939
-1.13370.6309
-0.45191.2594
-0.81760.9562
-0.69311.0525
0.00331.4664
-0.81240.8547
-1.3523-0.3224
stats=0.9613298.36310.0000
预测：
x=9.2;
y=3.1003+1.4481*x;
x,y
结果为：
x=9.2000
y=16.4228
该公司计划新建一个设备能力为9.2千瓦/人的企业，估计劳动生产率将为16.4228。
2、测定某矿脉的金属含量
一矿脉有13个相邻样本点，人为地设定一个原点，现测得各样本点与原点的距离x，与该样本点处某种金属含量y的一组数据如下，试建立合适的回归模型。
程序为：
x=[23457810111415161819];
y=[106.42109.2109.58109.5110109.93110.49110.59110.6110.9110.76111111.2];
[p,S]=polyfit（x,y,2）;
p
[Y,delta]=polyconf（p,x,S）;
Y
p=-0.01740.5334106.9328
Y=107.9303108.3769108.7889109.1661109.8165110.0896110.5317110.7007110.9995111.0296111.0251110.9119110.8033
升级会员 