单缝衍射的相对光强分布matlab代码

disp("单缝衍射的相对光强分布")
disp("数据处理与分析")
disp("1.画出单缝衍射相对光强分布曲线")
X = X';
Y = Y';
disp("出光口的位置l0，狭缝位置l1，光探头位置l2,狭缝到探头的距离Z1(mm)：")
l0 = 170
l1 = 250
l2 = 950
Z1 = l2 - l1
disp("把x化成sinθ")
sin_ = X./sqrt(Z^2 + X.^2);
disp("求出相对光强：")
Y_ = Y./Y(1);
plot(sin_,Y_,'o');
title("激光单缝衍射光强相对分布曲线")
xlabel("sinθ")
ylabel("Iθ/I0")
grid on;
hold on;
plot(sin_,Y_);
hold off;

disp("2.计算各次极大光强与主极大光强的比值，与理论值比较")
disp("理论比值：")
t1_std = 0.0469
t2_std = 0.0166
disp("第二次极大光强与主极大值的比值：")
t1 = 0.049
disp("第三次极大光强与主极大值的比值：")
t2 = 0.0189
disp("差异的原因：1.测量误差，2.其他组做实验带来的影响")

disp("3.绘制白光光强P与距离平方（1/Z^2)的关系曲线")
P = double(P);
P = P';
Z2 = 200:50:900;
Z_ = 1./Z2.^2;
plot(Z_,P,'o');
title("白光光强P与距离平方（1/Z^2)的关系曲线");
xlabel("1/Z²(mm)");
ylabel("P(uW)");
grid on;
hold on;
plot(Z_,P);
hold off;
disp("在误差允许的范围内，可认为探头在线性工作区。")