Лабораторные работы / msu_lab4

function pdemodel
[pde_fig,ax]=pdeinit;
pdetool('appl_cb',9);
set(ax,'DataAspectRatio',[1 0.88421052631578934 1]);
set(ax,'PlotBoxAspectRatio',[4.75 3.166666666666667 1]);
set(ax,'XLim',[0 9.5]);
set(ax,'YLim',[-1.6000000000000001 4]);
set(ax,'XTick',[ 0,...
 0.40000000000000002,...
 0.80000000000000004,...
 1.2000000000000002,...
 1.6000000000000001,...
 2,...
 2.4000000000000004,...
 2.8000000000000003,...
 3.2000000000000002,...
 3.6000000000000001,...
 4,...
 4.4000000000000004,...
 4.8000000000000007,...
 5.1999999999999993,...
 5.5999999999999996,...
 6,...
 6.4000000000000004,...
 6.7999999999999998,...
 7.1999999999999993,...
 7.5999999999999996,...
 8,...
 8.4000000000000004,...
 8.8000000000000007,...
 9.1999999999999993,...
 9.5999999999999996,...
 10,...
]);
set(ax,'YTick',[ -2,...
 -1.8,...
 -1.6000000000000001,...
 -1.3999999999999999,...
 -1.2,...
 -1,...
 -0.79999999999999982,...
 -0.59999999999999987,...
 -0.39999999999999991,...
 -0.19999999999999996,...
 0,...
 0.20000000000000018,...
 0.40000000000000036,...
 0.60000000000000009,...
 0.80000000000000027,...
 1,...
 1.2000000000000002,...
 1.4000000000000004,...
 1.6000000000000001,...
 1.8000000000000003,...
 2,...
 2.2000000000000002,...
 2.4000000000000004,...
 2.6000000000000005,...
 2.8000000000000007,...
 3,...
 3.1999999999999993,...
 3.3999999999999995,...
 3.5999999999999996,...
 3.7999999999999998,...
 4,...
 4.1999999999999993,...
 4.4000000000000004,...
 4.5999999999999996,...
 4.7999999999999998,...
 5,...
 5.1999999999999993,...
 5.4000000000000004,...
 5.5999999999999996,...
 5.7999999999999998,...
 6,...
 6.2000000000000002,...
 6.4000000000000004,...
 6.5999999999999996,...
 6.7999999999999998,...
 7,...
 7.2000000000000002,...
 7.4000000000000004,...
 7.5999999999999996,...
 7.7999999999999998,...
 8,...
]);
pdetool('gridon','on');

% Geometry description:
pderect([1.5 7.5 2.3999999999999999 2.6000000000000001],'R2');
pdepoly([ 2,...
 4.5,...
 7,...
],...
[ 2.6000000000000001,...
 3.4300000000000002,...
 2.6000000000000001,...
],...
 'P2');
pderect([4 4.7999999999999998 0 1],'R3');
pderect([2 2.2000000000000002 2.3999999999999999 0],'R4');
pderect([6.7999999999999998 7 2.3999999999999999 0],'R5');
pdepoly([ 1.5,...
 2,...
 4.5,...
 7,...
 7.5,...
 4.5,...
],...
[ 2.6000000000000001,...
 2.6000000000000001,...
 3.4300000000000002,...
 2.6000000000000001,...
 2.6000000000000001,...
 3.6000000000000001,...
],...
 'P1');
pderect([1.2 7.7999999999999998 0 -0.40000000000000002],'R6');
pdepoly([ 1.2000000000000002,...
 1.2000000000000002,...
 7.7999999999999998,...
 7.7999999999999998,...
 8.8000000000000007,...
 8.8000000000000007,...
 0.40000000000000002,...
 0.40000000000000002,...
],...
[ 0,...
 -0.39999999999999991,...
 -0.39999999999999991,...
 0,...
 0,...
 -1,...
 -1,...
 0,...
],...
 'P3');
pdepoly([ 4.5,...
 7.5,...
 7.5,...
 4.5,...
 1.5,...
 1.5,...
],...
[ 3.7999999999999998,...
 2.8000000000000007,...
 2.6000000000000005,...
 3.6000000000000001,...
 2.6000000000000005,...
 2.8000000000000007,...
],...
 'P4');
