24小时热门版块排行榜    

查看: 4318  |  回复: 12

yz457694

铁杆木虫 (正式写手)

【答案】应助回帖

星——星: 回帖置顶 2016-08-20 16:14:54
代码没缩进了,重发一下
函数代码:
CODE:
function F=myfun(t)     
  ubound=10000;     
  x=0:0.1:ubound;     
  y=sqrt(x)./(1+exp(x-t));     
  F=trapz(x,y);%梯形积分

求解代码:
CODE:
clear,clc
F=1.4893 ;   %给定的F值
ta=0;tb=100;
fa=F-myfun(ta);
fb=F-myfun(tb); %fa*fb<0即可,不满足的话可以适当放大区间
error=inf;
while error>0.001   %自定义精度     
    tc=(ta+tb)/2;     
    fc=F-myfun(tc);     
    if fa*fc<0         
        tb=tc;     
    else         
        ta=tc;     
    end     
    error=tb-ta;
end
t=(ta+tb)/2;
t

11楼2016-08-20 00:37:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

星——星

新虫 (初入文坛)

引用回帖:
11楼: Originally posted by yz457694 at 2016-08-20 00:37:51
代码没缩进了,重发一下
函数代码:

function F=myfun(t)     
  ubound=10000;     
  x=0:0.1:ubound;     
  y=sqrt(x)./(1+exp(x-t));     
  F=trapz(x,y);%梯形积分

求解代码:

clear,clc
F= ...

嗯嗯,非常感谢~我还有一个问题想请教下您,那要是正常的从0到无穷这个范围内对x进行积分,是不是也要用这种梯形积分的方法?
12楼2016-08-20 12:30:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yz457694

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
星——星: 金币+20, ★★★★★最佳答案, 十分感谢大侠的帮助 抱拳抱拳 2016-08-20 16:16:15
引用回帖:
12楼: Originally posted by 星——星 at 2016-08-20 12:30:18
嗯嗯,非常感谢~我还有一个问题想请教下您,那要是正常的从0到无穷这个范围内对x进行积分,是不是也要用这种梯形积分的方法?...

辛普森积分、梯形积分这些数值积分函数一般都要求确定的上下界才能积。对于简单的被积函数用int的时候能得到无穷积分的解析解,被积函数复杂的话int会报错,像你这个就是,就只能用上面这种做法得数值解。可以参考这篇:
http://wenku.baidu.com/link?url= ... gjWgepWTsFzHq0SFPHG

发自小木虫Android客户端
13楼2016-08-20 14:12:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 星——星 的主题更新
信息提示
请填处理意见