找回密码
 欢迎注册
查看: 20918|回复: 5

[求助] 差分方程的z变换解

[复制链接]
发表于 2010-5-25 15:58:21 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?欢迎注册

×
本帖最后由 wayne 于 2010-5-26 10:43 编辑 $p_{n+1}=14p_n+9q_n$ $q_{n+1}=3p_n+2q_n$ 求上述方程的解,其中$p_0=1,q_0=0$。希望写出详细过程,我推导了一下,总感觉不正确,原因是delta为一平方数,可能是符号有问题。谢谢了!
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2010-5-28 16:40:20 | 显示全部楼层
wayne能否给出上述方程的解,我解出 q_n=(11^n-5^n)/2,但是不对啊!其中得到 Q(z)=3z/(z^2-16z+55),不知错在哪里?
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2010-5-28 18:06:19 | 显示全部楼层
$q(n)= k/(4k-6)*(((14 * k - 9) / k) ^ n - ((2 * k + 9) / k) ^ n)$ $p(n)=(14+3k)^n-k*q(n)$ 其中$k=sqrt(7)-2$
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2010-5-28 18:17:10 | 显示全部楼层
将k代入,化简得: $q(n)=( (8+sqrt(63))^n- (8-sqrt(63) )^n)/(2*sqrt(7))$
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2010-5-28 18:41:54 | 显示全部楼层
2# wsc810 设Z{p(n)}= p(z),Z{q(n)}=q(z) 那么,原方程经过z变换后为 z(p(z)-p(0))=14p(z)+9q(z) z(q(z)-q(0))=3p(z)+2q(z) 代入初始值,解此方程得 p(z)=$\frac{(-2+z) z}{1-16 z+z^2}$ q(z)=$\frac{3 z}{1-16 z+z^2}$ 再进行z反变换即可 ========================================== 更简单的方法是,把方程组转化成二阶差分方程要: $q_{n+2}-2q_{n+1}=14(q_{n+1}-2q_n)+27q_n$ 即:$q_{n+2}-16q_{n+1}+q_n=0$ 接下来,就不说了
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2010-5-28 19:16:55 | 显示全部楼层
本帖最后由 056254628 于 2010-5-28 19:57 编辑 对于以下类型的双递推公式求解,一般方法是: $p(n+1)=a_1*p(n)+b_1*q(n)$ 式1 $ q(n+1)=a_2*p(n)+b_2*q(n)$ 式2 ----------------------------------------------------- 式1+式2*k得: $p(n+1)+k*q(n+1)=(a_1+k*a_2)*p(n)+(b_1+k*b_2)*q(n)$ 若 $1/k=(a_1+k*a_2)/(b_1+k*b_2)$ 即 $a_2*k^2+(a_1-b_2)*k-b_1=0$ 那么 $p(n+1)+k*q(n+1)$就是以 $a_1+k*a_2$ 为等比的数列 那么 $p(n)+k*q(n) =(p(0)+k*q(0))*(a_1+k*a_2)^n$ 即 $p(n)=(p(0)+k*q(0))*(a_1+k*a_2)^n-k*q(n)$ 式3 将式3代入式2,消去数列p(n),得到以下形的数列递推公式: $a(n+1)=e*a(n)+f*g^n$ e、f、g都是常数 而该递推公式是非常容易解的。 $a(n)=e^n*a(0)+(f/(e-g))*(e^n-g^n)$
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

小黑屋|手机版|数学研发网 ( 苏ICP备07505100号 )

GMT+8, 2024-11-22 18:45 , Processed in 0.024446 second(s), 16 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表