该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
clear all;
%求贝塞尔函数零点
maxv = 2;
maxs =44;
MAXS =100;
j= zeros(maxv, MAXS);
incr =4;
for v=0:maxv-1
if (v==0)
j(v+1,1) = fzero(@(x)besselj(v,x),2);
else
j(v+1,1) = fzero(@(x)besselj(v,x),3.8);
end
for s=2:maxs
j(v+1,s) = fzero(@(x)besselj(v,x),j(v+1,s-1)+incr);
end
end
for p=maxs+1:MAXS
j(1,p)=j(1,maxs)+pi*(p-maxs);
j(2,p)=j(2,maxs)+pi*(p-maxs);
end
T0 = 22;
k=17;%导热系数w/cm*k
r=0;%求解位置
t=1;%求解时间
b=0.5;%材料半径cm
a=0.025;%激光光斑半径cm
A=11.9;%材料热扩散率
d=1 ;%吸收率
h=0.01;%材料厚度cm
%画网格
Q=0.01;%间隔
xx=-0.3:Q:0.3