24小时热门版块排行榜    

查看: 3114  |  回复: 4

shali1165

木虫 (小有名气)

[交流] 【求助】请教Matlab高手求定积分近似值问题 已有2人参与

现有一积分 ,精确值不能求出,请高手帮忙求一下近似值,过程也贴出来给看一下,谢谢啊,
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijinfeng042

木虫 (小有名气)

Matlab


robert2020(金币+1):鼓励应助! 2010-05-25 22:23:30
shali1165(金币+1): 2010-05-26 08:40:22
我尝试一下数值积分吧~呵呵 估计梯形法是不太有精确一些的结果了 ....;找一下再说 呵呵 找不到就只有看书帮你再写一下了
工作了,偶尔会上来~可以关注新浪微博 @云是风的梦_Matlab
2楼2010-05-25 21:24:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gitan

木虫 (正式写手)

Visiting Scholar


robert2020(金币+1):鼓励交流经验!多谢应助! 2010-05-25 22:33:06
shali1165(金币+1):8 2010-05-26 08:40:30
shali1165(金币+8): 2010-05-26 08:40:54
该函数在积分区间内的 x=1.4304e-006 存在奇异点,因此在matlab中调用数值求解函数quad无法求出结果。
关于奇异点的求解,只需得到分母部分的零点即可,求解代码如下:
f1 =  @(x)(x.^0.9595 - 10*x.^1.1306);
x0 = fzero(f1,[3.57e-7 3.57e-6])
得到
x0 =

  1.4304e-006
3楼2010-05-25 22:23:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gitan

木虫 (正式写手)

Visiting Scholar


robert2020(金币+1):多谢应助! 2010-05-25 22:33:20
shali1165(金币+1):10 2010-05-26 08:42:02
shali1165(金币+9): 2010-05-26 08:42:58
用Mathematica可以得到楼主给出的积分的符号运算结果为:
163.156 x^0.0405 +  x^0.0405 (-163.156 +   24.6914 Hypergeometric2F1[0.236704, 1., 1.2367, 10. x^0.1711])

在积分区间求定积分的结果为
1.76457 - 10.6462 I
4楼2010-05-25 22:26:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijinfeng042

木虫 (小有名气)

Matlab


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by gitan at 2010-05-25 22:26:38:
用Mathematica可以得到楼主给出的积分的符号运算结果为:
163.156 x^0.0405 +  x^0.0405 (-163.156 +   24.6914 Hypergeometric2F1[0.236704, 1., 1.2367, 10. x^0.1711])

在积分区间求定积分的结果为
1.764 ...

那个结果换成数值是....?交流一下 呵呵 看一下我算得的怎样 我自己也不确定 呵呵
工作了,偶尔会上来~可以关注新浪微博 @云是风的梦_Matlab
5楼2010-05-26 17:34:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 shali1165 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见