pdepoly([ 2.2000000000000002,...
 4,...
 4,...
 4.8000000000000007,...
 4.8000000000000007,...
 6.7999999999999998,...
 6.7999999999999998,...
 2.2000000000000002,...
],...
[ 0,...
 0,...
 1,...
 1,...
 0,...
 0,...
 2.4000000000000004,...
 2.4000000000000004,...
],...
 'P5');
set(findobj(get(pde_fig,'Children'),'Tag','PDEEval'),'String','R2+P2+R3+R4+R5+P1+R6+P3+P4+P5')

% Boundary conditions:
pdetool('changemode',0)
pdesetbd(39,...
'dir',...
1,...
'1',...
'253')
pdesetbd(35,...
'dir',...
1,...
'1',...
'253')
pdesetbd(33,...
'dir',...
1,...
'1',...
'253')
pdesetbd(32,...
'dir',...
1,...
'1',...
'253')
pdesetbd(26,...
'dir',...
1,...
'1',...
'253')
pdesetbd(25,...
'dir',...
1,...
'1',...
'253')
pdesetbd(21,...
'dir',...
1,...
'1',...
'253')
pdesetbd(20,...
'dir',...
1,...
'1',...
'253')
pdesetbd(17,...
'dir',...
1,...
'1',...
'253')
pdesetbd(16,...
'dir',...
1,...
'1',...
'253')
pdesetbd(7,...
'dir',...
1,...
'1',...
'253')
pdesetbd(6,...
'dir',...
1,...
'1',...
'253')
pdesetbd(5,...
'dir',...
1,...
'1',...
'273')
pdesetbd(4,...
'dir',...
1,...
'1',...
'273')
pdesetbd(3,...
'dir',...
1,...
'1',...
'273')
pdesetbd(2,...
'dir',...
1,...
'1',...
'253')
pdesetbd(1,...
'dir',...
1,...
'1',...
'253')

% Mesh generation:
setappdata(pde_fig,'Hgrad',1.3);
setappdata(pde_fig,'refinemethod','regular');
setappdata(pde_fig,'jiggle',char('on','mean',''));
setappdata(pde_fig,'MesherVersion','preR2013a');
pdetool('initmesh')
pdetool('refine')
pdetool('refine')

% PDE coefficients:
pdeseteq(2,...
'0.2!0.2!0.4!0.15!0.026!0.15!1.7!47!47!0.026',...
'0!0!0!0!0!0!0!0!15!0',...
'(0)+(0).*(253)!(0)+(0).*(253)!(0)+(0).*(273)!(0)+(0).*(253)!(0)+(0).*(253)!(0)+(0).*(253)!(0)+(0).*(273)!(0)+(0).*(253)!(30000)+(15).*(273)!(0)+(0).*(273)',...
'(500).*(840)!(500).*(840)!(1500).*(850)!(200).*(2090)!(1.2).*(1005)!(450).*(2700)!(2500).*(840)!(7900).*(462)!(7900).*(462)!(1.2).*(1005)',...
'0:5000:300000',...
'253',...
'0.0',...
'[0 100]')
setappdata(pde_fig,'currparam',...
['500!500!1500!200!1.2!450!2500!7900!7900!1.2';...
'840!840!850!2090!1005!2700!840!462!462!1005';...
'0.2!0.2!0.4!0.15!0.026!0.15!1.7!47!47!0.026';...
'0!0!0!0!0!0!0!0!30000!0                    ';...
'0!0!0!0!0!0!0!0!15!0                       ';...
'253!253!273!253!253!253!273!253!273!273    '])

% Solve parameters:
setappdata(pde_fig,'solveparam',...
char('0','28056','10','pdeadworst',...
'0.5','longest','0','1E-4','','fixed','Inf'))

% Plotflags and user data strings:
setappdata(pde_fig,'plotflags',[1 1 1 1 2 1 7 1 1 0 0 61 1 1 0 1 0 1]);
setappdata(pde_fig,'colstring','');
setappdata(pde_fig,'arrowstring','');
setappdata(pde_fig,'deformstring','');
setappdata(pde_fig,'heightstring','');

% Solve PDE:
pdetool('solve')