《《MATLAB及应用》实验指导书》由会员分享，可在线阅读，更多相关《《MATLAB及应用》实验指导书（35页珍藏版）》请在金锄头文库上搜索。
1、MATLAB及应用实验指导书班 级： T943-3 姓 名： 尤敏 学 号： 20090430326 总评成绩： 汽车工程系电测与汽车数字应用中心目录实验04051001 MATLAB语言基础1实验04051002 MATLAB科学计算及绘图3实验04051003 MATLAB综合实例编程5实验04051001 MATLAB语言基础操作成绩报告成绩1实验目的1) 熟悉MATLAB的运行环境2) 掌握MATLAB的矩阵和数组的运算3) 掌握MATLAB符号表达式的创建4) 熟悉符号方程的求解2实验内容1. 课本第一章第2题.假设x=2和y=5。使用MATLAB计算一下公式的值。a. b. c.
2、d.解：x=3;y=5; a=y*x3/(x-y) b=3*x/2/y c=3/2*x*y d=x5/(x5-1) a =-67.5000b =0.9000c =22.5000d =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
3、8 11 5 4 1C = -5 9 10 2 6 13 8 11 15 5 4 1D = 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+Cdisp(A-B+C);A-B+Cdisp(验证结合律：(A+B)+C=A+(B+C);disp(A+B)+C);(A+B)+Cdisp(A+(B+C);A+(B+C)disp(验证交换律：A+B+C=B+C+A=A+C+B);disp(A+B+C);A+B+Cdisp(B+C+A);B+C+Adisp(A+C+B);A+C+B A+B+C
4、ans = -4 2 22 15A-B+Cans = -16 12 -2 19验证结合律：(A+B)+C=A+(B+C)(A+B)+Cans = -4 2 22 15A+(B+C)ans = -4 2 22 15验证交换律：A+B+C=B+C+A=A+C+BA+B+Cans = -4 2 22 15B+C+Aans = -4 2 22 15A+C+Bans = -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 0L =
5、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.46122.73248. 课本第三章第16题.解：function volume,area=torus(a,b)% a = inner diameter% b = outer diameter
6、volume=pi2/4*(a+b).*(b-a).2;area=pi2*(b.2-a.2);end在MATLAB中调用volume,area=torus(3,4)volume = 17.2718area = 69.08729. 课本第四章第2题.解：a=2:4;b=10,24,24;c=12,48,100;for k=1:3disp(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的数值分析判断是否显示实部、虚部if(delta0
7、)disp(complex value:real and image part:)disp(real(x1),imag(x2);elsedisp(real value:)disp(x1,x2);endend NO.1real value: -2 -3NO.2real value: -4 -4NO.3complex value:real and image part: -3 -4 10. 课本第四章第20题.解：a.t0=0;x0=-10;y0=144;for t = 0:0.01:4x=5*t-10;y=25*t2-120*t+144;if(x2+y2)(x02+y02)t0=t;x0=x;y
8、0=y;endenddisp(time:)disp(t0);disp(Minimium distance is:);d=sqrt(x02+y02);disp(d); time: 2.2300Minimium 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.2300Minimiom distance is: 1.35
9、81 3思考题1. MATLAB中，数组与矩阵在表示与应用上有哪些区别。答：在表示上一维数组相当于向量，二维数组相当于矩阵，所以矩阵是数组的子集。在应用上矩阵运算按线性变换定义，使用通常符号；数组运算按对应元素运算定义，使用点运算符。实验04051002 MATLAB科学计算及绘图操作成绩报告成绩1实验目的1) 熟悉MATLAB所提供的常用数值计算的函数(方程（组）的求解、插值、拟合)；2) 掌握MATLAB二维图形绘制命令及其图形控制（plot、loglog、contour、polar等）；3) 熟悉MATLAB三维图形绘制命令及其图形控制(mesh、surf等)。2实验内容1. 课本第五章
10、第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(A2+B2);d=abs(d);plot(t,d);xlabel(time(hour);ylabel(distance from boundary); 2. 课本第五章第22题.解：x=1990,1991,1992,1993,1994;
11、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=AbA*solution solution = 3.0000 5.0000 -2.0000b.A=6
12、-3 4;12 5 -7;-5 2 6;b=41;-26;14;solution=AbA*solution solution = 2 -3 54. 课本第六章第8题.解：a.A=-2,1;-2,1;b=-5;-5.00001;if(det(A)=0)if(rank(A)=rank(A b)disp(此方程组无解！);endend 此方程组无解！ 故此方程组无解！b.A=1,5,-1,6;2,-1,1,-2;-1,4,-1,3;3,-7,-2,1;det(A) ans = 80.0000 b=19;7;20;75;Ab ans = 17.0000 8.1250 -47.6250 -14.3750
