matlab任务书解读

1、MATLAB程设计任务书（-）一、名称： MATLAB 编程简介二、目的熟悉 MATLAB 编程环境，掌握 Help 命令、基本的变量类型、矩阵的基本运算、基本的 绘图函数和 M-file 的建立 。三、内容例题1、Help 命令helpcoshelpplothelpabshelpexphelp+例 2、 变量和矩阵运算(1) Matrix The basic variable type M=3M=1 2 6M=1 2 6; 4 6 7M13=M(1,3)size(M)(2) The Colon Operator ( : ) %Creating Array and Vector % v = s

2、tart: skip: end x1=0:2:10x2=0:1:10 (or x=0:10) t=-1:0.2:1%Accessing MatrixA=1 2 3;4 5 6;7 8 9A(2:3,1:2) x2(4:8)(3) Matrix Operations (A B)A=2 3 4; 6 9 8B=1 2 3; 5 8 7C1=A+BC2=A-BC3=A-4(4) Matrix Operations (A*B A.*B)% A*BA=2 3 4; 6 9 8B=1 2; 3 5; 8 7A*B% A.*BA=2 3 4; 6 9 8B=1 2 3 ;5 8 7A.*B(5) Matri

3、x Operations (B/A ,AC, B./A ,A.B) % B/A B*inv(A)% AC inv(A)*C%B./A B(i,j)/A(i,j)%A.B B(i,j)/A(i,j) Matrix Operatio ns ( A and .A) % a Operati onA=1 2 3; 4 5 6; 7 8 9b=AA2% .A OperationA=1 2 3; 4 5 6; 7 8 9 b=A.A2(7) Matrix Operations ( A and A.) % A共轭转置a=1+2i3+4i; 3+2i5+5ia% A .非共轭转置a.例 3 、 绘图函数 plo

4、t(x,y) ， stem(k,y)% plot(x,y)x=0:0.01:2;y=sin(2*pi*x);plot(x,y)% stem(k,y)k=0:50;y=exp(-0.1*k);stem(k,y)例 4 、M file% y(t)=sin(2t) + sin(5t)-2pi t （s+5 Is2 十5S + 25）12、已知系统函数为H s二- 1，试分别画出:=0,14,1,2时系统的零极点图。s+2o（s+14如果系统是稳定的，画出系统的幅度响应曲线。s 23、 已知F s = 3 笃，画出该系统的零极点分布图，求出系统的冲激响应、s3+2s2+2s + 1阶跃响应和频率响应。

5、4、已知系统函数H（s）的零极点分布如下图所示，讨论它们分别是哪种滤波器（低通、 根据零极点的位置设零极点的数值）带通、高通、带阻）。（提示:（2）/+ -3例2、已知一离散因果LTI系统的系统函数为： Wo0爲二十。讨，该系统的零极点说明：先将系统改写为：z2亠2z亠1H（zz3 -0,z2 -0.005z 0.3，然后用血p函数求系统的零极点。MATLAB 程序:% zeros and poles of H(z)b=1 2 1;a=1 -0.5 -0.005 0.3;r,p k=tf2zp(b,a)% r,p为列向量，求零极点zpla ne(b,a);运行结果:r =-1-10.5198

6、+ 0.5346i0.5198 - 0.5346i-0.5396k =1求出的零点为：z=1(2 阶)，z = 0 ;极点为：z = 0.5198 _0.5346i , z =-0.5396零极点图如下：mReal Part例3、利用MATLAB画出二_0.9时F(ej )：，的幅度频谱。e jMATLAB 程序:% amplitude freque ncy spectrum of a seque neeb=1;a 1=1 -0.9;a2=1 0.9;w=linspace(0,2*pi,512);%线性均分0到2pi的间隔，共512个点h1=freqz(b,a1,w);h2=freqz(b,a

7、2,w);plot(w/pi,abs(h1),w/pi,abs(h2),k:);xlabel(Omega八pi);lege nd(alpha=0.9,alpha=-0.9);运行结果:121086420IIIIiii00.20.40.60.811.21.41.61.82例4、已知一离散因果LTI系统的系统函数为:H(z)二z2 2z 1z3 - 0.5z2 - 0.005z0.3利用MATLAB画出该系统的零极点分布图，求系统的单位样值响应和幅频响应，并判断系统的稳 定性。MATLAB 程序:% Impulse response,amplitude frequency responseand stability analysis of %LTI H(z)b=0 1 2 1;a=1 -0.5 -0.005 0.3;figure(1);zpla ne(b,a);nu m=0 1 2 1;den=1 -0.5 -0.005 0.3;h=impz( nu m,de n)figure;stem(h);xlabel(k);ylabel(h(k);title(Impulse Resp on se);