a='rmgcbwyk-'
for lamda=0.0015509-0.0016*4:0.0016:0.0015509+0.0016*4
x=-300:300
I=(sin((M*pi*ns*d./lamda).*(x./R-nw*dertaL/ns*d))./sin((pi*ns*d./lamda).*(x./R-nw*dertaL/ns*d))).^2
[y,ps] = mapminmax(I)
hold on
plot(x,y,a())
end
(注:参数定义没写上去),倒数第二行a()里应该填上什么啊,怎么让图像每回都不一样啊