 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 8 | 1/1 | 返回列表 |
| 查看: 1110 | 回復: 7 | ||
[求助]
5點差分格式解邊值問題。。。求大神幫忙 已有1人參與
|
|
有沒有類似的程序，哪位大神給我發(fā)一下，或者幫我寫以下程序。 IMG_20151229_193059.jpg |
回復此樓» 猜你喜歡
307求調劑
已經有20人回復
-
285求調劑
已經有8人回復
-
材料與化工（0856）304求B區(qū)調劑
已經有8人回復
-
370求調劑
已經有3人回復
-
070300化學279求調劑
已經有5人回復
-
福建理工大學材料學院先進合金團隊招收考研調劑學生
已經有4人回復
-
英一數一408，總分284，二戰(zhàn)真誠求調劑
已經有5人回復
-
083000環(huán)境科學與工程調劑，總分281
已經有4人回復
-
求調劑
已經有3人回復
-
食品工程專碩一志愿中海洋309求調劑
已經有10人回復
米米米BO
金蟲 (初入文壇)
- 應助: 6 (幼兒園)
- 金幣: 2065.3
- 紅花: 4
- 帖子: 48
- 在線: 12.5小時
- 蟲號: 2814649
- 注冊: 2013-11-20
- 性別: GG
- 專業(yè): 土壤學
2樓2015-12-29 19:55:20
|
下面別人的程序，我該如何利用我的問題里面的邊值條件呢？ function [ u x y k]=wudianchafenfa1(h,m,n,kmax,ep) % g-s迭代法解五點差分法問題 %kmax為最大迭代次數 %m,n為x,y方向的網格數，例如（2-0）/0.01=200;aa syms temp; u=zeros(n+1,m+1); x=-1+(0:m)*h; y=-1+(0:n)*h; for(i=1:n+1) u(i,1)=sin(3*pi/4+pi*y(i)/4); u(i,m+1)=sin(3*pi/4+pi*y(i)/4); end for(i=1:n) for(j=1:m) f(i,j)=32*u(i,j); end end t=zeros(n-1,m-1); for(k=1:kmax) for(i=2:n) for(j=2:m) temp=h*h*f(i,j)/4+(u(i,j+1)+u(i,j-1)+u(i+1,j)+u(i-1,j))/4; t(i,j)=(temp-u(i,j))*(temp-u(i,j)); u(i,j)=temp; end end t(i,j)=sqrt(t(i,j)); if(k>kmax) break; end if(max(max(t))<ep) break; end end |
3樓2015-12-29 21:42:43
米米米BO
金蟲 (初入文壇)
- 應助: 6 (幼兒園)
- 金幣: 2065.3
- 紅花: 4
- 帖子: 48
- 在線: 12.5小時
- 蟲號: 2814649
- 注冊: 2013-11-20
- 性別: GG
- 專業(yè): 土壤學
4樓2015-12-30 00:13:43
5樓2015-12-30 07:36:59
米米米BO
金蟲 (初入文壇)
- 應助: 6 (幼兒園)
- 金幣: 2065.3
- 紅花: 4
- 帖子: 48
- 在線: 12.5小時
- 蟲號: 2814649
- 注冊: 2013-11-20
- 性別: GG
- 專業(yè): 土壤學
【答案】應助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
柱子234: 金幣+20, ★★★★★最佳答案 2015-12-31 08:57:13
柱子234: 金幣+20, ★★★★★最佳答案 2015-12-31 08:57:13
|
clear,clc; n=9; % grid number h=0.25; % step length u=zeros(n,n); u(:,1)=1; % top boundary u(:,n)=0; % bottom boundary for k=1:1000 % maximum iteration u0=u; % left boundary u(2,2:n-1)=(u(3,2:n-1)+u(2,1:n-2)+u(2,3:n)-0.5*h)/5.0; u(1,2:n-1)=u(2,2:n-1)-0.5*h; % right boundary u(n-1,2:n-1)=(u(n-2,2:n-1)+u(n-1,1:n-2)+u(n-1,3:n)-0.5*h)/5.0; u(n,2:n-1)=u(n-1,2:n-1)-0.5*h; % middle grid u(3:n-2,2:n-1)=(u(2:n-3,2:n-1)+u(4:n-1,2:n-1)+u(3:n-2,1:n-2)+u(3:n-2,3:n))/6.0; delta=max(max(abs(u-u0))); if(delta<1e-8) % precision disp(['Number of iteration: ',num2str(k)]); disp(u'); break; end end 運行下這個試試 |
6樓2015-12-30 10:13:51
米米米BO
金蟲 (初入文壇)
- 應助: 6 (幼兒園)
- 金幣: 2065.3
- 紅花: 4
- 帖子: 48
- 在線: 12.5小時
- 蟲號: 2814649
- 注冊: 2013-11-20
- 性別: GG
- 專業(yè): 土壤學
7樓2015-12-30 10:16:32
8樓2015-12-31 08:56:58
相關版塊跳轉
我要訂閱樓主 柱子234 的主題更新
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 英一數一408，總分284，二戰(zhàn)真誠求調劑 +3 | 12.27 2026-03-30 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿北化085600材料專碩275|有文章專利｜求調劑 +10 | Micky11223 2026-03-25 | 10/500 |
|
|
|
[考研] 材料求調劑 一志愿哈工大總分298分，前三科223分 +10 | dongfang59 2026-03-27 | 10/500 |
|
|
|
[考研] 297 地理學070500 復試求調劑 +3 | 小圓圈圈ooo 2026-03-30 | 3/150 |
|
|
|
[考研] 085600 286分 材料求調劑 +11 | 麻辣魷魚 2026-03-27 | 12/600 |
|
|
|
[考研] 085600 材料與化工 329分求調劑 +18 | Mr. Z 2026-03-25 | 19/950 |
|
|
|
[考研] 071010 323 分求調劑 +3 | Baekzhy 2026-03-27 | 3/150 |
|
|
|
[考研] 求調劑，一志愿 南京航空航天大學大學 ，080500材料科學與工程學碩 +6 | @taotao 2026-03-26 | 7/350 |
|
|
|
[考研] 材料與化工304求B區(qū)調劑 +4 | 邱gl 2026-03-26 | 7/350 |
|
|
|
[考研] 318一志愿吉林大學生物與醫(yī)藥 求調劑 +5 | 篤行致遠. 2026-03-28 | 5/250 |
|
|
|
[考研] 275求調劑 +15 | Micky11223 2026-03-25 | 20/1000 |
|
|
|
[考研] 一志愿北京理工大學本科211材料工程294求調劑 +8 | mikasa的圍巾 2026-03-28 | 8/400 |
|
|
|
[考研] 330分求調劑 +5 | qzenlc 2026-03-29 | 5/250 |
|
|
|
[考研] 數一英一271專碩（085401）求調劑，可跨 +7 | 前行必有光 2026-03-28 | 8/400 |
|
|
|
[考研] 299求調劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
|
[考研] 265求調劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
|
[考研] 一志愿鄭大085600，310分求調劑 +5 | 李瀟可 2026-03-26 | 5/250 |
|
|
|
[考研] 292求調劑 +4 | 求求了收下我吧?/a> 2026-03-26 | 4/200 |
|
|
|
[考研] 調劑 +4 | 柚柚yoyo 2026-03-26 | 4/200 |
|
|
|
[考研] 材料專碩331求調劑 +4 | 鮮當牛 2026-03-24 | 4/200 |
|