第一部分 MATLAB基础.doc

《第一部分 MATLAB基础.doc》由会员分享，可在线阅读，更多相关《第一部分 MATLAB基础.doc（14页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第一部分 MATLAB基础一.运算符矩阵运算：+ - * / 一般运算：.+ ._ .* ./二.矩阵1 2 3 4 5 ;1 2 3 4 5;1 2 3 4 5; 1 2 3 4 5;第一个数:间隔:最后一个数;1:1:5-1 2 3 4 5Zeros(n,n)Ones()Eyes()三.画图连续：Plot(x,y);离散：Stem(x,y);stairsXlabelYlabelTitle()Axis()Figure()Subplot四.函数Sin,cos,tanExp,real,imag,abs,length连续信号：Sign(符号信号),sinc(抽样信号),rectpuls(矩形脉冲信

2、号),tripuls(三角波脉冲信号),square(周期性矩形脉冲),sawtooth(周期性三角波信号)五.离散时间信号1正弦序列离散正弦序列的MATLAB表示与连续信号类似，只不过是用stem函数而不是用plot函数来画出序列的波形。下面就是正弦序列的MATLAB源程序。程序运行结果如图1.19所示。%正弦序列实现程序k=0:39;fk=sin(pi/6*k);stem(k,fk)图1.19 正弦序列波形2指数序列离散指数序列的一般形式为，可用MATLAB中的数组幂运算（即点幂运算）c*来实现。下面为用MATLAB编写绘制离散时间实指数序列波形的函数。function dszsu(c,a

3、,k1,k2)%c：指数序列的幅度%a：指数序列的底数%k1：绘制序列的起始序号%k2：绘制序列的终止序号k=k1:k2;x=c*(a.k);stem(k,x,filled)hold onplot(k1,k2,0,0)hold off利用上述函数，实现实指数波形MATLAB程序如下（其中值分别为）。%离散时间实指数序列实现程序subplot 221;dszsu(1,5/4,0,20);xlabel(k);title(f1k);subplot 222dszsu(1,3/4,0,20);xlabel(k);title(f2k);subplot 223;dszsu(1,-5/4,0,20);xlab

4、el(k);title(f3k);subplot 224;dszsu(1,-3/4,0,20);xlabel(k);title(f4k);程序运行结果如图1.20所示。如图可知，对于离散时间实指数序列，当的绝对值大于1时，序列为随时间发散的序列，当的绝对值小于1时，序列为随时间收敛的序列。同时可见，当的值小于零时，其波形在增长或衰减的同时，还交替地改变序列值的符号。图1.20 不同底数的实指数序列对于离散时间虚指数序列，可用通过调用下列绘制虚指数序列时域波形的MATLAB函数。function=dxzsu(n1,n2,w)%n1:绘制波形的虚指数序列的起始时间序号%n2：绘制波形的虚指数序列的

5、终止时间序号%w：虚指数序列的角频率k=n1:n2;f=exp(i*w*k);Xr=real(f)Xi=imag(f)Xa=abs(f)Xn=angle(f)subplot(2,2,1), stem(k,Xr,filled),title(实部);subplot(2,2,3), stem(k,Xi,filled),title(虚部);subplot(2,2,2), stem(k,Xa,filled),title(模);subplot(2,2,4), stem(k,Xn,filled),title(相角);利用上述函数，实现虚指数波形MATLAB程序如下（其中虚指数分别为）%离散时间虚指数实现程序

6、figure(1);dxzsu(0,20,pi/4);figure(2);dxzsu(0,20,2);程序运行结果如图1.21（a）、（b）所示。由图可见，只有当虚指数序列的角频率满足为有理数时，信号的实部和虚部和相角都为周期序列，否则为非周期序列。 （a）波形 (b) 波形图1.21 虚指数序列波形对于复指数序列，其一般形式为可以通过调用下面绘制复指数序列时域波形的MATLAB函数。function dfzsu(n1,n2,r,w)%n1:绘制波形的虚指数序列的起始时间序号%n2：绘制波形的虚指数序列的终止时间序号%w：虚指数序列的角频率%r: 指数序列的底数k=n1:n2;f=(r*exp

7、(i*w).k;Xr=real(f);Xi=imag(f);Xa=abs(f);Xn=angle(f);subplot(2,2,1), stem(k,Xr,filled),title(实部);subplot(2,2,3), stem(k,Xi,filled),title(虚部);subplot(2,2,2), stem(k,Xa,filled),title(模);subplot(2,2,4), stem(k,Xn,filled),title(相角);利用上述函数，实现复指数序列波形MATLAB程序如下。%复指数序列实现程序(r1)figure(1);dfzsu(0,20,1.2,pi/4);%

8、复指数序列实现程序(0r1时，复指数序列的实部和虚部分别为幅度按指数增长的正弦序列;当0r1 (b) 0r=0; %序列的起点为ns，终点为nf，在n=n0点处生成单位阶跃。（3）两个信号相加的生成函数sigadd.mfunction y,n=sigadd(x1,n1,x2,n2)n=min(min(n1),min(n2):max(max(n1),max(n2)y1=zeros(1,length(n);y2=y1;y1=(find(n=min(n1)&(n=min(n2)&(n=min(n1)&(n=min(n2)&(n=max(n2)=1)=x2;y=y1.*y2;（5）序列移位的生成函数s

9、igshift.mfunction y,n=sigshift(x,m,n0)n=m+n0;y=x（6）序列翻折的生成函数sigfold.mfunction y,n=sigfold(x,n)y=fliplr(x);n=-fliplr(n)（7）奇偶综合函数evenodd.m此函数可以将任一给定的序列x(n)分解为xe(n) 和xo(n)两部分。function xe,xo,m=evenodd(x,n)if (imag(x)=0) error(x is not a real sequence)； endm=-fliplr(n);m1=min(m,n);m2=max(m,n);m=m1:m2;nm=

10、n(1)-m(1);n1=1:length(n);x1=zeros(1,length(m);x1(n1+nm)=x;x=x1;xe=0.5*(x+fliplr(x); xo=0.5*(x-fliplr(x);(8）求卷积和求卷积和直接采用MATLAB中的函数conv，即y=conv(x,h)；它默认序列从n=0开始。但是如果序列是从一负值开始，即其中nx10或nh1=min(k1)&(k=min(k2)&(k=min(k1)&(k=min(k2)&(k=0;subplot 221stem(n,x)xlabel(n);ylabel(x(n);title(Step Sequence);grid on;%Decomposition of the Sequencexeven,xodd,m=evenodd(x,n);subplot 223stem(m,xeven);xlabel(m);ylabel(x even(n);title(Even Part);grid on;subplot 224stem(m,xodd);xlabel(m);ylabel(x odd(n);title(Odd Part);grid on;程序运行结果如图1.30所示。