[信息与通信]第八周 MATLAB在通信系统中的应用2.pptVIP

二、信号基本运算的MATLAB实现 2. 信号的相加与相乘 相加用算术运算符“+”实现 相乘用数组运算符“.*”实现 例:画信号Aeatcos(w0t+f)的波形 t=0:0.001:8; A=1; a=-0.4; w0=2*pi;phi=0; ft1=A*exp(a*t).*sin(w0*t+phi); plot(t,ft1) 二、信号基本运算的MATLAB实现 3. 离散序列的差分与求和  连续信号的微分与积分 差分 y=diff(f); 求和 y=sum(f(k1:k2)); 微分 y=diff(f)/h; h为数值计算所取时间间隔 定积分 quad(‘function_name’,a,b); function_name为被积函数名，a和b指定积分区间。 二、信号基本运算的MATLAB实现 3. 离散序列的差分与求和 连续信号的微分与积分 例:已知三角波f(t)，画出其微分与积分的波形 %differentiation h=0.001;t= -3:h:3; y1=diff(ft2)*1/h; plot(t(1:length(t)-1),y1) %integration t= -3:0.1:3; for x=1:length(t) y2(x)=quad(ft2, -3,t(x)); end plot(t,y2) 三角波f(t)微分与积分的波形 利用MATLAB进行系统的时域分析 连续时间系统零状态响应的求解 连续时间系统冲激响应和阶跃响应的求解 离散时间系统零状态响应的求解 离散时间系统单位脉冲响应的求解 离散卷积的计算 一、连续时间系统零状态响应的求解 t 表示计算系统响应的抽样点向量 f是系统输入信号向量 sys是LTI系统模型，借助tf函数获得 a=[a3, a2, a1, a0]; b=[b3, b2, b1, b0]; sys=tf(b,a) y=lsim(sys,f,t) sys=tf(b,a) b和a分别为微分方程右端和左端各项的系 二、连续时间系统冲激响应和阶跃响应的求解 连续时间系统冲激响应可用impulse函数直接求出，其调用形式为 y=impulse(sys,t) 连续时间系统阶跃响应可用step函数直接求出，其调用形式为 y=step(sys,t) t 表示计算系统响应的抽样点向量 sys 是LTI系统模型 三、离散时间系统零状态响应的求解 b , a 分别是差分方程左、右端的系数向量 b=[b0,b1,b2,,bM]; a=[a0,a1,a2, ,aN]; 可用MATLAB表示为 y=filter(b,a,f) f 表示输入序列, y 表示输出序列 四、离散时间系统单位脉冲响应的求解 b, a 分别是差分方程左、右端的系数向量 k 表示输出序列的取值范围 h 就是单位脉冲响应 h=impz(b,a,k) 五、离散卷积的计算 例：(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)=sin(2pt) u(t)。 %微分方程求解 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) 。 %连续时间系统的冲激响应 ts=0;te=5;dt=0.01; sys=tf([10],[1 2 100]); t=ts:dt:te; y=impulse(sys,t); plot(t,y); xlabel(Time(sec)) ylabel(h(t)) 例3 计算x[k]* y[k]并画出卷积结果，已知x[k]={1,2,3,4; k=0,1,2,3}，y[k]={1,1,1,1,1; k=0,1,2,3,4} 。 x=[1,2,3,4]; y=[1,1,1,1,1]; z=conv(x,y); N=length(z); stem(