 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 443 | 回復(fù): 0 | |||
zwb565055403新蟲(chóng) (小有名氣)
|
[交流]
編程為什么結(jié)果出來(lái)不對(duì)啊
|
|
我寫(xiě)的程序如下（雖然能運(yùn)行出來(lái)，結(jié)果卻不對(duì) 大家?guī)臀铱聪鲁绦蚴欠裼袉?wèn)題）： % Example 1 clear;clc %%%%%%%%%%%%%%%%%%%% 方程里的參數(shù)設(shè)置 %%%%%%%%%%%%%%%%%%%%%%%%%%%% alpha=0.5; beta=0.5; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% L=1; %區(qū)間x的長(zhǎng)度[0,1] Mx=L/0.0025; %區(qū)間x的等分的份數(shù) h=L/Mx; %區(qū)間x的劃分小區(qū)間長(zhǎng)度 x=[0:Mx]*h; % x的點(diǎn)值 %上面給出了x的點(diǎn)值，結(jié)合t點(diǎn)的點(diǎn)值和函數(shù) u(x,t)的值就可以畫(huà)出 曲面u(x,t),三維數(shù)組 N=1000; %時(shí)間t的等分份數(shù) tau=0.0025; % 時(shí)間t的劃分小區(qū)間長(zhǎng)度 t=[0:N]*tau; %時(shí)間t的點(diǎn)值 0點(diǎn)和1點(diǎn) 101個(gè) %%%%%%%%%%%%%%%%%%%%%%%%%%% 系數(shù)矩陣?yán)锩娴某?shù)定義 r R(i) %%%%%%%%%%%%%%%%%% R=(tau^alpha)*(gamma(2-alpha))/(h*h); % r常數(shù) %計(jì)算 r(j) for j=1:Mx-1 r(j)=beta*R/j; % 變數(shù) r_j end r; %%%%%%%%%%%%%%%%%%%%%%%%%%%%分?jǐn)?shù)階系數(shù)矩陣%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% v0=1+2*R-r; %主對(duì)角線上元素 A=diag(v0); v1=r-R ; v1=v1(2:end);%下次對(duì)角線上的元素 B=diag(v1,-1); v3=-R*ones(1,Mx-2); %上次對(duì)角線上的元素 C=diag(v3,1); D=A+B+C; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %非齊次項(xiàng)f(x,t)=2\sqrt(t/pi) sin(pi x)-pi t ((cos pi x)/(2x)-pi sin(pi x)) for k=1:size(t,2) %size(a,2)指行向量元素個(gè)數(shù) F(k,  =2*(t(k)^(1/2)/pi^(1/2)).*sin(pi*x) + pi*t(k).*(pi*sin(pi.*x) - cos(pi*x)./(2*x));%第k行向量 時(shí)間不變而空間位置變化 end F=F(2:end-1,2:end-1); F=F'; %去掉第一行和最后一行 去掉第一列和最后一列 對(duì)應(yīng) t=0 x=0 %%%%%%%%%%%%%%%%%%下面用遞推關(guān)系式得到關(guān)于 u 的數(shù)表 %%%%%%%%%%%%%%%%%%%%%%%%%%% u(:,1)=inv(D)*((tau^alpha)*gamma(2-alpha)*F(:,1)); % D*u1=u0+ tau^alpha)F1 例一中 u0=0 即初值函數(shù)為0函數(shù) %定義 w(k)=(1+k)^(1-alpha)-k^(1-alpha);d(k)=w(k)-w(k+1) for k=1:N w(k)=(1+k)^(1-tau)-k^(1-tau); end for k=1:N-1 d(k)=w(k)-w(k+1); %系數(shù) d(k) end d; %%%%%%%%%%%%%%%%%%%%%%%%%% 單獨(dú)算出 u2 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%% u(:,2)=inv(D)*((1-w(1))*u(:,1)+(tau^alpha)*gamma(2-alpha)*F(:,1)); %%%%%求和 \sum d(k)*u^(n-k) 列和 sum(A(:,a:b);2) for n=2:N-2 for k=1:n-1 s(:,k)=d(n-k)*u(:,k); end u(:,n+1)=inv(D)*((tau^alpha)*gamma(2-alpha)*F(:,n+1)+(1-w(1))*u(:,n)+sum(s(:,1:n-1),2)); % u0=0 end h=x(2:end-1); g=t(2:end-1); [h,g]=meshgrid(h,g);%生成網(wǎng)格 mesh(h,g,u');%畫(huà)出曲面圖 xlabel('x') %添加坐標(biāo)軸 ylabel('t') zlabel('u(x,t)') for k=1:N-1 for j=1:Mx-1 plot3(x(j),t(k),u(j,k)) e(j,k)=u(j,k)-t(k)*sin(pi*x(j)); %誤差函數(shù) end end e=max(max(abs(e))) %最大誤差 |
回復(fù)此樓» 猜你喜歡
314求調(diào)劑
已經(jīng)有8人回復(fù)
-
08工科 320總分 求調(diào)劑
已經(jīng)有5人回復(fù)
-
085601專碩，總分342求調(diào)劑，地區(qū)不限
已經(jīng)有5人回復(fù)
-
085600材料與化工調(diào)劑 324分
已經(jīng)有4人回復(fù)
-
311求調(diào)劑
已經(jīng)有6人回復(fù)
-
化工學(xué)碩306求調(diào)劑
已經(jīng)有10人回復(fù)
-
286求調(diào)劑
已經(jīng)有8人回復(fù)
-
材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦
已經(jīng)有9人回復(fù)
-
311求調(diào)劑
已經(jīng)有12人回復(fù)
-
271材料工程求調(diào)劑
已經(jīng)有3人回復(fù)
找到一些相關(guān)的精華帖子，希望有用哦~
在下編程菜鳥(niǎo)，這是我編的一個(gè)c++的程序，運(yùn)行不了總是有一個(gè)錯(cuò)誤，求大神！
已經(jīng)有6人回復(fù)
- 提示矩陣維數(shù)必須一致，不清楚在哪里編程出現(xiàn)了錯(cuò)誤，求指教
已經(jīng)有5人回復(fù)
- 研究結(jié)果表明不好，能否發(fā)出去文章？
已經(jīng)有7人回復(fù)
- 畢業(yè)論文急需要ansys模型，算出來(lái)結(jié)果不對(duì)啊，
已經(jīng)有3人回復(fù)
- 計(jì)算圓周率得到的結(jié)果錯(cuò)誤
已經(jīng)有10人回復(fù)
- 程序錯(cuò)誤求教�。。�！
已經(jīng)有3人回復(fù)
- 請(qǐng)大俠指教matlab的這個(gè)錯(cuò)誤提示到底是個(gè)神馬意思？
已經(jīng)有12人回復(fù)
- Fortran編程過(guò)程中遇到錯(cuò)誤，求大俠幫忙看下
已經(jīng)有5人回復(fù)
- Fortran提示如下錯(cuò)誤：
已經(jīng)有10人回復(fù)
- 怎么做MOFs連粉末都沒(méi)出來(lái)，出來(lái)的晶體顏色都不對(duì)。。。是什么問(wèn)題�。�
已經(jīng)有15人回復(fù)
- fortran 程序有write和沒(méi)有write兩種情況計(jì)算的結(jié)果為什么會(huì)不一樣?
已經(jīng)有14人回復(fù)
- 大家愛(ài)幫我看看這段MATLAB程序?yàn)樯冻鲥e(cuò)了
已經(jīng)有5人回復(fù)
- 兩次調(diào)用同一子程序 第二次調(diào)用得出的結(jié)果是錯(cuò)誤的。。
已經(jīng)有10人回復(fù)
- 急求fortran運(yùn)行錯(cuò)誤原因，在線等
已經(jīng)有7人回復(fù)
- 求助一下程序調(diào)試過(guò)程出現(xiàn)的錯(cuò)誤
已經(jīng)有4人回復(fù)
- 請(qǐng)教 fortran 運(yùn)行錯(cuò)誤的原因
已經(jīng)有13人回復(fù)
- 譚浩強(qiáng)C語(yǔ)言教材，P298頁(yè)照著輸入都不對(duì)！求原因
已經(jīng)有12人回復(fù)
- 【求助】DLPOLY計(jì)算出來(lái)的結(jié)果，是否要編程才能分析？
已經(jīng)有15人回復(fù)
- 【求助】有很多錯(cuò)誤的程序，找不出來(lái)了
已經(jīng)有18人回復(fù)
- 【求助】程序運(yùn)行后出現(xiàn)的錯(cuò)誤 access violation
已經(jīng)有15人回復(fù)
點(diǎn)擊這里搜索更多相關(guān)資源科研從小木蟲(chóng)開(kāi)始，人人為我，我為人人
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 學(xué)員J46msc 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 070300化學(xué)319求調(diào)劑 +6 | 錦鯉0909 2026-03-17 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 | 我愛(ài)生物生物愛(ài)?/a> 2026-03-17 | 3/150 |
|
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
|
[考研] 268求調(diào)劑 +7 | 好運(yùn)連綿不絕 2026-03-12 | 8/400 |
|
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
|
[考研] 211本，11408一志愿中科院277分，曾在中科院自動(dòng)化所實(shí)習(xí) +6 | Losir 2026-03-12 | 7/350 |
|
|
|
[論文投稿] 有沒(méi)有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[基金申請(qǐng)]
今年的國(guó)基金是打分制嗎？
 50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
|
[考研] 085600調(diào)劑 +5 | 漾漾123sun 2026-03-12 | 6/300 |
|
|
|
[考研] 一志愿華中師范071000，325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對(duì)?/a> 2026-03-13 | 5/250 |
|
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
|
[考研]
 機(jī)械專碩調(diào)劑
+3
|
笨笨兔子 2026-03-12 | 3/150 |
|
|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請(qǐng) +6 | zlingli 2026-03-13 | 6/300 |
|
|
|
[考研] 294求調(diào)劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過(guò)A區(qū)線 +6 | 步川酷紫123 2026-03-13 | 6/300 |
|
|
|
[考研] 281求調(diào)劑 +9 | Koxui 2026-03-12 | 11/550 |
|
|
|
[考研] 277求調(diào)劑 +4 | anchor17 2026-03-12 | 4/200 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-12 | 4/200 |
|