《MATLAB及应用》实验指导书
班 级: T943-3
姓 名: 尤敏
学 号: 20090430326
总评成绩:
汽车工程系
电测与汽车数字应用中心
目录
实验04051001 MATLAB语言基础 1
实验04051002 MATLAB科学计算及绘图 3
实验04051003 MATLAB综合实例编程 5
实验04051001 MATLAB语言基础
操作成绩
报告
软件系统测试报告下载sgs报告如何下载关于路面塌陷情况报告535n,sgs报告怎么下载竣工报告下载
成绩
1实验目的
1) 熟悉MATLAB的运行环境
2) 掌握MATLAB的矩阵和数组的运算
3) 掌握MATLAB符号表达式的创建
4) 熟悉符号方程的求解
2实验内容
1. 课本第一章第2题.假设x=2和y=5。使用MATLAB计算一下公式的值。
a.
b.
c.
d.
解:x=3;y=5;
a=y*x^3/(x-y)
b=3*x/2/y
c=3/2*x*y
d=x^5/(x^5-1)
a =-67.5000
b =0.9000
c =22.5000
d =1.0041
2. 课本第一章第18题.
解:
a=[13,182,-184,2503];
roots(a)
ans =
-15.6850
0.8425 + 3.4008i
0.8425 - 3.4008i
3. 课本第二章第6题.
解:
A=[3 7 -4 12;-5 9 10 2;6 13 8 11;15 5 4 1];
B=A(:,2:4)
C=A(2:4,:)
D=A(1:2,2:4)
B =
7 -4 12
9 10 2
13 8 11
5 4 1
C =
-5 9 10 2
6 13 8 11
15 5 4 1
D =
7 -4 12
9 10 2
4. 课本第二章第12题
解:
A=[-7,16;4,9];B=[6,-5;12,-2];C=[-3,-9;6,8];
disp('A+B+C');A+B+C
disp('A-B+C');A-B+C
disp('验证结合律:(A+B)+C=A+(B+C)');
disp('(A+B)+C');(A+B)+C
disp('A+(B+C)');A+(B+C)
disp('验证交换律:A+B+C=B+C+A=A+C+B');
disp('A+B+C');A+B+C
disp('B+C+A');B+C+A
disp('A+C+B');A+C+B
A+B+C
ans =
-4 2
22 15
A-B+C
ans =
-16 12
-2 19
验证结合律:(A+B)+C=A+(B+C)
(A+B)+C
ans =
-4 2
22 15
A+(B+C)
ans =
-4 2
22 15
验证交换律:A+B+C=B+C+A=A+C+B
A+B+C
ans =
-4 2
22 15
B+C+A
ans =
-4 2
22 15
A+C+B
ans =
-4 2
22 15
5. 课本第二章第32题.
解:
t=0:0.5:5;
P=zeros(11,3);
P(:,1)=2;
P(:,2)=10*t+3;
r=P(11,:)
v=[0 10 0];
L=cross(r,v)
r =
2 53 0
L =
0 0 20
6. 课本第三章第4题.
解:x=1:10;
y1=sinh(x);
y2=(exp(x)-exp(-x))/2;
y=y1-y2
y =
1.0e-15 *
-0.2220 0 0 0 0 0 0 0 0 0
7. 课本第三章第10题.
解:function t=ht(h,v0,g)
a=[-1/2*g,v0,-h];
t=roots(a);
计算时输入如下指令:
ht(100,50,9.81)
ans =
7.4612
2.7324
8. 课本第三章第16题.
解:
function [volume,area]=torus(a,b)
% a = inner diameter
% b = outer diameter
volume=pi^2/4*(a+b).*(b-a).^2;
area=pi^2*(b.^2-a.^2);
end
在MATLAB中调用
[volume,area]=torus(3,4)
volume =
17.2718
area =
69.0872
9. 课本第四章第2题.
解:
a=2:4;
b=[10,24,24];
c=[12,48,100];
for k=1:3
disp(['NO.',num2str(k)])
delta=b(k)^2-4*a(k)*c(k);
x1=(-b(k)+sqrt(delta))/(2*a(k));
x2=(-b(k)-sqrt(delta))/(2*a(k));
%根据delta的数值
分析
定性数据统计分析pdf销售业绩分析模板建筑结构震害分析销售进度分析表京东商城竞争战略分析
判断是否显示实部、虚部
if(delta<0)
disp('complex value:real and image part:')
disp([real(x1),imag(x2)]);
else
disp('real value:')
disp([x1,x2]);
end
end
NO.1
real value:
-2 -3
NO.2
real value:
-4 -4
NO.3
complex value:real and image part:
-3 -4
10. 课本第四章第20题.
解:a.
t0=0;
x0=-10;
y0=144;
for t = 0:0.01:4
x=5*t-10;
y=25*t^2-120*t+144;
if((x^2+y^2)<(x0^2+y0^2))
t0=t;
x0=x;
y0=y;
end
end
disp('time:')
disp(t0);
disp('Minimium distance is:');
d=sqrt(x0^2+y0^2);
disp(d);
time:
2.2300
Minimium distance is:
1.3581
b.
t=0:0.01:4;
x=5*t-10;
y=25*t.^2-120*t+144;
d=sqrt(x.^2+y.^2);
[d_min,pos]=min(d);
disp('time:')
disp(t(pos));
disp('Minimiom distance is:');
disp(d_min);
time:
2.2300
Minimiom distance is:
1.3581
3思考题
1. MATLAB中,数组与矩阵在表示与应用上有哪些区别。
答:在表示上一维数组相当于向量,二维数组相当于矩阵,所以矩阵是数组的子集。
在应用上矩阵运算按线性变换定义,使用通常符号;数组运算按对应元素运算定义,使用点运算符。
实验04051002 MATLAB科学计算及绘图
操作成绩
报告成绩
1实验目的
1) 熟悉MATLAB所提供的常用数值计算的
函
关于工期滞后的函关于工程严重滞后的函关于工程进度滞后的回复函关于征求同志党风廉政意见的函关于征求廉洁自律情况的复函
数(方程(组)的求解、插值、拟合);
2) 掌握MATLAB二维图形绘制命令及其图形控制(plot、loglog、contour、polar等);
3) 熟悉MATLAB三维图形绘制命令及其图形控制(mesh、surf等)。
2实验内容
1. 课本第五章第8题.
解:
a.
t=0:10;
x1=t;
y1=0.5*t.^2+10;
y2=2*x1+6;
plot(x1,y1,x1,y2,'--')
xlabel('time(hour)');
legend('ship','boundary')
由图可以看出,该船并未越过边界线。
b.
t=0:10;
x1=t;
y1=0.5*t.^2+10;
A=2;B=-1;C=6;
d=(A*x1+B*y1+C)/sqrt(A^2+B^2);
d=abs(d);
plot(t,d);
xlabel('time(hour)');
ylabel('distance from boundary');
2. 课本第五章第22题.
解:
x=[1990,1991,1992,1993,1994];
y=[18,19,21,17,20];
subplot(3,1,1);stem(x,y)
xlabel('year');ylabel('temperature');
subplot(3,1,2);bar(x,y)
xlabel('year');ylabel('temperature');
subplot(3,1,3);stairs(x,y)
xlabel('year');ylabel('temperature');
3. 课本第六章第6题.
解:
a.
A=[12 -5 0;-3 4 7;6 2 3];
b=[11;-3;22];
solution=A\b
A*solution
solution =
3.0000
5.0000
-2.0000
b.
A=[6 -3 4;12 5 -7;-5 2 6];
b=[41;-26;14];
solution=A\b
A*solution
solution =
2
-3
5
4. 课本第六章第8题.
解:a.
A=[-2,1;-2,1];
b=[-5;-5.00001];
if(det(A)==0)