matlab 奇异点,MATLAB含有奇异点的边值问题求解

% 此问题的关键在于方程中含有奇异点 r=0，一般参考例子均是求解不含有奇异点的边

%值问题，这里特意阐述一下奇异值的处理方法。代码如下：

关于奇异点处理的完整代码请文末企鹅号

function bianzhi_example

EE=Const;R1=Const;g=Const;K=Con;A=Con;evi=Con;KT=Con;

C=Con;C1=Con;C2=Con;gex=Con;Tw=Con;Vex=evi;

% 上述为方程参数

T = linspace(0,R1,15);%为生成初始化网络准备，网格数量可以修改

function yinit =

mat4init(t)

%对y初值的估计函数,由于y1(R1) = Tw,y2(0) =

0;假设一个满足上述条件的函数，这里

%函数有很多种，随机取一种计算，效果不佳可更换，初始函数选择是否合适可能对结果

%有较大影响

yinit = [  t^2/R1^2+Tw-1;

2*t/R1^2 ];

%

yinit = [

Tw;

%

0 ];

%另一种初设函数

end