資源簡介
我積極編寫的matlab M文件用來求解菲克第二定律的擴散

代碼片段和文件信息
function?pdemodel
[pde_figax]=pdeinit;
pdetool(‘appl_cb‘1);
set(ax‘DataAspectRatio‘[1?0.029999999999999995?1]);
set(ax‘PlotBoxAspectRatio‘[83.333333333333329?33.333333333333336?166666.66666666666]);
set(ax‘XLimMode‘‘auto‘);
set(ax‘YLimMode‘‘auto‘);
set(ax‘XTickMode‘‘auto‘);
set(ax‘YTickMode‘‘auto‘);
%?Geometry?description:
pderect([0?0.001?1.0000000000000001e-005?0]‘R1‘);
set(findobj(get(pde_fig‘Children‘)‘Tag‘‘PDEEval‘)‘String‘‘R1‘)
%?Boundary?conditions:
pdetool(‘changemode‘0)
pdesetbd(4...
‘dir‘...
1...
‘1‘...
‘0.001‘)
pdesetbd(3...
‘neu‘...
1...
‘0‘...
‘0‘)
pdesetbd(2...
‘dir‘...
1...
‘1‘...
‘0‘)
pdesetbd(1...
‘neu‘...
1...
‘0‘...
‘0‘)
%?Mesh?generation:
setappdata(pde_fig‘Hgrad‘1.3);
setappdata(pde_fig‘refinemethod‘‘regular‘);
pdetool(‘initmesh‘)
%?PDE?coefficients:
pdeseteq(2...
‘0.000000016‘...
‘0.0‘...
‘0‘...
‘1.0‘...
‘0:100‘...
‘0‘...
‘0.0‘...
‘[0?100]‘)
setappdata(pde_fig‘currparam‘...
[‘0.000000016‘;...
‘0.0????????‘;...
‘0??????????‘;...
‘1.0????????‘])
%?Solve?parameters:
setappdata(pde_fig‘solveparam‘...
str2mat(‘0‘‘1000‘‘10‘‘pdeadworst‘...
‘0.5‘‘longest‘‘0‘‘1E-4‘‘‘‘fixed‘‘Inf‘))
%?Plotflags?and?user?data?strings:
setappdata(pde_fig‘plotflags‘[1?1?1?1?1?1?1?1?0?0?0?101?1?0?1?0?0?1]);
setappdata(pde_fig‘colstring‘‘‘);
setappdata(pde_fig‘arrowstring‘‘‘);
setappdata(pde_fig‘deformstring‘‘‘);
setappdata(pde_fig‘heightstring‘‘‘);
%?Solve?PDE:
pdetool(‘solve‘)
?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件???????1471??2011-05-23?16:51??for?fick?law.m
-----------?---------??----------?-----??----
?????????????????1471????????????????????1
評論
共有 條評論