| 24小時(shí)熱門版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 1056 | 回復(fù): 3 | ||
| 【懸賞金幣】回答本帖問(wèn)題,作者h(yuǎn)ylb3020將贈(zèng)送您 30 個(gè)金幣 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
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)行分類,程序一直運(yùn)行不出來(lái)。希望得到大家的指點(diǎn),不勝感激。 |

木蟲(chóng) (正式寫(xiě)手)
流體力學(xué)愛(ài)好者

版主 (知名作家)

送紅花一朵
![]() ![]() |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料考研339求調(diào)劑 +3 | Karry*^_^* 2026-03-04 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 332材料求調(diào)劑 +6 | zjy101327 2026-03-05 | 7/350 |
|
|
[考研] 306求調(diào)劑 +4 | Bahati 2026-03-05 | 4/200 |
|
|
[考研] 303求調(diào)劑 +3 | forgman95 2026-03-05 | 5/250 |
|
|
[考研] 282求調(diào)劑 +7 | 夕~日 2026-03-05 | 8/400 |
|
|
[考研] 271求調(diào)劑 +7 | 月色c 2026-03-05 | 8/400 |
|
|
[考研] 歡迎采礦、地質(zhì)、巖土、計(jì)算機(jī)、人工智能等專業(yè)的同學(xué)報(bào)考 +10 | pin8023 2026-02-28 | 13/650 |
|
|
[考研] 0703 東華大學(xué) 理學(xué)化學(xué)均過(guò) a 線,總分 281 求調(diào)劑 +4 | 香香 Lu 2026-03-03 | 8/400 |
|
|
[考研]
材料學(xué)學(xué)碩308分/本科東北大學(xué)/一志愿西北工業(yè)大學(xué)/
5+3
|
蘇堯幺幺 2026-03-03 | 8/400 |
|
|
[考研] 中國(guó)科學(xué)技術(shù)大學(xué)材料與化工281求調(diào)劑,有科研和獲獎(jiǎng)經(jīng)歷 +11 | wsxw 2026-03-02 | 12/600 |
|
|
[考研] 0856材料專碩274能調(diào)劑去哪里? +3 | 22735 2026-03-04 | 4/200 |
|
|
[考研] 304求調(diào)劑 +7 | 曼殊2266 2026-02-28 | 8/400 |
|
|
[考研] 調(diào)劑材料學(xué)碩 +6 | 詞凝Y 2026-03-02 | 6/300 |
|
|
[考研]
|
旅行中的紫葡萄 2026-03-03 | 4/200 |
|
|
[考研] 一志愿中科大能動(dòng)297求調(diào)劑,本科川大 +3 | 邵11 2026-03-03 | 3/150 |
|
|
[考研] 268求調(diào)劑 +6 | 好運(yùn)連綿不絕 2026-03-02 | 6/300 |
|
|
[考研] 307求調(diào)劑 +6 | wyyyqx 2026-03-01 | 6/300 |
|
|
[考研] 一志愿華南理工大學(xué)材料與化工326分,求調(diào)劑 +3 | wujinrui1 2026-02-28 | 3/150 |
|
|
[考研] 調(diào)劑 +3 | 13853210211 2026-03-02 | 4/200 |
|
|
[考研] 311求調(diào)劑 +9 | 南迦720 2026-02-28 | 10/500 |
|