Я пытаюсь решить уравнение теплопроводности u t + Δ u = e u в 1 измерении от 0 до 1. Я хочу использовать IC и BC как 0 надомен.
Чтобы решить PDE, я хочу с нуля написать код конечного элемента, чтобы помочь мне лучше понять методы MATLAB и конечного элемента.Тем не менее, я борюсь с созданием кусочных функций, которые мы можем использовать в качестве меша.
npoint=5 %points in decomposition - points on our line
nelement=npoint-1 %number of intervals on our line
x=linspace(0,1,npoint) %creates vertices
intervals(1:nelement,1)=1:nelement %creates e2p, part 1????
(1:nelement,2)=2:npoint %creates e2p, part 2
%e2p this descirbes the sections which will be approximated, e.g. area of
%bounday 1,2
figure=plot(x,0*x, 'b-o', 'MarkerFaceColor','r') %draw composition
figure
Если бы кто-то мог посоветовать, как мне продолжить работу с кодом конечных элементов, или указать мне направление хорошего прохода, это было бы очень полезно.Большое спасибо.