 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 1048 | 回復: 3 | ||
| 【懸賞金幣】回答本帖問題，作者hylb3020將贈送您 30 個金幣 | ||
hylb3020木蟲 (正式寫手)
流體力學愛好者
|
[求助]
matlab程序， 被調(diào)用函數(shù)中使用 if 語句出錯 已有1人參與
|
|
|
clear all; close all; l=45; alpha1=-0.25; lambda1=0.25; n=128; x=l/n*[-n/2:n/2-1]; y=x; kx=(2*pi/l)*[0:n/2-1 -n/2:-1]; kx(1)=1e-6; ky=kx; [x,y]=meshgrid(x,y); [kx,ky]=meshgrid(kx,ky); h=2*pi/n; column=[00.5*(-1).^(1:n-1).*cot((1:n-1)*h/2)]'; d=(2*pi/l)*toeplitz(column,column([1n:-1:2])); column=[-pi^2/(3*h^2)-1/6-0.5*(-1).^(1:n-1)./sin(h*(1:n-1)/2).^2]; d2=(2*pi/l)^2*toeplitz(column); column=[0 (-1).^(1:n-1).*cot((1:n-1)*h/2).*(-pi^2/(2*h^2)+3/4*csc((1:n-1)*h/2).^2)]'; w=1+0.1*cos(0.31*x)+0.1*cos(0.31*y);%初始條件 t=[0 17 17.5 18 19 20 21 22 23 25 30 40] [t,wsol]=ode45('advection_diffusion', t, w( : ),[],n,d,d2,alpha1,lambda1,x); advection_diffusion.m function dw=advection_diffusion(t,w,dummy,n,d,d2,alpha1,lambda1, x); w=reshape(w,n,n); wt=fft2(w); if (t<=20) dw=reshape(i*(alpha1*d2*w+lambda1*w*d2')+t,n^2,1); else dw=reshape(i*(alpha1*d2*w+lambda1*w*d2')-t,n^2,1); end 我想求解一個偏微分方程，變量t 取了12個值，在被調(diào)用的函數(shù) advection_diffusion 中需要用到一個 if 語句，對 t進行分類，程序一直運行不出來。希望得到大家的指點，不勝感激。  |
回復此樓» 猜你喜歡
西安電子科學大學杭州研究院劉麗香教授招收智能多模態(tài)傳感器和微型儲能器件方向博士
已經(jīng)有13人回復
-
重慶交通大學光子學微結(jié)構(gòu)與器件課題組招收2026年碩士研究生信息
已經(jīng)有1人回復
-
物理學I論文潤色/翻譯怎么收費?
已經(jīng)有190人回復
-
一志愿鄭大材料學碩298分，求調(diào)劑
已經(jīng)有6人回復
-
尋合作：應力腐蝕多尺度模擬
已經(jīng)有3人回復
-
考研交流
已經(jīng)有0人回復
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復
獨孤神宇
版主 (知名作家)
- 應助: 490 (碩士)
- 貴賓: 0.008
- 金幣: 31016.3
- 散金: 802
- 紅花: 122
- 沙發(fā): 1
- 帖子: 5600
- 在線: 856.5小時
- 蟲號: 3522474
- 注冊: 2014-11-06
- 性別: GG
- 專業(yè): 機械動力學
- 管轄: 計算模擬
2樓2021-06-29 19:41:51
hylb3020
木蟲 (正式寫手)
流體力學愛好者
- 應助: 14 (小學生)
- 金幣: 1976.9
- 散金: 498
- 紅花: 5
- 帖子: 317
- 在線: 402.4小時
- 蟲號: 2020234
- 注冊: 2012-09-22
- 性別: GG
- 專業(yè): 流體力學
3樓2021-06-29 22:06:54
送紅花一朵 |
|
4樓2021-06-29 22:08:41
相關版塊跳轉(zhuǎn)
我要訂閱樓主 hylb3020 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學求調(diào)劑，總分320分，一志愿中南民族大學 +3 | 19286712440 2026-03-04 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿武漢理工大學-085602-總分296分-求調(diào)劑 +5 | 紫川葡柚 2026-03-04 | 5/250 |
|
|
[考研] 求調(diào)劑 +7 | 博斯特525 2026-03-04 | 7/350 |
|
|
|
[考研] 085600材料調(diào)劑 總分330 +3 | 池池丶 2026-03-03 | 3/150 |
|
|
|
[考研] 266求調(diào)劑 +6 | 哇塞王帥 2026-03-03 | 6/300 |
|
|
|
[考研] 材料專碩290求調(diào)劑 +3 | 杰尼龜aaa 2026-03-04 | 3/150 |
|
|
|
[論文投稿]
 主編終審卡者不動
 100+3
|
Stray2021 2026-03-03 | 3/150 |
|
|
|
[考研] 一志愿西交大材料學碩 346 求調(diào)劑 +3 | zju51 2026-03-04 | 3/150 |
|
|
|
[考研] 0703化學調(diào)劑 +4 | G212 2026-03-03 | 5/250 |
|
|
|
[考研] 0703化學 學碩 理工科均可 不區(qū)分研究方向 總分279求調(diào)劑 +5 | 1一11 2026-03-03 | 5/250 |
|
|
|
[考研] 理學，工學，農(nóng)學調(diào)劑，少走彎路，這里歡迎您！ +8 | likeihood 2026-03-02 | 11/550 |
|
|
|
[考研] 085602化學工程350，調(diào)劑，有沒有211的 +5 | 利好利好. 2026-03-02 | 9/450 |
|
|
|
[考研] 283求調(diào)劑 +9 | 鹿沫笙 2026-03-02 | 9/450 |
|
|
|
[考研] 一志愿東北大學化學314分求調(diào)劑 +5 | lr1212.. 2026-03-02 | 5/250 |
|
|
|
[考研] 化工京區(qū)271求調(diào)劑 +7 | 11ing 2026-03-02 | 7/350 |
|
|
|
[考研] 298求調(diào)劑 +10 | 人間唯你是清歡 2026-02-28 | 14/700 |
|
|
|
[基金申請]
本子寫完了，給DS兄弟看了，得了92分
+3
|
Doma 2026-03-01 | 7/350 |
|
|
|
[基金申請] 成果系統(tǒng)訪問量大，請一小時后再嘗試。---NSFC啥時候好哦，已經(jīng)兩天這樣了 +4 | NSFC2026我來了 2026-02-28 | 4/200 |
|
|
|
[考研]
尋找調(diào)劑
+4
|
LYidhsjabdj 2026-02-28 | 4/200 |
|
|
|
[論文投稿]
Optics letters投稿被拒求助
30+3
|
luckyry 2026-02-26 | 4/200 |
|