MATLAB应用实践分析.doc

浦 江 学 院 实 习 报 告 （ 2013 — 2014 学年 第 一 学期） 课程名称： MATLAB应用实践 班 级： 学 号： 姓 名： 指导教师： 2013 年 12 月 实验1 练习1、 练习2、 1、源程序： clear all n=input(n=); s=0; for i=1:2:2*n-1 s=s+i; end s 结果： 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 结果： 练习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 结果： 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 结果： 源程序： 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]; 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) 结果： 例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) 结果： 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=[5 6 0 0 0 1 5 6 0 0 0 1 5 6 0 0 0 1 5 6 0 0 0 1 5]; B=[1 0 0 0 1]; R_A=rank(A) X=A\B 例5-2源程序： A=[1 1 -3 -1;3 -1 -3 4;1 5 -9 -8]; B=[1 4 0]; X=A\B 实验二 单位脉冲源程序： n1=-4; n2=4; n0=0; n=n1:n2; x=[n==n0]; stem(n,x,filled); title(单位脉冲序列); xlabel(时间(n)); ylabel(序列值x(n)); 单位阶跃源程序： n1=-4; n2=4; n0=0; n=n1:n2; x=[nn0]; stem(n,x,filled); title(单位阶跃序列); xlabel(时间(n)); ylabel(序列值x(n)); （3）矩形源程序： n=1:40; x=sign(sign(8-n)+1); stem(n,x,filled); title(矩形序列); xlabel(时间(n)); ylabel(序列值x(n)); 正弦源程序： fs=40;n=10; n=0:n-1; x=sin(pi*n/5+pi/3); stem(n,x,filled); title(正弦波序列); xlabel(时间(n)); yl