matlab实习报告Word格式.docx
- 文档编号:22644271
- 上传时间:2023-02-05
- 格式:DOCX
- 页数:20
- 大小:328.01KB
matlab实习报告Word格式.docx
《matlab实习报告Word格式.docx》由会员分享,可在线阅读,更多相关《matlab实习报告Word格式.docx(20页珍藏版)》请在冰豆网上搜索。
3,通过下拉菜单File:
Preferences改变设置指令窗中的字体、颜色等。
如:
数据显示格式的设置:
Formatshort小数点后四位(缺省情况),如显示p为3.1416。
Formatlong小数点后十四位,如显示p为3.14159265358979。
Formantbank小数点后两位,如显示p为3.14。
Formatshorte小数点后四位科学记数法,如显示p为3.1416e+000。
Formatlonge小数点后十四五位科学记数法,如显示p为3.141592653589793e+000。
在指令窗中键入y1=2*sin(0.3*pi)/(1+sqrt(5))
然后利用回调指令计算y2=2*y1*cos(0.3*pi)/(1+sqrt(5))
y1=
0.5000
y2=
0.1816
4,了解软件中的一些功能和命令。
包括熟悉历史指令窗(CommandHistory)、熟悉当前目录浏览器(CurrentDiretory)、熟悉工作空间浏览器(WorkspaceBrowser),以及一些指令如:
clear清除当前工作区中的所有变量
clc清除指令窗中内容(未清除当前工作区中的变量)
clf清除图形窗口
cd设置当前工作目录
exit,quit退出Matlab
实习内容二数值数组及其运算
1,实习目的
掌握一维数组的创建和寻访;
掌握二维数组的创建和寻访;
掌握的Matlab矩阵和数组的运算;
熟悉Matlab关系操作和逻辑操作;
2,实习内容和步骤
(1)常量与变量
变量名区分字母大小写;
变量名必须以字母打头,其后可以是任意字母,数字,或下划线的组合。
(2)数值型向量(矩阵)的输入
任何矩阵(向量),可以直接按行方式输入每个元素;
同一行中的元素用逗号或者空格符来分隔行行之间用分号分隔。
所有元素处于[]中
(3)符号向量(矩阵)的输入
a.用函数sym定义符号矩阵:
函数sym实际是在定义一个符号表达式,这时的符号矩阵中的元素可以是任意的符号和表达式,没有长度限制,只需将方括号置于单引号中。
b.用函数sym定义符号矩阵
先定义矩阵中的每一个元素为一个符号变量,然后像普通矩阵一样输入符号矩阵。
(4)矩阵的运算
运算符:
+(加)、-(减)、*(乘)、./(右除)、.\(左除)、^(乘方)、’(转置)等;
常用函数:
det(行列式)、inv(逆矩阵)、rank(秩)、eig(特征值,特征向量)、rref(化矩阵为行最简形)
A=[20-1;
132];
B=[17-1;
423;
201];
M=A*B%矩阵A与B按矩阵运算相乘
det_B=det(B)%矩阵A的行列式
rank_A=rank(A)%矩阵A的秩
inv_B=inv(B)%矩阵B的逆矩阵
[V,D]=eig(B)%矩阵B的特征值矩阵V与特征向量构成的矩阵D
X=A/B%A/B=A*B-1,即XB=A,求X
M=
014-3
171310
det_B=20
rank_A=2
inv_B=
0.1000-0.35001.1500
0.10000.1500-0.3500
-0.20000.7000-1.3000
V=
-0.70940.74440.7444
-0.6675-0.3599+0.0218i-0.3599-0.0218i
-0.2263-0.5587-0.0607i-0.5587+0.0607i
D=
7.268000
0-1.6340+0.2861i0
00-1.6340-0.2861i
X=
0.4000-1.40003.6000
0.00001.5000-2.5000
(5)要在闭区间[0,2]上产生50个等距采样的一维数组,试用两种不同的指令实现,要寻访1到5个元素如何实现;
寻访7到最后一个元素如何实现;
要寻访2,6,8个元素如何实现;
寻访大于2的元素如何实现。
A=[0:
pi/25:
pi*2]
B=[A
(1)A
(2)A(3)A(4)A(5)]
C=[A(7):
A(51)]
A(find(A>
2))
A=
00.12570.25130.37700.5027
0.62830.75400.87961.00531.1310
1.25661.38231.50801.63361.7593
1.88502.01062.13632.26192.3876
2.51332.63892.76462.89033.0159
3.14163.26733.39293.51863.6442
3.76993.89564.02124.14694.2726
4.39824.52394.64964.77524.9009
5.02655.15225.27795.40355.5292
5.65495.78055.90626.03196.1575
6.2832
B=
C=
0.75400.87961.00531.13101.2566
1.38231.50801.63361.75931.8850
2.01062.13632.26192.38762.5133
2.63892.76462.89033.01593.1416
3.26733.39293.51863.64423.7699
3.89564.02124.14694.27264.3982
4.52394.64964.77524.90095.0265
5.15225.27795.40355.52925.6549
5.78055.90626.03196.15756.2832
ans=
5.78055.90626.03196.15756.2832
(6)在指令窗中输入
x=[23pi/29];
x=[2,3,pi/2,9]观察结果是否一样?
x=[23pi/29]
y=[2,3,pi/2,9]
x=
2.00003.00001.57089.0000
y=
(7)在指令窗中输入A=[123;
456]观察结果;
分别用“全下标”和“单下标”来寻访第2行第1列的元素,并说明二维数组元素的存储顺序。
>
a=[123;
456]
a=
123
456
a(2,1)
4
a
(2)
(8)在指令窗中输入B=【1,2,3;
6,7,8】,并计算A.*B;
和A*B并说明矩阵相乘和数组相乘的区别。
a=[123;
456];
b=[123;
678];
m=a.*b
n=a*b
m=
149
243548
n=
1444
32107
(9)在区间[0,10]中,绘出
的曲线。
1/20:
10;
y=1-exp(-0.5*t).*cos(2*t);
plot(y)
(10)字符串数组
clear
a=2;
w=3;
t=0:
0.01:
y=exp(-a*t).*sin(w*t);
[y_max,i_max]=max(y);
t_text=['
t='
num2str(t(i_max))];
y_text=['
y='
num2str(y_max)];
max_text=char('
maximum'
t_text,y_text);
tit=['
y=exp(-'
num2str(a),'
t)*sin('
num2str(w),'
t)'
];
plot(t,zeros(size(t)),'
k'
)
holdon
plot(t,y,'
b'
plot(t(i_max),y_max,'
r.'
'
MarkerSize'
20)
text(t(i_max)+0.3,y_max+.05,max_text)
title('
tit'
),xlabel('
t'
),ylabel('
y'
),holdoff
实习内容三Matlab图形绘制基础
1,实习目的
(1)进一步熟悉M文件的调试过程
(2)熟练掌握Matlab二维曲线的绘制
(3)掌握图形的修饰
2,实习内容
(1)用plot绘制y=sinx/x的图形
1.1用plot绘制y=sinx/x的图形
x=-2*pi:
pi/10:
2*pi;
y=sin(x)./x;
%准备数据
plot(x,y)%用plot作图
1.2
x=x+(x==0)*eps;
%(准备数据,处理分母零点))
plot(x,y)
1.3画出
包罗线
t=(0:
pi/100:
pi);
y1=sin(t);
y=-sin(t);
y2=sin(t).*sin(9*t);
t3=pi*(0:
9)/9;
y3=sin(t3).*sin(9*t3);
plot(t,y1,'
r:
'
t,y,'
t,y2,'
t3,y3,'
bo'
),axis([0,pi,-1,1])
gtext('
该点数值x'
三维曲线绘图:
(2)
输入以下程序,并得出结果。
t=-2*pi:
x=2*cos(t);
y=2*sin(t),z=3*t;
plot3(x,y,z,'
);
(3)三维曲面作图
x=-1:
0.05:
1;
y=x;
[x,y]=meshgrid(x,y);
z=x.^2+2*y.^2;
mesh(x,y,z)
(4)三维表面命令
surf(x,y,z)
(5)球面(sphere)及柱面(cylinder)的表达
5.1sphere(50)%绘制分格线条数为50的单位球面
5.2cylinder(8,50)
3,实习步骤
(1)plot指令的使用
X=1:
10
Y=[12345678]
Z=1:
-1:
-10
plot(x)
plot(x,z)
Columns1through9
123456789
Column10
10
Y=
1234567
Z=
10-1-2-3-4-5-6-7
Columns10through12
-8-9-10
(2)subplot指令的使用;
采样点数对绘图的影响的理解;
title、axis指令的使用
t1=(0:
11)/11*pi;
y1=sin(t1).*sin(9*t1);
t2=(0:
100)/100*pi;
y2=sin(t2).*sin(9*t2);
subplot(2,2,1),plot(t1,y1,'
),axis([0,pi,-1,1]),title('
子图
(1)'
),gtext('
不太完整的离散图'
subplot(2,2,2),plot(t2,y2,'
子图
(2)'
完整的离散图形'
subplot(2,2,3),plot(t1,y1,t1,y1,'
子图(3)'
是真的连续图形'
subplot(2,2,4),plot(t2,y2),axis([0,pi,-1,1]),title('
子图(4)'
正确的连续图形'
(3)其它绘图指令的使用及精细指令的使用
clf;
t=6*pi*(0:
100)/100;
y=1-exp(-0.3*t).*cos(0.7*t);
tt=t(find(abs(y-1)>
0.05));
ts=max(tt);
plot(t,y,'
r-'
linewidth'
3)
axis([-inf,6*pi,0.6,inf])
set(gca,'
xtick'
[2*pi,4*pi,6*pi],'
ytick'
[0.95,1,1.05,max(y)])
gridon
title('
\ity=1-e^{-\alphat}cos{\omegat}'
text(13.5,1.2,'
\fontsize{12}{\alpha}=0.3'
text(13.5,1.1,'
\fontsize{12}{\omega}=0.7'
holdon;
plot(ts,0.95,'
markersize'
10);
holdoff
cell_string{1}='
\fontsize{12}\uparrow'
;
cell_string{2}='
\fontsize{12}\fontname{隶书}镇定时间'
cell_string{3}='
\fontsize{6}'
cell_string{4}=['
\fontsize{14}\rmt_{s}=num2str(ts)'
text(ts,0.85,cell_string)
xlabel('
\fontsize{14}\bft\rightarrow'
ylabel('
实习总结:
通过一周的matlab上机实习,让我感受颇多。
人类的智慧是无穷的,能利用各个程序或软件来更方便的解决生活和学习中问题。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- matlab 实习 报告