MATLAB实训报告
项目一 MATLAB基本应用
1、 分别建立3×3的单位阵、3×2的零矩阵、每行、每列和两对
角线上的元素之和等于65的魔方矩阵。
2、 计算向量a=(1,2,3),b=(4,5,6)和c=(-3,-6,-3)的混合积a?(b?c)。
3、
?2x1?x2?x3?5用左除求下列方程组的根。??3x1?2x2?2x3?5
?5x?3x?x?1623?1
4、 求出函数s4+2s3-12s2-s+7在s=3处的值。
5、 分别计算
lim1x?0x、x?0x。
lim?1>> syms x; >> limit(1/x,x,0) ans = NaN
>> limit(1/x,x,0,'left') ans = -Inf 6、 求积分?0??x(1?x)2dx。
7、f(x)=ax2+bx+c,对x求二次微分,对a求二次微分。 >> syms x a b c; >> y=a*x^2+b*x+c; >> z=diff(y,'x',2) z = 2*a
>> k=diff(y,'a',2) k = 0
项目二 MATLAB图形系统
1、 画出x0.8×e-x在x∈[0, 4]之间的针状图形。
>> x=0:0.1:4; >> y=(x.^0.8).*exp(-x); >> stem(x,y)
0.40.350.30.250.20.150.10.05000.511.522.533.54
2、 绘制函数y(t)=1-2e-tsin(t) (0≤t≤8)的图形,且在x轴上标注
“Time”,y轴上标注 “Amplitude”,图形的标题为“Decaying Oscillating Exponential”。 >> x=0:1:8; >> figure
>> y=1-2*exp(-x).*sin(x); >> plot(x,y,'r-.')
>> line([0,5],[0,0])
>> xlabel('Time');ylabel('Amplitude'); >> title('Decaying Oscillating Exponential'); >> legend('y=1-2*exp(-x).*sin(x)') >> grid on
Decaying Oscillating Exponential1.4y=1-2*exp(-x).*sin(x)1.21Amplitude0.80.60.40.2001234Time5678
3、在一个图形窗口下绘制两个子图,分别显示下列曲线。
(1)y=sin2xcos3x;(2)y=0.4 x
要求给x轴、y轴加标注,每个子图加标题且两条曲线用不同颜色和线型表示。
>> x=linspace(0,2*pi,100);
>> subplot(2,2,1);plot(x,sin(2*x).*cos(3*x),'g-+'); >> xlabel('x');ylabel('y');title('sin(2*x).*cos(3*x)') >> subplot(2,2,2);plot(x,0.4*x,'r-.*'); >> xlabel('x');ylabel('y');title('0.4*x')