MATLAB应用实践解读.docxVIP

浦江学院 实习报告 2013 — 2014 学年 第 一 学期） 课程名称： MATLAB 应用实践 班 级： 学 号： 姓 名： 指导教师： 2013 年 12 月 实 习 报 告 实验 1 练习 1、 - 1 - 实 习 报 告 - 2 - 实 习 报 告 - 3 - 实 习 报 告 练习 2、 1、源程序： clear all n=input(n=); s=0; for i=1:2:2*n-1 s=s+i; end s - 4 - 实 习 报 告 结果： 2、源程序： function y=f(x) x=input(input:n=); n=length(x); if x0 x1 y=x; else if x=1 x=2 y=2-x; else y=0; end end - 5 - 实 习 报 告 结果： 练习 3、 LX0701 源程序： A=[1,1,1;1,2,3;1,3,6] B=[8,1,6;3,5,7;4,9,2] C=A+B D=A-B - 6 - 实 习 报 告 结果： LX0702 源程序： X=[2 3 4 5 1 2 2 1]； Y= [0 1 1 1 1 0 0 0 1 1 0 0]； Z=X*Y - 7 - 实 习 报 告 结果： 源程序： X=[-1 0 2]; a=2*X ; 结果： 源程序： X=[-1 0 2]; Y=[-2,-1,1]; Z=dot(X,Y) 结果： LX0703 源程序： a=[1 2 3]; b=[4 5 6]; c=cross(a,b) 结果： LX0704 源程序： a=[1,2,3]; b=[4 5 6]; c=[-3,6,-3]; - 8 - 实 习 报 告 x=dot(a,cross(b,c)) 结果： LX0751 源程序： A=[1 2 3; 2 2 1; 3 4 3]; inv(A) 结果： LX0752 源程序： B=[1, 2, 3, 1, 0, 0; 2, 2, 1, 0, 1, 0; 3, 4, 3, 0, 0, 1]; C=rref(B) % 化行最简形 X=C(:, 4:6) 结果： LX0706 源程序： A=[1 2 3; 2 2 1; 3 4 3]; D=det(A) 结果： - 9 - 实 习 报 告 例 3-7 A=sym( ‘[1/x,1/(x+1);1/(x+2),1/(x+3)] ’); B=sym( ‘[x,1;x+2,0] ’) C=B-A D=A\B 结 果 ： LX0710 源程序： syms x y p=expand((x+1)^3) q=expand(sin(x+y)) 结果： LX0714 源程序： A=[1 -2 2 3;-2 4 -1 3;-1 2 0 3;0 6 2 3;2 -6 3 4]; B=rank(A) -10- 实 习 报 告 结果： LX0715 源程序： a1=[1 -2 2 3]; a2=[-2 4 -1 3]; a3=[-1 2 0 3]; a4=[0 6 2 3]; a5=[2 -6 3 4]; A=[a1 a2 a3 a4 a5] format rat % 以有理格式输出 B=rref(A) % 求 A 的行最简形 LX0716 源程序： A=[56000 15600 01560 00156 00015]; B=[1 0 0 0 1]; R_A=rank(A) X=A\B -11- 实 习 报 告 例 5-2 源程序： A=[1 1 -3 -1;3 -1 -3 4;1 5 -9 -8]; B=[1 4 0]; X=A\B 实验二 （ 1） 单位脉冲源程序： n1=-4; n2=4; n0=0; n=n1:n2; x=[n==n0]; stem(n,x,filled); title( 单位脉冲序列 ); xlabel( 时间 (n)); ylabel( 序列值 x(n)); -12- 实 习 报 告 （ 2） 单位阶跃源程序： n1=-4; n2=4; n0=0; n=n1:n2; x=[nn0]; stem(n,x,filled); title( 单位阶跃序列 ); xlabel( 时间 (n)); ylabel( 序列值 x(n)); -13- 实 习 报 告 （ 3）矩形源程序： n=1:40; x=sign(sign(8-n)+1); stem(n,x,filled); title( 矩形序列 ); xlabel( 时间 (n)); ylabel( 序列值 x(n)); -14- 实 习 报 告 （ 3） 正弦源程序： fs=40;n=10; n=0:n-1; x=sin(pi*n/5+pi/3); stem(n,x,filled); title( 正弦