1、pi 表示圆周率，inf 表示无穷大，NaN（Not a Number）表示不定量，如0/0。3表达式数值 MATLAB用常用常规的十的十进制表示数制表示数值 用用i或或j作作为后后缀来表示复数的虚部来表示复数的虚部 例1.235e5表示1.235105，x=2+3j abs（x）求复数x的模 angle（x）求复数x的相角（弧度）real（x）求复数x的实部 imag（x）求复数x的虚部 conj（x）求复数x的共轭3表达式运算符号算数运算符算数运算符+加 -减 *乘 /除 乘方 矩阵的复共轭转置3表达式运算符号逻辑运算符逻辑运算符 A&B逻辑与（and）A|B逻辑或（or）A 逻辑非（no
2、t）值为0时表示逻辑假（F），其它任何非零值表示逻辑真。3表达式运算符号关系运算符关系运算符 A B 大于 A=B 大于等于 A=B 等于 A=B不等于3表达式运算符号冒号运算符冒号运算符表达式 1:10 表示产生一个行向量，它的值为 1 2 3 4 5 6 7 8 9 10 表达式 10:-2:1 表示产生一个递减的行向量，它的值为 10 8 6 4 2 4数组数组的构造用冒号用冒号:产生数生数组例 x=2:5 产生一个数组，它的值为 x（1）=2,x（2）=3,x（3）=4,x（4）=5x=linspace（0,2,11）将区间0，2均匀抽样11点作为数组x给2维数数组赋值时，用分号分号表
3、示一行的结束，如：z=1 2;3 4。用用linspace产生数生数组4数组MATLAB 提供了一些产生基本矩阵的函数 zeros 产生矩阵元素全为0的矩阵 ones 产生矩阵元素全为1的矩阵 rand 产生（0,1）均匀分布随机数矩阵 randn 产生正态分布随机数矩阵4数组数组的运算数数组和一个和一个标量相加或相乘量相加或相乘 例 y=x-1 z=3*x 2个数个数组的的对应元素相乘除元素相乘除 .*./例 z=x.*y 确定数确定数组大小的函数大小的函数 size（A）返回值数组A的行数和列数（二维）length（B）确定数组B的元素个数（一维）5.函数文件M文件的第一行包含functi
4、on功能:建立一个函数，可以同MATLAB的库函数一样使用.例：编一个绘制图示波形的函数。function y=tri（t）y=abs（t）eps s=s+1/（n*n）;n=n+1;end例:计算 的值，且误差小于10-6fprintf（s=%.5fn,s）8.plot函数-绘图函数（continuous）t=linspace（0,4*pi,512）;plot（t,sin（t）,t,cos（t）,-.）;title（my figure）;legend（sin（t）,cos（t）;9.stem函数-绘图函数（discrete）k=0:39;stem（k,cos（0.9*pi*k）;title（
5、cos（0.9pik）;Cos（0.9k）波形基本信号的基本信号的MATLAB表示表示 （指数信号Aeat,指数序列ak,抽样函数Sa（t）,正弦型信号,矩形脉冲信号,三角脉冲信号 信号基本运算的信号基本运算的MATLAB实现（尺度变换、翻转、时移、相加、相乘、差分与求和、微分与积分）信号的信号的MATLABMATLAB表示表示基本信号的MATLAB表示指数信号Aeaty=A*exp（a*t）;指数序列ak幂运算a.k实现正弦型信号内部函数cos（）和sin（）抽样函数Sa（t）sinc（t）矩形脉冲信号y=rectpuls（t,width）三角波脉冲信号y=tripuls（t,width,s
6、kew）%decaying exponential t=0:001:10;A=1;a=-0.4;ft=A*exp（a*t）;plot（t,ft）t=0:0.1:stem（t,ft）%rectpulst=0:0.001:4;T=1;ft=rectpuls（t-2*T,T）;plot（t,ft）axis（0,4,-0.5,1.5）%tripuls t=-3:3;ft=tripuls（t,4,0.5）;plot（t,ft）ft=tripuls（t,4,1）;%unit impuls sequence k=-50:50;delta=zeros（1,50）,1,zeros（1,50）;stem（k,de
7、lta）function f,k=impseq（k0,k1,k2）%产生 fk=delta（k-k0）；k1=k=k2 k=k1:k2;f=（k-k0）=0;k0=0;k1=-50;k2=50;f,k=impseq（k0,k1,k2）;stem（k,f）%unit step sequence k=-50:uk=zeros（1,50）,ones（1,51）;stem（k,uk）function f,k=stepseq（k0,k1,k2）%产生 fk=u（k-k0）;k1=k=0;f,k=stepseq（k0,k1,k2）;stem（k,f）信号基本运算的MATLAB实现 t=-3:ft1=tri
8、puls（2*t,4,0.5）;subplot（2,1,1）plot（t,ft1）title（f（2t）ft2=tripuls（2-2*t）,4,0.5）;subplot（2,1,2）plot（t,ft2）title（f（2-2t）信号的尺度变换、翻转、时移（平移）信号的尺度变换、翻转、时移（平移）已知三角波f（t），用MATLAB画出的f（2t）和f（2-2t）波形信号的相加与相乘信号的相加与相乘相加用算术运算符“+”实现 相乘用数组运算符“.*”实现 例画信号Aeatsin（w0t+f）的波形 t=0:8;w0=2*pi;phi=0;ft1=A*exp（a*t）.*sin（w0*t+phi
9、）;plot（t,ft1）离散序列的差分与求和离散序列的差分与求和差分y=diff（f）;求和y=sum（f（k1:k2）;连续信号的微分与积分连续信号的微分与积分微分 y=diff（f）/h;h为数值计算所取时间间隔定积分 quad（function_name,a,b）;function_name为被积函数名，a和b指定积分区间。例已知三角波f（t），画出其微分与积分的波形%differentiation h=0.001;t=-3:h:y1=diff（f2_2（t）/h;plot（t（1:length（t）-1）,y1）%integrationt=-3:for x=1:length（t）y2
10、（x）=quad（f2_2,-3,t（x）;endplot（t,y2）function yt=f2_2（t）yt=tripuls（t,4,0.5）;三角波f（t）微分与积分的波形利用利用MATLABMATLAB进行系统的时域分析进行系统的时域分析连续时间系系统零状零状态响响应的求解的求解 连续时间系系统冲激响冲激响应和和阶跃响响应的求解的求解 离散离散时间系系统零状零状态响响应的求解的求解 离散离散时间系系统单位脉冲响位脉冲响应的求解的求解 离散卷离散卷积的的计算算t 表示计算系统响应的抽样点向量a=a3,a2,a1,a0;b=b3,b2,b1,b0;sys=tf（b,a）1.连续时间系统零状
11、态响应的求解y=lsim（sys,f,t）sys=tf（b,a）b和a分别为微分方程右端和左端各项的系数向量f 是系统输入信号向量，sys 是LTI系统模型，借助tf函数获得连续时间系统冲激响应可用impulse函数直接求出，其调用形式为2.连续时间系统冲激响应和阶跃响应的求解y=impulse（sys,t）连续时间系统阶跃响应可用step函数直接求出，其调用形式为y=step（sys,t）t 表示计算系统响应的抽样点向量 sys 是LTI系统模型3.离散时间系统零状态响应的求解b,a 分别是差分方程左、右端的系数向量b=b0,b1,b2,bM;a=a0,a1,a2,aN;可用MATLAB表示
12、为y=filter（b,a,f）f 表示输入序列,y 表示输出序列4.离散时间系统单位脉冲响应的求解b,a 分别是差分方程左、右端的系数向量 k 表示输出序列的取值范围 h 就是单位脉冲响应h=impz（b,a,k）5.离散卷积的计算例如：（s3+2s+3）（s2+3s+2）可用下面MATLAB语句求出 a=1,0,2,3;b=1,3,2;c=conv（a,b）c=conv（a,b）式中a,b为待卷积两序列的向量表示，c是卷积结果。conv函数也可用于计算两个多项式的积例1 求系统y”（t）+2y（t）+100y（t）=10f（t）的零状态响应，已知f（t）=（sin2t）u（t）。%prog
13、ram3_1微分方程求解ts=0;te=5;dt=0.01;sys=tf（1,1 2 100）;t=ts:dt:te;f=10*sin（2*pi*t）;y=lsim（sys,f,t）;plot（t,y）;xlabel（Time（sec）ylabel（y（t）例2 求系统y”（t）+2y（t）+100y（t）=10f（t）的零状态响应，已知f（t）=d（t）。%program3_2连续时间系统的冲激响应ts=0;sys=tf（10,1 2 100）;y=impulse（sys,t）;xlabel（Time（sec）ylabel（h（t）例3 分析噪声干扰的信号fk=sk+dk通过M点滑动平均系统
14、的响应，其中sk=（2k）0.9k是原始信号，dk是噪声。R=51;d=rand（1,R）-0.5;k=0:R-1;s=2*k.*（0.9.k）;f=s+d;figure（1）;plot（k,d,r-.,k,s,b-,k,f,g-）;M=5;b=ones（M,1）/M;a=1;y=filter（b,a,f）;figure（2）;plot（k,s,b-,k,y,r-）;噪声干扰信号fk=sk+dk通过M点滑动平均系统的响应例4 求系统yk+3yk-1+2yk-1=10f k的单位脉冲响应。%program 3_4 离散系统的单位脉冲响应 k=0:a=1 3 2;b=1;h=impz（b,a,k）;s
