 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 956 | 回復: 1 | ||
crackpotop銀蟲 (小有名氣)
|
[求助]
matlab 求解積分問題
|
|
我想求一個含有復數(shù)的積分， 并作出積分結果的幅值與變量a的圖像，matlab程序如下： b=10*10^-6; ef=0.606; es=1.1; af=0.012; as=0.065; w=10; qf=; qs=; k=0; for a=-20*10^-6:1*10^-6:20*10^-6 k=k+1; e(k)= a; f=@(n)1./(ef*sqrt(n.^2+2*1i*w/af)+es*sqrt(n.^2+2*1i*w/as)).*sin(n*b).*cos(n*a)./(n*b); y= quad(f,0,inf); d(k)= abs(y); end plot(e,d) 然而最后的積分結果都是NaN 請問怎么辦？ |
回復此樓» 猜你喜歡
805600專碩材料與化工348分求調(diào)劑
已經(jīng)有6人回復
-
重慶大學材料與化工085600，初試370+，求求調(diào)劑建議
已經(jīng)有5人回復
-
309求調(diào)劑
已經(jīng)有3人回復
-
一志愿北交大材料工程，總分358
已經(jīng)有4人回復
-
302求調(diào)劑一志愿北航070300，本科鄭大化學
已經(jīng)有11人回復
-
化學工程專碩324分，一志愿中國礦業(yè)大學求調(diào)劑
已經(jīng)有7人回復
-
材料求調(diào)劑
已經(jīng)有11人回復
-
【求調(diào)劑】新能源材料本科，一志愿211，初試321
已經(jīng)有4人回復
-
336求調(diào)劑
已經(jīng)有5人回復
-
材料求調(diào)劑
已經(jīng)有8人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 求助MATLAB中程序看不懂的簡單問題！謝謝！
已經(jīng)有3人回復
- 關于matlab求導和積分的兩個問題
已經(jīng)有5人回復
- ＭＡＴＬＡＢ計算編程請教
已經(jīng)有15人回復
- 用含有定積分的隱函數(shù)擬合數(shù)據(jù)，matlab或者origin
已經(jīng)有18人回復
- 求助matlab一段代碼問題
已經(jīng)有3人回復
- matlab可以求解差分方程系統(tǒng)的解嗎？
已經(jīng)有6人回復
- 高斯方程積分編程matlab或C語言
已經(jīng)有15人回復
- matlab求函數(shù)積分出錯
已經(jīng)有12人回復
- matlab不定積分的問題
已經(jīng)有10人回復
- 請教MATLAB中的LMI求解問題，急�。�！
已經(jīng)有7人回復
- Matlab數(shù)值積分結果是虛數(shù)
已經(jīng)有6人回復
- 求高手解答Matlab求解微分方程的初始值和邊界值問題
已經(jīng)有5人回復
- matlab 積分求助
已經(jīng)有8人回復
- matlab 積分計算 算不出結果，求幫助�。�！
已經(jīng)有10人回復
- 關于Matlab算符號積分
已經(jīng)有6人回復
- matlab求積分，出現(xiàn)奇點怎么辦？
已經(jīng)有8人回復
- 【求助】積分微分方程matlab求解
已經(jīng)有6人回復
crackpotop
銀蟲 (小有名氣)
- 應助: 0 (幼兒園)
- 金幣: 409.5
- 散金: 50
- 帖子: 51
- 在線: 27.8小時
- 蟲號: 3294178
- 注冊: 2014-06-26
- 專業(yè): 微/納機械系統(tǒng)
|
看推薦的帖子，積分用 quadgk，出結果了~~ 【http://blog.sina.com.cn/s/blog_8138d87601013wrk.html】 四、[q,errbnd] = quadgk(fun,a,b,param1,val1,param2,val2,...) 自適應Gauss-Kronrod數(shù)值積分，適用于高精度和震蕩數(shù)值積分，支持無窮區(qū)間，并且能夠處理端點包含奇點的情況，同時還支持沿著不連續(xù)函數(shù)積分，復數(shù)域線性路徑的圍道積分法。 注意事項： 1.積分限[a,b]可以是[-inf,inf]，但必須快速衰減； 2.被積函數(shù)在端點可以有奇點，如果區(qū)間內(nèi)部有奇點，將以奇點區(qū)間劃分成多個，也就是說奇點只能出現(xiàn)在端點上； 3.被積函數(shù)可以劇烈震蕩； 4.可以計算不連續(xù)積分，此時需要用到'Waypoints'參數(shù)，'Waypoints'中的點必須嚴格單調(diào)； 5.可以計算圍道積分，此時需要用到'Waypoints'參數(shù)，并且為復數(shù)，各點之間使用直線連接； 6.param,val為函數(shù)的其它控制參數(shù)，比如上面的'waypoints'就是，具體看幫助。 出現(xiàn)錯誤： 1.'Reached the limit on the maximum number of intervals in use' 2.'Infinite or Not-a-Number function value encountered' 例4 計算有奇點積分int(exp(x)*log(x),0,1)。 >>F=@(x)exp(x).*log(x);%奇點必須在端點上，否則請先進行區(qū)間劃分。 >>Q = quadgk(F,0,1) Q = -1.3179 例5 計算半無限震蕩積分int(x^5*exp(-x)*sin(x),0,inf)。 >>F=@(x)x.^5.*exp(-x).*sin(x); >>fplot(F,[0,100])%繪圖，看看函數(shù)的圖形 >>[q,errbnd] = quadgk(F,0,inf,'RelTol',1e-8,'AbsTol',1e-12)%積分限中可以有inf，但必須快速收斂 q = -15.0000 errbnd = 9.4386e-009 例6 計算不連續(xù)積分，積分函數(shù)為f(x)=x^5*exp(-x)*sin(x)，但是人為定義f(2)=1000，f(5)=-100，積分區(qū)間為[1 10]。 >>F=@(x)x.^5.*exp(-x).*sin(x); >>[q,errbnd] = quadgk(F,1,10,'Waypoints',[2 5])%顯然2，5為間斷點 q = -10.9408 errbnd = 3.2296e-014 例7 計算圍道積分，在復數(shù)域內(nèi)，積分函數(shù)1/(2*z-1)，積分路徑為由[-1-i 1-i 1+i -1+i -1-i]圍成的矩形邊框。 >>Waypoints=[-1-i 1-i 1+i -1+i -1-i]; >>plot(Waypoints);%繪制積分路徑 >>xlabel('Real axis');ylabel('Image axis');axis([-1.5 1.5 -1.5 1.5]);grid on; >>Q = quadgk(@(z)1./(2*z - 1),-1-i,-1-i,'Waypoints',[1-i,1+i,-1+i])%注意各點間使用直線連接 ans = 0.0000 + 3.1416i >> quadgk(@(z)1./(2*z - 1),-1-i,-1-i,'Waypoints',Waypoints)%使用這個的效果也是一樣的，就是說始末點可以隨便包不包含在Waypoints中。 ans = 0.0000 + 3.1416i |
2樓2015-06-12 18:50:30
相關版塊跳轉(zhuǎn)
我要訂閱樓主 crackpotop 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 生物學327，求調(diào)劑 +5 | 書上的梅子 2026-04-01 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 289求調(diào)劑 +21 | 新時代材料 2026-03-27 | 23/1150 |
|
|
|
[考研] 295材料工程專碩求調(diào)劑 +19 | 1428151015 2026-03-27 | 19/950 |
|
|
|
[考研] 0856，269分求調(diào)劑 +8 | 有學上就行求求?/a> 2026-03-30 | 11/550 |
|
|
|
[考研] 349求調(diào)劑 +6 | 吃的不少 2026-04-01 | 6/300 |
|
|
|
[考研] 283求調(diào)劑 +9 | A child 2026-03-28 | 9/450 |
|
|
|
[考研] 材料與化工調(diào)劑一志愿大連海事085600，349 +9 | 吃的不少 2026-03-30 | 9/450 |
|
|
|
[考博] 26申博 +4 | 加油沖��！ 2026-03-26 | 4/200 |
|
|
|
[考研] 一志愿：西北大學，英一數(shù)一408-284分求調(diào)劑 +7 | 12.27 2026-03-27 | 7/350 |
|
|
|
[考研] 329求調(diào)劑，一志愿西北工業(yè)大學，材料工程（085601） +6 | 小小機靈蟲 2026-03-29 | 12/600 |
|
|
|
[考研] 263求調(diào)劑 +3 | DDDDuu 2026-03-27 | 3/150 |
|
|
|
[基金申請] 面上5B能上會嗎？ +8 | redcom 2026-03-29 | 8/400 |
|
|
|
[考研] 282求調(diào)劑 不挑專業(yè) 求收留 +4 | Yam. 2026-03-30 | 5/250 |
|
|
|
[考研] 本科211安全工程，初試290分，求調(diào)劑 +3 | 2719846834 2026-03-28 | 3/150 |
|
|
|
[考博] 材料專業(yè)申博 +5 | 杜雨婷dyt 2026-03-29 | 5/250 |
|
|
|
[考研] 08工科，295，接受跨專業(yè)調(diào)劑 +6 | lmnlzy 2026-03-30 | 6/300 |
|
|
|
[考研] 08工科求調(diào)劑286 +5 | tgs_001 2026-03-28 | 5/250 |
|
|
|
[考研] 調(diào)劑 +4 | GK72 2026-03-30 | 4/200 |
|
|
|
[考研] 086000生物與醫(yī)藥調(diào)劑 +5 | Feisty。 2026-03-28 | 9/450 |
|
|
|
[考研] 356求調(diào)劑 +3 | gysy?s?a 2026-03-28 | 3/150 |
|