1樓:匿名使用者
參考函式pdepe的用法
研究生課題做的就是球座標下的擴散方程
用matlab寫了一個解擴散方程的程式(用的是老師給的程式),可是執行顯示出錯,哪位高手幫幫忙!!
2樓:匿名使用者
%只修改了一點語
法錯誤,修正了向量維數,沒對資料做任何修改:
function demo_mol
clc;clear all;format long;
n=15;dpi=pi/n;
x=dpi:dpi:pi;
size(x)
u=sin(x);
t=0.:0.4:40.;
[t u]=ode45(@myfun,t,u);
uu(:,2:n)=u(:,1:n-1);
uu(:,1)=0;uu(:,n+1)=0;
x=[0. x];[xx yy]=meshgrid(x,t);
surf(xx,yy,uu);
xlabel('x');ylabel('t');zlabel('u(x,t)');
endfunction y=myfun(~,u)n=15;
x=zeros(1,n);
x(1)=-2*u(1)+u(2);
x(n-1)=u(n-2)-2*u(n-1);
for i=2:n-2
x(i)=u(i-1)-2*u(i)+u(i+1);
endy=x';end
怎麼用matlab實現反應擴散方程的畫圖問題?
3樓:匿名使用者
你要畫的是解析解?你可以通過某種演算法解得方程的解析解形式。然後畫圖。例如傅立葉方法等。
也可以畫出它的數值解。時間項和空間項用差分法(或者有限元和譜方法)離散。然後利用初值,迭代出解的形式。
如何用matlab求解時滯反應擴散方程組,附帶畫圖時程式的
4樓:匿名使用者
1將**儲存為m檔案2將下面的**複製到命令視窗,按回車鍵,執行即可fd1d_predator_prey
MATLAB實現基於熱擴散方程的影象去噪
提出了一種利用規則化各向異性熱擴散方程sar影象分割的橋樑檢測演算法。該演算法在perona和malik提出的各向異性熱擴散方程的基礎上構造了一個新的擴散函式,利用數值逼近理論,得到一個新的規則化擴散模型,用此模型對影象初始分割的最大後驗概率矩陣進行多尺度各向異性平滑,得到影象中河流的精確分割結果,...
如何用matlab求方程的整數解
這裡提供一個利用網格矩陣解二維不定方程整數解的方法 已知6整數解。matlab程式如下 x 5 29 y 14 35 x,y meshgrid x,y z 2 x 5 y k find z 126 x k y k 輸出為 ans 8 13 18 23 28 ans 22 20 18 16 14 即方...
Matlab中BP神經網路訓練結果求指導,萌新求各位大神給分
p 1 2 3 4 5 月 p p 50 t 2 3 4 5 6 月訓練樣本 t t 50 threshold 0 1 0 1 0 1 0 1 0 1 0 1 0 1 newff threshold,15,7 trainlm trainparam.epochs 2000 trainparam.goa...