信息工程学院电子科学与技术系徐利民.pptVIP

二、信号与系统的数值化处理根据此方程求出差分方程的初始条件：求解差分方程需要初始条件:y(-1),y(-2),怎么办？上例解得y(-1)=1；y(-2)=0.9。编程如下：二、信号与系统的数值化处理Ts=0.1;%时间间隔n=1:4/Ts+2;%坐标1,2代表初始条件坐标-2,-1F(n)=1;F(1)=0;F(2)=0;%激励信号初始条件Y(1)=0.9;Y(2)=1;%输入初始条件y(-2)=0.9;y(-1)=1fork=3:length(n)%迭代求解Y(k)=(13*F(k)-10*F(k-1)-100*Y(k-2)+240*Y(k-1))/143;endn1=(n-3)*Ts;subplot(2,1,1);plot(n1,Y);xlabel(t);ylabel(y(t));title(全响应(数值计算));gridon;%下图用实际响应进行对比subplot(2,1,2);ezplot(-exp(-3*t)+exp(-t)+1,[0,4]);xlabel(t);ylabel(y(t));title(全响应);gridon;二、信号与系统的数值化处理运行结果：二、信号与系统的数值化处理思考：1、如果连续系统给出的是初始条件y(0+)；y‘(0+)，那么求离散系统初始条件的方程组应该是什么？2、基于状态方程则只用一阶微分式就足够，请思考原理。3、如果计算微分用起点和终点斜率的平均值计算，则可以是计算精度由二阶提高到三阶(梯形预估-校正法)。二、信号与系统的数值化处理4、如果在每一时间区间内多取几个点的斜率值，将它们进行线性组合，则可能构造出精度更高的算法（龙格-库塔法）这些知识很重要，在MATLAB和SIMULINK中多有用到。请一定理解其原理！！！二、信号与系统的数值化处理例9?3：求解微分方程初值问题的数值解，求解范围为区间[0,0.5]。，程序如下：fun=inline(-2*y+2*x^2+2*x,x,y);[x,y]=ode23(fun,[0,0.5],1);plot(x,y,o-)运行结果：二、信号与系统的数值化处理SolverODE类型特点说明ode45非刚性单步算法；4、5阶Runge-Kutta方程；累计截断误差达(Δx)3大部分场合的首选算法ode23非刚性单步算法；2、3阶Runge-Kutta方程；累计截断误差达(Δx)3使用于精度较低的情形ode113非刚性多步法；Adams算法；高低精度均可到10-3~10-6计算时间比ode45短ode23t适度刚性采用梯形算法适度刚性情形ode15s刚性多步法；Gears反向数值微分；精度中等若ode45失效时，可尝试使用ode23s刚性单步法；2阶Rosebrock算法；低精度当精度较低时，计算时间比ode15s短ode23tb刚性梯形算法；低精度当精度较低时，计算时间比ode15s短二、信号与系统的数值化处理SIMULINK中首先要设置的参数：二、信号与系统的数值化处理用MATLAB进行信号与系统实践1234MATLAB的使用信号与系统的数值化处理工程应用基石：函数与工具箱应用课题的研究、分析与设计函数functiony=linspace(d1,d2,n)%LINSPACELinearlyspacedvector.%LINSPACE(X1,X2)generatesarowvectorof100linearly……ifnargin==2n=100;endy=[d1+(0:n-2)*(d2-d1)/(floor(n)-1)d2];函数头输入变量输出变量说明行函数代码，一般要对输出变量赋值H1行函数名，需要与.m文件同名函数1、大量数学函数1）三角函数和双曲函数sin正弦cos余弦tan正切asin反正弦acos反余弦atan反正切sinh双曲正弦cosh双曲余弦tanh双曲正切asinh反双曲正弦acosh反双曲余弦atanh反双曲正切cot余切sec正割csc余割acot反余切asec反正割acsc反余割coth双曲余切sech双曲正割csch双曲余割acoth反双曲余切asech反双曲正割acsch反双曲余割atan2四象限反正切函数2）指数函数exp指数log10常用对数pow22的幂log自然对数l