数模实验课报告重邮Word文档下载推荐.docx
- 文档编号:1460550
- 上传时间:2023-04-30
- 格式:DOCX
- 页数:50
- 大小:1.60MB
数模实验课报告重邮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
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数模 实验 报告