liuzhilin 发表于 2011-8-30 09:22:49

matlab积分

函数f=sqrt(x^2-1.655^2)/((2*x.^2-1.655^2)^2-4*x^2*sqrt(x^2-1.655^2)*sqrt(x^2-1)),请问,怎么用matlab实现函数从0到无穷的积分。急用啊~拜托各位大侠了!!

scott198510 发表于 2011-8-30 17:01:10

syms x
int (sqrt(x^2-1.655^2)/((2*x.^2-1.655^2)^2-4*x^2*sqrt(x^2-1.655^2)*sqrt(x^2-1)));
ans =

562949953421312/3915936170993989*log((4503599627370496*x^2-12335471969358474)^(1/2)-67108864*x)+42830321400457905640628044137938718153775251456/3915936170993989*sum(_R*log(((4503599627370496*x^2-12335471969358474)^(1/2)-67108864*x)^2+3260214573553528817903373391145451420397166265407105790348388825821655078913899152710995923430738527067941955289696535948426311255238659514488293905393471663852076268743403641381137070452291191072787727053892342930258903682209507482014868092905413320666734406227636197887583086843251968/13035127478251002887775378327167029453174844673395391012237696797715999414860768196024213746366202305818453997*_R^5-14790103234984823768589093099110182628101109263512010146524031516578444624354283878850683680691246511996770293402489912546346682070970599370234669290158204794610175761680988480440436714068133333750921998287573036930772103176081987677296960/3328738495485301395856623239597110280750810251981756352448564753470433350838669471994339355273*_R^4-115569359217190508579537004643013400307346048956350139732960004100080998279393681841208825089818912789960647490105116911516359008408268953312838728291011494823559145248744770497071465315070029872574922612519697777912450976/13035127478251002887775378327167029453174844673395391012237696797715999414860768196024213746366202305818453997*_R^3-1754585867135477983050029488035747960228970250447288367150304763437097882326175373192960364505270012722352877506735715681982554**25291559307905070867813992029988806289901048/3328738495485301395856623239597110280750810251981756352448564753470433350838669471994339355273*_R^2-2106316747340252127611392819870898749310218364487107919771487372040734909901536276169443706441297107962099335363942939512350627/850049221982175930040193216191363166673366223255693091024224134280604432791157*_R-2820113475822031919147355505623938617883982994448339362514841055789461084894837/868297321369675372923622146401960993407491135675649600191974852),_R = RootOf(4276757067287896547462864692377923325998681441949282408258244123390859189285846129446384848946198488081939387088135146672727412218729378318944572437576430786445651752410837102166404261076983657240108692535434151400420672233071421170231819359447273863090176*_Z^6+4509115671922**2268940072802402250002823832672051295311357691420947992828550981648331212237838142082628425368681758561400460548262262685830011386881819155440259408882873947517355211536303360*_Z^4-1700365964946755872561873273762178534367856542297923483755174012577684388897873178586770289431889057878185430567781273017615824*_Z^2-235148610639873671987531378709830125311756416501313603763958641))-13888506741160505489612158795776*sum(_R*log(((4503599627370496*x^2-12335471969358474)^(1/2)-67108864*x)^2-44386883729072930024713830524730772497919763301061537888059960386009887062503464631437650962666920551322963973167655289332627002581854166364412486660125015650756674139682604535837069336438251896060736392470529171976569704960/3845952103088485469590960270521627893219423549359*_R^5-738841855047926755737327299946655023294883903633087031657524204929821379958186025493096935940973980309596263952470181979026223346152022907665538901025556559880192312486670409229815439552226304/3845952103088485469590960270521627893219423549359*_R^4-55156538234210955489255926895443426863086214307610766848993508207952162471993097043474759401824187473142269709141694962610549118950397474372070954986635901220864/3845952103088485469590960270521627893219423549359*_R^3-517191884188827080769239224704145686090302312232230364925456281305245159519992771979297278985052128798216686076598071025650308144/3845952103088485469590960270521627893219423549359*_R^2-14405974093361738043639163056652178969394015912103651935808783262414024019013666179035031594818856/3845952103088485469590960270521627893219423549359*_R-61299486932437823258825465178854851172043075771278696853068331504/3845952103088485469590960270521627893219423549359),_R = RootOf(738841855047926755737327299946655023294883903633087031657524204929821379958186025493096935940973980309596263952470181979026223346152022907665538901025556559880192312486670409229815439552226304*_Z^6+568831013446011654766120941711213029312571055289489031823481772223982535230533000389942117193917024839291543821164421434290033664*_Z^4+98187299987246420497788356788416232522730664601457493786166480560*_Z^2-1))+171321285601831622562512176551754872615101005824*sum(_R*log(((4503599627370496*x^2-12335471969358474)^(1/2)-67108864*x)^2-4560008670862654963632763316996101709746244854194150688329874974098302336767590218505462816718936927422055419685667710401258073245926984466328458637295452087224857645796381479791122173606350934381902772634555269558925860806919703805323095557930460598927350960588485048899085835177920629488464594824963249349944909842048031826561157465679593825857710847096667134741282425270182641401856/57867226817067256486220417296618467964412442631674050290235424582138832452871682251034150599160162668764048310495977568821445363*_R^5-54361753801505565224611061665394375923849402576017247113969956658598590489086606590974259602000339096943219731173025328093751678865258503248269314547347022110018207698584581145958826482323621985098275694171047252969774652375502711086859901705173619327326878446786701015702564479744232247735256693648249827681109776252928/14777367222096144641319905012838229408781804484792823686960324115385370759438066365425733292955335239448672522967*_R^4-15952890380442631497812022111844086085141135534169589639877417982036446655384135533952363732297347261619956895264930714363191160475820960472484077745446021049551219076586771398983736496900782788526313697845170206902715001508055331414545268864877259018843615496306912653495658354472790501376/57867226817067256486220417296618467964412442631674050290235424582138832452871682251034150599160162668764048310495977568821445363*_R^3-279590542287090045872505792889052882601738491393832184745195238733413711878723377087015271598572452623254313388675132949467618461073599945055740535073869686962619690165864147650012958741467196451828930732995480301616321948928/14777367222096144641319905012838229408781804484792823686960324115385370759438066365425733292955335239448672522967*_R^2-13260032014914389907892391101131983659942297887973879894404201261668563378944417292007494675566437768025887105375902033646123602292084211559748825788621785184254058375852724197362204131640879216/57867226817067256486220417296618467964412442631674050290235424582138832452871682251034150599160162668764048310495977568821445363*_R-254181685415198777964554164616576002646221214099942320389579067996048465507016005305758910539447885779341692969026670032912413888/14777367222096144641319905012838229408781804484792823686960324115385370759438066365425733292955335239448672522967),_R = RootOf(1584373782844600978969856444497563869554721546311006147857810588100293342525303358387238330413637218094860551335848021353436055981194357118695231999384578726419188372363862854423101914347964963560197785251401521334245371891867986776715173820180330337913239391238895727947456014111274889831163303828844486462828591281340416*_Z^6+6030974991865406393662354946696765416819650062230190693519974191442769866811390591266102030602157328016446662879193971865049409550729257809654530223514557699142185030364709656317477913874867454192112001744321271894153368764416*_Z^4+5554152949526860135259549318128409366217482650273033364517926540098948511036090209330124313825951426667957142595761483475482393088*_Z^2-15334556095299063856352274132121))+281474976710656/3915936170993989*sum(_R*log(((x^2-1)^(1/2)-x)^2-12574835683793014910070648907409708879801505505874564075716639074085105656687018060594938169626314252288584705614276933316289491804008087409946452479508103756265042190738331105224184242285622083479830313628008066329841577/79459584170435581027306200981991796938852518345612973788208606964110400113753025666067310801469481495208893622703368092394995864341043161191809073619632313459438735011006307786170848061809225984349768622286989607043072*_R^5-682764944537281722125113507169189950255496870669394997554776481290370231671541652294705947735503262292305972759997776577010408974697602644716915294269464423581281392469498468583091286506134563591879469997/20291337933188582674582771577164380184889089560930011927295879630874957804876389809120926751850293983230061627027618881185042400793474641919279294016922090759935027985563219549063499935816151657714024448*_R^4-7447719632676217167791052985594196812664620716241178679966734396566362980546341993285307114821900838267042151680561690922702589872813475977177288505825025759705274869161648585153742527237251398217495563017798377678745821/39729792085217790513653100490995898469426259172806486894104303482055200056876512833033655400734740747604446811351684046197497932170521580595904536809816156729719367505503153893085424030904612992174884311143494803521536*_R^3-407526801668408706774235294181334865306878380675714230222427857540693897401073992778628773832620490282770265875399798863652376832431372165387086930828387975964452660004392180379101786051085717135081028213/10145668966594291337291385788582190092444544780465005963647939815437478902438194904560463375925146991615030813513809440592521200396737320959639647008461045379967513992781609774531749967908075828857012224*_R^2-55005801793736573459611724270435108720845401720078475414207470687330992573821411318320109823782210529930750815521073741710739636550963321726074127868675099194416945564535573680125748783525/5181733574589392885969904747408597452934776816672378511743851735551960838551367860311699402328850576236182565663523572940087598018823149114850687377624127837644963140538202460527611871232*_R-4030353802010581630438415033162016913793129096123294909858763263399844781998720275066335647673709427194517863871759203980746056445512877242783052017745599065701870174623825/1323242603638788191387090283404679440733489510391365241992554998675839777023885595200471406928251642808308758383090692334660210724473218601273134973842856503488379394981888),_R = RootOf(3075205963250293602311820482617296954645690125492223723936467653282965738496568498151859529872401639178770095336728998306112706815787808744957*_Z^6+3424576817463407947600142252511634482827526979626362153713235011071967084787465441202231459974888330391790861428512620069091742362980451690913*_Z^4-30058068470491955307378863854784701983710997169603915242405184811289261226686249295214745888914233025632118366194488675632181216787312022009*_Z^2-11123048427929556214249701755735699542572547802584909922675582942966284044517902700697566798317640536511556040018760733713775566784883382437))+562949953421312/3915936170993989*log((x^2-1)^(1/2)-x)

根本不能在0-1上面积分

liuzhilin 发表于 2011-9-2 10:20:31

scott198510 发表于 2011-8-30 17:01 static/image/common/back.gif
syms x
int (sqrt(x^2-1.655^2)/((2*x.^2-1.655^2)^2-4*x^2*sqrt(x^2-1.655^2)*sqrt(x^2-1)));
ans =


根本不能在0~1上面积分式什么意思呢?其实我要的就是个结果的虚数部分。

scott198510 发表于 2011-9-2 16:43:40

log((x^2-1)^(1/2)-x)

这个里面你说x能取值是多少?还有积分出来的函数的定义域?你没学过高数?

jyliao 发表于 2011-9-4 13:20:13

(0,1)f=i*Im
(1,1.655) f = Re + i*Im
(x>1.655) f=Re
Integrating in 3-segments
页: [1]
查看完整版本: matlab积分