| 6 | 1/1 | 返回列表 |
| 查看: 3540 | 回復(fù): 5 | ||
goleo木蟲(chóng) (正式寫(xiě)手)
|
[求助]
MATLAB求解偏微分方程組
|

木蟲(chóng) (正式寫(xiě)手)

木蟲(chóng) (正式寫(xiě)手)

銅蟲(chóng) (正式寫(xiě)手)
木蟲(chóng) (正式寫(xiě)手)
|
能否幫我看一下如下解偏微分方程的代碼,這是采用的什么數(shù)值解法,有限差分?但看著公式不像 NN=20; hk=1/NN; RT=[(1-1/sqrt(3))/2 (1+1/sqrt(3))/2]'; H=[(1-RT).^2.*(1+2*RT) RT.*(1-RT).^2 RT.^2.*(3-2*RT) RT.^2.*(RT-1)]; A=[6*RT.*(RT-1) (1-RT).*(1-3*RT) -6*RT.*(RT-1) RT.*(3*RT-2)]; B=[12*RT-6 6*RT-4 6-12*RT 6*RT-2]; E=diag([1,hk,1,hk]); H=H*E; A=A*E; B=B*E; MatrixM=zeros(2*NN,2*NN); MatrixH=zeros(2*NN+2,2*NN+2); MatrixA=zeros(2*NN+2,2*NN+2); MatrixB=zeros(2*NN+2,2*NN+2); for i=1:NN MatrixM(2*i-1:2*i,2*i-1:2*i+2)=H(:, ;MatrixH(2*i:2*i+1,2*i-1:2*i+2)=H(:, ;MatrixA(2*i:2*i+1,2*i-1:2*i+2)=A(:, ;MatrixB(2*i:2*i+1,2*i-1:2*i+2)=B(:, ;end %MatrixH(2*NN+1:2*NN+2,2*NN:2*NN+2)=H(:,[1 2 3]); %MatrixA(2*NN+1:2*NN+2,2*NN:2*NN+2)=A(:,[1 2 3]); %MatrixB(2*NN+1:2*NN+2,2*NN:2*NN+2)=B(:,[1 2 3]); MatrixH(1,1)=1; MatrixH(1,2)=-1/Peclet; %Peclet就是一個(gè)具體數(shù) MatrixH(2*NN+2,2*NN+2)=1; MatrixH; MatrixA; MatrixB; %-------------------------------------------------------------------------- A=MatrixA/hk; |

木蟲(chóng) (正式寫(xiě)手)

| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學(xué) 305求調(diào)劑 +3 | FY_yy 2026-03-14 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 材料專業(yè)求調(diào)劑 +3 | hanamiko 2026-03-18 | 3/150 |
|
|
[考研] 344求調(diào)劑 +5 | knight344 2026-03-16 | 6/300 |
|
|
[考研] 0854可跨調(diào)劑,一作一項(xiàng)核心論文五項(xiàng)專利,省、國(guó)級(jí)證書(shū)40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
[考研] 0703化學(xué)調(diào)劑 ,六級(jí)已過(guò),有科研經(jīng)歷 +10 | 曦熙兮 2026-03-15 | 10/500 |
|
|
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +6 | 薛云鵬 2026-03-13 | 6/300 |
|
|
[考研] 302求調(diào)劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
[考研] 307求調(diào)劑 +3 | 冷笙123 2026-03-17 | 3/150 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
[考研] 有沒(méi)有道鐵/土木的想調(diào)劑南林,給自己招師弟中~ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|
|
[考研] 招收0805(材料)調(diào)劑 +3 | 18595523086 2026-03-13 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | 程雨杭 2026-03-12 | 3/150 |
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級(jí)已過(guò) 本科山東雙非 求調(diào)劑! +3 | 不想理你 2026-03-12 | 3/150 |
|