发布时间 : 星期日 文章Matlab上机实习报告更新完毕开始阅读ece59ac50408763231126edb6f1aff00bed5700f
《MATLAB课程》课程实习报告（一） 实习日期：2016-9-13
1.先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量； (1)代码:
z1=2*sin(85*pi/180)/(1+exp(2)) 截图：
(2)代码： a=-3.0:0.1:3.0;
z2=(exp(0.3*a)-exp(-0.3*a))/2.*sin(a+0.3)+log((0.3+a)/2) 截图：
(3)代码： t=0:0.5:2.5;
z3=t.^2.*(t>=0&t<1)+(t.^2-1).*(t>=1&t<2)+(t.^2-2*t+1).*(t>=2&t<3) 截图：
2.设有矩阵A和B (1)求它们的乘积C
(2)将矩阵的右下角3*2子矩阵赋给D 代码：
A=[1 2 3 4 5;6 7 8 9 10;11 12 13 14 15;16 17 18 19 20;21 22 23 24 25] B=[3 0 16;17 -6 9;0 23 -4;9 7 0;4 13 11] C=A*B F=size(C)
1 / 4
D=C(F(1)-2:F(1),F(2)-1:F(2)) 截图：
3. （1）求之间能被21整除的数的个数。 代码： A=100:999; B=rem(A,21);
C=length(find(B==0)) 截图：
（2）建立一个字符串向量，删除其中大写字母。 代码:
> A='qWeRTyUiopASDfGHjklz'; k=find(A>='A'&A<='Z'); A(k)=[] 截图：
2 / 4
4.解线性方程组 （1）求方程的解；（2）将方程右边向量元素改成0.53，再求解，并比较的变化和解的相对变化。 代码:
A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6]; b=[0.95,0.67,0.52]'; x=inv(A)*b
B=[0.95,0.67,0.53]'; x=inv(A)*B 截图：
5.（1）取出A的前3行构成矩阵B，前两列构成矩阵C，右下角3*2子矩阵构成矩阵D，B与C的乘积构成矩阵E。 代码:
A=[23,10,-0.778,0;41,-45,65,5;32,5,0,32;6,-9.54,54,3.14] B=A(1:3,:) C=A(:,1:2) D=A(2:4,3:4) E=B*C 截图:
3 / 4
（2）分别求E?D、E＆D、E∣D和find（A>=10＆A<25）。 代码：
E
4 / 4