第四章作业讲解佳.pptxVIP

第四章作业讲解王文佳习题4.11.问题描述：利用MATLAB编程作出一个朗伯光源的辐射方向图和一个由 给出的辐射方向图。假定两个光源有相同的峰值辐射强度 ，并且两个辐射方向图都已归一化。2.输入输出描述：输入：输出：分别作出它们的辐射方向图。习题4.13.过程分析朗伯光源的发射方向图关系表达式:另一个光源发射方向图关系表达式:在MATLAB中即可画出其辐射方向图。习题4.14.matlab实现及结果展示clcclearclose all N=1000;theta=linspace(-pi/2,pi/2,N);B1=cos(theta);B2=cos(theta).^5;polar(theta,B1)hold onpolar(theta,B2,--r)习题4.14.matlab实现图像如下：习题4.21.问题描述：辐射方向图由 给出，利用MATLAB编程作出B( ),作为m的函数曲线，m的取值范围是[1,21]的奇数，视角为10°、20°、45°、60°，假设所有光源都有同样的峰值辐射强度 .2.输入输出描述：输入：m为1到21的奇数，视角为10°、20°、45°、60°，输出：利用MATLAB编程作出作为m的函数曲线。3.过程分析把 =10°、20°、45°、60°时，分别代入式进而画出 作为m的四条函数曲线。习题4.24.matlab实现及结果展示clc clear close all m=1:2:21;theta=[10,20,45,60]a=theta*pi/180;B=zeros(length(m),length(a)); for i = 1:1:length(a) B(:,i) = (cos(a(i))).^m;end plot(m,B,-o)xlabel(m)ylabel(B(\theta))习题4.24.matlab实现及结果展示图像如下：习题4.24.matlab实现及结果展示结果显示：图中蓝线、绿线、红线、青线分别为=10°、20°、45°、60°时， 作为m的四条函数曲线。 较小时， 随m的变化幅度较小。随着 的增大， 随m的变化幅度增大。习题4.31.问题描述：圆形发射区域半径为25μm LED有个朗伯辐射方向图，在80mA的驱动电流下有80轴向辐射强度，试用MATLAB编程计算有多少光功率能够耦合进纤芯直径为100 μm 、NA=0.22的阶跃折射率光纤中？ 2.输入输出描述：输入： rs=0.0025cm; B0 = 80 a = 0.0050cm; NA = 0.22;输出：3、过程分析习题4.34.matlab实现及结果展示clcclear close allrs = 0.0025; B0 = 80;a = 0.0050; NA = 0.22; Pstep = pi^2*rs^2*B0*(NA)^2;Pstep习题4.34.matlab实现及结果展示结果显示：有0.23884mW的光功率能够耦合进阶跃折射率光纤中。习题4.44、MATLAB实现：clc clear close all t = 0:0.1:2;eta1 = 2/pi*acos(t/2)-t/pi.*sqrt(1-(t/2).^2); for i = 1:length(t) eta2(i) = quad(@(x)(sqrt(1-x.^2)),t(i)/2,1)/pi*4;end a = abs(eta2-eta1);err = a./eta1;plot(t,err,r)xlabel(t=d/a)ylabel(err)图像如下：习题4.5习题4.5习题4.5习题4.5习题4.54、MATLAB实现：clc clear close all t = 0:0.1:2;eta1 = 2/pi*acos(t/2)-t/pi.*sqrt(1-(t/2).^2); for i = 1:length(t) eta3(i) = dblquad(@(x,y)(sqrt(x.^2+y.^2)=1 ... sqrt((x-t(i)).^2+y.^2)=1),-1,1,-1,1)/pi;end a = abs(eta3-eta1);err = a./eta1;plot(t,err,r)xlabel(t=d/a)ylabel(err)图像如下：matlab知识quad() 函数的使用法：quad(fun ,a, b)。其中，fun是一个函数的句柄，参数a,b分别是积分的下限与上限。dblquad() 函数的使用法：dblquad(fun ,xmin, xmax, ymin, ymax)。其中，fun是一个函数的句柄，xmin, xmax分别是变量x的下限与上限。ymin, ymax分别是变量y的下限与上限。阶跃型光纤（Step