§74常系数线性差分方程的求解.doc

§7.4 常系数线性差分方程 的求解 北京邮电大学电子工程学院 陈智娇 解法 第 页 1.迭代法 2.时域经典法：齐次解+特解（类比连续系统） 3.零输入响应+零状态响应（类比连续系统） 利用卷积和求系统的零状态响应 4. z变换法?反变换?y(n) （第八章） BUPT 尹霄丽 一．迭代法 第 页 解差分方程的基础方法,差分方程本身是一种递推关系， 例7-4-1 ? ? 0.9 ? 1? ? ? ? 1? 0, 已知 ，且 求解方程。 y n ? y n ? ? u n y ? ? ? ? ? ? n ? y ? y ? ? ? 0 0 0.9 1 1 1 ? ? ? ? n ? y ? y ? ? 1 1 0.9 0 1 1.9 ? ? ? ? n ? y ? y ? ? 2 2 0.9 1 1 2.71 ? ? ? ? n ? y ? y ? ? 3 3 0.9 2 1 3.439 ? ? ? ? 由递推关系,可得输出值：? ? y n 1, 1.9, 2.71, 3.439, ? ? ? ? ? ? ? n 0 ? ? ? 但得不到输出序列y n 的解析式,适合利用计算机求解。 BUPT 尹霄丽 编程实现 第 页 % iteration_demo.m % zero state clear N=50;a=0.9; n=1:N; x=ones(1,N);y(1)=x(1); for ii=2:N y(ii)=a*y(ii-1)+x(ii); end hh=stem(n-1,y);set(hh,'LineWidth',2); xlabel('n');ylabel('y(n)'); BUPT 尹霄丽 结果显示 第 页 10 9 8 7 6 5 4 3 2 1 0 0 5 10 15 20 25 30 35 40 45 50 n BUPT 尹霄丽 例7-4-2（迭代法示例） 第 页 某人因购房向银行贷款，当前的贷款月利率为I，还款 期限为N个月，采用等额均还方式，每月还款金额为R。 (1)请按照上述规定建立差分方程式描述还款过程。 (2)若贷款金额为70万元，当前贷款年利率为6.55% ? ? ? 1 （5年以上） I= 6.55% ，月还款额为5000。请计算 ? ? ? 12 ? 还款期限，并求出还款金额总数和利息总数。 解： 在7.3节已列出差分方程 y?n? ? ?1? I ? y?n ? 1?? x?n? ?n ? 1? ?0? 700000 y ? BUPT 尹霄丽 MATLAB程序 第 页 interest=6.55*0.8/12/100; y(1)=700000; x=5000; for ii=2:12*30 y(ii)=y(ii-1)*(1+interest)-x; if y(ii)<=0 x_last=y(ii-1); y(ii)=0; break; end end x_total=x*(ii-1)+x_last; BUPT 尹霄丽 还款曲线，年利息6.55%*0.8 第 页 70 218months,total number is:108万 元 60 50 40 30 20 10 0 0 50 100 150 200 250 month 贷款70万，每月还款5000元，还款218个月，本利共还108万 BUPT 尹霄丽 还款曲线，年利息6.55%*0.8 第 页 120 360months,total number is:90万 元 100 80 60 40 20 0 0 50 100 150 200 250 300 350 400 month 贷款70万，每月还款2500元，还款额小于利息 BUPT 尹霄丽 二．时域经典法 第 页 完全解=齐次解+特解 边界条件 BUPT 尹霄丽 齐次解 第 页 1.齐次解：齐次方程的解 ?n?? ay?n ?1?? 0 y ? ? ? ? ? ?不能全为零 但起始状态y ?1 , y ? 2 ,? y ? N ? ? ? ? ? ? ? ? y 0 y 1 y n y ? ? ? ? ? ? 1 0, ? ? ? ? ? ? ? y 1 y 0 y n 1 ? ? a ? ?是一个公比为 的几何级数 所以 说明 y n a , ? ? y ? n Ca n 或由特征方程 a 0, r r ? ? 可得 ? a 指数形式 ? ? y ? n ? Ca n Cr n BUPT 尹霄丽 求待定系数 C由边界条件决定 第 页 ?0? 2 给定 y ? ? ? Ca C y 0 ? 0 ? 齐次解 C ? 2 ? ? y ? 2 n a n 求差分方程齐次解步骤 差分方程? 特征方程?特征根? y(n)的解?由起始状态确定系数 BUPT 尹霄丽 根据特征根，解的三种情况