 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 429 | 回復(fù): 1 | ||
2008130170金蟲 (小有名氣)
|
[求助]
在求無窮積分時(shí)，上限及步長的選取
|
|
蟲友們，我算無窮積分時(shí)，發(fā)現(xiàn)積分上限和步長會(huì)影響結(jié)果，甚至出現(xiàn)NoN值，公式已在附件里，敬請多多指教！謝謝！如果能幫忙化簡下，更好，運(yùn)行太慢了 clear all lambda=1; tend=599.99;%積分上限 gamma=0.2; theta=0:0.001:pi/2; nn=1000;%步長決定積分的精確度 s= (tend-0)/nn; t=0:s:tend; m=length(t); d=sqrt(lambda.^2-2*lambda*gamma); for j=1:m k(j)=exp(-lambda*t(j)/2).*(cosh(d*t(j)/2)+lambda./d*sinh(d*t(j)/2)); kdot(j)=-0.5*lambda*k(j)+0.5*d*exp(-lambda*t(j)/2).*(sinh(d*t(j)/2)+lambda./d*cosh(d*t(j)/2)); p(j)=k(j)^(2); pdot(j)=2*k(j).*kdot(j); for q=1:length(theta) sigma1(j,q)=(cos(2*theta(q))^2*p(j)+0.5*sin(2*theta(q))^2)*pdot(j)./sqrt(cos(2*theta(q))^2*p(j)^2+sin(2*theta(q))^2*p(j)); sigma(j,q)=0.5*(sigma1(j,q)+abs(sigma1(j,q))); end end for q=1:length(theta) trap(q)=0; for i=1:nn-1 trap(q)=trap(q)+0.5*s*(sigma(i,q)+sigma(i+1,q)); end end [x,y]=max(trap)%積分后對theta求最大值 fig.png |
回復(fù)此樓» 猜你喜歡
一志愿中南化學(xué)337求調(diào)劑
已經(jīng)有4人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有23人回復(fù)
-
265求調(diào)劑
已經(jīng)有9人回復(fù)
-
279分求調(diào)劑 一志愿211
已經(jīng)有10人回復(fù)
-
招收調(diào)劑碩士
已經(jīng)有12人回復(fù)
-
296求調(diào)劑
已經(jīng)有8人回復(fù)
-
招收博士1-2人
已經(jīng)有3人回復(fù)
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有11人回復(fù)
-
274求調(diào)劑
已經(jīng)有8人回復(fù)
-
290求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 碳納米管MD模擬中積分時(shí)間步長選取時(shí)存在的問題
已經(jīng)有15人回復(fù)
2008130170
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1384.3
- 散金: 30
- 帖子: 65
- 在線: 27.9小時(shí)
- 蟲號(hào): 1884947
- 注冊: 2012-07-09
- 專業(yè): 光學(xué)
2樓2015-05-02 21:08:36
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 2008130170 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 307求調(diào)劑 +9 | 冷笙123 2026-03-17 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +3 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 5/250 |
|
|
|
[考研] 266求調(diào)劑 +5 | 陽陽哇塞 2026-03-14 | 10/500 |
|
|
|
[考研] 一志愿福大288有機(jī)化學(xué)，求調(diào)劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
|
[考研] 0703化學(xué) 305求調(diào)劑 +4 | FY_yy 2026-03-14 | 4/200 |
|
|
|
[考研] 354求調(diào)劑 +4 | Tyoumou 2026-03-18 | 7/350 |
|
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯！ 2026-03-15 | 18/900 |
|
|
|
[考研] 08工科 320總分 求調(diào)劑 +5 | 梨花珞晚風(fēng) 2026-03-17 | 5/250 |
|
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
|
[考研] 22408總分284求調(diào)劑 +3 | InAspic 2026-03-13 | 3/150 |
|
|
|
[考研] 266求調(diào)劑 +4 | 學(xué)員97LZgn 2026-03-13 | 4/200 |
|
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|