返回信息流本人小菜,做关于通信仿真时遇到这个方程解不出,希望大牛们帮帮指点:
其中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 阶:
有对这种方法熟悉的童鞋,帮忙讲解下哈~呷呷侬哩!
这是一条镜像帖。来源:北邮人论坛 / matlab / #8336同步于 2011/6/5
该镜像源已超过 30 天没有更新,可能在源站已被删除。
Matlab机器人发帖
求助:关于求含积分的方程 高斯—拉盖尔方法
blindfrog
2011/6/5镜像同步5 回复
订阅后,新回复会通过你的通知中心匿名送达。
5 条回复
代码是不是只计算了0到50之间的积分,但是我是0到+无穷的积分
为什么50到+无穷之间的积分可以忽略啊?
【 在 dajinsique 的大作中提到: 】
: lz看这个可以么,之间可能比较长,迭代初始值需要设定,最后的t就是求得的值
: 看懂原理后就把里面的fprintf()都去掉,速度会快一些。。
: 【 在 blindfrog 的大作中提到: 】
: ...................
只好做成这样了。。没办法积分到无穷了。。误差大。。
等牛人吧
【 在 blindfrog 的大作中提到: 】
: 代码是不是只计算了0到50之间的积分,但是我是0到+无穷的积分
: 为什么50到+无穷之间的积分可以忽略啊?
附件(3.6KB) diedai.m
谢谢啦~
【 在 dajinsique 的大作中提到: 】
: 只好做成这样了。。没办法积分到无穷了。。误差大。。
: 等牛人吧
: 【 在 blindfrog 的大作中提到: 】
: ...................