matlab拟合杨氏双缝干涉光强分布曲线

I_0为光强极值，d为双缝间距，lamda为波长，D1为双缝所在平面到光屏距离。

%Intensiy distribution
I_0=10;
d=0.02;D1=50;lamda=5e-7;
x=-67:67;
k=pi*d/(lamda*D1);

I=4*I_0*cosd(k*x).*cosd(k*x)

plot([0,0],[0,45],'k')
hold on
plot(x,I,'b')
hold off
text(40,42,'光强分布曲线I')
text(2,43,'y轴')