第9讲MATLAB在电路中的应用.pptVIP

利用MATLAB进行 信号与系统分析 ;信号的MATLAB表示 ;一、基本信号的MATLAB表示 ;%decaying exponential t=0:001:10; A=1; a=-0.4; ft=A*exp(a*t); plot(t,ft) t=0:0.1:10; A=1; a=-0.4; ft=A*exp(a*t); stem(t,ft);% rectpuls t=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:0.001:3; ft=tripuls(t,4,0.5); plot(t,ft) ft=tripuls(t,4,1); ;% unit impuls sequence ;% unit step sequence;二、信号基本运算的MATLAB实现;已知三角波f(t)，用MATLAB画出的f(2t)和f(2-2t) 波形;二、信号基本运算的MATLAB实现;二、信号基本运算的MATLAB实现;二、信号基本运算的MATLAB实现;三角波f(t)微分与积分的波形;利用MATLAB进行系统的时域分析;一、连续时间系统零状态响应的求解;二、连续时间系统冲激响应和阶跃响应的求解;三、离散时间系统零状态响应的求解;四、离散时间系统单位脉冲响应的求解;五、离散卷积的计算;例1 求系统 y(t)+2y(t)+100y(t)=10f(t) 的零状态响应，已知f(t)=sin(2pt) u(t)。;例2 求系统 y (t)+2y (t)+100y(t)=10f(t) 的零状态响应，已知f(t) =d (t) 。;例3 分析噪声干扰的信号f[k]=s[k]+d[k]通过M点滑动平均系统的响应，其中s[k]=(2k)0.9k是原始信号，d[k]是噪声。;例3 分析噪声干扰的信号f[k]=s[k]+d[k]通过M点滑动平均系统的响应，其中s[k]=(2k)0.9k是原始信号，d[k]是噪声。;例4 求系统y[k]+3y[k-1]+2y[k-1]=10f [k]的单位脉冲响应。;根据卷积公式： 因此编程的过程为： （1）写出h(t)＝e-0.1t的MATLAB表达式； （2）写出u(t)的MATLAB表达式； （3）利用MATLAB的卷积语句y=conv(u，h)求解 （4）画曲线plot(t,y)。;例5 计算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} 。;例6 方波分解为多次正弦波之和;利用MATLAB进行信号的频域分析;一、周期信号频谱的MATLAB实现;例1 试用MATLAB画出图示周期三角波信号的频谱。;画三角波信号频谱的MATLAB程序;程序运行结果;例2 求周期矩形脉冲的Fourier级数表示式。并用MATLAB求出由前N项Fourier级数系数得出的信号近似波形。;% Gibbs phenomenon;% Gibbs phenomenon;二、用数值积分分析非周期信号频谱;例3 试用数值方法近似计算三角波信号的频谱;例3 试用数值方法近似计算三角波信号的频谱;例3 试用数值方法近似计算三角波信号的频谱;利用MATLAB进行系统频域分析;一、连续系统频响特性的计算;一、连续系统频响特性的计算;一、连续系统频响特性的计算;二、周期信号通过系统的响应;二、周期信号通过系统的响应;二、周期信号通过系统的响应;三、离散系统频率响应的计算;b=[1];;三、离散系统频响特性的计算;利用MATLAB进行连续系统的复频域分析;一、部分分式展开的MATLAB实现;二、H(s)的零极点与系统特性的MATLAB计算;例1 用部分分式展开法求F(s)的反变换。;例2 用部分分式展开法求F(s)的反变换。;例2 用部分分式展开法求F(s)的反变换。;例3 试画出系统 的零极点分布图，求其单位冲激响应h(t)和频率响应H(j?)，并判断系统是否稳定。; 利用MATLAB进行离散系统的Z域分析;一、部分分式展开的MATLAB实现;二、H(z)的零极点与系统特性的MATLAB计算;例1 将F(z)用部分分式展开。;例2 试画出系统 的零极点分布图，求其单位冲激响应h[k]和频率响应H(ejW) 。;运行结果 ;利用MATLAB进行系统状态变量分析