BBYR Achieve
返回信息流
这是一条镜像帖。来源:北邮人论坛 / matlab / #8336同步于 2011/6/5
该镜像源已超过 30 天没有更新,可能在源站已被删除。
Matlab机器人发帖

求助:关于求含积分的方程 高斯—拉盖尔方法

blindfrog
2011/6/5镜像同步5 回复
本人小菜,做关于通信仿真时遇到这个方程解不出,希望大牛们帮帮指点: 其中a,b都是常数, α = 0.1,t是要求解的变量。 设 a = b = 1 我的代码如下: syms h t y = exp( -h - t / (a + b * h) ); f = int(y, h, 0, 'inf'); solve(f = 0.1, t) 程序执行到积分 int()出报错,错误如下: Warning: Explicit integral could not be found. 怎么办啊!懂的童鞋帮帮教教啊~谢谢啊~ PS: 网上看到说符号积分算不出可以用数值积分解,可是怎么做啊? %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%% 分割线 %%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 网上一个大牛说用 高斯—拉盖尔 积分解这个方程。 找到方法如下: 形如 的积分方程可以转换成以下形式: 其中xk 为 n+1 阶拉盖尔多项式 Ln = 0 的解 Ak 为系数 (怎么求?我不知道啊~~~ 看到了好几种,不知道那种好啊~~~懂得童鞋帮帮忙啊) 拉盖尔多项式为: 0 阶: L0 = 1 1 阶: L1 = -x+1 …… n 阶: 有对这种方法熟悉的童鞋,帮忙讲解下哈~呷呷侬哩!
订阅后,新回复会通过你的通知中心匿名送达。
5 条回复
blindfrog机器人#1 · 2011/6/6
代码是不是只计算了0到50之间的积分,但是我是0到+无穷的积分 为什么50到+无穷之间的积分可以忽略啊? 【 在 dajinsique 的大作中提到: 】 : lz看这个可以么,之间可能比较长,迭代初始值需要设定,最后的t就是求得的值 : 看懂原理后就把里面的fprintf()都去掉,速度会快一些。。 : 【 在 blindfrog 的大作中提到: 】 : ...................
dajinsique机器人#2 · 2011/6/6
哦,确实,是我弄错了。而且那个积分还有点毛病
dajinsique机器人#3 · 2011/6/6
只好做成这样了。。没办法积分到无穷了。。误差大。。 等牛人吧 【 在 blindfrog 的大作中提到: 】 : 代码是不是只计算了0到50之间的积分,但是我是0到+无穷的积分 : 为什么50到+无穷之间的积分可以忽略啊? 附件(3.6KB) diedai.m
blindfrog机器人#4 · 2011/6/6
谢谢啦~ 【 在 dajinsique 的大作中提到: 】 : 只好做成这样了。。没办法积分到无穷了。。误差大。。 : 等牛人吧 : 【 在 blindfrog 的大作中提到: 】 : ...................
blindfrog机器人#5 · 2011/6/7
更新了下 肿么大牛们都这么忙啊.....