單頻溫度場計算第11講（2學(xué)時）單頻溫度場計算基本方程其中：為溫度變化，為熱擴散率，為熱導(dǎo)率，

邊界條件為熱源其中：為樣品對入射光的吸收系數(shù)，為入射激發(fā)光的光強，為樣品材料的反射率，為光斑半徑。為空間熱源——————————解————————————基本參數(shù)：）計算的體系結(jié)構(gòu):functiony=Temain(U,f)（主程序，計算r,z向的溫度幅度分布）--------functiony=Temp(U,r,z)

（溫度計算程序，計算在（r,z）點的溫度幅度）----------------functiony=Gauss(U,r,z)(高斯積分程序，計算溫度場的高斯積分)------------------------functiony=Grxz(U,x,r,z)

（積分函數(shù)形式程序，同時實現(xiàn)積分區(qū)間轉(zhuǎn)換）--------------------------------functiony=Tqz(U,q,z)

（解公式的形式程序，計算具體公式）functiony=Temain(U,f)%U=[beta,P,kapa,alfa,L]hr=2e-4;hz=2e-4;forii=1:20forjj=1:20r=(ii-1)*hr;z=(jj-1)*hj;y(ii,jj)=Temp(U,r,z);endendplot(y)functiony=Temp(U,r,z)y=Gauss(U,r,z);主要程序結(jié)構(gòu)圖:functiony=Gauss(U,r,z)A=[-0.9681602395,-0.8360311073,-0.6133714327,-0.3242534234,0.,0.3242534234,0.6133714327,0.8360311073,0.9681602395];W=[0.0812743884,0.1806481607,0.2606106964,0.3123470770,0.3302393550,0.3123470770,0.2606106964,0.1806481607,0.0812743884];h=1/128;quad=0;forii=1:128a=(ii-1)*h;b=ii*h;N=length(A);T=zeros(1,N);T=((a+b)/2)+((b-a)/2)*A;quad=quad+((b-a)/2)*sum(W.*Grxz(U,x,r,z));endfunctiony=Grxz(U,x,r,z)q=-ln(x);y=Tqz(U,log(1./x),z).*besselj(0,r*log(1./x)).*log(1./x)./x;functiony=Tqz(U,q,z)beta=U(1);P=U(2);kapa=U(3);alfa=U(4);L=U(5);kesait=H=B5=B2=B1=y=B1*exp(-kesait*z)+B2*exp(kesait*(z-L))+B5*exp(-beta*z);基本方程其中：為溫度變化，為熱擴散率，為熱導(dǎo)率，

邊界條件為熱源其中：為樣品對入射光的吸收系數(shù)，為光斑半徑。為空間熱源測試題目要求:(1)寫出程序體系結(jié)構(gòu)或框架;(2)寫出程序的代碼;——————————解————————————基本參數(shù)：）A=[-0.9681602395,-0.8360311073,-0.6133714327,-0.3242534234,0.,0.3242534234,0.6133714327,0.8360311073,0.9681602395];W=[0.0812743884,0.1806481607,0.26061069