 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 1045 | 回復(fù): 3 | ||
| 【懸賞金幣】回答本帖問(wèn)題，作者h(yuǎn)ylb3020將贈(zèng)送您 30 個(gè)金幣 | ||
hylb3020木蟲(chóng) (正式寫(xiě)手)
流體力學(xué)愛(ài)好者
|
[求助]
matlab程序， 被調(diào)用函數(shù)中使用 if 語(yǔ)句出錯(cuò) 已有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 我想求解一個(gè)偏微分方程，變量t 取了12個(gè)值，在被調(diào)用的函數(shù) advection_diffusion 中需要用到一個(gè) if 語(yǔ)句，對(duì) t進(jìn)行分類(lèi)，程序一直運(yùn)行不出來(lái)。希望得到大家的指點(diǎn)，不勝感激。  |
回復(fù)此樓» 猜你喜歡
推薦一款可以AI輔助寫(xiě)作的Latex編輯器SmartLatexEditor，超級(jí)好用，推薦試試
已經(jīng)有11人回復(fù)
-
2026-CJ開(kāi)始申報(bào)了
已經(jīng)有1人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有101人回復(fù)
-
西安電子科學(xué)大學(xué)杭州研究院劉麗香教授招收智能多模態(tài)傳感器和微型儲(chǔ)能器件方向博士
已經(jīng)有13人回復(fù)
-
重慶交通大學(xué)光子學(xué)微結(jié)構(gòu)與器件課題組招收2026年碩士研究生信息
已經(jīng)有1人回復(fù)
-
一志愿鄭大材料學(xué)碩298分，求調(diào)劑
已經(jīng)有6人回復(fù)
-
尋合作：應(yīng)力腐蝕多尺度模擬
已經(jīng)有3人回復(fù)
-
考研交流
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
獨(dú)孤神宇
版主 (知名作家)
- 應(yīng)助: 490 (碩士)
- 貴賓: 0.008
- 金幣: 31016.3
- 散金: 802
- 紅花: 122
- 沙發(fā): 1
- 帖子: 5600
- 在線: 856.5小時(shí)
- 蟲(chóng)號(hào): 3522474
- 注冊(cè): 2014-11-06
- 性別: GG
- 專(zhuān)業(yè): 機(jī)械動(dòng)力學(xué)
- 管轄: 計(jì)算模擬
2樓2021-06-29 19:41:51
hylb3020
木蟲(chóng) (正式寫(xiě)手)
流體力學(xué)愛(ài)好者
- 應(yīng)助: 14 (小學(xué)生)
- 金幣: 1976.9
- 散金: 498
- 紅花: 5
- 帖子: 317
- 在線: 402.4小時(shí)
- 蟲(chóng)號(hào): 2020234
- 注冊(cè): 2012-09-22
- 性別: GG
- 專(zhuān)業(yè): 流體力學(xué)
3樓2021-06-29 22:06:54
送紅花一朵 |
|
4樓2021-06-29 22:08:41
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 hylb3020 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0857求調(diào)劑 +4 | CDDD1013 2026-03-04 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 材料328求調(diào)劑 +7 | 一個(gè)蘿卜02 2026-03-03 | 7/350 |
|
|
|
[考研] 312求調(diào)劑 +8 | 醋精華了一下發(fā)?/a> 2026-03-03 | 10/500 |
|
|
|
[考研] 347分電子信息專(zhuān)碩求調(diào)劑 +3 | 咕嘟咕嘟1 2026-03-03 | 3/150 |
|
|
|
[考研] 291求調(diào)劑 +3 | 咕嚕咕嚕123123 2026-03-02 | 4/200 |
|
|
|
[考研] 環(huán)境調(diào)劑 +8 | chenhanheng 2026-03-02 | 8/400 |
|
|
|
[考研] 化工專(zhuān)碩調(diào)劑 +4 | 利好利好. 2026-03-03 | 7/350 |
|
|
|
[考研] 298求調(diào)劑 +3 | 人間唯你是清歡 2026-03-03 | 4/200 |
|
|
|
[考研] 成績(jī)276，專(zhuān)業(yè)代碼0856求調(diào)劑 +7 | 小陳朵 2026-03-03 | 7/350 |
|
|
|
[考研] 化工270求調(diào)劑 +10 | 什么名字qwq 2026-03-02 | 10/500 |
|
|
|
[考研] 291求調(diào)劑 +3 | MuoLuo1312 2026-03-02 | 6/300 |
|
|
|
[考研] 0857調(diào)劑 +6 | 一ll半 2026-02-28 | 8/400 |
|
|
|
[考研] 0854復(fù)試調(diào)劑 276 +5 | wmm9 2026-03-01 | 7/350 |
|
|
|
[考研] 085600材料工程一志愿中科大總分312求調(diào)劑 +9 | 吃宵夜1 2026-02-28 | 11/550 |
|
|
|
[考研] 一志愿中石油（華東）本科齊魯工業(yè)大學(xué) +3 | 石能偉 2026-03-02 | 3/150 |
|
|
|
[基金申請(qǐng)] 此成果不能導(dǎo)入原因：元數(shù)據(jù)必填信息不完整，可 進(jìn)行補(bǔ)充。 +4 | Kittylucky 2026-03-02 | 5/250 |
|
|
|
[考研] 一志愿鄭大材料學(xué)碩298分，求調(diào)劑 +6 | wsl111 2026-03-01 | 6/300 |
|
|
|
[考研] 材料調(diào)劑 +6 | 愛(ài)擦汗的可樂(lè)冰 2026-02-28 | 7/350 |
|
|
|
[考研] 274求調(diào)劑 +3 | cgyzqwn 2026-03-01 | 7/350 |
|
|
|
[考研] 0856材料求調(diào)劑 +4 | 麻辣魷魚(yú) 2026-02-28 | 4/200 |
|