3Â¥: Originally posted by ÔÂÖ»À¶ at 2016-03-31 11:03:15
ÒÔÏ´úÂëÒªÍêÕû¸´ÖƵ½Ò»¸ömÎļþÖУ¬ÔÙÔËÐС£
function qiu_jifen
clear all;clc
format long
jifen=quadgk(@objfun,3.4,10.5)

function f=objfun(t)
T=150.92*t.^(-0.184);
f=exp(-755.8335./(T+70)).*exp ...

jifen=quadgk(@(t)exp(-755.8335./((150.92*t.^(-0.184))+70)).*exp(-70000*(143+(150.92*t.^(-0.184)))./2./((150.92*t.^(-0.184)).^2)./(143-(150.92*t.^(-0.184)))),3.4,10.5)

function qiu_jifen
clear all;clc
format long
jifen=quadgk(@objfun,3.4,10.5)

function f=objfun(t)
T=150.92*t.^(-0.184);
f=exp(-755.8335./(T+70)).*exp(-70000*(143+T)./2./(T.^2)./(143-T));